Fork me on GitHub

source: git/cards/delphes_card_ATLAS_PileUp.tcl@ aa3c287

ImprovedOutputFile Timing dual_readout llp 3.4.2pre16
Last change on this file since aa3c287 was 070a19b, checked in by Michele Selvaggi <michele.selvaggi@…>, 6 years ago

fix isolation input collection

  • Property mode set to 100644
File size: 21.0 KB
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# Energy flow merger
364####################
365
366module Merger EFlowMerger {
367# add InputArray InputArray
368 add InputArray Calorimeter/eflowTracks
369 add InputArray Calorimeter/eflowPhotons
370 add InputArray Calorimeter/eflowNeutralHadrons
371 set OutputArray eflow
372}
373
374######################
375# EFlowFilter
376######################
377
378module PdgCodeFilter EFlowFilter {
379 set InputArray EFlowMerger/eflow
380 set OutputArray eflow
381
382 add PdgCode {11}
383 add PdgCode {-11}
384 add PdgCode {13}
385 add PdgCode {-13}
386}
387
388
389#############
390# Rho pile-up
391#############
392
393module FastJetGridMedianEstimator Rho {
394
395 set InputArray Calorimeter/towers
396 set RhoOutputArray rho
397
398 # add GridRange rapmin rapmax drap dphi
399 # rapmin - the minimum rapidity extent of the grid
400 # rapmax - the maximum rapidity extent of the grid
401 # drap - the grid spacing in rapidity
402 # dphi - the grid spacing in azimuth
403
404 add GridRange -5.0 -2.5 1.0 1.0
405 add GridRange -2.5 2.5 0.5 0.5
406 add GridRange 2.5 5.0 1.0 1.0
407
408}
409
410
411#####################
412# Neutrino Filter
413#####################
414
415module PdgCodeFilter NeutrinoFilter {
416
417 set InputArray Delphes/stableParticles
418 set OutputArray filteredParticles
419
420 set PTMin 0.0
421
422 add PdgCode {12}
423 add PdgCode {14}
424 add PdgCode {16}
425 add PdgCode {-12}
426 add PdgCode {-14}
427 add PdgCode {-16}
428
429}
430
431#####################
432# MC truth jet finder
433#####################
434
435module FastJetFinder GenJetFinder {
436 set InputArray NeutrinoFilter/filteredParticles
437
438 set OutputArray jets
439
440 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
441 set JetAlgorithm 6
442 set ParameterR 0.6
443
444 set JetPTMin 20.0
445}
446
447#########################
448# Gen Missing ET merger
449########################
450
451module Merger GenMissingET {
452# add InputArray InputArray
453 add InputArray NeutrinoFilter/filteredParticles
454 set MomentumOutputArray momentum
455}
456
457
458############
459# Jet finder
460############
461
462module FastJetFinder FastJetFinder {
463 set InputArray Calorimeter/towers
464
465 set OutputArray jets
466
467 # 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
468 set AreaAlgorithm 5
469
470 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
471 set JetAlgorithm 6
472 set ParameterR 0.6
473
474 set JetPTMin 20.0
475}
476
477###########################
478# Jet Pile-Up Subtraction
479###########################
480
481module JetPileUpSubtractor JetPileUpSubtractor {
482 set JetInputArray FastJetFinder/jets
483 set RhoInputArray Rho/rho
484
485 set OutputArray jets
486
487 set JetPTMin 20.0
488}
489
490##################
491# Jet Energy Scale
492##################
493
494module EnergyScale JetEnergyScale {
495 set InputArray JetPileUpSubtractor/jets
496 set OutputArray jets
497
498 # scale formula for jets
499 set ScaleFormula {1.0}
500}
501
502###################
503# Photon efficiency
504###################
505
506module Efficiency PhotonEfficiency {
507 set InputArray Calorimeter/eflowPhotons
508 set OutputArray photons
509
510 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
511
512 # efficiency formula for photons
513 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
514 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
515 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
516 (abs(eta) > 2.5) * (0.00)}
517}
518
519##################
520# Photon isolation
521##################
522
523module Isolation PhotonIsolation {
524 set CandidateInputArray PhotonEfficiency/photons
525 set IsolationInputArray EFlowFilter/eflow
526 set RhoInputArray Rho/rho
527
528 set OutputArray photons
529
530 set DeltaRMax 0.5
531
532 set PTMin 0.5
533
534 set PTRatioMax 0.12
535}
536
537#####################
538# Electron efficiency
539#####################
540
541module Efficiency ElectronEfficiency {
542 set InputArray TrackPileUpSubtractor/electrons
543 set OutputArray electrons
544
545 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
546
547 # efficiency formula for electrons
548 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
549 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
550 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
551 (abs(eta) > 2.5) * (0.00)}
552}
553
554####################
555# Electron isolation
556####################
557
558module Isolation ElectronIsolation {
559 set CandidateInputArray ElectronEfficiency/electrons
560 set IsolationInputArray EFlowFilter/eflow
561 set RhoInputArray Rho/rho
562
563 set OutputArray electrons
564
565 set DeltaRMax 0.5
566
567 set PTMin 0.5
568
569 set PTRatioMax 0.12
570}
571
572#################
573# Muon efficiency
574#################
575
576module Efficiency MuonEfficiency {
577 set InputArray TrackPileUpSubtractor/muons
578 set OutputArray muons
579
580 # set EfficiencyFormula {efficiency as a function of eta and pt}
581
582 # efficiency formula for muons
583 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
584 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
585 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) +
586 (abs(eta) > 2.7) * (0.00)}
587}
588
589################
590# Muon isolation
591################
592
593module Isolation MuonIsolation {
594 set CandidateInputArray MuonEfficiency/muons
595 set IsolationInputArray EFlowFilter/eflow
596 set RhoInputArray Rho/rho
597
598 set OutputArray muons
599
600 set DeltaRMax 0.5
601
602 set PTMin 0.5
603
604 set PTRatioMax 0.25
605}
606
607###################
608# Missing ET merger
609###################
610
611module Merger MissingET {
612# add InputArray InputArray
613 add InputArray EFlowMergerAllTracks/eflow
614 set MomentumOutputArray momentum
615}
616
617
618##################
619# Scalar HT merger
620##################
621
622module Merger ScalarHT {
623# add InputArray InputArray
624 add InputArray UniqueObjectFinder/jets
625 add InputArray UniqueObjectFinder/electrons
626 add InputArray UniqueObjectFinder/photons
627 add InputArray UniqueObjectFinder/muons
628 set EnergyOutputArray energy
629}
630
631########################
632# Jet Flavor Association
633########################
634
635module JetFlavorAssociation JetFlavorAssociation {
636
637 set PartonInputArray Delphes/partons
638 set ParticleInputArray Delphes/allParticles
639 set ParticleLHEFInputArray Delphes/allParticlesLHEF
640 set JetInputArray JetEnergyScale/jets
641
642 set DeltaR 0.5
643 set PartonPTMin 1.0
644 set PartonEtaMax 2.5
645
646}
647
648###########
649# b-tagging
650###########
651
652module BTagging BTagging {
653 set JetInputArray JetEnergyScale/jets
654
655 set BitNumber 0
656
657 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
658 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
659 # gluon's PDG code has the lowest priority
660
661 # based on ATL-PHYS-PUB-2015-022
662
663 # default efficiency formula (misidentification rate)
664 add EfficiencyFormula {0} {0.002+7.3e-06*pt}
665
666 # efficiency formula for c-jets (misidentification rate)
667 add EfficiencyFormula {4} {0.20*tanh(0.02*pt)*(1/(1+0.0034*pt))}
668
669 # efficiency formula for b-jets
670 add EfficiencyFormula {5} {0.80*tanh(0.003*pt)*(30/(1+0.086*pt))}
671}
672
673#############
674# tau-tagging
675#############
676
677module TrackCountingTauTagging TauTagging {
678
679 set ParticleInputArray Delphes/allParticles
680 set PartonInputArray Delphes/partons
681 set TrackInputArray TrackMerger/tracks
682 set JetInputArray JetEnergyScale/jets
683
684 set DeltaR 0.2
685 set DeltaRTrack 0.2
686
687 set TrackPTMin 1.0
688
689 set TauPTMin 1.0
690 set TauEtaMax 2.5
691
692 # instructions: {n-prongs} {eff}
693
694 # 1 - one prong efficiency
695 # 2 - two or more efficiency
696 # -1 - one prong mistag rate
697 # -2 - two or more mistag rate
698
699 set BitNumber 0
700
701 # taken from ATL-PHYS-PUB-2015-045 (medium working point)
702 add EfficiencyFormula {1} {0.70}
703 add EfficiencyFormula {2} {0.60}
704 add EfficiencyFormula {-1} {0.02}
705 add EfficiencyFormula {-2} {0.01}
706
707}
708
709#####################################################
710# Find uniquely identified photons/electrons/tau/jets
711#####################################################
712
713module UniqueObjectFinder UniqueObjectFinder {
714# earlier arrays take precedence over later ones
715# add InputArray InputArray OutputArray
716 add InputArray PhotonIsolation/photons photons
717 add InputArray ElectronIsolation/electrons electrons
718 add InputArray MuonIsolation/muons muons
719 add InputArray JetEnergyScale/jets jets
720}
721
722##################
723# ROOT tree writer
724##################
725
726# tracks, towers and eflow objects are not stored by default in the output.
727# if needed (for jet constituent or other studies), uncomment the relevant
728# "add Branch ..." lines.
729
730module TreeWriter TreeWriter {
731# add Branch InputArray BranchName BranchClass
732 add Branch Delphes/allParticles Particle GenParticle
733
734# add Branch TrackMerger/tracks Track Track
735 add Branch Calorimeter/towers Tower Tower
736
737# add Branch Calorimeter/eflowTracks EFlowTrack Track
738# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
739# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
740
741 add Branch GenJetFinder/jets GenJet Jet
742 add Branch GenMissingET/momentum GenMissingET MissingET
743
744 add Branch UniqueObjectFinder/jets Jet Jet
745 add Branch UniqueObjectFinder/electrons Electron Electron
746 add Branch UniqueObjectFinder/photons Photon Photon
747 add Branch UniqueObjectFinder/muons Muon Muon
748 add Branch MissingET/momentum MissingET MissingET
749 add Branch ScalarHT/energy ScalarHT ScalarHT
750 add Branch Rho/rho Rho Rho
751 add Branch PileUpMerger/vertices Vertex Vertex
752
753}
Note: See TracBrowser for help on using the repository browser.