Changes in cards/FCC/FCChh.tcl [96c959a:eb332c8] in git
- File:
-
- 1 edited
-
cards/FCC/FCChh.tcl (modified) (22 diffs)
Legend:
- Unmodified
- Added
- Removed
-
cards/FCC/FCChh.tcl
r96c959a reb332c8 4 4 # Main authors: Michele Selvaggi (CERN) 5 5 # 6 # Released on: Dec. 1st, 20176 # Released on: Apr. 6th, 2017 7 7 # 8 8 # Configuration: FCC-hh baseline detector … … 13 13 14 14 set ExecutionPath { 15 15 16 16 ParticlePropagator 17 17 … … 24 24 MuonMomentumSmearing 25 25 26 DenseProp27 DenseMergeTracks28 TrackMergerProp29 26 TrackMerger 30 27 … … 51 48 GenMissingET 52 49 53 GenJetFinder02 54 GenJetFinder04 55 GenJetFinder08 56 57 FastJetFinder02 58 FastJetFinder04 59 FastJetFinder08 60 61 CaloJetFinder02 62 CaloJetFinder04 63 CaloJetFinder08 64 65 TrackJetFinder02 66 TrackJetFinder04 67 TrackJetFinder08 50 GenJetFinder 51 FastJetFinder 52 FatJetFinder 68 53 69 54 JetEnergyScale … … 82 67 } 83 68 84 #85 69 ################################# 86 70 # Propagate particles in cylinder … … 89 73 module ParticlePropagator ParticlePropagator { 90 74 set InputArray Delphes/stableParticles 75 91 76 set OutputArray stableParticles 92 77 set ChargedHadronOutputArray chargedHadrons … … 202 187 } 203 188 204 205 189 ############## 206 190 # Track merger 207 191 ############## 208 192 209 module Merger TrackMerger Prop{193 module Merger TrackMerger { 210 194 # add InputArray InputArray 211 195 add InputArray ChargedHadronMomentumSmearing/chargedHadrons … … 214 198 set OutputArray tracks 215 199 } 216 217 218 #################################219 # Dense Track propagation220 #################################221 222 module ParticlePropagator DenseProp {223 set InputArray Delphes/stableParticles224 225 set OutputArray stableParticles226 set NeutralOutputArray neutralParticles227 set ChargedHadronOutputArray chargedHadrons228 set ElectronOutputArray electrons229 set MuonOutputArray muons230 231 # radius of the magnetic field coverage, in m232 set Radius 0.45233 set RadiusMax 1.5234 # half-length of the magnetic field coverage, in m235 set HalfLength 0.8236 set HalfLengthMax 5237 238 # magnetic field239 set Bz 4.0240 }241 242 ####################243 # Dense Track merger244 ###################245 246 module Merger DenseMergeTracks {247 # add InputArray InputArray248 add InputArray DenseProp/chargedHadrons249 add InputArray DenseProp/electrons250 add InputArray DenseProp/muons251 set OutputArray tracks252 }253 254 255 ######################256 # Dense Track Filter257 ######################258 259 module DenseTrackFilter TrackMerger {260 set TrackInputArray TrackMergerProp/tracks261 set DenseChargedInputArray DenseMergeTracks/tracks262 263 set TrackOutputArray tracks264 265 set EtaPhiRes 0.001266 set EtaMax 6.0267 268 set pi [expr {acos(-1)}]269 270 set nbins_phi [expr {$pi/$EtaPhiRes} ]271 set nbins_phi [expr {int($nbins_phi)} ]272 273 set PhiBins {}274 for {set i -$nbins_phi} {$i <= $nbins_phi} {incr i} {275 add PhiBins [expr {$i * $pi/$nbins_phi}]276 }277 278 set nbins_eta [expr {$EtaMax/$EtaPhiRes} ]279 set nbins_eta [expr {int($nbins_eta)} ]280 281 for {set i -$nbins_eta} {$i <= $nbins_eta} {incr i} {282 set eta [expr {$i * $EtaPhiRes}]283 add EtaPhiBins $eta $PhiBins284 }285 }286 287 288 200 289 201 ############# … … 488 400 add InputArray ECal/ecalTowers 489 401 add InputArray HCal/hcalTowers 490 add InputArray MuonMomentumSmearing/muons491 402 set OutputArray towers 492 403 } … … 562 473 563 474 475 ##################### 476 # MC truth jet finder 477 ##################### 478 479 # TBC: is jet radius fine? 480 481 module FastJetFinder GenJetFinder { 482 # set InputArray NeutrinoFilter/filteredParticles 483 set InputArray Delphes/stableParticles 484 485 set OutputArray jets 486 487 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 488 set JetAlgorithm 6 489 set ParameterR 0.4 490 491 set JetPTMin 5.0 492 } 493 564 494 ######################### 565 495 # Gen Missing ET merger … … 574 504 575 505 576 ##################### 577 # MC truth jet finder 578 ##################### 579 580 # TBC: is jet radius fine? 581 582 module FastJetFinder GenJetFinder02 { 583 set InputArray NeutrinoFilter/filteredParticles 506 507 ############ 508 # Jet finder 509 ############ 510 511 # TBC need to include jet substructure variables 512 # TBC is jet radius fine? 513 514 module FastJetFinder FastJetFinder { 515 # set InputArray Calorimeter/towers 516 set InputArray EFlowMerger/eflow 517 518 set OutputArray jets 519 520 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 521 # 7: anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness 522 523 set JetAlgorithm 6 524 set ParameterR 0.4 525 526 set JetPTMin 30.0 527 } 528 529 ################## 530 # Fat Jet finder 531 ################## 532 533 module FastJetFinder FatJetFinder { 534 set InputArray EFlowMerger/eflow 584 535 585 536 set OutputArray jets … … 587 538 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 588 539 set JetAlgorithm 6 589 set ParameterR 0. 2540 set ParameterR 0.8 590 541 591 542 set ComputeNsubjettiness 1 … … 593 544 set AxisMode 4 594 545 595 set ComputeSoftDrop 1 596 set BetaSoftDrop 0.0 597 set SymmetryCutSoftDrop 0.1 598 set R0SoftDrop 0.2 599 600 set JetPTMin 25.0 601 } 602 603 604 ##################### 605 # MC truth jet finder 606 ##################### 607 608 # TBC: is jet radius fine? 609 610 module FastJetFinder GenJetFinder04 { 611 set InputArray NeutrinoFilter/filteredParticles 612 613 set OutputArray jets 614 615 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 616 set JetAlgorithm 6 617 set ParameterR 0.4 618 619 set ComputeNsubjettiness 1 620 set Beta 1.0 621 set AxisMode 4 622 623 set ComputeSoftDrop 1 624 set BetaSoftDrop 0.0 625 set SymmetryCutSoftDrop 0.1 626 set R0SoftDrop 0.4 627 628 set JetPTMin 25.0 629 } 630 ##################### 631 # MC truth jet finder 632 ##################### 633 634 # TBC: is jet radius fine? 635 636 module FastJetFinder GenJetFinder08 { 637 set InputArray NeutrinoFilter/filteredParticles 638 639 set OutputArray jets 640 641 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 642 set JetAlgorithm 6 643 set ParameterR 0.8 644 645 set ComputeNsubjettiness 1 646 set Beta 1.0 647 set AxisMode 4 546 set ComputeTrimming 1 547 set RTrim 0.2 548 set PtFracTrim 0.05 549 550 set ComputePruning 1 551 set ZcutPrun 0.1 552 set RcutPrun 0.5 553 set RPrun 0.8 648 554 649 555 set ComputeSoftDrop 1 … … 652 558 set R0SoftDrop 0.8 653 559 654 set JetPTMin 25.0 655 } 656 657 658 ################## 659 # Fast Jet finder 660 ################## 661 662 module FastJetFinder FastJetFinder02 { 663 set InputArray EFlowMerger/eflow 664 665 set OutputArray jets 666 667 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 668 set JetAlgorithm 6 669 set ParameterR 0.2 670 671 set ComputeNsubjettiness 1 672 set Beta 1.0 673 set AxisMode 4 674 675 set ComputeSoftDrop 1 676 set BetaSoftDrop 0.0 677 set SymmetryCutSoftDrop 0.1 678 set R0SoftDrop 0.2 679 680 set JetPTMin 25.0 681 } 682 683 ################## 684 # Fast Jet finder 685 ################## 686 687 module FastJetFinder FastJetFinder04 { 688 set InputArray EFlowMerger/eflow 689 690 set OutputArray jets 691 692 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 693 set JetAlgorithm 6 694 set ParameterR 0.4 695 696 set ComputeNsubjettiness 1 697 set Beta 1.0 698 set AxisMode 4 699 700 set ComputeSoftDrop 1 701 set BetaSoftDrop 0.0 702 set SymmetryCutSoftDrop 0.1 703 set R0SoftDrop 0.4 704 705 set JetPTMin 25.0 706 } 707 708 709 ################## 710 # Fast Jet finder 711 ################## 712 713 module FastJetFinder FastJetFinder08 { 714 set InputArray EFlowMerger/eflow 715 716 set OutputArray jets 717 718 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 719 set JetAlgorithm 6 720 set ParameterR 0.8 721 722 set ComputeNsubjettiness 1 723 set Beta 1.0 724 set AxisMode 4 725 726 set ComputeSoftDrop 1 727 set BetaSoftDrop 0.0 728 set SymmetryCutSoftDrop 0.1 729 set R0SoftDrop 0.8 730 731 set JetPTMin 25.0 732 } 733 734 735 736 ################## 737 # Fast Jet finder 738 ################## 739 740 module FastJetFinder CaloJetFinder02 { 741 set InputArray Calorimeter/towers 742 743 set OutputArray jets 744 745 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 746 set JetAlgorithm 6 747 set ParameterR 0.2 748 749 set ComputeNsubjettiness 1 750 set Beta 1.0 751 set AxisMode 4 752 753 set ComputeSoftDrop 1 754 set BetaSoftDrop 0.0 755 set SymmetryCutSoftDrop 0.1 756 set R0SoftDrop 0.2 757 758 set JetPTMin 25.0 759 } 760 761 ################## 762 # Fast Jet finder 763 ################## 764 765 module FastJetFinder CaloJetFinder04 { 766 set InputArray Calorimeter/towers 767 768 set OutputArray jets 769 770 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 771 set JetAlgorithm 6 772 set ParameterR 0.4 773 774 set ComputeNsubjettiness 1 775 set Beta 1.0 776 set AxisMode 4 777 778 set ComputeSoftDrop 1 779 set BetaSoftDrop 0.0 780 set SymmetryCutSoftDrop 0.1 781 set R0SoftDrop 0.4 782 783 set JetPTMin 25.0 784 } 785 786 787 ################## 788 # Fast Jet finder 789 ################## 790 791 module FastJetFinder CaloJetFinder08 { 792 set InputArray Calorimeter/towers 793 794 set OutputArray jets 795 796 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 797 set JetAlgorithm 6 798 set ParameterR 0.8 799 800 set ComputeNsubjettiness 1 801 set Beta 1.0 802 set AxisMode 4 803 804 set ComputeSoftDrop 1 805 set BetaSoftDrop 0.0 806 set SymmetryCutSoftDrop 0.1 807 set R0SoftDrop 0.8 808 809 set JetPTMin 25.0 810 } 811 812 813 ################## 814 # Fast Jet finder 815 ################## 816 817 module FastJetFinder TrackJetFinder02 { 818 set InputArray TrackMerger/tracks 819 820 set OutputArray jets 821 822 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 823 set JetAlgorithm 6 824 set ParameterR 0.2 825 826 set ComputeNsubjettiness 1 827 set Beta 1.0 828 set AxisMode 4 829 830 set ComputeSoftDrop 1 831 set BetaSoftDrop 0.0 832 set SymmetryCutSoftDrop 0.1 833 set R0SoftDrop 0.2 834 835 set JetPTMin 25.0 836 } 837 838 ################## 839 # Fast Jet finder 840 ################## 841 842 module FastJetFinder TrackJetFinder04 { 843 set InputArray TrackMerger/tracks 844 845 set OutputArray jets 846 847 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 848 set JetAlgorithm 6 849 set ParameterR 0.4 850 851 set ComputeNsubjettiness 1 852 set Beta 1.0 853 set AxisMode 4 854 855 set ComputeSoftDrop 1 856 set BetaSoftDrop 0.0 857 set SymmetryCutSoftDrop 0.1 858 set R0SoftDrop 0.4 859 860 set JetPTMin 25.0 861 } 862 863 864 ################## 865 # Fast Jet finder 866 ################## 867 868 module FastJetFinder TrackJetFinder08 { 869 set InputArray TrackMerger/tracks 870 871 set OutputArray jets 872 873 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 874 set JetAlgorithm 6 875 set ParameterR 0.8 876 877 set ComputeNsubjettiness 1 878 set Beta 1.0 879 set AxisMode 4 880 881 set ComputeSoftDrop 1 882 set BetaSoftDrop 0.0 883 set SymmetryCutSoftDrop 0.1 884 set R0SoftDrop 0.8 885 886 set JetPTMin 25.0 560 set JetPTMin 200.0 887 561 } 888 562 … … 894 568 895 569 module EnergyScale JetEnergyScale { 896 set InputArray FastJetFinder 04/jets570 set InputArray FastJetFinder/jets 897 571 set OutputArray jets 898 572 … … 1006 680 1007 681 1008 1009 682 ########### 1010 683 # b-tagging … … 1019 692 1020 693 (pt <= 10.0) * (0.00) + 1021 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.0 1) + \1022 (abs(eta) < 2.5) * (pt > 500.0 && pt < 15000.0) * (0.01)*(1.0 - pt/15000.) + \1023 (abs(eta) < 2.5) * (pt > 15000.0) * (0.00) + \1024 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.00 75) + \1025 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0.0075)*(1.0 - pt/15000.) + \1026 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 15000.0) * (0.000) + \694 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.001) + \ 695 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.001)*(1.0 - pt/20000.) + \ 696 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \ 697 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.00075) + \ 698 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.00075)*(1.0 - pt/20000.) + \ 699 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \ 1027 700 (abs(eta) > 4.0) * (0.00)} 1028 701 … … 1030 703 1031 704 (pt <= 10.0) * (0.00) + 1032 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.0 5) + \1033 (abs(eta) < 2.5) * (pt > 500.0 && pt < 15000.0) * (0.05)*(1.0 - pt/15000.) + \1034 (abs(eta) < 2.5) * (pt > 15000.0) * (0.000) + \705 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.04) + \ 706 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.04)*(1.0 - pt/20000.) + \ 707 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \ 1035 708 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.03) + \ 1036 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0.03)*(1.0 - pt/15000.) + \1037 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 15000.0) * (0.000) + \709 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.03)*(1.0 - pt/20000.) + \ 710 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \ 1038 711 (abs(eta) > 4.0) * (0.00)} 1039 712 … … 1042 715 (pt <= 10.0) * (0.00) + 1043 716 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.85) + 1044 (abs(eta) < 2.5) * (pt > 500.0 && pt < 15000.0) * (0.85)*(1.0 - pt/15000.) +1045 (abs(eta) < 2.5) * (pt > 15000.0) * (0.000) +717 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.85)*(1.0 - pt/20000.) + 718 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + 1046 719 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.64) + 1047 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0.64)*(1.0 - pt/15000.) +1048 (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 15000.0) * (0.000) +720 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.64)*(1.0 - pt/20000.) + 721 (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + 1049 722 (abs(eta) >= 4.0) * (0.00)} 1050 723 … … 1069 742 1070 743 (pt <= 10.0) * (0.00) + 744 (abs(eta) < 4.0) * (pt > 10.0) * (0.10) + \ 745 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)} 746 747 add EfficiencyFormula {5} { 748 749 (pt <= 10.0) * (0.00) + 1071 750 (abs(eta) < 4.0) * (pt > 10.0) * (0.25) + \ 1072 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}1073 1074 add EfficiencyFormula {5} {1075 1076 (pt <= 10.0) * (0.00) +1077 (abs(eta) < 4.0) * (pt > 10.0) * (0.03) + \1078 751 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)} 1079 752 … … 1156 829 add Branch Delphes/allParticles Particle GenParticle 1157 830 831 add Branch GenJetFinder/jets GenJet Jet 1158 832 add Branch GenMissingET/momentum GenMissingET MissingET 1159 add Branch GenJetFinder02/jets GenJet Jet1160 833 1161 834 add Branch TrackMerger/tracks Track Track … … 1171 844 add Branch UniqueObjectFinder/jets Jet Jet 1172 845 1173 add Branch GenJetFinder02/jets GenJet02 Jet 1174 add Branch GenJetFinder04/jets GenJet04 Jet 1175 add Branch GenJetFinder08/jets GenJet08 Jet 1176 1177 add Branch FastJetFinder02/jets ParticleFlowJet02 Jet 1178 add Branch FastJetFinder04/jets ParticleFlowJet04 Jet 1179 add Branch FastJetFinder08/jets ParticleFlowJet08 Jet 1180 1181 add Branch CaloJetFinder02/jets CaloJet02 Jet 1182 add Branch CaloJetFinder04/jets CaloJet04 Jet 1183 add Branch CaloJetFinder08/jets CaloJet08 Jet 1184 1185 add Branch TrackJetFinder02/jets TrackJet02 Jet 1186 add Branch TrackJetFinder04/jets TrackJet04 Jet 1187 add Branch TrackJetFinder08/jets TrackJet08 Jet 846 add Branch FatJetFinder/jets FatJet Jet 1188 847 1189 848 add Branch MissingET/momentum MissingET MissingET
Note:
See TracChangeset
for help on using the changeset viewer.
