MA5SandBox: delphes_card_atlas_2016_32.tcl

File delphes_card_atlas_2016_32.tcl, 21.0 KB (added by Benjamin Fuks, 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 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
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.6
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