MA5SandBox: delphes_card_atlas_2016_25.tcl

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