Fork me on GitHub

source: git/cards/delphes_card_HLLHC.tcl@ 97f1b2d

ImprovedOutputFile Timing dual_readout llp 3.4.2pre09
Last change on this file since 97f1b2d was 97f1b2d, checked in by Michele Selvaggi <michele.selvaggi@…>, 7 years ago

fix muon smearing
  • Property mode set to 100644
File size: 58.7 KB
Line 
1#
2# Beta card for HL-LHC and HE-LHC studies
3#
4# Main authors: Michele Selvaggi (CERN)
5#
6# Released on: Dec. 6th, 2017
7#
8#
9#######################################
10# Order of execution of various modules
11#######################################
12
13set ExecutionPath {
14
15 ParticlePropagator
16
17 ChargedHadronTrackingEfficiency
18 ElectronTrackingEfficiency
19 MuonTrackingEfficiency
20
21 ChargedHadronMomentumSmearing
22 ElectronMomentumSmearing
23 MuonMomentumSmearing
24
25 DenseProp
26 DenseMergeTracks
27 TrackMergerProp
28 TrackMerger
29
30 ECal
31 HCal
32
33 Calorimeter
34 EFlowMerger
35 EFlowFilter
36
37 PhotonIsolation
38 PhotonEfficiency
39
40 ElectronFilter
41 ElectronIsolation
42 ElectronEfficiency
43
44 ChargedHadronFilter
45
46 MuonIsolation
47 MuonEfficiency
48
49 MissingET
50
51 NeutrinoFilter
52 GenJetFinder
53 GenMissingET
54
55 FastJetFinder
56 FatJetFinder
57
58 JetEnergyScale
59
60 JetFlavorAssociation
61
62 BTagging
63 TauTagging
64
65 UniqueObjectFinder
66
67 ScalarHT
68
69 TreeWriter
70}
71
72
73#################################
74# Propagate particles in cylinder
75#################################
76
77module ParticlePropagator ParticlePropagator {
78 set InputArray Delphes/stableParticles
79
80 set OutputArray stableParticles
81 set ChargedHadronOutputArray chargedHadrons
82 set ElectronOutputArray electrons
83 set MuonOutputArray muons
84
85 # radius of the magnetic field coverage, in m
86 set Radius 1.2
87 # half-length of the magnetic field coverage, in m
88 set HalfLength 3.25
89
90 # magnetic field
91 set Bz 3.0
92}
93
94####################################
95# Charged hadron tracking efficiency
96####################################
97
98module Efficiency ChargedHadronTrackingEfficiency {
99 set InputArray ParticlePropagator/chargedHadrons
100 set OutputArray chargedHadrons
101
102 # TBC (which eta_max ? which pT min?)
103
104 # tracking efficiency formula for charged hadrons
105
106 set EfficiencyFormula {
107 (pt <= 0.5) * (0.00) + \
108 (abs(eta) <= 1.2) * (pt > 0.5 && pt <= 1.0) * (pt * 0.90) + \
109 (abs(eta) <= 1.2) * (pt > 1.0) * (0.95) + \
110 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1.0) * (pt*0.85) + \
111 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.90) + \
112 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.5 && pt <= 1.0) * (pt*0.80) + \
113 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.85) + \
114 (abs(eta) > 4.0) * (0.00)
115 }
116
117}
118
119##############################
120# Electron tracking efficiency
121##############################
122
123module Efficiency ElectronTrackingEfficiency {
124 set InputArray ParticlePropagator/electrons
125 set OutputArray electrons
126
127# TBC (which eta_max ?)
128# putting same as charged hadrons for now...
129
130 set EfficiencyFormula {
131 (pt <= 0.5) * (0.00) + \
132 (abs(eta) <= 1.2) * (pt > 0.5 && pt <= 1.0) * (pt * 0.90) + \
133 (abs(eta) <= 1.2) * (pt > 1.0) * (0.95) + \
134 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1.0) * (pt*0.85) + \
135 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.90)/(-0.033*log10(pt) + 1.033) + \
136 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.5 && pt <= 1.0) * (pt*0.80) + \
137 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.85)/(-0.033*log10(pt) + 1.033) + \
138 (abs(eta) > 4.0) * (0.00)
139 }
140
141}
142##########################
143# Muon tracking efficiency
144##########################
145
146module Efficiency MuonTrackingEfficiency {
147 set InputArray ParticlePropagator/muons
148 set OutputArray muons
149
150 set EfficiencyFormula {
151 (pt <= 0.5) * (0.00) + \
152 (abs(eta) <= 1.2) * (pt > 0.5 && pt <= 1.0) * (pt * 1.00) + \
153 (abs(eta) <= 1.2) * (pt > 1.0) * (1.00) + \
154 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 0.5 && pt <= 1.0) * (pt*1.00) + \
155 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 1.0) * (1.00) + \
156 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 0.5 && pt <= 1.0) * (pt*0.95) + \
157 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
158 (abs(eta) > 4.0) * (0.00)
159
160 }
161
162}
163
164########################################
165# Momentum resolution for charged tracks
166########################################
167
168module MomentumSmearing ChargedHadronMomentumSmearing {
169 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
170 set OutputArray chargedHadrons
171
172 # --- CMS resolution Phase II ---
173
174 set ResolutionFormula { 2*((abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00457888) + \
175 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004579 + (pt-1.000000)* 0.000045) + \
176 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.004983 + (pt-10.000000)* 0.000047) + \
177 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.009244*pt/100.000000) + \
178 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00505011) + \
179 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005050 + (pt-1.000000)* 0.000033) + \
180 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005343 + (pt-10.000000)* 0.000043) + \
181 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.009172*pt/100.000000) + \
182 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00510573) + \
183 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005106 + (pt-1.000000)* 0.000023) + \
184 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005317 + (pt-10.000000)* 0.000042) + \
185 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.009077*pt/100.000000) + \
186 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00578020) + \
187 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005780 + (pt-1.000000)* -0.000000) + \
188 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.005779 + (pt-10.000000)* 0.000038) + \
189 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.009177*pt/100.000000) + \
190 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00728723) + \
191 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.007287 + (pt-1.000000)* -0.000031) + \
192 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007011 + (pt-10.000000)* 0.000038) + \
193 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.010429*pt/100.000000) + \
194 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01045117) + \
195 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010451 + (pt-1.000000)* -0.000051) + \
196 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.009989 + (pt-10.000000)* 0.000043) + \
197 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.013867*pt/100.000000) + \
198 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01477199) + \
199 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014772 + (pt-1.000000)* -0.000128) + \
200 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.013616 + (pt-10.000000)* 0.000035) + \
201 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.016800*pt/100.000000) + \
202 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01731474) + \
203 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017315 + (pt-1.000000)* -0.000208) + \
204 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.015439 + (pt-10.000000)* 0.000030) + \
205 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.018161*pt/100.000000) + \
206 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01942025) + \
207 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.019420 + (pt-1.000000)* -0.000417) + \
208 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.015669 + (pt-10.000000)* 0.000026) + \
209 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.018039*pt/100.000000) + \
210 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.02201432) + \
211 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.022014 + (pt-1.000000)* -0.000667) + \
212 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.016012 + (pt-10.000000)* 0.000045) + \
213 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.020098*pt/100.000000) + \
214 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.02574300) + \
215 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.025743 + (pt-1.000000)* -0.001118) + \
216 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.015681 + (pt-10.000000)* 0.000051) + \
217 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.020289*pt/100.000000) + \
218 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.02885821) + \
219 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.028858 + (pt-1.000000)* -0.001345) + \
220 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.016753 + (pt-10.000000)* 0.000053) + \
221 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.021524*pt/100.000000) + \
222 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.03204812) + \
223 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.032048 + (pt-1.000000)* -0.001212) + \
224 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.021138 + (pt-10.000000)* 0.000037) + \
225 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.024477*pt/100.000000) + \
226 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.03950405) + \
227 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.039504 + (pt-1.000000)* -0.001386) + \
228 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.027026 + (pt-10.000000)* 0.000037) + \
229 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.030392*pt/100.000000) + \
230 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.04084751) + \
231 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.040848 + (pt-1.000000)* -0.001780) + \
232 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.024824 + (pt-10.000000)* 0.000029) + \
233 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.027445*pt/100.000000) + \
234 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04532425) + \
235 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045324 + (pt-1.000000)* -0.002497) + \
236 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.022851 + (pt-10.000000)* 0.000024) + \
237 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.025053*pt/100.000000) + \
238 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.06418925) + \
239 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.064189 + (pt-1.000000)* -0.004055) + \
240 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.027691 + (pt-10.000000)* 0.000034) + \
241 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.030710*pt/100.000000) + \
242 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.07682500) + \
243 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.076825 + (pt-1.000000)* -0.004510) + \
244 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.036234 + (pt-10.000000)* 0.000049) + \
245 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.040629*pt/100.000000) + \
246 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.09796358) + \
247 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.097964 + (pt-1.000000)* -0.005758) + \
248 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.046145 + (pt-10.000000)* 0.000069) + \
249 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.052345*pt/100.000000) + \
250 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.13415929) + \
251 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.134159 + (pt-1.000000)* -0.008283) + \
252 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.059612 + (pt-10.000000)* 0.000111) + \
253 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.069617*pt/100.000000))
254 }
255
256}
257
258
259###################################
260# Momentum resolution for electrons
261###################################
262
263module MomentumSmearing ElectronMomentumSmearing {
264 set InputArray ElectronTrackingEfficiency/electrons
265 set OutputArray electrons
266
267 set ResolutionFormula { 2*((abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00457888) + \
268 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004579 + (pt-1.000000)* 0.000045) + \
269 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.004983 + (pt-10.000000)* 0.000047) + \
270 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.009244*pt/100.000000) + \
271 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00505011) + \
272 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005050 + (pt-1.000000)* 0.000033) + \
273 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005343 + (pt-10.000000)* 0.000043) + \
274 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.009172*pt/100.000000) + \
275 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00510573) + \
276 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005106 + (pt-1.000000)* 0.000023) + \
277 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005317 + (pt-10.000000)* 0.000042) + \
278 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.009077*pt/100.000000) + \
279 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00578020) + \
280 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005780 + (pt-1.000000)* -0.000000) + \
281 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.005779 + (pt-10.000000)* 0.000038) + \
282 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.009177*pt/100.000000) + \
283 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00728723) + \
284 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.007287 + (pt-1.000000)* -0.000031) + \
285 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007011 + (pt-10.000000)* 0.000038) + \
286 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.010429*pt/100.000000) + \
287 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01045117) + \
288 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010451 + (pt-1.000000)* -0.000051) + \
289 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.009989 + (pt-10.000000)* 0.000043) + \
290 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.013867*pt/100.000000) + \
291 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01477199) + \
292 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014772 + (pt-1.000000)* -0.000128) + \
293 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.013616 + (pt-10.000000)* 0.000035) + \
294 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.016800*pt/100.000000) + \
295 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01731474) + \
296 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017315 + (pt-1.000000)* -0.000208) + \
297 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.015439 + (pt-10.000000)* 0.000030) + \
298 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.018161*pt/100.000000) + \
299 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01942025) + \
300 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.019420 + (pt-1.000000)* -0.000417) + \
301 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.015669 + (pt-10.000000)* 0.000026) + \
302 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.018039*pt/100.000000) + \
303 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.02201432) + \
304 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.022014 + (pt-1.000000)* -0.000667) + \
305 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.016012 + (pt-10.000000)* 0.000045) + \
306 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.020098*pt/100.000000) + \
307 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.02574300) + \
308 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.025743 + (pt-1.000000)* -0.001118) + \
309 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.015681 + (pt-10.000000)* 0.000051) + \
310 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.020289*pt/100.000000) + \
311 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.02885821) + \
312 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.028858 + (pt-1.000000)* -0.001345) + \
313 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.016753 + (pt-10.000000)* 0.000053) + \
314 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.021524*pt/100.000000) + \
315 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.03204812) + \
316 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.032048 + (pt-1.000000)* -0.001212) + \
317 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.021138 + (pt-10.000000)* 0.000037) + \
318 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.024477*pt/100.000000) + \
319 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.03950405) + \
320 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.039504 + (pt-1.000000)* -0.001386) + \
321 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.027026 + (pt-10.000000)* 0.000037) + \
322 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.030392*pt/100.000000) + \
323 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.04084751) + \
324 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.040848 + (pt-1.000000)* -0.001780) + \
325 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.024824 + (pt-10.000000)* 0.000029) + \
326 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.027445*pt/100.000000) + \
327 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04532425) + \
328 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045324 + (pt-1.000000)* -0.002497) + \
329 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.022851 + (pt-10.000000)* 0.000024) + \
330 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.025053*pt/100.000000) + \
331 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.06418925) + \
332 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.064189 + (pt-1.000000)* -0.004055) + \
333 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.027691 + (pt-10.000000)* 0.000034) + \
334 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.030710*pt/100.000000) + \
335 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.07682500) + \
336 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.076825 + (pt-1.000000)* -0.004510) + \
337 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.036234 + (pt-10.000000)* 0.000049) + \
338 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.040629*pt/100.000000) + \
339 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.09796358) + \
340 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.097964 + (pt-1.000000)* -0.005758) + \
341 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.046145 + (pt-10.000000)* 0.000069) + \
342 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.052345*pt/100.000000) + \
343 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.13415929) + \
344 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.134159 + (pt-1.000000)* -0.008283) + \
345 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.059612 + (pt-10.000000)* 0.000111) + \
346 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.069617*pt/100.000000))
347 }
348}
349
350
351###############################
352# Momentum resolution for muons
353###############################
354
355module MomentumSmearing MuonMomentumSmearing {
356 set InputArray MuonTrackingEfficiency/muons
357 set OutputArray muons
358
359 # Equal to Tk resolution up to pT = 100 GeV, then tk slope/5.
360
361 set ResolutionFormula { ((abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00457888) + \
362 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004579 + (pt-1.000000)* 0.000045) + \
363 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.004983 + (pt-10.000000)* 0.000047) + \
364 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.009244/5.*pt/100.000000 + 4./5.*0.009244) + \
365 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00505011) + \
366 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005050 + (pt-1.000000)* 0.000033) + \
367 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005343 + (pt-10.000000)* 0.000043) + \
368 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.009172/5.*pt/100.000000 + 4./5.*0.009172) + \
369 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00510573) + \
370 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005106 + (pt-1.000000)* 0.000023) + \
371 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005317 + (pt-10.000000)* 0.000042) + \
372 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.009077/5.*pt/100.000000 + 4./5.*0.009077) + \
373 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00578020) + \
374 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005780 + (pt-1.000000)* -0.000000) + \
375 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.005779 + (pt-10.000000)* 0.000038) + \
376 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.009177/5.*pt/100.000000 + 4./5.*0.009177) + \
377 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00728723) + \
378 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.007287 + (pt-1.000000)* -0.000031) + \
379 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007011 + (pt-10.000000)* 0.000038) + \
380 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.010429/5.*pt/100.000000 + 4./5.*0.010429) + \
381 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01045117) + \
382 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010451 + (pt-1.000000)* -0.000051) + \
383 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.009989 + (pt-10.000000)* 0.000043) + \
384 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.013867/5.*pt/100.000000 + 4./5.*0.013867) + \
385 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01477199) + \
386 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014772 + (pt-1.000000)* -0.000128) + \
387 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.013616 + (pt-10.000000)* 0.000035) + \
388 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.016800/5.*pt/100.000000 + 4./5.*0.016800) + \
389 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01731474) + \
390 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017315 + (pt-1.000000)* -0.000208) + \
391 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.015439 + (pt-10.000000)* 0.000030) + \
392 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.018161/5.*pt/100.000000 + 4./5.*0.018161) + \
393 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01942025) + \
394 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.019420 + (pt-1.000000)* -0.000417) + \
395 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.015669 + (pt-10.000000)* 0.000026) + \
396 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.018039/5.*pt/100.000000 + 4./5.*0.018039) + \
397 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.02201432) + \
398 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.022014 + (pt-1.000000)* -0.000667) + \
399 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.016012 + (pt-10.000000)* 0.000045) + \
400 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.020098/5.*pt/100.000000 + 4./5.*0.020098) + \
401 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.02574300) + \
402 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.025743 + (pt-1.000000)* -0.001118) + \
403 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.015681 + (pt-10.000000)* 0.000051) + \
404 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.020289/5.*pt/100.000000 + 4./5.*0.020289) + \
405 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.02885821) + \
406 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.028858 + (pt-1.000000)* -0.001345) + \
407 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.016753 + (pt-10.000000)* 0.000053) + \
408 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.021524/5.*pt/100.000000 + 4./5.*0.021524) + \
409 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.03204812) + \
410 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.032048 + (pt-1.000000)* -0.001212) + \
411 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.021138 + (pt-10.000000)* 0.000037) + \
412 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.024477/5.*pt/100.000000 + 4./5.*0.024477) + \
413 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.03950405) + \
414 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.039504 + (pt-1.000000)* -0.001386) + \
415 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.027026 + (pt-10.000000)* 0.000037) + \
416 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.030392/5.*pt/100.000000 + 4./5.*0.030392) + \
417 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.04084751) + \
418 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.040848 + (pt-1.000000)* -0.001780) + \
419 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.024824 + (pt-10.000000)* 0.000029) + \
420 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.027445/5.*pt/100.000000 + 4./5.*0.027445) + \
421 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04532425) + \
422 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045324 + (pt-1.000000)* -0.002497) + \
423 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.022851 + (pt-10.000000)* 0.000024) + \
424 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.025053/5.*pt/100.000000 + 4./5.*0.025053) + \
425 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.06418925) + \
426 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.064189 + (pt-1.000000)* -0.004055) + \
427 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.027691 + (pt-10.000000)* 0.000034) + \
428 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.030710/5.*pt/100.000000 + 4./5.*0.030710) + \
429 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.07682500) + \
430 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.076825 + (pt-1.000000)* -0.004510) + \
431 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.036234 + (pt-10.000000)* 0.000049) + \
432 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.040629/5.*pt/100.000000 + 4./5.*0.040629) + \
433 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.09796358) + \
434 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.097964 + (pt-1.000000)* -0.005758) + \
435 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.046145 + (pt-10.000000)* 0.000069) + \
436 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.052345/5.*pt/100.000000 + 4./5.*0.052345) + \
437 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.13415929) + \
438 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.134159 + (pt-1.000000)* -0.008283) + \
439 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.059612 + (pt-10.000000)* 0.000111) + \
440 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.069617/5.*pt/100.000000 + 4./5.*0.069617))
441 }
442
443}
444
445##############
446# Track merger
447##############
448
449module Merger TrackMergerProp {
450# add InputArray InputArray
451 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
452 add InputArray ElectronMomentumSmearing/electrons
453 add InputArray MuonMomentumSmearing/muons
454 set OutputArray tracks
455}
456
457
458#################################
459# Dense Track propagation
460#################################
461
462module ParticlePropagator DenseProp {
463 set InputArray Delphes/stableParticles
464
465 set OutputArray stableParticles
466 set NeutralOutputArray neutralParticles
467 set ChargedHadronOutputArray chargedHadrons
468 set ElectronOutputArray electrons
469 set MuonOutputArray muons
470
471 # radius of the magnetic field coverage, in m
472 set Radius 0.3
473 set RadiusMax 1.2
474 # half-length of the magnetic field coverage, in m
475 set HalfLength 0.7
476 set HalfLengthMax 3.25
477
478 # magnetic field
479 set Bz 3.0
480}
481
482
483####################
484# Dense Track merger
485###################
486
487module Merger DenseMergeTracks {
488# add InputArray InputArray
489 add InputArray DenseProp/chargedHadrons
490 add InputArray DenseProp/electrons
491 add InputArray DenseProp/muons
492 set OutputArray tracks
493}
494
495
496######################
497# Dense Track Filter
498######################
499
500module DenseTrackFilter TrackMerger {
501 set TrackInputArray TrackMergerProp/tracks
502 set DenseChargedInputArray DenseMergeTracks/tracks
503
504 set TrackOutputArray tracks
505
506 set EtaPhiRes 0.003
507 set EtaMax 4.0
508
509 set pi [expr {acos(-1)}]
510
511 set nbins_phi [expr {$pi/$EtaPhiRes} ]
512 set nbins_phi [expr {int($nbins_phi)} ]
513
514 set PhiBins {}
515 for {set i -$nbins_phi} {$i <= $nbins_phi} {incr i} {
516 add PhiBins [expr {$i * $pi/$nbins_phi}]
517 }
518
519 set nbins_eta [expr {$EtaMax/$EtaPhiRes} ]
520 set nbins_eta [expr {int($nbins_eta)} ]
521
522 for {set i -$nbins_eta} {$i <= $nbins_eta} {incr i} {
523 set eta [expr {$i * $EtaPhiRes}]
524 add EtaPhiBins $eta $PhiBins
525 }
526}
527
528
529#############
530# ECAL
531#############
532
533module SimpleCalorimeter ECal {
534 set ParticleInputArray ParticlePropagator/stableParticles
535 set TrackInputArray TrackMerger/tracks
536
537 set TowerOutputArray ecalTowers
538 set EFlowTrackOutputArray eflowTracks
539 set EFlowTowerOutputArray eflowPhotons
540
541 set IsEcal true
542
543 set EnergyMin 0.5
544 set EnergySignificanceMin 1.0
545
546 set SmearTowerCenter true
547
548 set pi [expr {acos(-1)}]
549
550 # lists of the edges of each tower in eta and phi
551 # each list starts with the lower edge of the first tower
552 # the list ends with the higher edged of the last tower
553
554 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
555
556 set PhiBins {}
557 for {set i -180} {$i <= 180} {incr i} {
558 add PhiBins [expr {$i * $pi/180.0}]
559 }
560
561 # 0.02 unit in eta up to eta = 1.5 (barrel)
562 for {set i -85} {$i <= 86} {incr i} {
563 set eta [expr {$i * 0.0174}]
564 add EtaPhiBins $eta $PhiBins
565 }
566
567 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
568
569 set PhiBins {}
570 for {set i -180} {$i <= 180} {incr i} {
571 add PhiBins [expr {$i * $pi/180.0}]
572 }
573
574 # 0.02 unit in eta up to eta = 3
575 for {set i 1} {$i <= 84} {incr i} {
576 set eta [expr { -2.958 + $i * 0.0174}]
577 add EtaPhiBins $eta $PhiBins
578 }
579
580 for {set i 1} {$i <= 84} {incr i} {
581 set eta [expr { 1.4964 + $i * 0.0174}]
582 add EtaPhiBins $eta $PhiBins
583 }
584
585 # take present CMS granularity for HF
586
587 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
588 set PhiBins {}
589 for {set i -18} {$i <= 18} {incr i} {
590 add PhiBins [expr {$i * $pi/18.0}]
591 }
592
593 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} {
594 add EtaPhiBins $eta $PhiBins
595 }
596
597
598 add EnergyFraction {0} {0.0}
599 # energy fractions for e, gamma and pi0
600 add EnergyFraction {11} {1.0}
601 add EnergyFraction {22} {1.0}
602 add EnergyFraction {111} {1.0}
603 # energy fractions for muon, neutrinos and neutralinos
604 add EnergyFraction {12} {0.0}
605 add EnergyFraction {13} {0.0}
606 add EnergyFraction {14} {0.0}
607 add EnergyFraction {16} {0.0}
608 add EnergyFraction {1000022} {0.0}
609 add EnergyFraction {1000023} {0.0}
610 add EnergyFraction {1000025} {0.0}
611 add EnergyFraction {1000035} {0.0}
612 add EnergyFraction {1000045} {0.0}
613 # energy fractions for K0short and Lambda
614 add EnergyFraction {310} {0.3}
615 add EnergyFraction {3122} {0.3}
616
617
618 # ECAL is CMS
619 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
620 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117
621
622 set ResolutionFormula { 2*((abs(eta) <= 1.50) * sqrt(energy^2*0.009^2 + energy*0.12^2 + 0.45^2) +
623 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
624 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
625 (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2)) + \
626 (abs(eta) >= 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
627
628
629}
630
631#############
632# HCAL
633#############
634module SimpleCalorimeter HCal {
635 set ParticleInputArray ParticlePropagator/stableParticles
636 set TrackInputArray ECal/eflowTracks
637
638 set TowerOutputArray hcalTowers
639 set EFlowTrackOutputArray eflowTracks
640 set EFlowTowerOutputArray eflowNeutralHadrons
641
642 set IsEcal false
643
644 set EnergyMin 1.0
645 set EnergySignificanceMin 2.0
646
647 set SmearTowerCenter true
648
649 set pi [expr {acos(-1)}]
650
651 # lists of the edges of each tower in eta and phi
652 # each list starts with the lower edge of the first tower
653 # the list ends with the higher edged of the last tower
654
655 # 10 degrees towers
656 set PhiBins {}
657 for {set i -18} {$i <= 18} {incr i} {
658 add PhiBins [expr {$i * $pi/18.0}]
659 }
660 foreach eta {-3.2 -2.5 -2.4 -2.3 -2.2 -2.1 -2 -1.9 -1.8 -1.7 -1.6 -1.5 -1.4 -1.3 -1.2 -1.1 -1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 2.6 3.3} {
661 add EtaPhiBins $eta $PhiBins
662 }
663
664 # 20 degrees towers
665 set PhiBins {}
666 for {set i -9} {$i <= 9} {incr i} {
667 add PhiBins [expr {$i * $pi/9.0}]
668 }
669 foreach eta {-4.9 -4.7 -4.5 -4.3 -4.1 -3.9 -3.7 -3.5 -3.3 -3 -2.8 -2.6 2.8 3 3.2 3.5 3.7 3.9 4.1 4.3 4.5 4.7 4.9} {
670 add EtaPhiBins $eta $PhiBins
671 }
672
673 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
674 add EnergyFraction {0} {1.0}
675 # energy fractions for e, gamma and pi0
676 add EnergyFraction {11} {0.0}
677 add EnergyFraction {22} {0.0}
678 add EnergyFraction {111} {0.0}
679 # energy fractions for muon, neutrinos and neutralinos
680 add EnergyFraction {12} {0.0}
681 add EnergyFraction {13} {0.0}
682 add EnergyFraction {14} {0.0}
683 add EnergyFraction {16} {0.0}
684 add EnergyFraction {1000022} {0.0}
685 add EnergyFraction {1000023} {0.0}
686 add EnergyFraction {1000025} {0.0}
687 add EnergyFraction {1000035} {0.0}
688 add EnergyFraction {1000045} {0.0}
689 # energy fractions for K0short and Lambda
690 add EnergyFraction {310} {0.7}
691 add EnergyFraction {3122} {0.7}
692
693 # HCAL is ALTAS
694 # http://arxiv.org/pdf/hep-ex/0004009v1
695 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
696 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
697 set ResolutionFormula { 2*((abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
698 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2)) +
699 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
700}
701
702
703#################
704# Electron filter
705#################
706
707module PdgCodeFilter ElectronFilter {
708 set InputArray HCal/eflowTracks
709 set OutputArray electrons
710 set Invert true
711 add PdgCode {11}
712 add PdgCode {-11}
713}
714
715
716######################
717# ChargedHadronFilter
718######################
719
720module PdgCodeFilter ChargedHadronFilter {
721 set InputArray HCal/eflowTracks
722 set OutputArray chargedHadrons
723
724 add PdgCode {11}
725 add PdgCode {-11}
726 add PdgCode {13}
727 add PdgCode {-13}
728}
729
730
731###################################################
732# Tower Merger (in case not using e-flow algorithm)
733###################################################
734
735module Merger Calorimeter {
736# add InputArray InputArray
737 add InputArray ECal/ecalTowers
738 add InputArray HCal/hcalTowers
739 set OutputArray towers
740}
741
742####################
743# Energy flow merger
744####################
745
746module Merger EFlowMerger {
747# add InputArray InputArray
748 add InputArray HCal/eflowTracks
749 add InputArray ECal/eflowPhotons
750 add InputArray HCal/eflowNeutralHadrons
751 set OutputArray eflow
752}
753
754######################
755# EFlowFilter
756######################
757
758module PdgCodeFilter EFlowFilter {
759 set InputArray EFlowMerger/eflow
760 set OutputArray eflow
761
762 add PdgCode {11}
763 add PdgCode {-11}
764 add PdgCode {13}
765 add PdgCode {-13}
766}
767
768
769###################
770# Missing ET merger
771###################
772
773module Merger MissingET {
774# add InputArray InputArray
775 add InputArray EFlowMerger/eflow
776 set MomentumOutputArray momentum
777}
778
779
780
781##################
782# Scalar HT merger
783##################
784
785module Merger ScalarHT {
786# add InputArray InputArray
787 add InputArray EFlowMerger/eflow
788 set EnergyOutputArray energy
789}
790
791#################
792# Neutrino Filter
793#################
794
795module PdgCodeFilter NeutrinoFilter {
796
797 set InputArray Delphes/stableParticles
798 set OutputArray filteredParticles
799
800 set PTMin 0.0
801
802 add PdgCode {12}
803 add PdgCode {14}
804 add PdgCode {16}
805 add PdgCode {-12}
806 add PdgCode {-14}
807 add PdgCode {-16}
808
809}
810
811
812#####################
813# MC truth jet finder
814#####################
815
816module FastJetFinder GenJetFinder {
817 set InputArray NeutrinoFilter/filteredParticles
818
819 set OutputArray jets
820
821 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
822 set JetAlgorithm 6
823 set ParameterR 0.4
824
825 set JetPTMin 15.0
826}
827
828#########################
829# Gen Missing ET merger
830########################
831
832module Merger GenMissingET {
833
834# add InputArray InputArray
835 add InputArray NeutrinoFilter/filteredParticles
836 set MomentumOutputArray momentum
837}
838
839
840############
841# Jet finder
842############
843
844module FastJetFinder FastJetFinder {
845# set InputArray Calorimeter/towers
846 set InputArray EFlowMerger/eflow
847
848 set OutputArray jets
849
850 set JetAlgorithm 6
851 set ParameterR 0.4
852
853 set ComputeNsubjettiness 1
854 set Beta 1.0
855 set AxisMode 4
856
857 set ComputeSoftDrop 1
858 set BetaSoftDrop 0.0
859 set SymmetryCutSoftDrop 0.1
860 set R0SoftDrop 0.4
861
862 set JetPTMin 30.0
863}
864
865##################
866# Fat Jet finder
867##################
868
869module FastJetFinder FatJetFinder {
870 set InputArray EFlowMerger/eflow
871
872 set OutputArray jets
873
874 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
875 set JetAlgorithm 6
876 set ParameterR 0.8
877
878 set ComputeNsubjettiness 1
879 set Beta 1.0
880 set AxisMode 4
881
882 set ComputeSoftDrop 1
883 set BetaSoftDrop 0.0
884 set SymmetryCutSoftDrop 0.1
885 set R0SoftDrop 0.8
886
887 set JetPTMin 30.0
888}
889
890
891
892##################
893# Jet Energy Scale
894##################
895
896module EnergyScale JetEnergyScale {
897 set InputArray FastJetFinder/jets
898 set OutputArray jets
899
900 # scale formula for jets
901 set ScaleFormula {1.00}
902}
903
904
905########################
906# Jet Flavor Association
907########################
908
909module JetFlavorAssociation JetFlavorAssociation {
910
911 set PartonInputArray Delphes/partons
912 set ParticleInputArray Delphes/allParticles
913 set ParticleLHEFInputArray Delphes/allParticlesLHEF
914 set JetInputArray JetEnergyScale/jets
915
916 set DeltaR 0.5
917 set PartonPTMin 5.0
918 set PartonEtaMax 4.0
919
920}
921
922##################
923# Photon isolation
924##################
925
926# TBC: check values for iso cuts
927
928module Isolation PhotonIsolation {
929 set CandidateInputArray ECal/eflowPhotons
930 set IsolationInputArray EFlowFilter/eflow
931
932 set OutputArray photons
933
934 set DeltaRMax 0.3
935 set PTMin 0.0
936 set PTRatioMax 0.1
937}
938
939
940
941###################
942# Photon efficiency
943###################
944
945module Efficiency PhotonEfficiency {
946 set InputArray PhotonIsolation/photons
947 set OutputArray photons
948
949 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
950
951 set EfficiencyFormula {
952 (pt <= 1.0) * (0.00) + \
953 (abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0) * (0.70) +
954 (abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) +
955 (abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
956
957 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0) * (0.60) +
958 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 5.0 && pt < 10.0) * (0.80) +
959 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) +
960
961 (abs(eta) > 4.0) * (0.00)}
962
963}
964
965
966
967####################
968# Electron isolation
969####################
970
971# TBC: check values for iso cuts
972
973module Isolation ElectronIsolation {
974 set CandidateInputArray ElectronFilter/electrons
975 set IsolationInputArray EFlowFilter/eflow
976
977 set OutputArray electrons
978
979 set DeltaRMax 0.3
980
981 set PTMin 0.0
982
983 set PTRatioMax 0.1
984}
985
986
987#######################
988# Electron efficiency #
989#######################
990
991module Efficiency ElectronEfficiency {
992
993 set InputArray ElectronIsolation/electrons
994 set OutputArray electrons
995
996 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
997 # efficiency formula for electrons (taken from CMS PhaseII card)
998 # efficiency for low pT leptons is set artifically to 100%. Analyzers should rescale according to proper lepton Id
999 set EfficiencyFormula {
1000 (pt <= 2.0) * (0.00) +
1001 (abs(eta) <= 1.45 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1002 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
1003 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
1004 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
1005 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
1006 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
1007 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
1008 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
1009 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1010 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
1011 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
1012 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
1013 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
1014 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1015 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
1016 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
1017 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
1018 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
1019 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
1020 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1021 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1022 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1023 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1024 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1025 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
1026 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1027 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1028 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1029 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1030 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1031 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1032 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
1033 (abs(eta) > 4.0) * (0.00)
1034
1035 }
1036
1037}
1038
1039
1040
1041################
1042# Muon isolation
1043################
1044
1045# TBC: check values for iso cuts
1046
1047module Isolation MuonIsolation {
1048 set CandidateInputArray MuonMomentumSmearing/muons
1049 set IsolationInputArray EFlowFilter/eflow
1050
1051 set OutputArray muons
1052
1053 set DeltaRMax 0.3
1054
1055 set PTMin 0.0
1056
1057 set PTRatioMax 0.2
1058}
1059
1060####################
1061# Muon Efficiency #
1062####################
1063
1064module Efficiency MuonEfficiency {
1065 set InputArray MuonIsolation/muons
1066 set OutputArray muons
1067 # tracking + TightID efficiency formula for muons
1068 set EfficiencyFormula {
1069 (pt <= 2.0) * (0.00) +
1070 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 2.0 && pt <= 4.0) * (0.04) +
1071 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 4.0 && pt <= 6.0) * (0.43) +
1072 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 6.0 && pt <= 8.0) * (0.53) +
1073 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 8.0 && pt <= 10.0) * (0.68) +
1074 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 10.0 && pt <= 20.0) * (0.81) +
1075 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 20.0 && pt <= 35.0) * (0.91) +
1076 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 35.0 && pt <= 50.0) * (0.96) +
1077 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1078 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 2.0 && pt <= 4.0) * (0.05) +
1079 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 4.0 && pt <= 6.0) * (0.47) +
1080 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 6.0 && pt <= 8.0) * (0.56) +
1081 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 8.0 && pt <= 10.0) * (0.69) +
1082 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 10.0 && pt <= 20.0) * (0.79) +
1083 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 20.0 && pt <= 35.0) * (0.93) +
1084 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 35.0 && pt <= 50.0) * (0.94) +
1085 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1086 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 2.0 && pt <= 4.0) * (0.16) +
1087 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 4.0 && pt <= 6.0) * (0.48) +
1088 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 6.0 && pt <= 8.0) * (0.53) +
1089 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 8.0 && pt <= 10.0) * (0.66) +
1090 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 10.0 && pt <= 20.0) * (0.79) +
1091 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 20.0 && pt <= 35.0) * (0.89) +
1092 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 35.0 && pt <= 50.0) * (0.95) +
1093 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1094 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 2.0 && pt <= 4.0) * (0.24) +
1095 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 4.0 && pt <= 6.0) * (0.44) +
1096 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 6.0 && pt <= 8.0) * (0.51) +
1097 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 8.0 && pt <= 10.0) * (0.71) +
1098 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.77) +
1099 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 35.0) * (0.91) +
1100 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 35.0 && pt <= 50.0) * (0.92) +
1101 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1102 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 2.0 && pt <= 4.0) * (0.23) +
1103 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 4.0 && pt <= 6.0) * (0.35) +
1104 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 6.0 && pt <= 8.0) * (0.43) +
1105 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 8.0 && pt <= 10.0) * (0.57) +
1106 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.63) +
1107 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 35.0) * (0.71) +
1108 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 35.0 && pt <= 50.0) * (0.76) +
1109 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 50.0 && pt <= 14000.0) * (0.81) +
1110 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 2.0 && pt <= 4.0) * (0.05) + \
1111 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1112 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1113 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1114 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1115 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
1116 (abs(eta) > 4.0) * (0.00)
1117 }
1118}
1119
1120
1121
1122###########
1123# b-tagging
1124###########
1125
1126module BTagging BTagging {
1127 set JetInputArray JetEnergyScale/jets
1128
1129 set BitNumber 0
1130
1131 add EfficiencyFormula {0} {
1132
1133 (pt <= 10.0) * (0.00) +
1134 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.01) + \
1135 (abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.01)*(1.0 - pt/5000.) + \
1136 (abs(eta) < 2.5) * (pt > 5000.0) * (0.00) + \
1137 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.0075) + \
1138 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.0075)*(1.0 - pt/5000.) + \
1139 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) + \
1140 (abs(eta) > 4.0) * (0.00)}
1141
1142 add EfficiencyFormula {4} {
1143
1144 (pt <= 10.0) * (0.00) +
1145 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.10) + \
1146 (abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.10)*(1.0 - pt/5000.) + \
1147 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
1148 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.06) + \
1149 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.06)*(1.0 - pt/5000.) + \
1150 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) + \
1151 (abs(eta) > 4.0) * (0.00)}
1152
1153 add EfficiencyFormula {5} {
1154
1155 (pt <= 10.0) * (0.00) +
1156 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.75) +
1157 (abs(eta) < 2.5) * (pt > 500.0 && pt < 5000.0) * (0.75)*(1.0 - pt/5000.) +
1158 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) +
1159 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.60) +
1160 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 5000.0) * (0.60)*(1.0 - pt/5000.) +
1161 (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.000) +
1162 (abs(eta) >= 4.0) * (0.00)}
1163
1164}
1165
1166
1167#############
1168# tau-tagging
1169#############
1170
1171
1172module TauTagging TauTagging {
1173 set ParticleInputArray Delphes/allParticles
1174 set PartonInputArray Delphes/partons
1175 set JetInputArray JetEnergyScale/jets
1176
1177 set DeltaR 0.5
1178
1179 set TauPTMin 1.0
1180
1181 set TauEtaMax 4.0
1182
1183 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1184 add EfficiencyFormula {0} {
1185
1186 (pt <= 10.0) * (0.00) +
1187 (abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.01) + \
1188 (abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.01) *(1. - pt/5000.) + \
1189 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
1190 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.0075) + \
1191 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.0075)*(1. - pt/5000.) + \
1192 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
1193 (abs(eta) > 4.0) * (0.00)}
1194
1195 add EfficiencyFormula {11} {
1196
1197 (pt <= 10.0) * (0.00) +
1198 (abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.005) + \
1199 (abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.005) *(1. - pt/5000.) + \
1200 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
1201 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.00375) + \
1202 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.00375)*(1. - pt/5000.) + \
1203 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
1204 (abs(eta) > 4.0) * (0.00)}
1205
1206 add EfficiencyFormula {15} {
1207
1208 (pt <= 10.0) * (0.00) +
1209 (abs(eta) < 2.5) * (pt > 10.0 && pt < 1000.0) * (0.6) + \
1210 (abs(eta) < 2.5) * (pt > 1000.0 && pt < 5000.0) * (0.6) *(1. - pt/5000.) + \
1211 (abs(eta) < 2.5) * (pt > 5000.0) * (0.000) + \
1212 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 1000.0) * (0.45) + \
1213 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 1000.0 && pt < 5000.0) * (0.45)*(1. - pt/5000.) + \
1214 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0) * (0.00) + \
1215 (abs(eta) > 4.0) * (0.00)}
1216
1217}
1218
1219#####################################################
1220# Find uniquely identified photons/electrons/tau/jets
1221#####################################################
1222
1223module UniqueObjectFinder UniqueObjectFinder {
1224# earlier arrays take precedence over later ones
1225# add InputArray InputArray OutputArray
1226 add InputArray PhotonEfficiency/photons photons
1227 add InputArray ElectronEfficiency/electrons electrons
1228 add InputArray MuonEfficiency/muons muons
1229 add InputArray JetEnergyScale/jets jets
1230
1231}
1232
1233
1234##################
1235# ROOT tree writer
1236##################
1237
1238module TreeWriter TreeWriter {
1239# add Branch InputArray BranchName BranchClass
1240 add Branch Delphes/allParticles Particle GenParticle
1241
1242 add Branch GenJetFinder/jets GenJet Jet
1243 add Branch GenMissingET/momentum GenMissingET MissingET
1244
1245 add Branch TrackMerger/tracks Track Track
1246 add Branch Calorimeter/towers Tower Tower
1247
1248 add Branch HCal/eflowTracks EFlowTrack Track
1249 add Branch ECal/eflowPhotons EFlowPhoton Tower
1250 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1251
1252 add Branch UniqueObjectFinder/photons Photon Photon
1253 add Branch UniqueObjectFinder/electrons Electron Electron
1254 add Branch UniqueObjectFinder/muons Muon Muon
1255 add Branch UniqueObjectFinder/jets Jet Jet
1256
1257 add Branch FatJetFinder/jets FatJet Jet
1258
1259 add Branch MissingET/momentum MissingET MissingET
1260 add Branch ScalarHT/energy ScalarHT ScalarHT
1261}
Note: See TracBrowser for help on using the repository browser.