Fork me on GitHub

source: git/examples/delphes_card_ATLAS_PileUp.tcl@ fc4e460

ImprovedOutputFile Timing dual_readout llp
Last change on this file since fc4e460 was be2222c, checked in by Michele <michele.selvaggi@…>, 10 years ago

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