Changeset 4a0d9d5 in git for cards/CMS_PhaseII/CMS_PhaseII_200PU.tcl
- Timestamp:
- Jun 6, 2016, 4:47:16 PM (8 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- bd0a1aa
- Parents:
- 4406bf8 (diff), 740d430 (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the(diff)
links above to see all the changes relative to each parent. - File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
cards/CMS_PhaseII/CMS_PhaseII_200PU.tcl
r4406bf8 r4a0d9d5 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 … … 40 41 EFlowMerger 41 42 43 LeptonFilterNoLep 44 LeptonFilterLep 45 RunPUPPIBase 42 46 RunPUPPI 47 48 PhotonFilter 43 49 44 50 PhotonIsolation 45 51 PhotonEfficiency 46 52 47 53 ElectronIsolation 48 54 ElectronEfficiency 49 55 50 56 MuonIsolation 51 57 MuonLooseIdEfficiency 52 58 MuonTightIdEfficiency 53 59 54 60 NeutrinoFilter 55 61 56 62 MissingET 63 PuppiMissingET 57 64 GenMissingET 58 65 GenPileUpMissingET … … 62 69 63 70 ScalarHT 64 71 65 72 JetEnergyScale 66 73 … … 72 79 73 80 TauTagging 74 81 75 82 GenParticleFilter 76 83 77 84 TreeWriter 78 85 } … … 91 98 92 99 # pre-generated minbias input file 93 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 94 101 95 102 # average expected pile up 96 103 set MeanPileUp 200 97 104 98 105 # maximum spread in the beam direction in m 99 set ZVertexSpread 0. 15106 set ZVertexSpread 0.25 100 107 101 108 # maximum spread in time in s 102 set TVertexSpread 1.5E-09109 set TVertexSpread 800E-12 103 110 104 111 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s) 105 106 #set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))} 107 108 109 set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) + 110 (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) + 111 (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) + 112 (abs(t) > 1.0e-09) * (abs(z) > 0.15) * (0.00)} 112 set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))} 113 113 114 114 } … … 167 167 set OutputArray electrons 168 168 # tracking efficiency formula for electrons 169 set EfficiencyFormula { 169 set EfficiencyFormula { 170 170 (pt <= 0.2) * (0.00) + \ 171 171 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \ … … 199 199 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \ 200 200 (abs(eta) > 4.0) * (0.00) 201 201 202 202 } 203 203 } … … 212 212 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons 213 213 set OutputArray chargedHadrons 214 # resolution formula for charged hadrons , 215 214 # resolution formula for charged hadrons , 215 216 216 # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html 217 217 source trackMomentumResolution.tcl … … 229 229 230 230 # resolution formula for electrons 231 231 232 232 # taking something flat in energy for now, ECAL will take over at high energy anyway. 233 233 # inferred from hep-ex/1306.2016 and 1502.02701 234 set ResolutionFormula { 235 234 set ResolutionFormula { 235 236 236 (abs(eta) <= 1.5) * (1+0.64*abs(eta)^2)*(energy*0.028) + 237 237 (abs(eta) > 1.5 && abs(eta) <= 1.75) * (energy*0.037) + 238 238 (abs(eta) > 1.75 && abs(eta) <= 2.15) * (energy*0.038) + 239 239 (abs(eta) > 2.15 && abs(eta) <= 3.00) * (energy*0.044) + 240 (abs(eta) > 3.00 && abs(eta) <= 4.00) * (energy*0.10)} 240 (abs(eta) > 3.00 && abs(eta) <= 4.00) * (energy*0.10)} 241 241 242 242 } … … 250 250 set OutputArray muons 251 251 # resolution formula for muons 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 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 255 255 # http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html 256 source muonMomentumResolution.tcl 256 source muonMomentumResolution.tcl 257 257 } 258 258 … … 284 284 set EFlowTowerOutputArray eflowPhotons 285 285 286 set IsEcal true 287 286 set IsEcal true 287 288 288 set EnergyMin 0.5 289 289 set EnergySignificanceMin 1.0 … … 298 298 299 299 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5 300 300 301 301 set PhiBins {} 302 302 for {set i -180} {$i <= 180} {incr i} { … … 311 311 312 312 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL) 313 313 314 314 set PhiBins {} 315 315 for {set i -180} {$i <= 180} {incr i} { … … 328 328 } 329 329 330 # take present CMS granularity for HF 331 330 # take present CMS granularity for HF 331 332 332 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0 333 333 set PhiBins {} … … 335 335 add PhiBins [expr {$i * $pi/18.0}] 336 336 } 337 337 338 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} { 339 339 add EtaPhiBins $eta $PhiBins … … 361 361 362 362 # set ResolutionFormula {resolution formula as a function of eta and energy} 363 363 364 364 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701 365 365 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117 366 366 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) + 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) + 368 368 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \ 369 369 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \ … … 385 385 set EFlowTowerOutputArray eflowNeutralHadrons 386 386 387 set IsEcal false 388 387 set IsEcal false 388 389 389 set EnergyMin 1.0 390 390 set EnergySignificanceMin 1.0 … … 397 397 # each list starts with the lower edge of the first tower 398 398 # the list ends with the higher edged of the last tower 399 399 400 400 # assume 0.087 x 0.087 resolution in eta,phi in the barrel |eta| < 1.5 401 401 … … 409 409 410 410 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- HCAL) 411 411 412 412 set PhiBins {} 413 413 for {set i -180} {$i <= 180} {incr i} { … … 426 426 } 427 427 428 # take present CMS granularity for HF 429 428 # take present CMS granularity for HF 429 430 430 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0 431 431 set PhiBins {} … … 433 433 add PhiBins [expr {$i * $pi/18.0}] 434 434 } 435 435 436 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} { 437 437 add EtaPhiBins $eta $PhiBins … … 478 478 } 479 479 480 480 481 ########################## 481 482 # Track pile-up subtractor … … 547 548 ######################################### 548 549 549 module RunPUPPI RunPUPPI { 550 module PdgCodeFilter LeptonFilterNoLep { 551 set InputArray HCal/eflowTracks 552 set OutputArray eflowTracksNoLeptons 553 set Invert false 554 add PdgCode {13} 555 add PdgCode {-13} 556 add PdgCode {11} 557 add PdgCode {-11} 558 } 559 560 module PdgCodeFilter LeptonFilterLep { 561 set InputArray HCal/eflowTracks 562 set OutputArray eflowTracksLeptons 563 set Invert true 564 add PdgCode {11} 565 add PdgCode {-11} 566 add PdgCode {13} 567 add PdgCode {-13} 568 } 569 570 module RunPUPPI RunPUPPIBase { 550 571 ## input information 551 set TrackInputArray TrackMerger/tracks572 set TrackInputArray LeptonFilterNoLep/eflowTracksNoLeptons 552 573 set NeutralInputArray NeutralEFlowMerger/eflowTowers 553 574 set PVInputArray PileUpMerger/vertices 554 575 set MinPuppiWeight 0.05 555 576 set UseExp false 556 557 ## define puppi algorithm parameters (more than one for the same eta region is possible) 558 add EtaMinBin 0. 2.5 2.5 3.0 3.0 559 add EtaMaxBin 2.5 3.0 3.0 10.0 10.0 560 add PtMinBin 0. 0.5 0.5 0.5 0.5 561 add ConeSizeBin 0.25 0.25 0.25 0.25 0.25 562 add RMSPtMinBin 0.1 0.5 0.5 0.5 0.5 563 add RMSScaleFactorBin 1.0 1.0 1.0 1.0 1.0 564 add NeutralMinEBin 0.2 1.0 1.0 1.5 1.5 565 add NeutralPtSlope 0.02 0.02 0.02 0.02 0.02 566 add ApplyCHS true true true true true 567 add UseCharged true false false false false 568 add ApplyLowPUCorr true true true true true 569 add MetricId 5 5 1 5 1 577 set UseNoLep false 578 579 ## define puppi algorithm parameters (more than one for the same eta region is possible) 580 add EtaMinBin 0.0 1.5 4.0 581 add EtaMaxBin 1.5 4.0 10.0 582 add PtMinBin 0.0 0.0 0.0 583 add ConeSizeBin 0.2 0.2 0.2 584 add RMSPtMinBin 0.1 0.5 0.5 585 add RMSScaleFactorBin 1.0 1.0 1.0 586 add NeutralMinEBin 0.2 0.2 0.5 587 add NeutralPtSlope 0.006 0.013 0.067 588 add ApplyCHS true true true 589 add UseCharged true true false 590 add ApplyLowPUCorr true true true 591 add MetricId 5 5 5 592 add CombId 0 0 0 570 593 571 594 ## output name … … 573 596 set OutputArrayTracks puppiTracks 574 597 set OutputArrayNeutrals puppiNeutrals 575 } 576 577 598 } 599 600 module Merger RunPUPPI { 601 add InputArray RunPUPPIBase/PuppiParticles 602 add InputArray LeptonFilterLep/eflowTracksLeptons 603 set OutputArray PuppiParticles 604 } 578 605 579 606 ################### … … 588 615 } 589 616 617 module Merger PuppiMissingET { 618 #add InputArray InputArray 619 add InputArray RunPUPPI/PuppiParticles 620 #add InputArray EFlowMerger/eflow 621 set MomentumOutputArray momentum 622 } 623 590 624 ################### 591 # Ger PileUp Missing ET 625 # Ger PileUp Missing ET 592 626 ################### 593 627 … … 636 670 module FastJetFinder GenJetFinder { 637 671 set InputArray NeutrinoFilter/filteredParticles 638 672 639 673 set OutputArray jets 640 674 … … 689 723 690 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 691 738 #################### 692 739 # Photon isolation # … … 694 741 695 742 module Isolation PhotonIsolation { 696 743 697 744 # particle for which calculate the isolation 698 set CandidateInputArray ECal/eflowPhotons699 745 set CandidateInputArray PhotonFilter/photons 746 700 747 # isolation collection 701 set IsolationInputArray EFlowMerger/eflow702 748 set IsolationInputArray RunPUPPI/PuppiParticles 749 703 750 # output array 704 751 set OutputArray photons 705 752 753 # veto isolation cand. based on proximity to input cand. 754 set DeltaRMin 0.01 755 set UseMiniCone true 756 706 757 # isolation cone 707 758 set DeltaRMax 0.3 708 709 # minimum pT 759 760 # minimum pT 710 761 set PTMin 1.0 711 762 712 763 # iso ratio to cut 713 764 set PTRatioMax 9999. 765 714 766 } 715 767 … … 721 773 722 774 module Efficiency PhotonEfficiency { 723 775 724 776 ## input particles 725 set InputArray PhotonIsolation/photons 777 set InputArray PhotonIsolation/photons 726 778 ## output particles 727 779 set OutputArray photons … … 741 793 742 794 module Isolation ElectronIsolation { 743 795 744 796 set CandidateInputArray ElectronFilter/electrons 745 797 746 798 # isolation collection 747 set IsolationInputArray EFlowMerger/eflow 748 799 set IsolationInputArray RunPUPPI/PuppiParticles 800 #set IsolationInputArray EFlowMerger/eflow 801 749 802 set OutputArray electrons 750 803 751 804 set DeltaRMax 0.3 752 805 set PTMin 1.0 … … 762 815 763 816 module Efficiency ElectronEfficiency { 764 817 765 818 set InputArray ElectronIsolation/electrons 766 819 set OutputArray electrons 767 820 768 821 # set EfficiencyFormula {efficiency formula as a function of eta and pt} 769 822 # efficiency formula for electrons … … 773 826 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \ 774 827 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \ 775 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \ 776 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \ 777 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \ 778 (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) + \ 779 832 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \ 780 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \ 781 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \ 782 (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) + \ 783 836 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \ 784 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \ 785 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \ 786 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \ 787 (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) + \ 788 841 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \ 789 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \ 790 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \ 791 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \ 792 (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) + \ 793 846 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \ 794 847 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \ … … 797 850 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \ 798 851 (abs(eta) > 4.0) * (0.00) 799 852 800 853 } 801 854 } … … 807 860 module Isolation MuonIsolation { 808 861 set CandidateInputArray MuonMomentumSmearing/muons 809 862 810 863 # isolation collection 811 set IsolationInputArray EFlowMerger/eflow812 864 set IsolationInputArray RunPUPPI/PuppiParticles 865 813 866 set OutputArray muons 814 867 815 868 set DeltaRMax 0.3 816 869 set PTMin 1.0 … … 885 938 886 939 set BitNumber 1 887 940 888 941 source btagMedium.tcl 889 942 } 890 943 891 892 944 893 945 ############# … … 900 952 set BitNumber 2 901 953 902 source btagTight.tcl 903 } 904 905 906 954 source btagTight.tcl 955 } 907 956 908 957 … … 930 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) ) + \ 931 980 (abs(eta) > 2.3) * (0.000) 932 } 981 } 933 982 } 934 983 … … 939 988 940 989 module StatusPidFilter GenParticleFilter { 941 990 942 991 set InputArray Delphes/allParticles 943 992 set OutputArray filteredParticles 944 993 set PTMin 5.0 945 994 946 995 } 947 996 … … 955 1004 add Branch GenParticleFilter/filteredParticles Particle GenParticle 956 1005 add Branch PileUpMerger/vertices Vertex Vertex 957 1006 958 1007 add Branch GenJetFinder/jets GenJet Jet 959 1008 add Branch GenMissingET/momentum GenMissingET MissingET … … 962 1011 # add Branch ECal/eflowPhotons EFlowPhoton Tower 963 1012 # add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower 964 1013 965 1014 add Branch PhotonEfficiency/photons Photon Photon 966 1015 add Branch ElectronEfficiency/electrons Electron Electron 967 1016 add Branch MuonLooseIdEfficiency/muons MuonLoose Muon 968 1017 add Branch MuonTightIdEfficiency/muons MuonTight Muon 969 1018 970 1019 add Branch JetEnergyScale/jets Jet Jet 971 1020 972 1021 add Branch MissingET/momentum MissingET MissingET 1022 add Branch PuppiMissingET/momentum PuppiMissingET MissingET 973 1023 add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET 974 1024 add Branch ScalarHT/energy ScalarHT ScalarHT
Note:
See TracChangeset
for help on using the changeset viewer.