Fork me on GitHub

Ticket #1115: delphes_card_ATLAS_0424.tcl

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