Fork me on GitHub

source: git/cards/delphes_card_ATLAS_PileUp.tcl@ 679a87e

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 679a87e was 934d037, checked in by Pavel Demin <pavel.demin@…>, 9 years ago

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