Fork me on GitHub

source: svn/trunk/examples/delphes_card_ATLAS_PileUp.tcl

Last change on this file was 1397, checked in by Michele Selvaggi, 10 years ago

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