Fork me on GitHub

source: svn/trunk/examples/delphes_card_ATLAS_PileUp.tcl@ 1392

Last change on this file since 1392 was 1378, checked in by Pavel Demin, 11 years ago

fix formatting
File size: 21.0 KB
RevLine 
[1260]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 ElectronEnergySmearing
16 MuonMomentumSmearing
17
18 TrackMerger
19 Calorimeter
20 TrackPileUpSubtractor
[1357]21 NeutralTowerMerger
[1260]22 EFlowMerger
23
24 GenJetFinder
25
26 Rho
27 FastJetFinder
28 JetPileUpSubtractor
29
30 JetEnergyScale
31
32 PhotonEfficiency
33 PhotonIsolation
34
35 ElectronEfficiency
36 ElectronIsolation
37
38 MuonEfficiency
39 MuonIsolation
40
41 MissingET
42
43 BTagging
44 TauTagging
45
46 UniqueObjectFinder
47
48 ScalarHT
49
50 TreeWriter
51}
52
53###############
54# PileUp Merger
55###############
56
57module PileUpMerger PileUpMerger {
58 set InputArray Delphes/stableParticles
59
[1326]60 set ParticleOutputArray stableParticles
61 set VertexOutputArray vertices
[1260]62
63 # pre-generated minbias input file
64 set PileUpFile MinBias.pileup
65
66 # average expected pile up
67 set MeanPileUp 50
[1378]68
69 # maximum spread in the beam direction in m
[1345]70 set ZVertexSpread 0.10
[1378]71
[1345]72 # maximum spread in time in s
73 set TVertexSpread 1.5E-09
74
75 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
[1378]76
[1345]77 set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
78
79 #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) + \
80 # (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) + \
81 # (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) + \
82 # (abs(t) > 1.0e-09) * (abs(z) > 0.15) * (0.00)}
83
84
[1260]85}
86
87#################################
88# Propagate particles in cylinder
89#################################
90
91module ParticlePropagator ParticlePropagator {
92 set InputArray PileUpMerger/stableParticles
93
94 set OutputArray stableParticles
95 set ChargedHadronOutputArray chargedHadrons
96 set ElectronOutputArray electrons
97 set MuonOutputArray muons
98
99 # radius of the magnetic field coverage, in m
100 set Radius 1.15
101 # half-length of the magnetic field coverage, in m
102 set HalfLength 3.51
103
104 # magnetic field
105 set Bz 2.0
106}
107
108####################################
109# Charged hadron tracking efficiency
110####################################
111
112module Efficiency ChargedHadronTrackingEfficiency {
113 set InputArray ParticlePropagator/chargedHadrons
114 set OutputArray chargedHadrons
115
116 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
117
118 # tracking efficiency formula for charged hadrons
119 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
120 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
121 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \
122 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + \
123 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + \
124 (abs(eta) > 2.5) * (0.00)}
125}
126
127##############################
128# Electron tracking efficiency
129##############################
130
131module Efficiency ElectronTrackingEfficiency {
132 set InputArray ParticlePropagator/electrons
133 set OutputArray electrons
134
135 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
136
137 # tracking efficiency formula for electrons
138 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
139 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \
140 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \
141 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \
142 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + \
143 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + \
144 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + \
145 (abs(eta) > 2.5) * (0.00)}
146}
147
148##########################
149# Muon tracking efficiency
150##########################
151
152module Efficiency MuonTrackingEfficiency {
153 set InputArray ParticlePropagator/muons
154 set OutputArray muons
155
156 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
157
158 # tracking efficiency formula for muons
159 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
160 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \
161 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \
162 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
163 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \
164 (abs(eta) > 2.5) * (0.00)}
165}
166
167########################################
168# Momentum resolution for charged tracks
169########################################
170
171module MomentumSmearing ChargedHadronMomentumSmearing {
172 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
173 set OutputArray chargedHadrons
174
175 # set ResolutionFormula {resolution formula as a function of eta and pt}
176
177 # resolution formula for charged hadrons
[1279]178 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
179 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
180 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
181 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
182 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
183 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \
184 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
185 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
[1260]186}
187
188#################################
189# Energy resolution for electrons
190#################################
191
192module EnergySmearing ElectronEnergySmearing {
193 set InputArray ElectronTrackingEfficiency/electrons
194 set OutputArray electrons
195
196 # set ResolutionFormula {resolution formula as a function of eta and energy}
197
198 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.015) + \
199 (abs(eta) <= 2.5) * (energy > 2.5e1) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
200 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
201 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
202
203}
204
205###############################
206# Momentum resolution for muons
207###############################
208
209module MomentumSmearing MuonMomentumSmearing {
210 set InputArray MuonTrackingEfficiency/muons
211 set OutputArray muons
212
213 # set ResolutionFormula {resolution formula as a function of eta and pt}
214
215 # resolution formula for muons
216 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
217 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 5.0e1) * (0.03) + \
218 (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) + \
219 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.07) + \
220 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.04) + \
221 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 5.0e1) * (0.04) + \
222 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) + \
223 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.10)}
224}
225
226##############
227# Track merger
228##############
229
230module Merger TrackMerger {
231# add InputArray InputArray
232 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
233 add InputArray ElectronEnergySmearing/electrons
[1322]234 add InputArray MuonMomentumSmearing/muons
[1260]235 set OutputArray tracks
236}
237
238#############
239# Calorimeter
240#############
241
242module Calorimeter Calorimeter {
243 set ParticleInputArray ParticlePropagator/stableParticles
244 set TrackInputArray TrackMerger/tracks
245
246 set TowerOutputArray towers
247 set PhotonOutputArray photons
248
249 set EFlowTrackOutputArray eflowTracks
[1357]250 set EFlowPhotonOutputArray eflowPhotons
251 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
[1378]252
[1260]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.9420^2 + energy*0.075^2)}
310}
311
312##########################
313# Track pile-up subtractor
314##########################
315
316module TrackPileUpSubtractor TrackPileUpSubtractor {
317# add InputArray InputArray OutputArray
318 add InputArray Calorimeter/eflowTracks eflowTracks
319 add InputArray ElectronEnergySmearing/electrons electrons
320 add InputArray MuonMomentumSmearing/muons muons
321
322 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
323 # Z vertex resolution in m
324 set ZVertexResolution 0.0001
325}
326
327####################
[1357]328# Neutral tower merger
329####################
330
331module Merger NeutralTowerMerger {
332# add InputArray InputArray
333 add InputArray Calorimeter/eflowPhotons
334 add InputArray Calorimeter/eflowNeutralHadrons
335 set OutputArray eflowTowers
[1378]336}
[1357]337
338
339####################
[1260]340# Energy flow merger
341####################
[1378]342
[1260]343module Merger EFlowMerger {
344# add InputArray InputArray
[1357]345 add InputArray Calorimeter/eflowTracks
346 add InputArray Calorimeter/eflowPhotons
347 add InputArray Calorimeter/eflowNeutralHadrons
[1260]348 set OutputArray eflow
[1378]349}
[1260]350
[1357]351
[1260]352#############
353# Rho pile-up
354#############
355
356module FastJetFinder Rho {
357 set InputArray Calorimeter/towers
358
359 set ComputeRho true
360 set RhoOutputArray rho
361
362 # 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
363 set AreaAlgorithm 5
364
365 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
366 set JetAlgorithm 4
367 set ParameterR 0.6
368 set GhostEtaMax 5.0
369 set RhoEtaMax 5.0
370
[1322]371 add RhoEtaRange 0.0 5.0
372
[1260]373 set JetPTMin 0.0
374}
375
376#####################
377# MC truth jet finder
378#####################
379
380module FastJetFinder GenJetFinder {
381 set InputArray Delphes/stableParticles
382
383 set OutputArray jets
384
385 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
386 set JetAlgorithm 6
387 set ParameterR 0.6
388
389 set JetPTMin 20.0
390}
391
392############
393# Jet finder
394############
395
396module FastJetFinder FastJetFinder {
397 set InputArray Calorimeter/towers
398
399 set OutputArray jets
400
401 # 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
402 set AreaAlgorithm 5
403
404 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
405 set JetAlgorithm 6
406 set ParameterR 0.6
407
408 set JetPTMin 20.0
409}
410
411###########################
412# Jet Pile-Up Subtraction
413###########################
414
415module JetPileUpSubtractor JetPileUpSubtractor {
416 set JetInputArray FastJetFinder/jets
417 set RhoInputArray Rho/rho
418
419 set OutputArray jets
420
421 set JetPTMin 20.0
422}
423
424##################
425# Jet Energy Scale
426##################
427
428module EnergyScale JetEnergyScale {
429 set InputArray JetPileUpSubtractor/jets
430 set OutputArray jets
431
432 # scale formula for jets
433 set ScaleFormula {1.0}
434}
435
436###################
437# Photon efficiency
438###################
439
440module Efficiency PhotonEfficiency {
441 set InputArray Calorimeter/photons
442 set OutputArray photons
443
444 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
445
446 # efficiency formula for photons
447 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
448 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
449 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
450 (abs(eta) > 2.5) * (0.00)}
451}
452
453##################
454# Photon isolation
455##################
456
457module Isolation PhotonIsolation {
458 set CandidateInputArray PhotonEfficiency/photons
459 set IsolationInputArray EFlowMerger/eflow
460 set RhoInputArray Rho/rho
461
462 set OutputArray photons
463
464 set DeltaRMax 0.5
465
466 set PTMin 0.5
467
468 set PTRatioMax 0.1
469}
470
471#####################
472# Electron efficiency
473#####################
474
475module Efficiency ElectronEfficiency {
476 set InputArray TrackPileUpSubtractor/electrons
477 set OutputArray electrons
478
479 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
480
481 # efficiency formula for electrons
482 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
483 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
484 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
485 (abs(eta) > 2.5) * (0.00)}
486}
487
488####################
489# Electron isolation
490####################
491
492module Isolation ElectronIsolation {
493 set CandidateInputArray ElectronEfficiency/electrons
494 set IsolationInputArray EFlowMerger/eflow
495 set RhoInputArray Rho/rho
496
497 set OutputArray electrons
498
499 set DeltaRMax 0.5
500
501 set PTMin 0.5
502
503 set PTRatioMax 0.1
504}
505
506#################
507# Muon efficiency
508#################
509
510module Efficiency MuonEfficiency {
511 set InputArray TrackPileUpSubtractor/muons
512 set OutputArray muons
513
514 # set EfficiencyFormula {efficiency as a function of eta and pt}
515
516 # efficiency formula for muons
517 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
518 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
519 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) + \
520 (abs(eta) > 2.7) * (0.00)}
521}
522
523################
524# Muon isolation
525################
526
527module Isolation MuonIsolation {
528 set CandidateInputArray MuonEfficiency/muons
529 set IsolationInputArray EFlowMerger/eflow
530 set RhoInputArray Rho/rho
531
532 set OutputArray muons
533
534 set DeltaRMax 0.5
535
536 set PTMin 0.5
537
538 set PTRatioMax 0.1
539}
540
541###################
542# Missing ET merger
543###################
544
545module Merger MissingET {
546# add InputArray InputArray
[1357]547 add InputArray EFlowMerger/eflow
[1260]548 set MomentumOutputArray momentum
549}
550
[1357]551
[1260]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# b-tagging
567###########
568
569module BTagging BTagging {
570 set PartonInputArray Delphes/partons
571 set JetInputArray JetEnergyScale/jets
572
573 set BitNumber 0
574
575 set DeltaR 0.5
576
577 set PartonPTMin 1.0
578
579 set PartonEtaMax 2.5
580
581 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
582 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
583 # gluon's PDG code has the lowest priority
584
585 # default efficiency formula (misidentification rate)
586 add EfficiencyFormula {0} {0.001}
587
588 # efficiency formula for c-jets (misidentification rate)
589 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
590 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
591 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
592 (abs(eta) > 2.5) * (0.000)}
593
594 # efficiency formula for b-jets
595 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \
596 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
597 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
598 (abs(eta) > 2.5) * (0.000)}
599}
600
601module TauTagging TauTagging {
602 set ParticleInputArray Delphes/allParticles
603 set PartonInputArray Delphes/partons
604 set JetInputArray JetEnergyScale/jets
605
606 set DeltaR 0.5
607
608 set TauPTMin 1.0
609
610 set TauEtaMax 2.5
611
612 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
613
614 # default efficiency formula (misidentification rate)
615 add EfficiencyFormula {0} {0.001}
616 # efficiency formula for tau-jets
617 add EfficiencyFormula {15} {0.4}
618}
619
620#####################################################
621# Find uniquely identified photons/electrons/tau/jets
622#####################################################
623
624module UniqueObjectFinder UniqueObjectFinder {
625# earlier arrays take precedence over later ones
626# add InputArray InputArray OutputArray
627 add InputArray PhotonIsolation/photons photons
628 add InputArray ElectronIsolation/electrons electrons
629 add InputArray MuonIsolation/muons muons
630 add InputArray JetEnergyScale/jets jets
631}
632
633##################
634# ROOT tree writer
635##################
636
[1357]637# tracks, towers and eflow objects are not stored by default in the output.
[1378]638# if needed (for jet constituent or other studies), uncomment the relevant
639# "add Branch ..." lines.
[1357]640
[1260]641module TreeWriter TreeWriter {
642# add Branch InputArray BranchName BranchClass
643 add Branch Delphes/allParticles Particle GenParticle
[1357]644
645# add Branch TrackMerger/tracks Track Track
646# add Branch Calorimeter/towers Tower Tower
[1378]647
[1260]648# add Branch Calorimeter/eflowTracks EFlowTrack Track
[1357]649# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
650# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
651
[1260]652 add Branch GenJetFinder/jets GenJet Jet
653 add Branch UniqueObjectFinder/jets Jet Jet
654 add Branch UniqueObjectFinder/electrons Electron Electron
655 add Branch UniqueObjectFinder/photons Photon Photon
656 add Branch UniqueObjectFinder/muons Muon Muon
657 add Branch MissingET/momentum MissingET MissingET
658 add Branch ScalarHT/energy ScalarHT ScalarHT
659 add Branch Rho/rho Rho Rho
[1326]660 add Branch PileUpMerger/vertices Vertex Vertex
[1357]661
[1260]662}
663
[1357]664
665
Note: See TracBrowser for help on using the repository browser.