Fork me on GitHub

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

Last change on this file since 1397 was 1397, checked in by mselvaggi, 6 years ago

added Set RandomSeed? 0 to all data cards

File size: 21.1 KB
Line 
1set RandomSeed 0
2
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
23  NeutralTowerMerger
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
62  set ParticleOutputArray stableParticles
63  set VertexOutputArray vertices
64
65  # pre-generated minbias input file
66  set PileUpFile MinBias.pileup
67
68  # average expected pile up
69  set MeanPileUp 50
70
71  # maximum spread in the beam direction in m
72  set ZVertexSpread 0.10
73
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)
78
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
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
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)}
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
236  add InputArray MuonMomentumSmearing/muons
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
252  set EFlowPhotonOutputArray eflowPhotons
253  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
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*0.075^2)}
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####################
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
338}
339
340
341####################
342# Energy flow merger
343####################
344
345module Merger EFlowMerger {
346# add InputArray InputArray
347  add InputArray Calorimeter/eflowTracks
348  add InputArray Calorimeter/eflowPhotons
349  add InputArray Calorimeter/eflowNeutralHadrons
350  set OutputArray eflow
351}
352
353
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
373  add RhoEtaRange 0.0 5.0
374
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
549  add InputArray EFlowMerger/eflow
550  set MomentumOutputArray momentum
551}
552
553
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
639# tracks, towers and eflow objects are not stored by default in the output.
640# if needed (for jet constituent or other studies), uncomment the relevant
641# "add Branch ..." lines.
642
643module TreeWriter TreeWriter {
644# add Branch InputArray BranchName BranchClass
645  add Branch Delphes/allParticles Particle GenParticle
646
647#  add Branch TrackMerger/tracks Track Track
648#  add Branch Calorimeter/towers Tower Tower
649
650#  add Branch Calorimeter/eflowTracks EFlowTrack Track
651#  add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
652#  add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
653
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
662  add Branch PileUpMerger/vertices Vertex Vertex
663
664}
665
666
667
Note: See TracBrowser for help on using the repository browser.