Fork me on GitHub

source: git/cards/delphes_card_HLLHC.tcl@ 8a58fff

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 8a58fff was 8a58fff, checked in by Michele Selvaggi <michele.selvaggi@…>, 6 years ago

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