Fork me on GitHub

Ticket #1308: delphes_card.dat

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