MA5SandBox: delphes_card_atlas_exot_2018_30.tcl

File delphes_card_atlas_exot_2018_30.tcl, 62.3 KB (added by Benjamin Fuks, 4 years ago)
Line 
1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6 ParticlePropagator
7
8 ChargedHadronTrackingEfficiency
9 ElectronTrackingEfficiency
10 MuonTrackingEfficiency
11
12 ChargedHadronMomentumSmearing
13 ElectronMomentumSmearing
14 MuonMomentumSmearing
15
16 TrackMerger
17
18 ECal
19 HCal
20
21 Calorimeter
22 EFlowMerger
23 EFlowFilter
24
25 PhotonEfficiency
26
27 ElectronFilter
28 ElectronEfficiency
29
30 ChargedHadronFilter
31
32 MuonEfficiency
33
34 MissingET
35
36 NeutrinoFilter
37 GenJetFinder
38 GenMissingET
39
40 FastJetFinder
41
42 JetEnergyScale
43
44 JetFlavorAssociation
45
46 BTagging
47 TauTagging
48
49 UniqueObjectFinder
50
51 ScalarHT
52
53 TreeWriter
54}
55
56#################################
57# Propagate particles in cylinder
58#################################
59
60module ParticlePropagator ParticlePropagator {
61 set InputArray Delphes/stableParticles
62
63 set OutputArray stableParticles
64 set ChargedHadronOutputArray chargedHadrons
65 set ElectronOutputArray electrons
66 set MuonOutputArray muons
67
68 # radius of the magnetic field coverage, in m
69 set Radius 1.15
70 # half-length of the magnetic field coverage, in m
71 set HalfLength 3.51
72
73 # magnetic field
74 set Bz 2.0
75}
76
77####################################
78# Charged hadron tracking efficiency
79####################################
80
81module Efficiency ChargedHadronTrackingEfficiency {
82 set InputArray ParticlePropagator/chargedHadrons
83 set OutputArray chargedHadrons
84
85 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
86
87 # tracking efficiency formula for charged hadrons
88 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
89 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
90 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
91 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
92 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
93 (abs(eta) > 2.5) * (0.00)}
94}
95
96##############################
97# Electron tracking efficiency
98##############################
99
100module Efficiency ElectronTrackingEfficiency {
101 set InputArray ParticlePropagator/electrons
102 set OutputArray electrons
103
104 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
105
106 # tracking efficiency formula for electrons
107 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
108 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
109 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
110 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
111 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
112 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
113 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
114 (abs(eta) > 2.5) * (0.00)}
115}
116
117##########################
118# Muon tracking efficiency
119##########################
120
121module Efficiency MuonTrackingEfficiency {
122 set InputArray ParticlePropagator/muons
123 set OutputArray muons
124
125 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
126
127 # tracking efficiency formula for muons
128 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
129 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
130 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) +
131 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
132 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) +
133 (abs(eta) > 2.5) * (0.00)}
134}
135
136########################################
137# Momentum resolution for charged tracks
138########################################
139
140module MomentumSmearing ChargedHadronMomentumSmearing {
141 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
142 set OutputArray chargedHadrons
143
144 # set ResolutionFormula {resolution formula as a function of eta and pt}
145
146 # resolution formula for charged hadrons
147 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
148 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
149 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
150}
151
152###################################
153# Momentum resolution for electrons
154###################################
155
156module MomentumSmearing ElectronMomentumSmearing {
157 set InputArray ElectronTrackingEfficiency/electrons
158 set OutputArray electrons
159
160 # set ResolutionFormula {resolution formula as a function of eta and energy}
161
162 # resolution formula for electrons
163 set ResolutionFormula {
164 (abs(eta) < 2.5) * (pt > 0.1) * 0.0009 }
165
166}
167
168###############################
169# Momentum resolution for muons
170###############################
171
172module MomentumSmearing MuonMomentumSmearing {
173 set InputArray MuonTrackingEfficiency/muons
174 set OutputArray muons
175
176 # set ResolutionFormula {resolution formula as a function of eta and pt}
177 # resolution formula for muons
178 set ResolutionFormula {
179 (abs(eta) <= 0.5) * (pt < 1000 && pt > 0.1) * sqrt(4.1e-3^2 + pt^2*1.7e-4^2) +
180 (abs(eta) > 0.5 && abs(eta) < 1.05) * (pt < 1000 && pt > 0.1) * sqrt(4.1e-3^2 + pt^2*1.7e-4^2) +
181 (abs(eta) >= 1.05 && abs(eta) <= 1.5) * (pt < 1000 && pt > 0.1) * sqrt(5.5e-3^2 + pt^2*3.4e-4^2) +
182 (abs(eta) > 1.5 && abs(eta) < 2.0) * (pt < 1000 && pt > 0.1) * sqrt(5.5e-3^2 + pt^2*3.4e-4^2) +
183 (abs(eta) >= 2.0) * (pt < 1000 && pt > 0.1) * sqrt(9e-3^2 + pt^2*5e-5^2) +
184
185 (abs(eta) < 1.05) * (pt >= 1000 ) * (4.1e-3+pt*1.7e-4) +
186 (abs(eta) >= 1.05 && abs(eta) < 2.0) * (pt >= 1000 ) * (5.5e-3+pt*3.4e-4) +
187 (abs(eta) >= 2.0) * (pt >= 1000 ) * (9e-3+pt*5e-5) }
188
189}
190
191
192
193
194##############
195# Track merger
196##############
197
198module Merger TrackMerger {
199# add InputArray InputArray
200 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
201 add InputArray ElectronMomentumSmearing/electrons
202 add InputArray MuonMomentumSmearing/muons
203 set OutputArray tracks
204}
205
206
207#############
208# ECAL
209#############
210
211module SimpleCalorimeter ECal {
212 set ParticleInputArray ParticlePropagator/stableParticles
213 set TrackInputArray TrackMerger/tracks
214
215 set TowerOutputArray ecalTowers
216 set EFlowTrackOutputArray eflowTracks
217 set EFlowTowerOutputArray eflowPhotons
218
219 set IsEcal true
220
221 set EnergyMin 0.5
222 set EnergySignificanceMin 2.0
223
224 set SmearTowerCenter true
225
226 set pi [expr {acos(-1)}]
227
228 # lists of the edges of each tower in eta and phi
229 # each list starts with the lower edge of the first tower
230 # the list ends with the higher edged of the last tower
231
232 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
233
234 set PhiBins {}
235 for {set i -180} {$i <= 180} {incr i} {
236 add PhiBins [expr {$i * $pi/180.0}]
237 }
238
239 # 0.02 unit in eta up to eta = 1.5 (barrel)
240 for {set i -85} {$i <= 86} {incr i} {
241 set eta [expr {$i * 0.0174}]
242 add EtaPhiBins $eta $PhiBins
243 }
244
245 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0
246 set PhiBins {}
247 for {set i -180} {$i <= 180} {incr i} {
248 add PhiBins [expr {$i * $pi/180.0}]
249 }
250
251 # 0.02 unit in eta up to eta = 3
252 for {set i 1} {$i <= 84} {incr i} {
253 set eta [expr { -2.958 + $i * 0.0174}]
254 add EtaPhiBins $eta $PhiBins
255 }
256
257 for {set i 1} {$i <= 84} {incr i} {
258 set eta [expr { 1.4964 + $i * 0.0174}]
259 add EtaPhiBins $eta $PhiBins
260 }
261
262 # take present CMS granularity for HF
263
264 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
265 set PhiBins {}
266 for {set i -18} {$i <= 18} {incr i} {
267 add PhiBins [expr {$i * $pi/18.0}]
268 }
269
270 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} {
271 add EtaPhiBins $eta $PhiBins
272 }
273
274
275 add EnergyFraction {0} {0.0}
276 # energy fractions for e, gamma and pi0
277 add EnergyFraction {11} {1.0}
278 add EnergyFraction {22} {1.0}
279 add EnergyFraction {111} {1.0}
280 # energy fractions for muon, neutrinos and neutralinos
281 add EnergyFraction {12} {0.0}
282 add EnergyFraction {13} {0.0}
283 add EnergyFraction {14} {0.0}
284 add EnergyFraction {16} {0.0}
285 add EnergyFraction {1000022} {0.0}
286 add EnergyFraction {1000023} {0.0}
287 add EnergyFraction {1000025} {0.0}
288 add EnergyFraction {1000035} {0.0}
289 add EnergyFraction {1000045} {0.0}
290 # energy fractions for K0short and Lambda
291 add EnergyFraction {310} {0.3}
292 add EnergyFraction {3122} {0.3}
293
294 # set ResolutionFormula {resolution formula as a function of eta and energy}
295
296 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
297 # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
298 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
299 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
300 set ResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) +
301 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
302
303
304}
305
306
307
308#############
309# HCAL
310#############
311
312module SimpleCalorimeter HCal {
313 set ParticleInputArray ParticlePropagator/stableParticles
314 set TrackInputArray ECal/eflowTracks
315
316 set TowerOutputArray hcalTowers
317 set EFlowTrackOutputArray eflowTracks
318 set EFlowTowerOutputArray eflowNeutralHadrons
319
320 set IsEcal false
321
322 set EnergyMin 1.0
323 set EnergySignificanceMin 2.0
324
325 set SmearTowerCenter true
326
327 set pi [expr {acos(-1)}]
328
329 # lists of the edges of each tower in eta and phi
330 # each list starts with the lower edge of the first tower
331 # the list ends with the higher edged of the last tower
332
333 # 10 degrees towers
334 set PhiBins {}
335 for {set i -18} {$i <= 18} {incr i} {
336 add PhiBins [expr {$i * $pi/18.0}]
337 }
338 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} {
339 add EtaPhiBins $eta $PhiBins
340 }
341
342 # 20 degrees towers
343 set PhiBins {}
344 for {set i -9} {$i <= 9} {incr i} {
345 add PhiBins [expr {$i * $pi/9.0}]
346 }
347 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} {
348 add EtaPhiBins $eta $PhiBins
349 }
350
351 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
352 add EnergyFraction {0} {1.0}
353 # energy fractions for e, gamma and pi0
354 add EnergyFraction {11} {0.0}
355 add EnergyFraction {22} {0.0}
356 add EnergyFraction {111} {0.0}
357 # energy fractions for muon, neutrinos and neutralinos
358 add EnergyFraction {12} {0.0}
359 add EnergyFraction {13} {0.0}
360 add EnergyFraction {14} {0.0}
361 add EnergyFraction {16} {0.0}
362 add EnergyFraction {1000022} {0.0}
363 add EnergyFraction {1000023} {0.0}
364 add EnergyFraction {1000025} {0.0}
365 add EnergyFraction {1000035} {0.0}
366 add EnergyFraction {1000045} {0.0}
367 # energy fractions for K0short and Lambda
368 add EnergyFraction {310} {0.7}
369 add EnergyFraction {3122} {0.7}
370
371 # http://arxiv.org/pdf/hep-ex/0004009v1
372 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
373 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
374 set ResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
375 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
376 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
377}
378
379
380#################
381# Electron filter
382#################
383
384module PdgCodeFilter ElectronFilter {
385 set InputArray HCal/eflowTracks
386 set OutputArray electrons
387 set Invert true
388 add PdgCode {11}
389 add PdgCode {-11}
390}
391
392######################
393# ChargedHadronFilter
394######################
395
396module PdgCodeFilter ChargedHadronFilter {
397 set InputArray HCal/eflowTracks
398 set OutputArray chargedHadrons
399
400 add PdgCode {11}
401 add PdgCode {-11}
402 add PdgCode {13}
403 add PdgCode {-13}
404}
405
406
407
408###################################################
409# Tower Merger (in case not using e-flow algorithm)
410###################################################
411
412module Merger Calorimeter {
413# add InputArray InputArray
414 add InputArray ECal/ecalTowers
415 add InputArray HCal/hcalTowers
416 add InputArray MuonMomentumSmearing/muons
417 set OutputArray towers
418}
419
420####################
421# Energy flow merger
422####################
423
424module Merger EFlowMerger {
425# add InputArray InputArray
426 add InputArray HCal/eflowTracks
427 add InputArray ECal/eflowPhotons
428 add InputArray HCal/eflowNeutralHadrons
429 set OutputArray eflow
430}
431
432######################
433# EFlowFilter
434######################
435
436module PdgCodeFilter EFlowFilter {
437 set InputArray EFlowMerger/eflow
438 set OutputArray eflow
439
440 add PdgCode {11}
441 add PdgCode {-11}
442 add PdgCode {13}
443 add PdgCode {-13}
444}
445
446###################
447# Photon efficiency
448###################
449
450module Efficiency PhotonEfficiency {
451 set InputArray ECal/eflowPhotons
452 set OutputArray photons
453
454 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
455
456 # efficiency formula for photons
457 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
458 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
459 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
460 (abs(eta) > 2.5) * (0.00)}
461}
462
463##################
464# Photon isolation
465##################
466
467module Isolation PhotonIsolation {
468# set CandidateInputArray PhotonEfficiency/photons
469# set IsolationInputArray EFlowFilter/eflow
470
471# set OutputArray photons
472
473# set DeltaRMax 0.5
474
475# set PTMin 0.5
476
477# set PTRatioMax 0.12
478}
479
480
481#####################
482# Electron efficiency
483#####################
484
485module Efficiency ElectronEfficiency {
486 set InputArray ElectronFilter/electrons
487 set OutputArray electrons
488
489 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
490
491 # efficiency formula for electrons
492 ### First, use eff vs. eta 1D plot (pT > 4.5)
493 ### Papers always give Data eff, and SF. We have to use Data eff/SF for MC eff
494 ### For eff vs. pT 1D plot, MC eff (pT : 35~50 averaged = 0.9401 )
495 ### For other pT region, we multiply (eff of pT / eff of pT(35~45)=0.7819)
496
497 # ### Tight-ID * pT eff ratio
498 set EfficiencyFormula {
499
500 (pt <= 4.5) * (0.00) +
501
502 (eta >= -2.47 && eta < -2.37 ) * (pt > 4.5 && pt <= 25.0) * (0.6532 / 0.8416 * 0.7328 / 0.8108) +
503 (eta >= -2.37 && eta < -2.01 ) * (pt > 4.5 && pt <= 25.0) * (0.7238 / 0.9567 * 0.7328 / 0.8108) +
504 (eta >= -2.01 && eta < -1.81 ) * (pt > 4.5 && pt <= 25.0) * (0.7402 / 0.9520 * 0.7328 / 0.8108) +
505 (eta >= -1.81 && eta < -1.52 ) * (pt > 4.5 && pt <= 25.0) * (0.7639 / 0.9346 * 0.7328 / 0.8108) +
506 (eta >= -1.52 && eta < -1.37 ) * (pt > 4.5 && pt <= 25.0) * (0.6310 / 0.9126 * 0.7328 / 0.8108) +
507 (eta >= -1.37 && eta < -1.15 ) * (pt > 4.5 && pt <= 25.0) * (0.7532 / 0.9220 * 0.7328 / 0.8108) +
508 (eta >= -1.15 && eta < -0.80 ) * (pt > 4.5 && pt <= 25.0) * (0.7808 / 0.9551 * 0.7328 / 0.8108) +
509 (eta >= -0.80 && eta < -0.60 ) * (pt > 4.5 && pt <= 25.0) * (0.8058 / 0.9756 * 0.7328 / 0.8108) +
510 (eta >= -0.60 && eta < -0.10 ) * (pt > 4.5 && pt <= 25.0) * (0.8136 / 0.9898 * 0.7328 / 0.8108) +
511 (eta >= -0.10 && eta < 0.00 ) * (pt > 4.5 && pt <= 25.0) * (0.7656 / 0.9803 * 0.7328 / 0.8108) +
512 (eta >= 0.00 && eta < 0.10 ) * (pt > 4.5 && pt <= 25.0) * (0.7466 / 0.9709 * 0.7328 / 0.8108) +
513 (eta >= 0.10 && eta < 0.60 ) * (pt > 4.5 && pt <= 25.0) * (0.8146 / 0.9866 * 0.7328 / 0.8108) +
514 (eta >= 0.60 && eta < 0.80 ) * (pt > 4.5 && pt <= 25.0) * (0.8112 / 0.9787 * 0.7328 / 0.8108) +
515 (eta >= 0.80 && eta < 1.15 ) * (pt > 4.5 && pt <= 25.0) * (0.7850 / 0.9598 * 0.7328 / 0.8108) +
516 (eta >= 1.15 && eta < 1.37 ) * (pt > 4.5 && pt <= 25.0) * (0.7584 / 0.9283 * 0.7328 / 0.8108) +
517 (eta >= 1.37 && eta < 1.52 ) * (pt > 4.5 && pt <= 25.0) * (0.6457 / 0.9268 * 0.7328 / 0.8108) +
518 (eta >= 1.52 && eta < 1.81 ) * (pt > 4.5 && pt <= 25.0) * (0.7731 / 0.9425 * 0.7328 / 0.8108) +
519 (eta >= 1.81 && eta < 2.01 ) * (pt > 4.5 && pt <= 25.0) * (0.7400 / 0.9583 * 0.7328 / 0.8108) +
520 (eta >= 2.01 && eta < 2.37 ) * (pt > 4.5 && pt <= 25.0) * (0.7271 / 0.9567 * 0.7328 / 0.8108) +
521 (eta >= 2.37 && eta < 2.47 ) * (pt > 4.5 && pt <= 25.0) * (0.6712 / 0.8559 * 0.7328 / 0.8108) +
522
523 (eta >= -2.47 && eta < -2.37 ) * (pt > 25.0 && pt <= 30.0) * (0.6532 / 0.8416 * 0.7613 / 0.8108) +
524 (eta >= -2.37 && eta < -2.01 ) * (pt > 25.0 && pt <= 30.0) * (0.7238 / 0.9567 * 0.7613 / 0.8108) +
525 (eta >= -2.01 && eta < -1.81 ) * (pt > 25.0 && pt <= 30.0) * (0.7402 / 0.9520 * 0.7613 / 0.8108) +
526 (eta >= -1.81 && eta < -1.52 ) * (pt > 25.0 && pt <= 30.0) * (0.7639 / 0.9346 * 0.7613 / 0.8108) +
527 (eta >= -1.52 && eta < -1.37 ) * (pt > 25.0 && pt <= 30.0) * (0.6310 / 0.9126 * 0.7613 / 0.8108) +
528 (eta >= -1.37 && eta < -1.15 ) * (pt > 25.0 && pt <= 30.0) * (0.7532 / 0.9220 * 0.7613 / 0.8108) +
529 (eta >= -1.15 && eta < -0.80 ) * (pt > 25.0 && pt <= 30.0) * (0.7808 / 0.9551 * 0.7613 / 0.8108) +
530 (eta >= -0.80 && eta < -0.60 ) * (pt > 25.0 && pt <= 30.0) * (0.8058 / 0.9756 * 0.7613 / 0.8108) +
531 (eta >= -0.60 && eta < -0.10 ) * (pt > 25.0 && pt <= 30.0) * (0.8136 / 0.9898 * 0.7613 / 0.8108) +
532 (eta >= -0.10 && eta < 0.00 ) * (pt > 25.0 && pt <= 30.0) * (0.7656 / 0.9803 * 0.7613 / 0.8108) +
533 (eta >= 0.00 && eta < 0.10 ) * (pt > 25.0 && pt <= 30.0) * (0.7466 / 0.9709 * 0.7613 / 0.8108) +
534 (eta >= 0.10 && eta < 0.60 ) * (pt > 25.0 && pt <= 30.0) * (0.8146 / 0.9866 * 0.7613 / 0.8108) +
535 (eta >= 0.60 && eta < 0.80 ) * (pt > 25.0 && pt <= 30.0) * (0.8112 / 0.9787 * 0.7613 / 0.8108) +
536 (eta >= 0.80 && eta < 1.15 ) * (pt > 25.0 && pt <= 30.0) * (0.7850 / 0.9598 * 0.7613 / 0.8108) +
537 (eta >= 1.15 && eta < 1.37 ) * (pt > 25.0 && pt <= 30.0) * (0.7584 / 0.9283 * 0.7613 / 0.8108) +
538 (eta >= 1.37 && eta < 1.52 ) * (pt > 25.0 && pt <= 30.0) * (0.6457 / 0.9268 * 0.7613 / 0.8108) +
539 (eta >= 1.52 && eta < 1.81 ) * (pt > 25.0 && pt <= 30.0) * (0.7731 / 0.9425 * 0.7613 / 0.8108) +
540 (eta >= 1.81 && eta < 2.01 ) * (pt > 25.0 && pt <= 30.0) * (0.7400 / 0.9583 * 0.7613 / 0.8108) +
541 (eta >= 2.01 && eta < 2.37 ) * (pt > 25.0 && pt <= 30.0) * (0.7271 / 0.9567 * 0.7613 / 0.8108) +
542 (eta >= 2.37 && eta < 2.47 ) * (pt > 25.0 && pt <= 30.0) * (0.6712 / 0.8559 * 0.7613 / 0.8108) +
543
544 (eta >= -2.47 && eta < -2.37 ) * (pt > 30.0 && pt <= 35.0) * (0.6532 / 0.8416 * 0.7811 / 0.8108) +
545 (eta >= -2.37 && eta < -2.01 ) * (pt > 30.0 && pt <= 35.0) * (0.7238 / 0.9567 * 0.7811 / 0.8108) +
546 (eta >= -2.01 && eta < -1.81 ) * (pt > 30.0 && pt <= 35.0) * (0.7402 / 0.9520 * 0.7811 / 0.8108) +
547 (eta >= -1.81 && eta < -1.52 ) * (pt > 30.0 && pt <= 35.0) * (0.7639 / 0.9346 * 0.7811 / 0.8108) +
548 (eta >= -1.52 && eta < -1.37 ) * (pt > 30.0 && pt <= 35.0) * (0.6310 / 0.9126 * 0.7811 / 0.8108) +
549 (eta >= -1.37 && eta < -1.15 ) * (pt > 30.0 && pt <= 35.0) * (0.7532 / 0.9220 * 0.7811 / 0.8108) +
550 (eta >= -1.15 && eta < -0.80 ) * (pt > 30.0 && pt <= 35.0) * (0.7808 / 0.9551 * 0.7811 / 0.8108) +
551 (eta >= -0.80 && eta < -0.60 ) * (pt > 30.0 && pt <= 35.0) * (0.8058 / 0.9756 * 0.7811 / 0.8108) +
552 (eta >= -0.60 && eta < -0.10 ) * (pt > 30.0 && pt <= 35.0) * (0.8136 / 0.9898 * 0.7811 / 0.8108) +
553 (eta >= -0.10 && eta < 0.00 ) * (pt > 30.0 && pt <= 35.0) * (0.7656 / 0.9803 * 0.7811 / 0.8108) +
554 (eta >= 0.00 && eta < 0.10 ) * (pt > 30.0 && pt <= 35.0) * (0.7466 / 0.9709 * 0.7811 / 0.8108) +
555 (eta >= 0.10 && eta < 0.60 ) * (pt > 30.0 && pt <= 35.0) * (0.8146 / 0.9866 * 0.7811 / 0.8108) +
556 (eta >= 0.60 && eta < 0.80 ) * (pt > 30.0 && pt <= 35.0) * (0.8112 / 0.9787 * 0.7811 / 0.8108) +
557 (eta >= 0.80 && eta < 1.15 ) * (pt > 30.0 && pt <= 35.0) * (0.7850 / 0.9598 * 0.7811 / 0.8108) +
558 (eta >= 1.15 && eta < 1.37 ) * (pt > 30.0 && pt <= 35.0) * (0.7584 / 0.9283 * 0.7811 / 0.8108) +
559 (eta >= 1.37 && eta < 1.52 ) * (pt > 30.0 && pt <= 35.0) * (0.6457 / 0.9268 * 0.7811 / 0.8108) +
560 (eta >= 1.52 && eta < 1.81 ) * (pt > 30.0 && pt <= 35.0) * (0.7731 / 0.9425 * 0.7811 / 0.8108) +
561 (eta >= 1.81 && eta < 2.01 ) * (pt > 30.0 && pt <= 35.0) * (0.7400 / 0.9583 * 0.7811 / 0.8108) +
562 (eta >= 2.01 && eta < 2.37 ) * (pt > 30.0 && pt <= 35.0) * (0.7271 / 0.9567 * 0.7811 / 0.8108) +
563 (eta >= 2.37 && eta < 2.47 ) * (pt > 30.0 && pt <= 35.0) * (0.6712 / 0.8559 * 0.7811 / 0.8108) +
564
565 (eta >= -2.47 && eta < -2.37 ) * (pt > 35.0 && pt <= 50.0) * (0.6532 / 0.8416) +
566 (eta >= -2.37 && eta < -2.01 ) * (pt > 35.0 && pt <= 50.0) * (0.7238 / 0.9567) +
567 (eta >= -2.01 && eta < -1.81 ) * (pt > 35.0 && pt <= 50.0) * (0.7402 / 0.9520) +
568 (eta >= -1.81 && eta < -1.52 ) * (pt > 35.0 && pt <= 50.0) * (0.7639 / 0.9346) +
569 (eta >= -1.52 && eta < -1.37 ) * (pt > 35.0 && pt <= 50.0) * (0.6310 / 0.9126) +
570 (eta >= -1.37 && eta < -1.15 ) * (pt > 35.0 && pt <= 50.0) * (0.7532 / 0.9220) +
571 (eta >= -1.15 && eta < -0.80 ) * (pt > 35.0 && pt <= 50.0) * (0.7808 / 0.9551) +
572 (eta >= -0.80 && eta < -0.60 ) * (pt > 35.0 && pt <= 50.0) * (0.8058 / 0.9756) +
573 (eta >= -0.60 && eta < -0.10 ) * (pt > 35.0 && pt <= 50.0) * (0.8136 / 0.9898) +
574 (eta >= -0.10 && eta < 0.00 ) * (pt > 35.0 && pt <= 50.0) * (0.7656 / 0.9803) +
575 (eta >= 0.00 && eta < 0.10 ) * (pt > 35.0 && pt <= 50.0) * (0.7466 / 0.9709) +
576 (eta >= 0.10 && eta < 0.60 ) * (pt > 35.0 && pt <= 50.0) * (0.8146 / 0.9866) +
577 (eta >= 0.60 && eta < 0.80 ) * (pt > 35.0 && pt <= 50.0) * (0.8112 / 0.9787) +
578 (eta >= 0.80 && eta < 1.15 ) * (pt > 35.0 && pt <= 50.0) * (0.7850 / 0.9598) +
579 (eta >= 1.15 && eta < 1.37 ) * (pt > 35.0 && pt <= 50.0) * (0.7584 / 0.9283) +
580 (eta >= 1.37 && eta < 1.52 ) * (pt > 35.0 && pt <= 50.0) * (0.6457 / 0.9268) +
581 (eta >= 1.52 && eta < 1.81 ) * (pt > 35.0 && pt <= 50.0) * (0.7731 / 0.9425) +
582 (eta >= 1.81 && eta < 2.01 ) * (pt > 35.0 && pt <= 50.0) * (0.7400 / 0.9583) +
583 (eta >= 2.01 && eta < 2.37 ) * (pt > 35.0 && pt <= 50.0) * (0.7271 / 0.9567) +
584 (eta >= 2.37 && eta < 2.47 ) * (pt > 35.0 && pt <= 50.0) * (0.6712 / 0.8559) +
585
586 (eta >= -2.47 && eta < -2.37 ) * (pt > 50.0 && pt <= 60.0) * (0.6487 / 0.8310) * (0.8170 / 0.9621) +
587 (eta >= -2.37 && eta < -2.01 ) * (pt > 50.0 && pt <= 60.0) * (0.7237 / 0.9551) * (0.8170 / 0.9621) +
588 (eta >= -2.01 && eta < -1.81 ) * (pt > 50.0 && pt <= 60.0) * (0.7402 / 0.9528) * (0.8170 / 0.9621) +
589 (eta >= -1.81 && eta < -1.52 ) * (pt > 50.0 && pt <= 60.0) * (0.7644 / 0.9339) * (0.8170 / 0.9621) +
590 (eta >= -1.52 && eta < -1.37 ) * (pt > 50.0 && pt <= 60.0) * (0.6303 / 0.9126) * (0.8170 / 0.9621) +
591 (eta >= -1.37 && eta < -1.15 ) * (pt > 50.0 && pt <= 60.0) * (0.7526 / 0.9222) * (0.8170 / 0.9621) +
592 (eta >= -1.15 && eta < -0.80 ) * (pt > 50.0 && pt <= 60.0) * (0.7803 / 0.9558) * (0.8170 / 0.9621) +
593 (eta >= -0.80 && eta < -0.60 ) * (pt > 50.0 && pt <= 60.0) * (0.8053 / 0.9749) * (0.8170 / 0.9621) +
594 (eta >= -0.60 && eta < -0.10 ) * (pt > 50.0 && pt <= 60.0) * (0.8132 / 0.9871) * (0.8170 / 0.9621) +
595 (eta >= -0.10 && eta < 0.00 ) * (pt > 50.0 && pt <= 60.0) * (0.7618 / 0.9753) * (0.8170 / 0.9621) +
596 (eta >= 0.00 && eta < 0.10 ) * (pt > 50.0 && pt <= 60.0) * (0.7421 / 0.9635) * (0.8170 / 0.9621) +
597 (eta >= 0.10 && eta < 0.60 ) * (pt > 50.0 && pt <= 60.0) * (0.8145 / 0.9852) * (0.8170 / 0.9621) +
598 (eta >= 0.60 && eta < 0.80 ) * (pt > 50.0 && pt <= 60.0) * (0.8105 / 0.9782) * (0.8170 / 0.9621) +
599 (eta >= 0.80 && eta < 1.15 ) * (pt > 50.0 && pt <= 60.0) * (0.7842 / 0.9617) * (0.8170 / 0.9621) +
600 (eta >= 1.15 && eta < 1.37 ) * (pt > 50.0 && pt <= 60.0) * (0.7579 / 0.9309) * (0.8170 / 0.9621) +
601 (eta >= 1.37 && eta < 1.52 ) * (pt > 50.0 && pt <= 60.0) * (0.6461 / 0.9286) * (0.8170 / 0.9621) +
602 (eta >= 1.52 && eta < 1.81 ) * (pt > 50.0 && pt <= 60.0) * (0.7724 / 0.9454) * (0.8170 / 0.9621) +
603 (eta >= 1.81 && eta < 2.01 ) * (pt > 50.0 && pt <= 60.0) * (0.7439 / 0.9598) * (0.8170 / 0.9621) +
604 (eta >= 2.01 && eta < 2.37 ) * (pt > 50.0 && pt <= 60.0) * (0.7263 / 0.9577) * (0.8170 / 0.9621) +
605 (eta >= 2.37 && eta < 2.47 ) * (pt > 50.0 && pt <= 60.0) * (0.6684 / 0.8530) * (0.8170 / 0.9621) +
606
607 (eta >= -2.47 && eta < -2.37 ) * (pt > 60.0 && pt <= 80.0) * (0.6487 / 0.8310) * (0.8170 / 0.9621) +
608 (eta >= -2.37 && eta < -2.01 ) * (pt > 60.0 && pt <= 80.0) * (0.7237 / 0.9551) * (0.8170 / 0.9621) +
609 (eta >= -2.01 && eta < -1.81 ) * (pt > 60.0 && pt <= 80.0) * (0.7402 / 0.9528) * (0.8170 / 0.9621) +
610 (eta >= -1.81 && eta < -1.52 ) * (pt > 60.0 && pt <= 80.0) * (0.7644 / 0.9339) * (0.8170 / 0.9621) +
611 (eta >= -1.52 && eta < -1.37 ) * (pt > 60.0 && pt <= 80.0) * (0.6303 / 0.9126) * (0.8170 / 0.9621) +
612 (eta >= -1.37 && eta < -1.15 ) * (pt > 60.0 && pt <= 80.0) * (0.7526 / 0.9222) * (0.8170 / 0.9621) +
613 (eta >= -1.15 && eta < -0.80 ) * (pt > 60.0 && pt <= 80.0) * (0.7803 / 0.9558) * (0.8170 / 0.9621) +
614 (eta >= -0.80 && eta < -0.60 ) * (pt > 60.0 && pt <= 80.0) * (0.8053 / 0.9749) * (0.8170 / 0.9621) +
615 (eta >= -0.60 && eta < -0.10 ) * (pt > 60.0 && pt <= 80.0) * (0.8132 / 0.9871) * (0.8170 / 0.9621) +
616 (eta >= -0.10 && eta < 0.00 ) * (pt > 60.0 && pt <= 80.0) * (0.7618 / 0.9753) * (0.8170 / 0.9621) +
617 (eta >= 0.00 && eta < 0.10 ) * (pt > 60.0 && pt <= 80.0) * (0.7421 / 0.9635) * (0.8170 / 0.9621) +
618 (eta >= 0.10 && eta < 0.60 ) * (pt > 60.0 && pt <= 80.0) * (0.8145 / 0.9852) * (0.8170 / 0.9621) +
619 (eta >= 0.60 && eta < 0.80 ) * (pt > 60.0 && pt <= 80.0) * (0.8105 / 0.9782) * (0.8170 / 0.9621) +
620 (eta >= 0.80 && eta < 1.15 ) * (pt > 60.0 && pt <= 80.0) * (0.7842 / 0.9617) * (0.8170 / 0.9621) +
621 (eta >= 1.15 && eta < 1.37 ) * (pt > 60.0 && pt <= 80.0) * (0.7579 / 0.9309) * (0.8170 / 0.9621) +
622 (eta >= 1.37 && eta < 1.52 ) * (pt > 60.0 && pt <= 80.0) * (0.6461 / 0.9286) * (0.8170 / 0.9621) +
623 (eta >= 1.52 && eta < 1.81 ) * (pt > 60.0 && pt <= 80.0) * (0.7724 / 0.9454) * (0.8170 / 0.9621) +
624 (eta >= 1.81 && eta < 2.01 ) * (pt > 60.0 && pt <= 80.0) * (0.7439 / 0.9598) * (0.8170 / 0.9621) +
625 (eta >= 2.01 && eta < 2.37 ) * (pt > 60.0 && pt <= 80.0) * (0.7263 / 0.9577) * (0.8170 / 0.9621) +
626 (eta >= 2.37 && eta < 2.47 ) * (pt > 60.0 && pt <= 80.0) * (0.6684 / 0.8530) * (0.8170 / 0.9621) +
627
628 (eta >= -2.47 && eta < -2.37 ) * (pt > 80.0) * (0.6487 / 0.8310) * (0.8170 / 0.9621) +
629 (eta >= -2.37 && eta < -2.01 ) * (pt > 80.0) * (0.7237 / 0.9551) * (0.8170 / 0.9621) +
630 (eta >= -2.01 && eta < -1.81 ) * (pt > 80.0) * (0.7402 / 0.9528) * (0.8170 / 0.9621) +
631 (eta >= -1.81 && eta < -1.52 ) * (pt > 80.0) * (0.7644 / 0.9339) * (0.8170 / 0.9621) +
632 (eta >= -1.52 && eta < -1.37 ) * (pt > 80.0) * (0.6303 / 0.9126) * (0.8170 / 0.9621) +
633 (eta >= -1.37 && eta < -1.15 ) * (pt > 80.0) * (0.7526 / 0.9222) * (0.8170 / 0.9621) +
634 (eta >= -1.15 && eta < -0.80 ) * (pt > 80.0) * (0.7803 / 0.9558) * (0.8170 / 0.9621) +
635 (eta >= -0.80 && eta < -0.60 ) * (pt > 80.0) * (0.8053 / 0.9749) * (0.8170 / 0.9621) +
636 (eta >= -0.60 && eta < -0.10 ) * (pt > 80.0) * (0.8132 / 0.9871) * (0.8170 / 0.9621) +
637 (eta >= -0.10 && eta < 0.00 ) * (pt > 80.0) * (0.7618 / 0.9753) * (0.8170 / 0.9621) +
638 (eta >= 0.00 && eta < 0.10 ) * (pt > 80.0) * (0.7421 / 0.9635) * (0.8170 / 0.9621) +
639 (eta >= 0.10 && eta < 0.60 ) * (pt > 80.0) * (0.8145 / 0.9852) * (0.8170 / 0.9621) +
640 (eta >= 0.60 && eta < 0.80 ) * (pt > 80.0) * (0.8105 / 0.9782) * (0.8170 / 0.9621) +
641 (eta >= 0.80 && eta < 1.15 ) * (pt > 80.0) * (0.7842 / 0.9617) * (0.8170 / 0.9621) +
642 (eta >= 1.15 && eta < 1.37 ) * (pt > 80.0) * (0.7579 / 0.9309) * (0.8170 / 0.9621) +
643 (eta >= 1.37 && eta < 1.52 ) * (pt > 80.0) * (0.6461 / 0.9286) * (0.8170 / 0.9621) +
644 (eta >= 1.52 && eta < 1.81 ) * (pt > 80.0) * (0.7724 / 0.9454) * (0.8170 / 0.9621) +
645 (eta >= 1.81 && eta < 2.01 ) * (pt > 80.0) * (0.7439 / 0.9598) * (0.8170 / 0.9621) +
646 (eta >= 2.01 && eta < 2.37 ) * (pt > 80.0) * (0.7263 / 0.9577) * (0.8170 / 0.9621) +
647 (eta >= 2.37 && eta < 2.47 ) * (pt > 80.0) * (0.6684 / 0.8530) * (0.8170 / 0.9621) +
648
649 (abs(eta) > 2.47 ) * (0.00)}
650}
651
652####################
653# Electron isolation
654####################
655
656module Isolation ElectronIsolation {
657# set CandidateInputArray ElectronEfficiency/electrons
658# set IsolationInputArray EFlowFilter/eflow
659
660# set OutputArray electrons
661
662# set DeltaRMax 0.5
663
664# set PTMin 0.5
665
666# set PTRatioMax 0.12
667}
668
669#################
670# Muon efficiency
671#################
672
673module Efficiency MuonEfficiency {
674 set InputArray MuonMomentumSmearing/muons
675 set OutputArray muons
676
677 # set EfficiencyFormula {efficiency as a function of eta and pt}
678
679 # efficiency formula for muons
680 # efficiency formula for muons
681 #set EfficiencyFormula { (pt <= 10.0) * (0.00) +
682 # (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
683 # (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) +
684 # (abs(eta) > 2.7) * (0.00)}
685
686 # ### Reco & highpT-ID eff (* pT eff ratio) * track-based Iso eff
687 set EfficiencyFormula { (pt <= 4.0) * (0.0000) +
688
689 (eta >= -2.5 && eta < -2.4 ) * (pt > 4.0 && pt < 20.0) * (0.9835 * 0.781) +
690 (eta >= -2.4 && eta < -2.3 ) * (pt > 4.0 && pt < 20.0) * (0.9835 * 0.781) +
691 (eta >= -2.3 && eta < -2.2 ) * (pt > 4.0 && pt < 20.0) * (0.9236 * 0.781) +
692 (eta >= -2.2 && eta < -2.1 ) * (pt > 4.0 && pt < 20.0) * (0.9578 * 0.781) +
693 (eta >= -2.1 && eta < -2.0 ) * (pt > 4.0 && pt < 20.0) * (0.8979 * 0.781) +
694 (eta >= -2.0 && eta < -1.9 ) * (pt > 4.0 && pt < 20.0) * (0.9214 * 0.781) +
695 (eta >= -1.9 && eta < -1.8 ) * (pt > 4.0 && pt < 20.0) * (0.9814 * 0.781) +
696 (eta >= -1.8 && eta < -1.7 ) * (pt > 4.0 && pt < 20.0) * (0.9792 * 0.781) +
697 (eta >= -1.7 && eta < -1.6 ) * (pt > 4.0 && pt < 20.0) * (0.9835 * 0.781) +
698 (eta >= -1.6 && eta < -1.5 ) * (pt > 4.0 && pt < 20.0) * (0.9814 * 0.781) +
699 (eta >= -1.5 && eta < -1.4 ) * (pt > 4.0 && pt < 20.0) * (0.5554 * 0.781) +
700 (eta >= -1.4 && eta < -1.3 ) * (pt > 4.0 && pt < 20.0) * (0.4847 * 0.781) +
701 (eta >= -1.3 && eta < -1.2 ) * (pt > 4.0 && pt < 20.0) * (0.6945 * 0.781) +
702 (eta >= -1.2 && eta < -1.1 ) * (pt > 4.0 && pt < 20.0) * (0.8722 * 0.781) +
703 (eta >= -1.1 && eta < -1.0 ) * (pt > 4.0 && pt < 20.0) * (0.9193 * 0.781) +
704 (eta >= -1.0 && eta < -0.9 ) * (pt > 4.0 && pt < 20.0) * (0.8850 * 0.781) +
705 (eta >= -0.9 && eta < -0.8 ) * (pt > 4.0 && pt < 20.0) * (0.9856 * 0.781) +
706 (eta >= -0.8 && eta < -0.7 ) * (pt > 4.0 && pt < 20.0) * (0.7523 * 0.781) +
707 (eta >= -0.7 && eta < -0.6 ) * (pt > 4.0 && pt < 20.0) * (0.9128 * 0.781) +
708 (eta >= -0.6 && eta < -0.5 ) * (pt > 4.0 && pt < 20.0) * (0.9450 * 0.781) +
709 (eta >= -0.5 && eta < -0.4 ) * (pt > 4.0 && pt < 20.0) * (0.8208 * 0.781) +
710 (eta >= -0.4 && eta < -0.3 ) * (pt > 4.0 && pt < 20.0) * (0.9171 * 0.781) +
711 (eta >= -0.3 && eta < -0.2 ) * (pt > 4.0 && pt < 20.0) * (0.9428 * 0.781) +
712 (eta >= -0.2 && eta < -0.1 ) * (pt > 4.0 && pt < 20.0) * (0.9128 * 0.781) +
713 (eta >= -0.1 && eta < 0.0 ) * (pt > 4.0 && pt < 20.0) * (0.3434 * 0.781) +
714 (eta >= 0.0 && eta < 0.1 ) * (pt > 4.0 && pt < 20.0) * (0.3713 * 0.781) +
715 (eta >= 0.1 && eta < 0.2 ) * (pt > 4.0 && pt < 20.0) * (0.9557 * 0.781) +
716 (eta >= 0.2 && eta < 0.3 ) * (pt > 4.0 && pt < 20.0) * (0.9428 * 0.781) +
717 (eta >= 0.3 && eta < 0.4 ) * (pt > 4.0 && pt < 20.0) * (0.9128 * 0.781) +
718 (eta >= 0.4 && eta < 0.5 ) * (pt > 4.0 && pt < 20.0) * (0.8208 * 0.781) +
719 (eta >= 0.5 && eta < 0.6 ) * (pt > 4.0 && pt < 20.0) * (0.9450 * 0.781) +
720 (eta >= 0.6 && eta < 0.7 ) * (pt > 4.0 && pt < 20.0) * (0.9128 * 0.781) +
721 (eta >= 0.7 && eta < 0.8 ) * (pt > 4.0 && pt < 20.0) * (0.7502 * 0.781) +
722 (eta >= 0.8 && eta < 0.9 ) * (pt > 4.0 && pt < 20.0) * (0.9856 * 0.781) +
723 (eta >= 0.9 && eta < 1.0 ) * (pt > 4.0 && pt < 20.0) * (0.8807 * 0.781) +
724 (eta >= 1.0 && eta < 1.1 ) * (pt > 4.0 && pt < 20.0) * (0.9257 * 0.781) +
725 (eta >= 1.1 && eta < 1.2 ) * (pt > 4.0 && pt < 20.0) * (0.8743 * 0.781) +
726 (eta >= 1.2 && eta < 1.3 ) * (pt > 4.0 && pt < 20.0) * (0.6924 * 0.781) +
727 (eta >= 1.3 && eta < 1.4 ) * (pt > 4.0 && pt < 20.0) * (0.4804 * 0.781) +
728 (eta >= 1.4 && eta < 1.5 ) * (pt > 4.0 && pt < 20.0) * (0.5489 * 0.781) +
729 (eta >= 1.5 && eta < 1.6 ) * (pt > 4.0 && pt < 20.0) * (0.9749 * 0.781) +
730 (eta >= 1.6 && eta < 1.7 ) * (pt > 4.0 && pt < 20.0) * (0.9835 * 0.781) +
731 (eta >= 1.7 && eta < 1.8 ) * (pt > 4.0 && pt < 20.0) * (0.9749 * 0.781) +
732 (eta >= 1.8 && eta < 1.9 ) * (pt > 4.0 && pt < 20.0) * (0.9814 * 0.781) +
733 (eta >= 1.9 && eta < 2.0 ) * (pt > 4.0 && pt < 20.0) * (0.9193 * 0.781) +
734 (eta >= 2.0 && eta < 2.1 ) * (pt > 4.0 && pt < 20.0) * (0.8979 * 0.781) +
735 (eta >= 2.1 && eta < 2.2 ) * (pt > 4.0 && pt < 20.0) * (0.9557 * 0.781) +
736 (eta >= 2.2 && eta < 2.3 ) * (pt > 4.0 && pt < 20.0) * (0.9450 * 0.781) +
737 (eta >= 2.3 && eta < 2.4 ) * (pt > 4.0 && pt < 20.0) * (0.9814 * 0.781) +
738 (eta >= 2.4 && eta < 2.5 ) * (pt > 4.0 && pt < 20.0) * (0.9814 * 0.781) +
739
740 (eta >= -2.5 && eta < -2.4 ) * (pt >= 20.0 && pt < 55.0) * (0.9835 * 0.804) +
741 (eta >= -2.4 && eta < -2.3 ) * (pt >= 20.0 && pt < 55.0) * (0.9835 * 0.804) +
742 (eta >= -2.3 && eta < -2.2 ) * (pt >= 20.0 && pt < 55.0) * (0.9236 * 0.804) +
743 (eta >= -2.2 && eta < -2.1 ) * (pt >= 20.0 && pt < 55.0) * (0.9578 * 0.804) +
744 (eta >= -2.1 && eta < -2.0 ) * (pt >= 20.0 && pt < 55.0) * (0.8979 * 0.804) +
745 (eta >= -2.0 && eta < -1.9 ) * (pt >= 20.0 && pt < 55.0) * (0.9214 * 0.804) +
746 (eta >= -1.9 && eta < -1.8 ) * (pt >= 20.0 && pt < 55.0) * (0.9814 * 0.804) +
747 (eta >= -1.8 && eta < -1.7 ) * (pt >= 20.0 && pt < 55.0) * (0.9792 * 0.804) +
748 (eta >= -1.7 && eta < -1.6 ) * (pt >= 20.0 && pt < 55.0) * (0.9835 * 0.804) +
749 (eta >= -1.6 && eta < -1.5 ) * (pt >= 20.0 && pt < 55.0) * (0.9814 * 0.804) +
750 (eta >= -1.5 && eta < -1.4 ) * (pt >= 20.0 && pt < 55.0) * (0.5554 * 0.804) +
751 (eta >= -1.4 && eta < -1.3 ) * (pt >= 20.0 && pt < 55.0) * (0.4847 * 0.804) +
752 (eta >= -1.3 && eta < -1.2 ) * (pt >= 20.0 && pt < 55.0) * (0.6945 * 0.804) +
753 (eta >= -1.2 && eta < -1.1 ) * (pt >= 20.0 && pt < 55.0) * (0.8722 * 0.804) +
754 (eta >= -1.1 && eta < -1.0 ) * (pt >= 20.0 && pt < 55.0) * (0.9193 * 0.804) +
755 (eta >= -1.0 && eta < -0.9 ) * (pt >= 20.0 && pt < 55.0) * (0.8850 * 0.804) +
756 (eta >= -0.9 && eta < -0.8 ) * (pt >= 20.0 && pt < 55.0) * (0.9856 * 0.804) +
757 (eta >= -0.8 && eta < -0.7 ) * (pt >= 20.0 && pt < 55.0) * (0.7523 * 0.804) +
758 (eta >= -0.7 && eta < -0.6 ) * (pt >= 20.0 && pt < 55.0) * (0.9128 * 0.804) +
759 (eta >= -0.6 && eta < -0.5 ) * (pt >= 20.0 && pt < 55.0) * (0.9450 * 0.804) +
760 (eta >= -0.5 && eta < -0.4 ) * (pt >= 20.0 && pt < 55.0) * (0.8208 * 0.804) +
761 (eta >= -0.4 && eta < -0.3 ) * (pt >= 20.0 && pt < 55.0) * (0.9171 * 0.804) +
762 (eta >= -0.3 && eta < -0.2 ) * (pt >= 20.0 && pt < 55.0) * (0.9428 * 0.804) +
763 (eta >= -0.2 && eta < -0.1 ) * (pt >= 20.0 && pt < 55.0) * (0.9128 * 0.804) +
764 (eta >= -0.1 && eta < 0.0 ) * (pt >= 20.0 && pt < 55.0) * (0.3434 * 0.804) +
765 (eta >= 0.0 && eta < 0.1 ) * (pt >= 20.0 && pt < 55.0) * (0.3713 * 0.804) +
766 (eta >= 0.1 && eta < 0.2 ) * (pt >= 20.0 && pt < 55.0) * (0.9557 * 0.804) +
767 (eta >= 0.2 && eta < 0.3 ) * (pt >= 20.0 && pt < 55.0) * (0.9428 * 0.804) +
768 (eta >= 0.3 && eta < 0.4 ) * (pt >= 20.0 && pt < 55.0) * (0.9128 * 0.804) +
769 (eta >= 0.4 && eta < 0.5 ) * (pt >= 20.0 && pt < 55.0) * (0.8208 * 0.804) +
770 (eta >= 0.5 && eta < 0.6 ) * (pt >= 20.0 && pt < 55.0) * (0.9450 * 0.804) +
771 (eta >= 0.6 && eta < 0.7 ) * (pt >= 20.0 && pt < 55.0) * (0.9128 * 0.804) +
772 (eta >= 0.7 && eta < 0.8 ) * (pt >= 20.0 && pt < 55.0) * (0.7502 * 0.804) +
773 (eta >= 0.8 && eta < 0.9 ) * (pt >= 20.0 && pt < 55.0) * (0.9856 * 0.804) +
774 (eta >= 0.9 && eta < 1.0 ) * (pt >= 20.0 && pt < 55.0) * (0.8807 * 0.804) +
775 (eta >= 1.0 && eta < 1.1 ) * (pt >= 20.0 && pt < 55.0) * (0.9257 * 0.804) +
776 (eta >= 1.1 && eta < 1.2 ) * (pt >= 20.0 && pt < 55.0) * (0.8743 * 0.804) +
777 (eta >= 1.2 && eta < 1.3 ) * (pt >= 20.0 && pt < 55.0) * (0.6924 * 0.804) +
778 (eta >= 1.3 && eta < 1.4 ) * (pt >= 20.0 && pt < 55.0) * (0.4804 * 0.804) +
779 (eta >= 1.4 && eta < 1.5 ) * (pt >= 20.0 && pt < 55.0) * (0.5489 * 0.804) +
780 (eta >= 1.5 && eta < 1.6 ) * (pt >= 20.0 && pt < 55.0) * (0.9749 * 0.804) +
781 (eta >= 1.6 && eta < 1.7 ) * (pt >= 20.0 && pt < 55.0) * (0.9835 * 0.804) +
782 (eta >= 1.7 && eta < 1.8 ) * (pt >= 20.0 && pt < 55.0) * (0.9749 * 0.804) +
783 (eta >= 1.8 && eta < 1.9 ) * (pt >= 20.0 && pt < 55.0) * (0.9814 * 0.804) +
784 (eta >= 1.9 && eta < 2.0 ) * (pt >= 20.0 && pt < 55.0) * (0.9193 * 0.804) +
785 (eta >= 2.0 && eta < 2.1 ) * (pt >= 20.0 && pt < 55.0) * (0.8979 * 0.804) +
786 (eta >= 2.1 && eta < 2.2 ) * (pt >= 20.0 && pt < 55.0) * (0.9557 * 0.804) +
787 (eta >= 2.2 && eta < 2.3 ) * (pt >= 20.0 && pt < 55.0) * (0.9450 * 0.804) +
788 (eta >= 2.3 && eta < 2.4 ) * (pt >= 20.0 && pt < 55.0) * (0.9814 * 0.804) +
789 (eta >= 2.4 && eta < 2.5 ) * (pt >= 20.0 && pt < 55.0) * (0.9814 * 0.804) +
790
791 (eta >= -2.5 && eta < -2.4 ) * (pt >= 55.0 && pt < 3000.0) * (0.9835) * (0.77 - pt * 0.00008) +
792 (eta >= -2.4 && eta < -2.3 ) * (pt >= 55.0 && pt < 3000.0) * (0.9835) * (0.77 - pt * 0.00008) +
793 (eta >= -2.3 && eta < -2.2 ) * (pt >= 55.0 && pt < 3000.0) * (0.9236) * (0.77 - pt * 0.00008) +
794 (eta >= -2.2 && eta < -2.1 ) * (pt >= 55.0 && pt < 3000.0) * (0.9578) * (0.77 - pt * 0.00008) +
795 (eta >= -2.1 && eta < -2.0 ) * (pt >= 55.0 && pt < 3000.0) * (0.8979) * (0.77 - pt * 0.00008) +
796 (eta >= -2.0 && eta < -1.9 ) * (pt >= 55.0 && pt < 3000.0) * (0.9214) * (0.77 - pt * 0.00008) +
797 (eta >= -1.9 && eta < -1.8 ) * (pt >= 55.0 && pt < 3000.0) * (0.9814) * (0.77 - pt * 0.00008) +
798 (eta >= -1.8 && eta < -1.7 ) * (pt >= 55.0 && pt < 3000.0) * (0.9792) * (0.77 - pt * 0.00008) +
799 (eta >= -1.7 && eta < -1.6 ) * (pt >= 55.0 && pt < 3000.0) * (0.9835) * (0.77 - pt * 0.00008) +
800 (eta >= -1.6 && eta < -1.5 ) * (pt >= 55.0 && pt < 3000.0) * (0.9814) * (0.77 - pt * 0.00008) +
801 (eta >= -1.5 && eta < -1.4 ) * (pt >= 55.0 && pt < 3000.0) * (0.5554) * (0.77 - pt * 0.00008) +
802 (eta >= -1.4 && eta < -1.3 ) * (pt >= 55.0 && pt < 3000.0) * (0.4847) * (0.77 - pt * 0.00008) +
803 (eta >= -1.3 && eta < -1.2 ) * (pt >= 55.0 && pt < 3000.0) * (0.6945) * (0.77 - pt * 0.00008) +
804 (eta >= -1.2 && eta < -1.1 ) * (pt >= 55.0 && pt < 3000.0) * (0.8722) * (0.77 - pt * 0.00008) +
805 (eta >= -1.1 && eta < -1.0 ) * (pt >= 55.0 && pt < 3000.0) * (0.9193) * (0.77 - pt * 0.00008) +
806 (eta >= -1.0 && eta < -0.9 ) * (pt >= 55.0 && pt < 3000.0) * (0.8850) * (0.77 - pt * 0.00008) +
807 (eta >= -0.9 && eta < -0.8 ) * (pt >= 55.0 && pt < 3000.0) * (0.9856) * (0.77 - pt * 0.00008) +
808 (eta >= -0.8 && eta < -0.7 ) * (pt >= 55.0 && pt < 3000.0) * (0.7523) * (0.77 - pt * 0.00008) +
809 (eta >= -0.7 && eta < -0.6 ) * (pt >= 55.0 && pt < 3000.0) * (0.9128) * (0.77 - pt * 0.00008) +
810 (eta >= -0.6 && eta < -0.5 ) * (pt >= 55.0 && pt < 3000.0) * (0.9450) * (0.77 - pt * 0.00008) +
811 (eta >= -0.5 && eta < -0.4 ) * (pt >= 55.0 && pt < 3000.0) * (0.8208) * (0.77 - pt * 0.00008) +
812 (eta >= -0.4 && eta < -0.3 ) * (pt >= 55.0 && pt < 3000.0) * (0.9171) * (0.77 - pt * 0.00008) +
813 (eta >= -0.3 && eta < -0.2 ) * (pt >= 55.0 && pt < 3000.0) * (0.9428) * (0.77 - pt * 0.00008) +
814 (eta >= -0.2 && eta < -0.1 ) * (pt >= 55.0 && pt < 3000.0) * (0.9128) * (0.77 - pt * 0.00008) +
815 (eta >= -0.1 && eta < 0.0 ) * (pt >= 55.0 && pt < 3000.0) * (0.3434) * (0.77 - pt * 0.00008) +
816 (eta >= 0.0 && eta < 0.1 ) * (pt >= 55.0 && pt < 3000.0) * (0.3713) * (0.77 - pt * 0.00008) +
817 (eta >= 0.1 && eta < 0.2 ) * (pt >= 55.0 && pt < 3000.0) * (0.9557) * (0.77 - pt * 0.00008) +
818 (eta >= 0.2 && eta < 0.3 ) * (pt >= 55.0 && pt < 3000.0) * (0.9428) * (0.77 - pt * 0.00008) +
819 (eta >= 0.3 && eta < 0.4 ) * (pt >= 55.0 && pt < 3000.0) * (0.9128) * (0.77 - pt * 0.00008) +
820 (eta >= 0.4 && eta < 0.5 ) * (pt >= 55.0 && pt < 3000.0) * (0.8208) * (0.77 - pt * 0.00008) +
821 (eta >= 0.5 && eta < 0.6 ) * (pt >= 55.0 && pt < 3000.0) * (0.9450) * (0.77 - pt * 0.00008) +
822 (eta >= 0.6 && eta < 0.7 ) * (pt >= 55.0 && pt < 3000.0) * (0.9128) * (0.77 - pt * 0.00008) +
823 (eta >= 0.7 && eta < 0.8 ) * (pt >= 55.0 && pt < 3000.0) * (0.7502) * (0.77 - pt * 0.00008) +
824 (eta >= 0.8 && eta < 0.9 ) * (pt >= 55.0 && pt < 3000.0) * (0.9856) * (0.77 - pt * 0.00008) +
825 (eta >= 0.9 && eta < 1.0 ) * (pt >= 55.0 && pt < 3000.0) * (0.8807) * (0.77 - pt * 0.00008) +
826 (eta >= 1.0 && eta < 1.1 ) * (pt >= 55.0 && pt < 3000.0) * (0.9257) * (0.77 - pt * 0.00008) +
827 (eta >= 1.1 && eta < 1.2 ) * (pt >= 55.0 && pt < 3000.0) * (0.8743) * (0.77 - pt * 0.00008) +
828 (eta >= 1.2 && eta < 1.3 ) * (pt >= 55.0 && pt < 3000.0) * (0.6924) * (0.77 - pt * 0.00008) +
829 (eta >= 1.3 && eta < 1.4 ) * (pt >= 55.0 && pt < 3000.0) * (0.4804) * (0.77 - pt * 0.00008) +
830 (eta >= 1.4 && eta < 1.5 ) * (pt >= 55.0 && pt < 3000.0) * (0.5489) * (0.77 - pt * 0.00008) +
831 (eta >= 1.5 && eta < 1.6 ) * (pt >= 55.0 && pt < 3000.0) * (0.9749) * (0.77 - pt * 0.00008) +
832 (eta >= 1.6 && eta < 1.7 ) * (pt >= 55.0 && pt < 3000.0) * (0.9835) * (0.77 - pt * 0.00008) +
833 (eta >= 1.7 && eta < 1.8 ) * (pt >= 55.0 && pt < 3000.0) * (0.9749) * (0.77 - pt * 0.00008) +
834 (eta >= 1.8 && eta < 1.9 ) * (pt >= 55.0 && pt < 3000.0) * (0.9814) * (0.77 - pt * 0.00008) +
835 (eta >= 1.9 && eta < 2.0 ) * (pt >= 55.0 && pt < 3000.0) * (0.9193) * (0.77 - pt * 0.00008) +
836 (eta >= 2.0 && eta < 2.1 ) * (pt >= 55.0 && pt < 3000.0) * (0.8979) * (0.77 - pt * 0.00008) +
837 (eta >= 2.1 && eta < 2.2 ) * (pt >= 55.0 && pt < 3000.0) * (0.9557) * (0.77 - pt * 0.00008) +
838 (eta >= 2.2 && eta < 2.3 ) * (pt >= 55.0 && pt < 3000.0) * (0.9450) * (0.77 - pt * 0.00008) +
839 (eta >= 2.3 && eta < 2.4 ) * (pt >= 55.0 && pt < 3000.0) * (0.9814) * (0.77 - pt * 0.00008) +
840 (eta >= 2.4 && eta < 2.5 ) * (pt >= 55.0 && pt < 3000.0) * (0.9814) * (0.77 - pt * 0.00008) +
841
842 (eta >= -2.5 && eta < -2.4 ) * (pt >= 3000.0) * (0.9835 * 0.53) +
843 (eta >= -2.4 && eta < -2.3 ) * (pt >= 3000.0) * (0.9835 * 0.53) +
844 (eta >= -2.3 && eta < -2.2 ) * (pt >= 3000.0) * (0.9236 * 0.53) +
845 (eta >= -2.2 && eta < -2.1 ) * (pt >= 3000.0) * (0.9578 * 0.53) +
846 (eta >= -2.1 && eta < -2.0 ) * (pt >= 3000.0) * (0.8979 * 0.53) +
847 (eta >= -2.0 && eta < -1.9 ) * (pt >= 3000.0) * (0.9214 * 0.53) +
848 (eta >= -1.9 && eta < -1.8 ) * (pt >= 3000.0) * (0.9814 * 0.53) +
849 (eta >= -1.8 && eta < -1.7 ) * (pt >= 3000.0) * (0.9792 * 0.53) +
850 (eta >= -1.7 && eta < -1.6 ) * (pt >= 3000.0) * (0.9835 * 0.53) +
851 (eta >= -1.6 && eta < -1.5 ) * (pt >= 3000.0) * (0.9814 * 0.53) +
852 (eta >= -1.5 && eta < -1.4 ) * (pt >= 3000.0) * (0.5554 * 0.53) +
853 (eta >= -1.4 && eta < -1.3 ) * (pt >= 3000.0) * (0.4847 * 0.53) +
854 (eta >= -1.3 && eta < -1.2 ) * (pt >= 3000.0) * (0.6945 * 0.53) +
855 (eta >= -1.2 && eta < -1.1 ) * (pt >= 3000.0) * (0.8722 * 0.53) +
856 (eta >= -1.1 && eta < -1.0 ) * (pt >= 3000.0) * (0.9193 * 0.53) +
857 (eta >= -1.0 && eta < -0.9 ) * (pt >= 3000.0) * (0.8850 * 0.53) +
858 (eta >= -0.9 && eta < -0.8 ) * (pt >= 3000.0) * (0.9856 * 0.53) +
859 (eta >= -0.8 && eta < -0.7 ) * (pt >= 3000.0) * (0.7523 * 0.53) +
860 (eta >= -0.7 && eta < -0.6 ) * (pt >= 3000.0) * (0.9128 * 0.53) +
861 (eta >= -0.6 && eta < -0.5 ) * (pt >= 3000.0) * (0.9450 * 0.53) +
862 (eta >= -0.5 && eta < -0.4 ) * (pt >= 3000.0) * (0.8208 * 0.53) +
863 (eta >= -0.4 && eta < -0.3 ) * (pt >= 3000.0) * (0.9171 * 0.53) +
864 (eta >= -0.3 && eta < -0.2 ) * (pt >= 3000.0) * (0.9428 * 0.53) +
865 (eta >= -0.2 && eta < -0.1 ) * (pt >= 3000.0) * (0.9128 * 0.53) +
866 (eta >= -0.1 && eta < 0.0 ) * (pt >= 3000.0) * (0.3434 * 0.53) +
867 (eta >= 0.0 && eta < 0.1 ) * (pt >= 3000.0) * (0.3713 * 0.53) +
868 (eta >= 0.1 && eta < 0.2 ) * (pt >= 3000.0) * (0.9557 * 0.53) +
869 (eta >= 0.2 && eta < 0.3 ) * (pt >= 3000.0) * (0.9428 * 0.53) +
870 (eta >= 0.3 && eta < 0.4 ) * (pt >= 3000.0) * (0.9128 * 0.53) +
871 (eta >= 0.4 && eta < 0.5 ) * (pt >= 3000.0) * (0.8208 * 0.53) +
872 (eta >= 0.5 && eta < 0.6 ) * (pt >= 3000.0) * (0.9450 * 0.53) +
873 (eta >= 0.6 && eta < 0.7 ) * (pt >= 3000.0) * (0.9128 * 0.53) +
874 (eta >= 0.7 && eta < 0.8 ) * (pt >= 3000.0) * (0.7502 * 0.53) +
875 (eta >= 0.8 && eta < 0.9 ) * (pt >= 3000.0) * (0.9856 * 0.53) +
876 (eta >= 0.9 && eta < 1.0 ) * (pt >= 3000.0) * (0.8807 * 0.53) +
877 (eta >= 1.0 && eta < 1.1 ) * (pt >= 3000.0) * (0.9257 * 0.53) +
878 (eta >= 1.1 && eta < 1.2 ) * (pt >= 3000.0) * (0.8743 * 0.53) +
879 (eta >= 1.2 && eta < 1.3 ) * (pt >= 3000.0) * (0.6924 * 0.53) +
880 (eta >= 1.3 && eta < 1.4 ) * (pt >= 3000.0) * (0.4804 * 0.53) +
881 (eta >= 1.4 && eta < 1.5 ) * (pt >= 3000.0) * (0.5489 * 0.53) +
882 (eta >= 1.5 && eta < 1.6 ) * (pt >= 3000.0) * (0.9749 * 0.53) +
883 (eta >= 1.6 && eta < 1.7 ) * (pt >= 3000.0) * (0.9835 * 0.53) +
884 (eta >= 1.7 && eta < 1.8 ) * (pt >= 3000.0) * (0.9749 * 0.53) +
885 (eta >= 1.8 && eta < 1.9 ) * (pt >= 3000.0) * (0.9814 * 0.53) +
886 (eta >= 1.9 && eta < 2.0 ) * (pt >= 3000.0) * (0.9193 * 0.53) +
887 (eta >= 2.0 && eta < 2.1 ) * (pt >= 3000.0) * (0.8979 * 0.53) +
888 (eta >= 2.1 && eta < 2.2 ) * (pt >= 3000.0) * (0.9557 * 0.53) +
889 (eta >= 2.2 && eta < 2.3 ) * (pt >= 3000.0) * (0.9450 * 0.53) +
890 (eta >= 2.3 && eta < 2.4 ) * (pt >= 3000.0) * (0.9814 * 0.53) +
891 (eta >= 2.4 && eta < 2.5 ) * (pt >= 3000.0) * (0.9814 * 0.53) +
892
893 (abs(eta) > 2.5) * (0.00)}
894}
895
896
897################
898# Muon isolation
899################
900
901module Isolation MuonIsolation {
902# set CandidateInputArray MuonEfficiency/muons
903# set IsolationInputArray EFlowFilter/eflow
904
905# set OutputArray muons
906
907# set DeltaRMax 0.5
908
909# set PTMin 0.5
910
911# set PTRatioMax 0.25
912}
913
914###################
915# Missing ET merger
916###################
917
918module Merger MissingET {
919# add InputArray InputArray
920 add InputArray Calorimeter/towers
921 set MomentumOutputArray momentum
922}
923
924##################
925# Scalar HT merger
926##################
927
928module Merger ScalarHT {
929# add InputArray InputArray
930 add InputArray JetEnergyScale/jets
931 add InputArray ElectronEfficiency/electrons
932 add InputArray PhotonEfficiency/photons
933 add InputArray MuonEfficiency/muons
934 set EnergyOutputArray energy
935}
936
937
938#####################
939# Neutrino Filter
940#####################
941
942module PdgCodeFilter NeutrinoFilter {
943
944 set InputArray Delphes/stableParticles
945 set OutputArray filteredParticles
946
947 set PTMin 0.0
948
949 add PdgCode {12}
950 add PdgCode {14}
951 add PdgCode {16}
952 add PdgCode {-12}
953 add PdgCode {-14}
954 add PdgCode {-16}
955
956}
957
958#####################
959# MC truth jet finder
960#####################
961
962module FastJetFinder GenJetFinder {
963 set InputArray NeutrinoFilter/filteredParticles
964
965 set OutputArray jets
966
967 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
968 set JetAlgorithm 6
969 set ParameterR 0.6
970
971 set JetPTMin 20.0
972}
973
974
975#########################
976# Gen Missing ET merger
977########################
978
979module Merger GenMissingET {
980# add InputArray InputArray
981 add InputArray NeutrinoFilter/filteredParticles
982 set MomentumOutputArray momentum
983}
984
985
986
987############
988# Jet finder
989############
990
991module FastJetFinder FastJetFinder {
992 set InputArray Calorimeter/towers
993
994 set OutputArray jets
995
996 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
997 set JetAlgorithm 6
998 set ParameterR 0.6
999
1000 set JetPTMin 20.0
1001}
1002
1003##################
1004# Jet Energy Scale
1005##################
1006
1007module EnergyScale JetEnergyScale {
1008 set InputArray FastJetFinder/jets
1009 set OutputArray jets
1010
1011 # scale formula for jets
1012 set ScaleFormula { sqrt( (3.0 - 0.2*(abs(eta)))^2 / pt + 1.0 ) }
1013}
1014
1015########################
1016# Jet Flavor Association
1017########################
1018
1019module JetFlavorAssociation JetFlavorAssociation {
1020
1021 set PartonInputArray Delphes/partons
1022 set ParticleInputArray Delphes/allParticles
1023 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1024 set JetInputArray JetEnergyScale/jets
1025
1026 set DeltaR 0.5
1027 set PartonPTMin 1.0
1028 set PartonEtaMax 2.5
1029
1030}
1031
1032###########
1033# b-tagging
1034###########
1035
1036module BTagging BTagging {
1037 set JetInputArray JetEnergyScale/jets
1038
1039 set BitNumber 0
1040
1041 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1042 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
1043 # gluon's PDG code has the lowest priority
1044
1045 # based on ATL-PHYS-PUB-2015-022
1046
1047 # default efficiency formula (misidentification rate)
1048 add EfficiencyFormula {0} {0.002+7.3e-06*pt}
1049
1050 # efficiency formula for c-jets (misidentification rate)
1051 add EfficiencyFormula {4} {0.20*tanh(0.02*pt)*(1/(1+0.0034*pt))}
1052
1053 # efficiency formula for b-jets
1054 add EfficiencyFormula {5} {0.80*tanh(0.003*pt)*(30/(1+0.086*pt))}
1055}
1056
1057#############
1058# tau-tagging
1059#############
1060
1061module TrackCountingTauTagging TauTagging {
1062
1063 set ParticleInputArray Delphes/allParticles
1064 set PartonInputArray Delphes/partons
1065 set TrackInputArray TrackMerger/tracks
1066 set JetInputArray JetEnergyScale/jets
1067
1068 set DeltaR 0.2
1069 set DeltaRTrack 0.2
1070
1071 set TrackPTMin 1.0
1072
1073 set TauPTMin 1.0
1074 set TauEtaMax 2.5
1075
1076 # instructions: {n-prongs} {eff}
1077
1078 # 1 - one prong efficiency
1079 # 2 - two or more efficiency
1080 # -1 - one prong mistag rate
1081 # -2 - two or more mistag rate
1082
1083 set BitNumber 0
1084
1085 # taken from ATL-PHYS-PUB-2015-045 (medium working point)
1086 add EfficiencyFormula {1} {0.70}
1087 add EfficiencyFormula {2} {0.60}
1088 add EfficiencyFormula {-1} {0.02}
1089 add EfficiencyFormula {-2} {0.01}
1090
1091}
1092
1093#####################################################
1094# Find uniquely identified photons/electrons/tau/jets
1095#####################################################
1096
1097module UniqueObjectFinder UniqueObjectFinder {
1098# earlier arrays take precedence over later ones
1099# add InputArray InputArray OutputArray
1100 add InputArray PhotonEfficiency/photons photons
1101 add InputArray ElectronEfficiency/electrons electrons
1102 add InputArray MuonEfficiency/muons muons
1103 add InputArray JetEnergyScale/jets jets
1104
1105}
1106
1107
1108##################
1109# ROOT tree writer
1110##################
1111
1112# tracks, towers and eflow objects are not stored by default in the output.
1113# if needed (for jet constituent or other studies), uncomment the relevant
1114# "add Branch ..." lines.
1115
1116module TreeWriter TreeWriter {
1117# add Branch InputArray BranchName BranchClass
1118 add Branch Delphes/allParticles Particle GenParticle
1119
1120 add Branch TrackMerger/tracks Track Track
1121 add Branch Calorimeter/towers Tower Tower
1122
1123 add Branch HCal/eflowTracks EFlowTrack Track
1124 add Branch ECal/eflowPhotons EFlowPhoton Tower
1125 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1126
1127 add Branch GenJetFinder/jets GenJet Jet
1128 add Branch GenMissingET/momentum GenMissingET MissingET
1129
1130 add Branch UniqueObjectFinder/jets Jet Jet
1131 add Branch UniqueObjectFinder/electrons Electron Electron
1132 add Branch UniqueObjectFinder/photons Photon Photon
1133 add Branch UniqueObjectFinder/muons Muon Muon
1134 add Branch MissingET/momentum MissingET MissingET
1135 add Branch ScalarHT/energy ScalarHT ScalarHT
1136}
1137