PublicAnalysisDatabase: delphes_card_ATLAS_CONF_2016_086.tcl

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