Fork me on GitHub

Ticket #1372: card_ATLAS_subjet.tcl

File card_ATLAS_subjet.tcl, 21.2 KB (added by andrea.diluca, 6 years ago)

card

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.5
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.5
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.5
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.6
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# Jet finder
615############
616
617module FastJetFinder FastJetFinder {
618
619 set InputArray EFlowMerger/eflow
620 set OutputArray jets
621
622 set JetAlgorithm 6
623 set ParameterR 1.0
624
625 set ComputeNsubjettiness 1
626 set Beta 1.0
627 set AxisMode 1
628
629 set ComputeTrimming 1
630 set RTrim 0.2
631 set PtFracTrim 0.05
632
633 set ComputePruning 1
634 set ZcutPrun 0.1
635 set RcutPrun 0.5
636 set RPrun 0.8
637
638 set ComputeSoftDrop 1
639 set BetaSoftDrop 0.0
640 set SymmetryCutSoftDrop 0.1
641 set R0SoftDrop 0.8
642
643 set JetPTMin 200.0
644}
645
646
647##################
648# Jet Energy Scale
649##################
650
651module EnergyScale JetEnergyScale {
652 set InputArray FastJetFinder/jets
653 set OutputArray jets
654
655 # scale formula for jets
656 set ScaleFormula { sqrt( (3.0 - 0.2*(abs(eta)))^2 / pt + 1.0 ) }
657}
658
659########################
660# Jet Flavor Association
661########################
662
663module JetFlavorAssociation JetFlavorAssociation {
664
665 set PartonInputArray Delphes/partons
666 set ParticleInputArray Delphes/allParticles
667 set ParticleLHEFInputArray Delphes/allParticlesLHEF
668 set JetInputArray JetEnergyScale/jets
669
670 set DeltaR 0.5
671 set PartonPTMin 1.0
672 set PartonEtaMax 2.5
673
674}
675
676###########
677# b-tagging
678###########
679
680module BTagging BTagging {
681 set JetInputArray JetEnergyScale/jets
682
683 set BitNumber 0
684
685 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
686 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
687 # gluon's PDG code has the lowest priority
688
689 # based on ATL-PHYS-PUB-2015-022
690
691 # default efficiency formula (misidentification rate)
692 add EfficiencyFormula {0} {0.002+7.3e-06*pt}
693
694 # efficiency formula for c-jets (misidentification rate)
695 add EfficiencyFormula {4} {0.20*tanh(0.02*pt)*(1/(1+0.0034*pt))}
696
697 # efficiency formula for b-jets
698 add EfficiencyFormula {5} {0.80*tanh(0.003*pt)*(30/(1+0.086*pt))}
699}
700
701#############
702# tau-tagging
703#############
704
705module TrackCountingTauTagging TauTagging {
706
707 set ParticleInputArray Delphes/allParticles
708 set PartonInputArray Delphes/partons
709 set TrackInputArray TrackMerger/tracks
710 set JetInputArray JetEnergyScale/jets
711
712 set DeltaR 0.2
713 set DeltaRTrack 0.2
714
715 set TrackPTMin 1.0
716
717 set TauPTMin 1.0
718 set TauEtaMax 2.5
719
720 # instructions: {n-prongs} {eff}
721
722 # 1 - one prong efficiency
723 # 2 - two or more efficiency
724 # -1 - one prong mistag rate
725 # -2 - two or more mistag rate
726
727 set BitNumber 0
728
729 # taken from ATL-PHYS-PUB-2015-045 (medium working point)
730 add EfficiencyFormula {1} {0.70}
731 add EfficiencyFormula {2} {0.60}
732 add EfficiencyFormula {-1} {0.02}
733 add EfficiencyFormula {-2} {0.01}
734
735}
736
737#####################################################
738# Find uniquely identified photons/electrons/tau/jets
739#####################################################
740
741module UniqueObjectFinder UniqueObjectFinder {
742# earlier arrays take precedence over later ones
743# add InputArray InputArray OutputArray
744 add InputArray PhotonIsolation/photons photons
745 add InputArray ElectronIsolation/electrons electrons
746 add InputArray MuonIsolation/muons muons
747 add InputArray JetEnergyScale/jets jets
748}
749
750##################
751# ROOT tree writer
752##################
753
754# tracks, towers and eflow objects are not stored by default in the output.
755# if needed (for jet constituent or other studies), uncomment the relevant
756# "add Branch ..." lines.
757
758module TreeWriter TreeWriter {
759# add Branch InputArray BranchName BranchClass
760 add Branch Delphes/allParticles Particle GenParticle
761
762 add Branch TrackMerger/tracks Track Track
763 add Branch Calorimeter/towers Tower Tower
764
765 add Branch HCal/eflowTracks EFlowTrack Track
766 add Branch ECal/eflowPhotons EFlowPhoton Tower
767 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
768
769 add Branch GenJetFinder/jets GenJet Jet
770 add Branch GenMissingET/momentum GenMissingET MissingET
771
772 add Branch UniqueObjectFinder/jets Jet Jet
773 add Branch UniqueObjectFinder/electrons Electron Electron
774 add Branch UniqueObjectFinder/photons Photon Photon
775 add Branch UniqueObjectFinder/muons Muon Muon
776 add Branch MissingET/momentum MissingET MissingET
777 add Branch ScalarHT/energy ScalarHT ScalarHT
778}
779