Fork me on GitHub

source: git/cards/delphes_card_ATLAS_PileUp.tcl@ 61f2714

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 61f2714 was 59d6164, checked in by Michele Selvaggi <michele.selvaggi@…>, 9 years ago

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