Changeset 740d430 in git
- Timestamp:
- Jun 6, 2016, 4:43:16 PM (8 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 479f808, 4a0d9d5
- Parents:
- e2339af
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
cards/CMS_PhaseII/CMS_PhaseII_200PU.tcl
re2339af r740d430 1 set MaxEvents 10 2 1 3 # 2 4 # Phase II - Pile-Up 3 5 # 4 6 # Main authors: Michele Selvaggi (UCL) 5 # 6 # Released on: 7 # 8 # Released on: 7 9 # 8 10 # Version: v01 … … 14 16 15 17 set ExecutionPath { 16 17 PileUpMerger 18 18 19 PileUpMerger 19 20 ParticlePropagator 20 21 … … 28 29 29 30 TrackMerger 30 31 31 32 ECal 32 33 HCal … … 34 35 ElectronFilter 35 36 TrackPileUpSubtractor 36 37 37 38 TowerMerger 38 39 NeutralEFlowMerger … … 45 46 RunPUPPI 46 47 48 PhotonFilter 47 49 48 50 PhotonIsolation 49 51 PhotonEfficiency 50 52 51 53 ElectronIsolation 52 54 ElectronEfficiency 53 55 54 56 MuonIsolation 55 57 MuonLooseIdEfficiency 56 58 MuonTightIdEfficiency 57 59 58 60 NeutrinoFilter 59 61 … … 67 69 68 70 ScalarHT 69 71 70 72 JetEnergyScale 71 73 … … 77 79 78 80 TauTagging 79 81 80 82 GenParticleFilter 81 83 82 84 TreeWriter 83 85 } … … 96 98 97 99 # pre-generated minbias input file 98 set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup 100 set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup 99 101 100 102 # average expected pile up 101 103 set MeanPileUp 200 102 104 103 105 # maximum spread in the beam direction in m 104 106 set ZVertexSpread 0.25 … … 109 111 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s) 110 112 set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))} 111 113 112 114 } 113 115 … … 165 167 set OutputArray electrons 166 168 # tracking efficiency formula for electrons 167 set EfficiencyFormula { 169 set EfficiencyFormula { 168 170 (pt <= 0.2) * (0.00) + \ 169 171 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \ … … 197 199 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \ 198 200 (abs(eta) > 4.0) * (0.00) 199 201 200 202 } 201 203 } … … 210 212 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons 211 213 set OutputArray chargedHadrons 212 # resolution formula for charged hadrons , 213 214 # resolution formula for charged hadrons , 215 214 216 # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html 215 217 source trackMomentumResolution.tcl … … 227 229 228 230 # resolution formula for electrons 229 231 230 232 # taking something flat in energy for now, ECAL will take over at high energy anyway. 231 233 # inferred from hep-ex/1306.2016 and 1502.02701 232 set ResolutionFormula { 233 234 set ResolutionFormula { 235 234 236 (abs(eta) <= 1.5) * (1+0.64*abs(eta)^2)*(energy*0.028) + 235 237 (abs(eta) > 1.5 && abs(eta) <= 1.75) * (energy*0.037) + 236 238 (abs(eta) > 1.75 && abs(eta) <= 2.15) * (energy*0.038) + 237 239 (abs(eta) > 2.15 && abs(eta) <= 3.00) * (energy*0.044) + 238 (abs(eta) > 3.00 && abs(eta) <= 4.00) * (energy*0.10)} 240 (abs(eta) > 3.00 && abs(eta) <= 4.00) * (energy*0.10)} 239 241 240 242 } … … 248 250 set OutputArray muons 249 251 # resolution formula for muons 250 251 # up to |eta| < 2.8 take measurement from tracking + muon chambers 252 # for |eta| > 2.8 and pT < 5.0 take measurement from tracking alone taken from 252 253 # up to |eta| < 2.8 take measurement from tracking + muon chambers 254 # for |eta| > 2.8 and pT < 5.0 take measurement from tracking alone taken from 253 255 # http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html 254 source muonMomentumResolution.tcl 256 source muonMomentumResolution.tcl 255 257 } 256 258 … … 282 284 set EFlowTowerOutputArray eflowPhotons 283 285 284 set IsEcal true 285 286 set IsEcal true 287 286 288 set EnergyMin 0.5 287 289 set EnergySignificanceMin 1.0 … … 296 298 297 299 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5 298 300 299 301 set PhiBins {} 300 302 for {set i -180} {$i <= 180} {incr i} { … … 309 311 310 312 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL) 311 313 312 314 set PhiBins {} 313 315 for {set i -180} {$i <= 180} {incr i} { … … 326 328 } 327 329 328 # take present CMS granularity for HF 329 330 # take present CMS granularity for HF 331 330 332 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0 331 333 set PhiBins {} … … 333 335 add PhiBins [expr {$i * $pi/18.0}] 334 336 } 335 337 336 338 foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} { 337 339 add EtaPhiBins $eta $PhiBins … … 359 361 360 362 # set ResolutionFormula {resolution formula as a function of eta and energy} 361 363 362 364 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701 363 365 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117 364 366 365 set ResolutionFormula { (abs(eta) <= 1.50) * (1+0.64*abs(eta)^2)*sqrt(energy^2*0.009^2 + energy*0.12^2 + 0.45^2) + 367 set ResolutionFormula { (abs(eta) <= 1.50) * (1+0.64*abs(eta)^2)*sqrt(energy^2*0.009^2 + energy*0.12^2 + 0.45^2) + 366 368 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \ 367 369 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \ … … 383 385 set EFlowTowerOutputArray eflowNeutralHadrons 384 386 385 set IsEcal false 386 387 set IsEcal false 388 387 389 set EnergyMin 1.0 388 390 set EnergySignificanceMin 1.0 … … 395 397 # each list starts with the lower edge of the first tower 396 398 # the list ends with the higher edged of the last tower 397 399 398 400 # assume 0.087 x 0.087 resolution in eta,phi in the barrel |eta| < 1.5 399 401 … … 407 409 408 410 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- HCAL) 409 411 410 412 set PhiBins {} 411 413 for {set i -180} {$i <= 180} {incr i} { … … 424 426 } 425 427 426 # take present CMS granularity for HF 427 428 # take present CMS granularity for HF 429 428 430 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0 429 431 set PhiBins {} … … 431 433 add PhiBins [expr {$i * $pi/18.0}] 432 434 } 433 435 434 436 foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} { 435 437 add EtaPhiBins $eta $PhiBins … … 476 478 } 477 479 480 478 481 ########################## 479 482 # Track pile-up subtractor … … 573 576 set UseExp false 574 577 set UseNoLep false 575 576 ## define puppi algorithm parameters (more than one for the same eta region is possible) 578 579 ## define puppi algorithm parameters (more than one for the same eta region is possible) 577 580 add EtaMinBin 0.0 1.5 4.0 578 581 add EtaMaxBin 1.5 4.0 10.0 … … 585 588 add ApplyCHS true true true 586 589 add UseCharged true true false 587 add ApplyLowPUCorr true true true 590 add ApplyLowPUCorr true true true 588 591 add MetricId 5 5 5 589 592 add CombId 0 0 0 … … 593 596 set OutputArrayTracks puppiTracks 594 597 set OutputArrayNeutrals puppiNeutrals 595 } 598 } 596 599 597 600 module Merger RunPUPPI { … … 620 623 621 624 ################### 622 # Ger PileUp Missing ET 625 # Ger PileUp Missing ET 623 626 ################### 624 627 … … 667 670 module FastJetFinder GenJetFinder { 668 671 set InputArray NeutrinoFilter/filteredParticles 669 672 670 673 set OutputArray jets 671 674 … … 720 723 721 724 725 ################# 726 # Photon filter 727 ################# 728 729 module PdgCodeFilter PhotonFilter { 730 set InputArray ECal/eflowPhotons 731 set OutputArray photons 732 set Invert true 733 set PTMin 5.0 734 add PdgCode {22} 735 } 736 737 722 738 #################### 723 739 # Photon isolation # … … 725 741 726 742 module Isolation PhotonIsolation { 727 743 728 744 # particle for which calculate the isolation 729 set CandidateInputArray ECal/eflowPhotons730 745 set CandidateInputArray PhotonFilter/photons 746 731 747 # isolation collection 732 set IsolationInputArray EFlowMerger/eflow733 748 set IsolationInputArray RunPUPPI/PuppiParticles 749 734 750 # output array 735 751 set OutputArray photons 736 752 753 # veto isolation cand. based on proximity to input cand. 754 set DeltaRMin 0.01 755 set UseMiniCone true 756 737 757 # isolation cone 738 758 set DeltaRMax 0.3 739 740 # minimum pT 759 760 # minimum pT 741 761 set PTMin 1.0 742 762 743 763 # iso ratio to cut 744 764 set PTRatioMax 9999. 765 745 766 } 746 767 … … 752 773 753 774 module Efficiency PhotonEfficiency { 754 775 755 776 ## input particles 756 set InputArray PhotonIsolation/photons 777 set InputArray PhotonIsolation/photons 757 778 ## output particles 758 779 set OutputArray photons … … 772 793 773 794 module Isolation ElectronIsolation { 774 795 775 796 set CandidateInputArray ElectronFilter/electrons 776 797 777 798 # isolation collection 778 set IsolationInputArray EFlowMerger/eflow 779 799 set IsolationInputArray RunPUPPI/PuppiParticles 800 #set IsolationInputArray EFlowMerger/eflow 801 780 802 set OutputArray electrons 781 803 782 804 set DeltaRMax 0.3 783 805 set PTMin 1.0 … … 793 815 794 816 module Efficiency ElectronEfficiency { 795 817 796 818 set InputArray ElectronIsolation/electrons 797 819 set OutputArray electrons 798 820 799 821 # set EfficiencyFormula {efficiency formula as a function of eta and pt} 800 822 # efficiency formula for electrons … … 804 826 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \ 805 827 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \ 806 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \ 807 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \ 808 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \ 809 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \ 828 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \ 829 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \ 830 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \ 831 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \ 810 832 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \ 811 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \ 812 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \ 813 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \ 833 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \ 834 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \ 835 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \ 814 836 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \ 815 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \ 816 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \ 817 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \ 818 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \ 837 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \ 838 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \ 839 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \ 840 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \ 819 841 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \ 820 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \ 821 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \ 822 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \ 823 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \ 842 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \ 843 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \ 844 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \ 845 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \ 824 846 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \ 825 847 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \ … … 828 850 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \ 829 851 (abs(eta) > 4.0) * (0.00) 830 852 831 853 } 832 854 } … … 838 860 module Isolation MuonIsolation { 839 861 set CandidateInputArray MuonMomentumSmearing/muons 840 862 841 863 # isolation collection 842 set IsolationInputArray EFlowMerger/eflow843 864 set IsolationInputArray RunPUPPI/PuppiParticles 865 844 866 set OutputArray muons 845 867 846 868 set DeltaRMax 0.3 847 869 set PTMin 1.0 … … 916 938 917 939 set BitNumber 1 918 940 919 941 source btagMedium.tcl 920 942 } 921 943 922 923 944 924 945 ############# … … 931 952 set BitNumber 2 932 953 933 source btagTight.tcl 934 } 935 936 937 954 source btagTight.tcl 955 } 938 956 939 957 … … 961 979 add EfficiencyFormula {15} { (abs(eta) < 2.3) * 0.97*0.77*( (0.32 + 0.01*pt - 0.000054*pt*pt )*(pt<100)+0.78*(pt>100) ) + \ 962 980 (abs(eta) > 2.3) * (0.000) 963 } 981 } 964 982 } 965 983 … … 970 988 971 989 module StatusPidFilter GenParticleFilter { 972 990 973 991 set InputArray Delphes/allParticles 974 992 set OutputArray filteredParticles 975 993 set PTMin 5.0 976 994 977 995 } 978 996 … … 986 1004 add Branch GenParticleFilter/filteredParticles Particle GenParticle 987 1005 add Branch PileUpMerger/vertices Vertex Vertex 988 1006 989 1007 add Branch GenJetFinder/jets GenJet Jet 990 1008 add Branch GenMissingET/momentum GenMissingET MissingET … … 993 1011 # add Branch ECal/eflowPhotons EFlowPhoton Tower 994 1012 # add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower 995 1013 996 1014 add Branch PhotonEfficiency/photons Photon Photon 997 1015 add Branch ElectronEfficiency/electrons Electron Electron 998 1016 add Branch MuonLooseIdEfficiency/muons MuonLoose Muon 999 1017 add Branch MuonTightIdEfficiency/muons MuonTight Muon 1000 1018 1001 1019 add Branch JetEnergyScale/jets Jet Jet 1002 1020 1003 1021 add Branch MissingET/momentum MissingET MissingET 1004 1022 add Branch PuppiMissingET/momentum PuppiMissingET MissingET
Note:
See TracChangeset
for help on using the changeset viewer.