MA5SandBox: delphes_card_atlas_susy_2018_06.tcl

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