- Timestamp:
- Aug 22, 2016, 7:27:26 PM (8 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 1270747
- Parents:
- 8f4a953
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
cards/CMS_PhaseII/CMS_PhaseII_Substructure_PIX4022_0PU.tcl
r8f4a953 r8cbe720 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 … … 121 139 122 140 # maximum spread in the beam direction in m 123 set ZVertexSpread 0 141 set ZVertexSpread 0.25 124 142 125 143 # maximum spread in time in s 126 set TVertexSpread 0144 set TVertexSpread 800E-12 127 145 128 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 set VertexDistributionFormula 1147 set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))} 130 148 131 149 } … … 273 291 source muonMomentumResolution.tcl 274 292 } 275 276 277 278 293 279 294 ############## … … 527 542 } 528 543 544 ######################## 545 # Reco PU filter 546 ######################## 547 548 module RecoPuFilter RecoPuFilter { 549 set InputArray HCal/eflowTracks 550 set OutputArray eflowTracks 551 } 529 552 530 553 ################################################### … … 539 562 } 540 563 541 542 564 #################### 543 565 # Neutral eflow erger … … 551 573 } 552 574 553 554 #################### 575 ##################### 555 576 # Energy flow merger 556 #################### 577 ##################### 557 578 558 579 module Merger EFlowMerger { … … 564 585 } 565 586 566 ############################ ######567 # Energy flow merger (all tracks)568 ############################ ######569 570 module Merger EFlowMerger AllTracks{587 ############################ 588 # Energy flow merger no PU 589 ############################ 590 591 module Merger EFlowMergerCHS { 571 592 # add InputArray InputArray 572 add InputArray TrackMerger/tracks593 add InputArray RecoPuFilter/eflowTracks 573 594 add InputArray PhotonEnergySmearing/eflowPhotons 574 595 add InputArray HCal/eflowNeutralHadrons … … 695 716 } 696 717 697 698 718 ##################### 699 719 # MC truth jet finder … … 735 755 } 736 756 737 ############ 757 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 ############## 738 788 # Jet finder 739 ############ 789 ############## 740 790 741 791 module FastJetFinder FastJetFinder { 742 792 # set InputArray TowerMerger/towers 743 set InputArray EFlowMerger /eflow793 set InputArray EFlowMergerCHS/eflow 744 794 745 795 set OutputArray jets 796 797 set AreaAlgorithm 5 746 798 747 799 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt … … 755 807 module FastJetFinder FastJetFinderAK8 { 756 808 # set InputArray TowerMerger/towers 757 set InputArray EFlowMerger /eflow809 set InputArray EFlowMergerCHS/eflow 758 810 759 811 set OutputArray jets 812 813 set AreaAlgorithm 5 760 814 761 815 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt … … 784 838 } 785 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 786 866 module FastJetFinder FastJetFinderPUPPI { 787 867 # set InputArray TowerMerger/towers … … 833 913 834 914 module EnergyScale JetEnergyScale { 835 set InputArray FastJetFinder/jets915 set InputArray JetPileUpSubtractor/jets 836 916 set OutputArray jets 837 917 … … 841 921 842 922 module EnergyScale JetEnergyScaleAK8 { 843 set InputArray FastJetFinderAK8/jets923 set InputArray JetPileUpSubtractorAK8/jets 844 924 set OutputArray jets 845 925 … … 907 987 } 908 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 } 909 1015 910 1016 … … 929 1035 930 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 931 1056 ###################### 932 1057 # Electron isolation # … … 949 1074 } 950 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 } 951 1095 952 1096 … … 995 1139 } 996 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 997 1186 ################## 998 1187 # Muon isolation # … … 1013 1202 } 1014 1203 1015 1016 ################## 1017 # Muon Loose Id # 1018 ################## 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 ##################### 1019 1226 1020 1227 module Efficiency MuonLooseIdEfficiency { … … 1038 1245 } 1039 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 } 1040 1271 1041 1272 … … 1250 1481 module StatusPidFilter GenParticleFilter { 1251 1482 1252 set InputArray 1483 set InputArray Delphes/allParticles 1253 1484 set OutputArray filteredParticles 1254 1485 set PTMin 5.0 … … 1278 1509 add Branch MuonLooseIdEfficiency/muons MuonLoose Muon 1279 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 1280 1516 1281 1517 add Branch JetEnergyScale/jets Jet Jet … … 1289 1525 add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET 1290 1526 add Branch ScalarHT/energy ScalarHT ScalarHT 1291 } 1527 1528 }
Note:
See TracChangeset
for help on using the changeset viewer.