MA5SandBox: delphes_card_atlas_susy_2019_08.tcl

File delphes_card_atlas_susy_2019_08.tcl, 28.9 KB (added by Benjamin Fuks, 4 years ago)
Line 
1#set MaxEvents 10000
2
3#######################################
4# Order of execution of various modules
5#######################################
6
7set ExecutionPath {
8 ParticlePropagator
9
10 ChargedHadronTrackingEfficiency
11 ElectronTrackingEfficiency
12 MuonTrackingEfficiency
13
14 ChargedHadronMomentumSmearing
15 ElectronMomentumSmearing
16 MuonMomentumSmearing
17
18 TrackMerger
19
20 ECal
21 HCal
22
23 Calorimeter
24 EFlowMerger
25 EFlowFilter
26
27 PhotonEfficiency
28
29 ElectronFilter
30 ElectronEfficiency
31
32
33 ChargedHadronFilter
34
35 MuonEfficiency
36
37
38 MissingET
39
40 NeutrinoFilter
41 GenJetFinder
42 GenMissingET
43
44 FastJetFinder
45
46 JetEnergyScale
47
48 JetFlavorAssociation
49
50 BTagging
51 TauTagging
52
53
54 TreeWriter
55}
56
57#################################
58# Propagate particles in cylinder
59#################################
60
61module ParticlePropagator ParticlePropagator {
62 set InputArray Delphes/stableParticles
63
64 set OutputArray stableParticles
65 set ChargedHadronOutputArray chargedHadrons
66 set ElectronOutputArray electrons
67 set MuonOutputArray muons
68
69 # radius of the magnetic field coverage, in m
70 set Radius 1.15
71 # half-length of the magnetic field coverage, in m
72 set HalfLength 3.51
73
74 # magnetic field
75 set Bz 2.0
76}
77
78####################################
79# Charged hadron tracking efficiency
80####################################
81
82module Efficiency ChargedHadronTrackingEfficiency {
83 set InputArray ParticlePropagator/chargedHadrons
84 set OutputArray chargedHadrons
85
86 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
87
88 # tracking efficiency formula for charged hadrons
89 #set EfficiencyFormula { (pt <= 0.1) * (0.00) +
90 # (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (1.0) +
91 # (abs(eta) <= 1.5) * (pt > 1.0) * (1.0) +
92 # (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (1.0) +
93 # (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (1.0) +
94 # (abs(eta) > 2.5) * (0.00)}
95 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
96 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
97 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
98 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
99 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
100 (abs(eta) > 2.5) * (0.00)}
101
102}
103
104##############################
105# Electron tracking efficiency
106##############################
107
108module Efficiency ElectronTrackingEfficiency {
109 set InputArray ParticlePropagator/electrons
110 set OutputArray electrons
111
112 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
113
114 # tracking efficiency formula for electrons
115 # set EfficiencyFormula { (pt <= 0.1) * (0.00) +
116 # (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (1.0) +
117 # (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (1.0) +
118 # (abs(eta) <= 1.5) * (pt > 1.0e2) * (1.0) +
119 # (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (1.0) +
120 # (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (1.0) +
121 # (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (1.0) +
122 # (abs(eta) > 2.5) * (0.00)}
123
124 # 1908.00005 figure 5
125set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
126 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.6) + \
127 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \
128 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \
129 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + \
130 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \
131 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.95) + \
132 (abs(eta) > 2.5) * (0.00)}
133# set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
134# (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \
135# (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \
136# (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \
137# (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + \
138# (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + \
139# (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + \
140# (abs(eta) > 2.5) * (0.00)}
141
142}
143
144##########################
145# Muon tracking efficiency
146##########################
147
148module Efficiency MuonTrackingEfficiency {
149 set InputArray ParticlePropagator/muons
150 set OutputArray muons
151
152 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
153
154 # tracking efficiency formula for muons
155 # set EfficiencyFormula { (pt <= 0.1) * (0.00) +
156 # (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.99) +
157 # (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) +
158 # (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 0.1 && pt <= 1.0) * (0.99) +
159 # (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 1.0) * (0.99) +
160 # (abs(eta) > 2.7) * (0.00)}
161
162set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
163 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \
164 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \
165 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
166 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \
167 (abs(eta) > 2.5) * (0.00)}
168
169
170}
171
172########################################
173# Momentum resolution for charged tracks
174########################################
175
176module MomentumSmearing ChargedHadronMomentumSmearing {
177 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
178 set OutputArray chargedHadrons
179
180 # set ResolutionFormula {resolution formula as a function of eta and pt}
181
182 # resolution formula for charged hadrons
183 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
184 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
185 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
186}
187
188###################################
189# Momentum resolution for electrons
190###################################
191
192module MomentumSmearing ElectronMomentumSmearing {
193 set InputArray ElectronTrackingEfficiency/electrons
194 set OutputArray electrons
195
196 # set ResolutionFormula {resolution formula as a function of eta and energy}
197
198 # resolution formula for electrons
199 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
200 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
201 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)
202 }
203}
204
205###############################
206# Momentum resolution for muons
207###############################
208
209module MomentumSmearing MuonMomentumSmearing {
210 set InputArray MuonTrackingEfficiency/muons
211 set OutputArray muons
212
213 # set ResolutionFormula {resolution formula as a function of eta and pt}
214 # resolution formula for muons
215 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
216 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
217 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
218}
219
220##############
221# Track merger
222##############
223
224module Merger TrackMerger {
225# add InputArray InputArray
226 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
227 add InputArray ElectronMomentumSmearing/electrons
228 add InputArray MuonMomentumSmearing/muons
229 set OutputArray tracks
230}
231
232
233#############
234# ECAL
235#############
236
237module SimpleCalorimeter ECal {
238 set ParticleInputArray ParticlePropagator/stableParticles
239 set TrackInputArray TrackMerger/tracks
240
241 set TowerOutputArray ecalTowers
242 set EFlowTrackOutputArray eflowTracks
243 set EFlowTowerOutputArray eflowPhotons
244
245 set IsEcal true
246
247 set EnergyMin 0.5
248 #set EnergyMin 0.2
249 set EnergySignificanceMin 2.0
250
251# set EnergyMin 0.25
252# set EnergySignificanceMin 1.0
253
254
255
256
257 set SmearTowerCenter true
258
259 set pi [expr {acos(-1)}]
260
261 # lists of the edges of each tower in eta and phi
262 # each list starts with the lower edge of the first tower
263 # the list ends with the higher edged of the last tower
264
265 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
266
267 set PhiBins {}
268 for {set i -180} {$i <= 180} {incr i} {
269 add PhiBins [expr {$i * $pi/180.0}]
270 }
271
272 # 0.02 unit in eta up to eta = 1.5 (barrel)
273 for {set i -85} {$i <= 86} {incr i} {
274 set eta [expr {$i * 0.0174}]
275 add EtaPhiBins $eta $PhiBins
276 }
277
278 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0
279 set PhiBins {}
280 for {set i -180} {$i <= 180} {incr i} {
281 add PhiBins [expr {$i * $pi/180.0}]
282 }
283
284 # 0.02 unit in eta up to eta = 3
285 for {set i 1} {$i <= 84} {incr i} {
286 set eta [expr { -2.958 + $i * 0.0174}]
287 add EtaPhiBins $eta $PhiBins
288 }
289
290 for {set i 1} {$i <= 84} {incr i} {
291 set eta [expr { 1.4964 + $i * 0.0174}]
292 add EtaPhiBins $eta $PhiBins
293 }
294
295 # take present CMS granularity for HF
296
297 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
298 set PhiBins {}
299 for {set i -18} {$i <= 18} {incr i} {
300 add PhiBins [expr {$i * $pi/18.0}]
301 }
302
303 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} {
304 add EtaPhiBins $eta $PhiBins
305 }
306
307
308 add EnergyFraction {0} {0.0}
309 # energy fractions for e, gamma and pi0
310 add EnergyFraction {11} {1.0}
311 add EnergyFraction {22} {1.0}
312 add EnergyFraction {111} {1.0}
313 # energy fractions for muon, neutrinos and neutralinos
314 add EnergyFraction {12} {0.0}
315 add EnergyFraction {13} {0.0}
316 add EnergyFraction {14} {0.0}
317 add EnergyFraction {16} {0.0}
318 add EnergyFraction {1000022} {0.0}
319 add EnergyFraction {1000023} {0.0}
320 add EnergyFraction {1000025} {0.0}
321 add EnergyFraction {1000035} {0.0}
322 add EnergyFraction {1000045} {0.0}
323 # energy fractions for K0short and Lambda
324 add EnergyFraction {310} {0.3}
325 add EnergyFraction {3122} {0.3}
326
327 # set ResolutionFormula {resolution formula as a function of eta and energy}
328
329 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
330 # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
331 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
332 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
333 set ResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) +
334 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
335 # set ResolutionFormula {1.0}
336
337 # https://cds.cern.ch/record/2037117/files/ATL-CAL-SLIDE-2015-388.pdf
338 #
339 #set ResolutionFormula { 0.17 oplus 0.1*sqrt(energy) oplus 0.007 * energy}
340
341
342 #set ResolutionFormula {0.01}
343 #set ResolutionFormula { sqrt(0.17^2 + energy*0.1^2 + energy^2*0.007^2)}
344
345}
346
347
348
349#############
350# HCAL
351#############
352
353module SimpleCalorimeter HCal {
354 set ParticleInputArray ParticlePropagator/stableParticles
355 set TrackInputArray ECal/eflowTracks
356
357 set TowerOutputArray hcalTowers
358 set EFlowTrackOutputArray eflowTracks
359 set EFlowTowerOutputArray eflowNeutralHadrons
360
361 set IsEcal false
362
363# set EnergyMin 1.0
364# set EnergySignificanceMin 2.0
365 #set EnergyMin 0.2
366 #set EnergySignificanceMin 1.0
367 set EnergyMin 0.5
368 set EnergySignificanceMin 1.0
369 set SmearTowerCenter true
370
371 set pi [expr {acos(-1)}]
372
373 # lists of the edges of each tower in eta and phi
374 # each list starts with the lower edge of the first tower
375 # the list ends with the higher edged of the last tower
376
377 # 10 degrees towers
378 set PhiBins {}
379 for {set i -18} {$i <= 18} {incr i} {
380 add PhiBins [expr {$i * $pi/18.0}]
381 }
382 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} {
383 add EtaPhiBins $eta $PhiBins
384 }
385
386 # 20 degrees towers
387 set PhiBins {}
388 for {set i -9} {$i <= 9} {incr i} {
389 add PhiBins [expr {$i * $pi/9.0}]
390 }
391 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} {
392 add EtaPhiBins $eta $PhiBins
393 }
394
395 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
396 add EnergyFraction {0} {1.0}
397 # energy fractions for e, gamma and pi0
398 add EnergyFraction {11} {0.0}
399 add EnergyFraction {22} {0.0}
400 add EnergyFraction {111} {0.0}
401 # energy fractions for muon, neutrinos and neutralinos
402 add EnergyFraction {12} {0.0}
403 add EnergyFraction {13} {0.0}
404 add EnergyFraction {14} {0.0}
405 add EnergyFraction {16} {0.0}
406 add EnergyFraction {1000022} {0.0}
407 add EnergyFraction {1000023} {0.0}
408 add EnergyFraction {1000025} {0.0}
409 add EnergyFraction {1000035} {0.0}
410 add EnergyFraction {1000045} {0.0}
411 # energy fractions for K0short and Lambda
412 add EnergyFraction {310} {0.7}
413 add EnergyFraction {3122} {0.7}
414
415 # http://arxiv.org/pdf/hep-ex/0004009v1
416 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
417 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
418 set ResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
419 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
420 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
421
422# https://cds.cern.ch/record/2037117/files/ATL-CAL-SLIDE-2015-388.pdf
423 #set ResolutionFormula { 0.057 * energy oplus 0.529 * sqrt(energy)}
424 #set ResolutionFormula { sqrt(energy^2 * 0.057^2 + energy * 0.529^2)}
425 #set ResolutionFormula {0.01}
426}
427
428
429#################
430# Electron filter
431#################
432
433module PdgCodeFilter ElectronFilter {
434 set InputArray HCal/eflowTracks
435 set OutputArray electrons
436 set Invert true
437 add PdgCode {11}
438 add PdgCode {-11}
439}
440
441######################
442# ChargedHadronFilter
443######################
444
445module PdgCodeFilter ChargedHadronFilter {
446 set InputArray HCal/eflowTracks
447 set OutputArray chargedHadrons
448
449 add PdgCode {11}
450 add PdgCode {-11}
451 add PdgCode {13}
452 add PdgCode {-13}
453}
454
455
456
457###################################################
458# Tower Merger (in case not using e-flow algorithm)
459###################################################
460
461module Merger Calorimeter {
462# add InputArray InputArray
463 add InputArray ECal/ecalTowers
464 add InputArray HCal/hcalTowers
465 add InputArray MuonMomentumSmearing/muons
466 set OutputArray towers
467}
468
469####################
470# Energy flow merger
471####################
472
473## Calorimeter/towers -> EFlowMerger/eflow ??
474
475module Merger EFlowMerger {
476# add InputArray InputArray
477 add InputArray HCal/eflowTracks
478 add InputArray ECal/eflowPhotons
479 add InputArray HCal/eflowNeutralHadrons
480 set OutputArray eflow
481}
482
483######################
484# EFlowFilter
485######################
486
487module PdgCodeFilter EFlowFilter {
488 set InputArray EFlowMerger/eflow
489 set OutputArray eflow
490
491 add PdgCode {11}
492 add PdgCode {-11}
493 add PdgCode {13}
494 add PdgCode {-13}
495}
496
497###################
498# Photon efficiency
499###################
500
501module Efficiency PhotonEfficiency {
502 set InputArray ECal/eflowPhotons
503 set OutputArray photons
504
505 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
506## don't really care about this
507 # efficiency formula for photons
508# set EfficiencyFormula { (pt <= 10.0) * (0.00) +
509# (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
510# (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
511# (abs(eta) > 2.5) * (0.00)}
512
513
514 set EfficiencyFormula { (pt <= 5.0) * (0.00) +
515 (abs(eta) <= 1.5) * (pt > 5.0) * (0.99) +
516 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 5.0) * (0.99) +
517 (abs(eta) > 2.7) * (0.00)}
518}
519
520##################
521# Photon isolation
522##################
523
524# module Isolation PhotonIsolation {
525# set CandidateInputArray PhotonEfficiency/photons
526# set IsolationInputArray EFlowFilter/eflow
527#
528# set OutputArray photons
529#
530# set DeltaRMax 0.5
531#
532# set PTMin 0.5
533#
534# set PTRatioMax 0.12
535# }
536
537
538#####################
539# Electron efficiency
540#####################
541
542module Efficiency ElectronEfficiency {
543 set InputArray ElectronFilter/electrons
544 set OutputArray electrons
545
546 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
547
548 set EfficiencyFormula {
549 (pt <= 5.0) * (0.00) +
550 (abs(eta) <= 2.5) *(pt > 5.0 && pt < 6.0) * (0.85) +
551 (abs(eta) <= 2.5) *(pt >= 6.0) * (0.98) +
552 (abs(eta) > 2.5) * 0
553 }
554
555 # set EfficiencyFormula {
556 # (pt <= 10.0) * (0.00) +
557 # (abs(eta) > 0.1 && abs(eta) < 2.47) *(pt > 10.0 && pt <= 15.0) * (0.80) +
558 # (abs(eta) > 0.1 && abs(eta) < 2.47) *(pt > 15.0 && pt <= 35.0) * (0.82) +
559 # (abs(eta) > 0.1 && abs(eta) < 2.47) *(pt > 35.0 && pt <= 60.0) * (0.89) +
560 # (abs(eta) > 0.1 && abs(eta) < 2.47) *(pt > 60) * (0.93) +
561 # (abs(eta) > 2.47) * 0}
562}
563
564####################
565# Electron isolation
566####################
567
568
569# module Isolation ElectronIsolation {
570# set CandidateInputArray ElectronEfficiency/electrons
571# set IsolationInputArray EFlowFilter/eflow
572#
573# set OutputArray electrons
574#
575# set DeltaRMax 0.5
576#
577# set PTMin 0.5
578#
579# set PTRatioMax 0.12
580# }
581
582
583#################
584# Muon efficiency
585#################
586
587module Efficiency MuonEfficiency {
588 set InputArray MuonMomentumSmearing/muons
589 set OutputArray muons
590
591 # set EfficiencyFormula {efficiency as a function of eta and pt}
592
593
594 # set EfficiencyFormula {
595 # (pt <= 10.0) * (0.00) +
596 # (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
597 # (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
598 # (abs(eta) > 2.5) * (0.00)
599 # }
600 # set EfficiencyFormula {
601 # (pt <= 5.0) * (0.00) +
602 # (abs(eta) < 2.7) *(pt > 5.0) * (0.99) +
603 # (abs(eta) > 2.7) * 0
604 # }
605
606 ## 1603.05598
607 set EfficiencyFormula {
608 (pt <= 4.0) * (0.00) +
609 (pt > 4.0) * (abs(eta) <= 0.1) * (0.9) +
610 (pt > 4.0 && pt <= 6.0) *(abs(eta) >0.1 && abs(eta) <= 2.5) *(0.98) +
611 (pt > 6.0) *(abs(eta) >0.1 && abs(eta) <= 2.5) *(0.99) +
612 (abs(eta) > 2.5 && abs(eta) <= 2.7) *(pt <= 6.0) * (0.00) +
613 (abs(eta) > 2.5 && abs(eta) <= 2.7) *(pt > 6.0) * (0.9) +
614 (abs(eta) > 2.7) * (0.00)
615 }
616
617
618 # set EfficiencyFormula {
619 # (pt <= 6.0) * (0.00) + \
620 # (abs(eta) <= 0.1) * (pt > 6.0) * (0.573) + \
621 # (abs(eta) > 0.1 && abs(eta) <=1.05) * (pt > 6.0) * (0.947) + \
622 # (abs(eta) > 1.05 && abs(eta) <=1.5) * (pt > 6.0) * (0.943) + \
623 # (abs(eta) > 1.5 && abs(eta) <=2.) * (pt > 6.0) * (0.949) + \
624 # (abs(eta) > 2. && abs(eta) <=2.5) * (pt > 6.0) * (0.949) + \
625 # (abs(eta) > 2.5) * (0.00)}
626
627
628 # set EfficiencyFormula {
629 # (pt <= 5.0) * (0.00) +
630 # (abs(eta) < 0.1) * (pt > 5.0)* (0.9) +
631 # (abs(eta) >= 0.1 && abs(eta) < 2.5) * (pt > 5.0 && pt < 10.0)* (0.98) +
632 # (abs(eta) >= 2.5 && abs(eta) < 2.7) * (pt > 5.0 && pt < 10.0)* (0.0) +
633 # (abs(eta) >= 0.1 && abs(eta) < 2.7) * (pt > 10.0)* (0.99)
634 # }
635
636}
637
638################
639# Muon isolation
640################
641
642# module Isolation MuonIsolation {
643# set CandidateInputArray MuonEfficiency/muons
644# set IsolationInputArray EFlowFilter/eflow
645#
646# set OutputArray muons
647#
648# set DeltaRMax 0.5
649#
650# set PTMin 0.5
651#
652# set PTRatioMax 0.25
653# }
654
655###################
656# Missing ET merger
657###################
658
659module Merger MissingET {
660# add InputArray InputArray
661 add InputArray Calorimeter/towers
662 set MomentumOutputArray momentum
663}
664
665##################
666# Scalar HT merger
667##################
668
669#module Merger ScalarHT {
670## add InputArray InputArray
671# add InputArray UniqueObjectFinder/jets
672# add InputArray UniqueObjectFinder/electrons
673# add InputArray UniqueObjectFinder/photons
674# add InputArray UniqueObjectFinder/muons
675# set EnergyOutputArray energy
676#}
677
678
679#####################
680# Neutrino Filter
681#####################
682
683module PdgCodeFilter NeutrinoFilter {
684
685 set InputArray Delphes/stableParticles
686 set OutputArray filteredParticles
687
688 set PTMin 0.0
689
690 add PdgCode {12}
691 add PdgCode {14}
692 add PdgCode {16}
693 add PdgCode {-12}
694 add PdgCode {-14}
695 add PdgCode {-16}
696
697 add PdgCode {1000022}
698}
699
700#####################
701# MC truth jet finder
702#####################
703
704module FastJetFinder GenJetFinder {
705 set InputArray NeutrinoFilter/filteredParticles
706
707 set OutputArray jets
708
709 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
710 set JetAlgorithm 6
711 set ParameterR 0.4
712
713 set JetPTMin 20.0
714}
715
716
717#########################
718# Gen Missing ET merger
719########################
720
721module Merger GenMissingET {
722# add InputArray InputArray
723 add InputArray NeutrinoFilter/filteredParticles
724 set MomentumOutputArray momentum
725}
726
727
728
729############
730# Jet finder
731############
732
733module FastJetFinder FastJetFinder {
734 set InputArray Calorimeter/towers
735
736 set OutputArray jets
737
738 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
739 set JetAlgorithm 6
740 set ParameterR 0.4
741
742 set JetPTMin 20.0
743}
744
745##################
746# Jet Energy Scale
747##################
748
749module EnergyScale JetEnergyScale {
750 set InputArray FastJetFinder/jets
751 set OutputArray jets
752
753 # scale formula for jets
754 set ScaleFormula { sqrt( (3.0 - 0.2*(abs(eta)))^2 / pt + 1.0 ) }
755}
756
757########################
758# Jet Flavor Association
759########################
760
761module JetFlavorAssociation JetFlavorAssociation {
762
763 set PartonInputArray Delphes/partons
764 set ParticleInputArray Delphes/allParticles
765 set ParticleLHEFInputArray Delphes/allParticlesLHEF
766 set JetInputArray JetEnergyScale/jets
767
768 set DeltaR 0.5
769 set PartonPTMin 1.0
770 set PartonEtaMax 2.5
771
772}
773
774###########
775# b-tagging
776###########
777
778module BTagging BTagging {
779 set JetInputArray JetEnergyScale/jets
780
781 set BitNumber 0
782
783 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
784 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
785 # gluon's PDG code has the lowest priority
786
787 # based on ATL-PHYS-PUB-2015-022
788
789 # default efficiency formula (misidentification rate)
790 add EfficiencyFormula {0} {0.002+7.3e-06*pt}
791
792 # efficiency formula for c-jets (misidentification rate)
793 add EfficiencyFormula {4} {0.20*tanh(0.02*pt)*(1/(1+0.0034*pt))}
794
795 # efficiency formula for b-jets
796 # add EfficiencyFormula {5} {0.85*tanh(0.003*pt)*(30/(1+0.086*pt))}
797
798 # add EfficiencyFormula {5} { 0.77}
799
800 # add EfficiencyFormula {0} {
801 # (0<=pt && pt<=20 ) * (1.0/16.0 ) +
802 # (pt>20 && pt<=40 ) * (1.0/32.0 ) +
803# (pt>40 && pt<=60 ) * (1.0/120.0) +
804# (pt>60 && pt<=90 ) * (1.0/170.0) +
805# (pt>90 && pt<=150) * (1.0/195.0) +
806# (pt>150 && pt<=200) * (1.0/180.0) +
807# (pt>200 && pt<=300) * (1.0/85.0 ) +
808# (pt>300 && pt<=500) * (1.0/30.0 ) +
809# (pt>500 && pt<=800) * (1.0/10.0 ) +
810# (pt>800 ) * (1.0/10.0)
811#
812# }
813
814
815 # efficiency formula for c-jets (misidentification rate)
816# add EfficiencyFormula {4} {
817# (0<=pt && pt<=20 ) * (1.0/1.65) +
818# (pt>20 && pt<=40 ) * (1.0/3.3 ) +
819# (pt>40 && pt<=60 ) * (1.0/4.8 ) +
820# (pt>60 && pt<=90 ) * (1.0/6.0 ) +
821# (pt>90 && pt<=150) * (1.0/7.0 ) +
822# (pt>150 && pt<=200) * (1.0/6.2 ) +
823# (pt>200 && pt<=300) * (1.0/5.5 ) +
824# (pt>300 && pt<=500) * (1.0/4.2 ) +
825# (pt>500 && pt<=800) * (1.0/3.2 ) +
826# (pt>800 ) * (1.0/3.2 )
827# }
828
829 # efficiency formula for b-jets
830# add EfficiencyFormula {5} {
831# (0<=pt && pt<=20 ) * (0.34 ) +
832# (pt>20 && pt<=40 ) * (0.68 ) +
833# (pt>40 && pt<=60 ) * (0.77 ) +
834# (pt>60 && pt<=90 ) * (0.8 ) +
835# (pt>90 && pt<=150) * (0.82 ) +
836# (pt>150 && pt<=200) * (0.83 ) +
837# (pt>200 && pt<=300) * (0.78 ) +
838# (pt>300 && pt<=500) * (0.74 ) +
839# (pt>500 && pt<=800) * (0.67 ) +
840# (pt>800 ) * (0.67)
841# }
842
843 # efficiency formula for b-jets
844# add EfficiencyFormula {5} {
845# (0<=pt && pt<=20 ) * (0.34 ) +
846# (pt>20 && pt<=350 ) * (0.9 ) +
847# (pt>350) * (0.77 )
848 # }
849
850# add EfficiencyFormula {5} {
851# (0<=pt && pt<=20 ) * (0.34 ) +
852# (pt>20 && pt<=40 ) * (0.68 ) +
853# (pt>40 && pt<=60 ) * (0.77 ) +
854# (pt>60 && pt<=90 ) * (0.8 ) +
855# (pt>90 && pt<=150) * (0.82 ) +
856# (pt>150 && pt<=200) * (0.83 ) +
857# (pt>200 && pt<=300) * (0.78 ) +
858# (pt>300) * (0.77 )
859# }
860
861add EfficiencyFormula {5} {
862(0<=pt && pt<=20 ) * (0.34 ) +
863(pt>20 && pt<=300 ) * (0.85 ) +
864(pt>300) * (0.77 )
865}
866
867
868
869
870}
871
872#############
873# tau-tagging
874#############
875
876module TrackCountingTauTagging TauTagging {
877
878 set ParticleInputArray Delphes/allParticles
879 set PartonInputArray Delphes/partons
880 set TrackInputArray TrackMerger/tracks
881 set JetInputArray JetEnergyScale/jets
882
883 set DeltaR 0.2
884 set DeltaRTrack 0.2
885
886 set TrackPTMin 1.0
887
888 set TauPTMin 1.0
889 set TauEtaMax 2.5
890
891 # instructions: {n-prongs} {eff}
892
893 # 1 - one prong efficiency
894 # 2 - two or more efficiency
895 # -1 - one prong mistag rate
896 # -2 - two or more mistag rate
897
898 set BitNumber 0
899
900 # taken from ATL-PHYS-PUB-2015-045 (medium working point)
901 add EfficiencyFormula {1} {0.70}
902 add EfficiencyFormula {2} {0.60}
903 add EfficiencyFormula {-1} {0.02}
904 add EfficiencyFormula {-2} {0.01}
905
906}
907
908#####################################################
909# Find uniquely identified photons/electrons/tau/jets
910#####################################################
911
912# module UniqueObjectFinder UniqueObjectFinder {
913# # earlier arrays take precedence over later ones
914# # add InputArray InputArray OutputArray
915# add InputArray PhotonIsolation/photons photons
916# add InputArray ElectronIsolation/electrons electrons
917# add InputArray MuonIsolation/muons muons
918# add InputArray JetEnergyScale/jets jets
919# }
920
921##################
922# ROOT tree writer
923##################
924
925# tracks, towers and eflow objects are not stored by default in the output.
926# if needed (for jet constituent or other studies), uncomment the relevant
927# "add Branch ..." lines.
928
929module TreeWriter TreeWriter {
930 # add Branch InputArray BranchName BranchClass
931 add Branch Delphes/allParticles Particle GenParticle
932
933 add Branch TrackMerger/tracks Track Track
934# add Branch Calorimeter/towers Tower Tower
935
936 add Branch HCal/eflowTracks EFlowTrack Track
937 add Branch ECal/eflowPhotons EFlowPhoton Tower
938 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
939
940# add Branch GenJetFinder/jets GenJet Jet
941# add Branch GenMissingET/momentum GenMissingET MissingET
942
943 add Branch JetEnergyScale/jets Jet Jet
944 add Branch ElectronEfficiency/electrons Electron Electron
945 add Branch PhotonEfficiency/photons Photon Photon
946 add Branch MuonEfficiency/muons Muon Muon
947 add Branch MissingET/momentum MissingET MissingET
948# add Branch ScalarHT/energy ScalarHT ScalarHT
949 }
950