Changeset 8f4a953 in git
- Timestamp:
- Aug 22, 2016, 7:27:03 PM (8 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 8cbe720
- Parents:
- 308ab1a
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
cards/CMS_PhaseII/CMS_PhaseII_Substructure_PIX4022_200PU.tcl
r308ab1a r8f4a953 1 set MaxEvents 100 2 1 3 # 2 4 # Phase II - No Pile-Up … … 33 35 PhotonEnergySmearing 34 36 ElectronFilter 37 35 38 TrackPileUpSubtractor 39 RecoPuFilter 36 40 37 41 TowerMerger 38 42 NeutralEFlowMerger 39 EFlowMergerAllTracks 43 40 44 EFlowMerger 45 EFlowMergerCHS 46 Rho 41 47 42 48 LeptonFilterNoLep … … 48 54 49 55 PhotonIsolation 56 PhotonIsolationCHS 50 57 PhotonEfficiency 58 PhotonEfficiencyCHS 51 59 52 60 ElectronIsolation 61 ElectronIsolationCHS 62 53 63 ElectronEfficiency 64 ElectronEfficiencyCHS 54 65 55 66 MuonIsolation 67 MuonIsolationCHS 68 56 69 MuonLooseIdEfficiency 57 70 MuonTightIdEfficiency 71 72 MuonLooseIdEfficiencyCHS 73 MuonTightIdEfficiencyCHS 58 74 59 75 NeutrinoFilter … … 68 84 FastJetFinder 69 85 FastJetFinderAK8 86 JetPileUpSubtractor 87 JetPileUpSubtractorAK8 70 88 FastJetFinderPUPPI 71 89 FastJetFinderPUPPIAK8 … … 114 132 115 133 # pre-generated minbias input file 116 set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup117 #set PileUpFile MinBias.pileup134 #set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup 135 set PileUpFile MinBias.pileup 118 136 119 137 # average expected pile up … … 126 144 set TVertexSpread 800E-12 127 145 128 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s) 146 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s) - {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))} 129 147 set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))} 130 148 … … 138 156 module ParticlePropagator ParticlePropagator { 139 157 set InputArray PileUpMerger/stableParticles 158 #set InputArray Delphes/stableParticles 140 159 141 160 set OutputArray stableParticles … … 272 291 source muonMomentumResolution.tcl 273 292 } 274 275 276 277 293 278 294 ############## … … 526 542 } 527 543 544 ######################## 545 # Reco PU filter 546 ######################## 547 548 module RecoPuFilter RecoPuFilter { 549 set InputArray HCal/eflowTracks 550 set OutputArray eflowTracks 551 } 528 552 529 553 ################################################### … … 538 562 } 539 563 540 541 564 #################### 542 565 # Neutral eflow erger … … 550 573 } 551 574 552 553 #################### 575 ##################### 554 576 # Energy flow merger 555 #################### 577 ##################### 556 578 557 579 module Merger EFlowMerger { … … 563 585 } 564 586 565 ############################ ######566 # Energy flow merger (all tracks)567 ############################ ######568 569 module Merger EFlowMerger AllTracks{587 ############################ 588 # Energy flow merger no PU 589 ############################ 590 591 module Merger EFlowMergerCHS { 570 592 # add InputArray InputArray 571 add InputArray TrackMerger/tracks593 add InputArray RecoPuFilter/eflowTracks 572 594 add InputArray PhotonEnergySmearing/eflowPhotons 573 595 add InputArray HCal/eflowNeutralHadrons … … 694 716 } 695 717 696 697 718 ##################### 698 719 # MC truth jet finder … … 735 756 736 757 737 738 ############ 758 ############# 759 # Rho pile-up 760 ############# 761 762 module FastJetFinder Rho { 763 # set InputArray Calorimeter/towers 764 set InputArray EFlowMergerCHS/eflow 765 766 set ComputeRho true 767 set RhoOutputArray rho 768 769 # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area 770 set AreaAlgorithm 5 771 772 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt 773 set JetAlgorithm 4 774 set ParameterR 0.4 775 set GhostEtaMax 5.0 776 777 add RhoEtaRange -5.0 -4.0 778 add RhoEtaRange -4.0 -1.5 779 add RhoEtaRange -1.5 1.5 780 add RhoEtaRange 1.5 4.0 781 add RhoEtaRange 4.0 5.0 782 783 set JetPTMin 0.0 784 } 785 786 787 ############## 739 788 # Jet finder 740 ############ 789 ############## 741 790 742 791 module FastJetFinder FastJetFinder { 743 792 # set InputArray TowerMerger/towers 744 set InputArray EFlowMerger /eflow793 set InputArray EFlowMergerCHS/eflow 745 794 746 795 set OutputArray jets 796 797 set AreaAlgorithm 5 747 798 748 799 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt … … 756 807 module FastJetFinder FastJetFinderAK8 { 757 808 # set InputArray TowerMerger/towers 758 set InputArray EFlowMerger /eflow809 set InputArray EFlowMergerCHS/eflow 759 810 760 811 set OutputArray jets 812 813 set AreaAlgorithm 5 761 814 762 815 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt … … 785 838 } 786 839 840 ########################### 841 # Jet Pile-Up Subtraction 842 ########################### 843 844 module JetPileUpSubtractor JetPileUpSubtractor { 845 set JetInputArray FastJetFinder/jets 846 set RhoInputArray Rho/rho 847 848 set OutputArray jets 849 850 set JetPTMin 15.0 851 } 852 853 ############################## 854 # Jet Pile-Up Subtraction AK8 855 ############################## 856 857 module JetPileUpSubtractor JetPileUpSubtractorAK8 { 858 set JetInputArray FastJetFinderAK8/jets 859 set RhoInputArray Rho/rho 860 861 set OutputArray jets 862 863 set JetPTMin 15.0 864 } 865 787 866 module FastJetFinder FastJetFinderPUPPI { 788 867 # set InputArray TowerMerger/towers … … 834 913 835 914 module EnergyScale JetEnergyScale { 836 set InputArray FastJetFinder/jets915 set InputArray JetPileUpSubtractor/jets 837 916 set OutputArray jets 838 917 … … 842 921 843 922 module EnergyScale JetEnergyScaleAK8 { 844 set InputArray FastJetFinderAK8/jets923 set InputArray JetPileUpSubtractorAK8/jets 845 924 set OutputArray jets 846 925 … … 908 987 } 909 988 989 990 ######################## 991 # Photon isolation CHS # 992 ######################## 993 994 module Isolation PhotonIsolationCHS { 995 996 # particle for which calculate the isolation 997 set CandidateInputArray PhotonFilter/photons 998 999 # isolation collection 1000 set IsolationInputArray EFlowMerger/eflow 1001 1002 # output array 1003 set OutputArray photons 1004 1005 # isolation cone 1006 set DeltaRMax 0.3 1007 1008 # minimum pT 1009 set PTMin 1.0 1010 1011 # iso ratio to cut 1012 set PTRatioMax 9999. 1013 1014 } 910 1015 911 1016 … … 930 1035 931 1036 1037 ##################### 1038 # Photon efficiency # 1039 ##################### 1040 1041 module Efficiency PhotonEfficiencyCHS { 1042 1043 ## input particles 1044 set InputArray PhotonIsolationCHS/photons 1045 ## output particles 1046 set OutputArray photons 1047 # set EfficiencyFormula {efficiency formula as a function of eta and pt} 1048 # efficiency formula for photons 1049 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \ 1050 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \ 1051 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \ 1052 (abs(eta) > 4.0) * (0.00)} 1053 1054 } 1055 932 1056 ###################### 933 1057 # Electron isolation # … … 950 1074 } 951 1075 1076 1077 ########################## 1078 # Electron isolation CHS # 1079 ########################## 1080 1081 module Isolation ElectronIsolationCHS { 1082 1083 set CandidateInputArray ElectronFilter/electrons 1084 1085 # isolation collection 1086 set IsolationInputArray EFlowMerger/eflow 1087 1088 set OutputArray electrons 1089 1090 set DeltaRMax 0.3 1091 set PTMin 1.0 1092 set PTRatioMax 9999. 1093 1094 } 952 1095 953 1096 … … 996 1139 } 997 1140 1141 ########################### 1142 # Electron efficiency CHS # 1143 ########################### 1144 1145 module Efficiency ElectronEfficiencyCHS { 1146 1147 set InputArray ElectronIsolationCHS/electrons 1148 set OutputArray electrons 1149 1150 # set EfficiencyFormula {efficiency formula as a function of eta and pt} 1151 # efficiency formula for electrons 1152 set EfficiencyFormula { 1153 (pt <= 4.0) * (0.00) + \ 1154 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \ 1155 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \ 1156 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \ 1157 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \ 1158 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \ 1159 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \ 1160 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \ 1161 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \ 1162 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \ 1163 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \ 1164 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \ 1165 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \ 1166 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \ 1167 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \ 1168 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \ 1169 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \ 1170 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \ 1171 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \ 1172 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \ 1173 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \ 1174 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \ 1175 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \ 1176 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \ 1177 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \ 1178 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \ 1179 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \ 1180 (abs(eta) > 4.0) * (0.00) 1181 1182 } 1183 } 1184 1185 998 1186 ################## 999 1187 # Muon isolation # … … 1014 1202 } 1015 1203 1016 1017 ################## 1018 # Muon Loose Id # 1019 ################## 1204 ###################### 1205 # Muon isolation CHS # 1206 ###################### 1207 1208 module Isolation MuonIsolationCHS { 1209 set CandidateInputArray MuonMomentumSmearing/muons 1210 1211 # isolation collection 1212 set IsolationInputArray EFlowMerger/eflow 1213 1214 set OutputArray muons 1215 1216 set DeltaRMax 0.3 1217 set PTMin 1.0 1218 set PTRatioMax 9999. 1219 1220 } 1221 1222 1223 ##################### 1224 # Muon Loose Id # 1225 ##################### 1020 1226 1021 1227 module Efficiency MuonLooseIdEfficiency { … … 1039 1245 } 1040 1246 1247 1248 ##################### 1249 # Muon Loose Id CHS # 1250 ##################### 1251 1252 module Efficiency MuonLooseIdEfficiencyCHS { 1253 set InputArray MuonIsolationCHS/muons 1254 set OutputArray muons 1255 # tracking + LooseID efficiency formula for muons 1256 source muonLooseId.tcl 1257 1258 } 1259 1260 1261 ###################### 1262 # Muon Tight Id CHS # 1263 ###################### 1264 1265 module Efficiency MuonTightIdEfficiencyCHS { 1266 set InputArray MuonIsolationCHS/muons 1267 set OutputArray muons 1268 # tracking + TightID efficiency formula for muons 1269 source muonTightId.tcl 1270 } 1041 1271 1042 1272 … … 1251 1481 module StatusPidFilter GenParticleFilter { 1252 1482 1253 set InputArray 1483 set InputArray Delphes/allParticles 1254 1484 set OutputArray filteredParticles 1255 1485 set PTMin 5.0 … … 1279 1509 add Branch MuonLooseIdEfficiency/muons MuonLoose Muon 1280 1510 add Branch MuonTightIdEfficiency/muons MuonTight Muon 1511 1512 add Branch PhotonEfficiencyCHS/photons PhotonCHS Photon 1513 add Branch ElectronEfficiencyCHS/electrons ElectronCHS Electron 1514 add Branch MuonLooseIdEfficiencyCHS/muons MuonLooseCHS Muon 1515 add Branch MuonTightIdEfficiencyCHS/muons MuonTightCHS Muon 1281 1516 1282 1517 add Branch JetEnergyScale/jets Jet Jet … … 1290 1525 add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET 1291 1526 add Branch ScalarHT/energy ScalarHT ScalarHT 1292 } 1527 1528 } -
modules/ModulesLinkDef.h
r308ab1a r8f4a953 58 58 #include "modules/StatusPidFilter.h" 59 59 #include "modules/PdgCodeFilter.h" 60 #include "modules/RecoPuFilter.h" 60 61 #include "modules/Cloner.h" 61 62 #include "modules/Weighter.h" … … 103 104 #pragma link C++ class StatusPidFilter+; 104 105 #pragma link C++ class PdgCodeFilter+; 106 #pragma link C++ class RecoPuFilter+; 105 107 #pragma link C++ class Cloner+; 106 108 #pragma link C++ class Weighter+;
Note:
See TracChangeset
for help on using the changeset viewer.