Fork me on GitHub

source: git/cards/delphes_card_ATLAS_PileUp.tcl@ 4a0d9d5

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 4a0d9d5 was df5084b, checked in by Alexandre Mertens <alexandre.mertens@…>, 9 years ago

updating tracking pt resolution

  • Property mode set to 100644
File size: 21.1 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 # based on arXiv:1405.6569
185 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.5e-4^2) +
186 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*2.5e-4^2) +
187 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*5.5e-4^2)}
188}
189
190###################################
191# Momentum resolution for electrons
192###################################
193
194module MomentumSmearing ElectronMomentumSmearing {
195 set InputArray ElectronTrackingEfficiency/electrons
196 set OutputArray electrons
197
198 # set ResolutionFormula {resolution formula as a function of eta and energy}
199
200 # resolution formula for electrons
201 # based on arXiv:1405.6569
202 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
203 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
204 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
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) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
219 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
220 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
221}
222
223##############
224# Track merger
225##############
226
227module Merger TrackMerger {
228# add InputArray InputArray
229 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
230 add InputArray ElectronMomentumSmearing/electrons
231 add InputArray MuonMomentumSmearing/muons
232 set OutputArray tracks
233}
234
235#############
236# Calorimeter
237#############
238
239module Calorimeter Calorimeter {
240 set ParticleInputArray ParticlePropagator/stableParticles
241 set TrackInputArray TrackMerger/tracks
242
243 set TowerOutputArray towers
244 set PhotonOutputArray photons
245
246 set EFlowTrackOutputArray eflowTracks
247 set EFlowPhotonOutputArray eflowPhotons
248 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
249
250 set ECalEnergyMin 0.5
251 set HCalEnergyMin 1.0
252
253 set ECalEnergySignificanceMin 1.0
254 set HCalEnergySignificanceMin 1.0
255
256 set SmearTowerCenter true
257
258 set pi [expr {acos(-1)}]
259
260 # lists of the edges of each tower in eta and phi
261 # each list starts with the lower edge of the first tower
262 # the list ends with the higher edged of the last tower
263
264 # 10 degrees towers
265 set PhiBins {}
266 for {set i -18} {$i <= 18} {incr i} {
267 add PhiBins [expr {$i * $pi/18.0}]
268 }
269 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} {
270 add EtaPhiBins $eta $PhiBins
271 }
272
273 # 20 degrees towers
274 set PhiBins {}
275 for {set i -9} {$i <= 9} {incr i} {
276 add PhiBins [expr {$i * $pi/9.0}]
277 }
278 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} {
279 add EtaPhiBins $eta $PhiBins
280 }
281
282 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
283 add EnergyFraction {0} {0.0 1.0}
284 # energy fractions for e, gamma and pi0
285 add EnergyFraction {11} {1.0 0.0}
286 add EnergyFraction {22} {1.0 0.0}
287 add EnergyFraction {111} {1.0 0.0}
288 # energy fractions for muon, neutrinos and neutralinos
289 add EnergyFraction {12} {0.0 0.0}
290 add EnergyFraction {13} {0.0 0.0}
291 add EnergyFraction {14} {0.0 0.0}
292 add EnergyFraction {16} {0.0 0.0}
293 add EnergyFraction {1000022} {0.0 0.0}
294 add EnergyFraction {1000023} {0.0 0.0}
295 add EnergyFraction {1000025} {0.0 0.0}
296 add EnergyFraction {1000035} {0.0 0.0}
297 add EnergyFraction {1000045} {0.0 0.0}
298 # energy fractions for K0short and Lambda
299 add EnergyFraction {310} {0.3 0.7}
300 add EnergyFraction {3122} {0.3 0.7}
301
302 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
303 # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
304 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
305 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
306 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) +
307 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
308
309 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
310 # http://arxiv.org/pdf/hep-ex/0004009v1
311 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
312 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
313 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
314 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
315}
316
317#################
318# Electron filter
319#################
320
321module PdgCodeFilter ElectronFilter {
322 set InputArray Calorimeter/eflowTracks
323 set OutputArray electrons
324 set Invert true
325 add PdgCode {11}
326 add PdgCode {-11}
327}
328
329##########################
330# Track pile-up subtractor
331##########################
332
333module TrackPileUpSubtractor TrackPileUpSubtractor {
334# add InputArray InputArray OutputArray
335 add InputArray Calorimeter/eflowTracks eflowTracks
336 add InputArray ElectronFilter/electrons electrons
337 add InputArray MuonMomentumSmearing/muons muons
338
339 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
340 # Z vertex resolution in m
341 set ZVertexResolution 0.0001
342}
343
344####################
345# Neutral tower merger
346####################
347
348module Merger NeutralTowerMerger {
349# add InputArray InputArray
350 add InputArray Calorimeter/eflowPhotons
351 add InputArray Calorimeter/eflowNeutralHadrons
352 set OutputArray eflowTowers
353}
354
355##################################
356# Energy flow merger (all tracks)
357##################################
358
359module Merger EFlowMergerAllTracks {
360# add InputArray InputArray
361 add InputArray TrackMerger/tracks
362 add InputArray Calorimeter/eflowPhotons
363 add InputArray Calorimeter/eflowNeutralHadrons
364 set OutputArray eflow
365}
366
367
368####################
369# Energy flow merger
370####################
371
372module Merger EFlowMerger {
373# add InputArray InputArray
374 add InputArray Calorimeter/eflowTracks
375 add InputArray Calorimeter/eflowPhotons
376 add InputArray Calorimeter/eflowNeutralHadrons
377 set OutputArray eflow
378}
379
380#############
381# Rho pile-up
382#############
383
384module FastJetGridMedianEstimator Rho {
385
386 set InputArray Calorimeter/towers
387 set RhoOutputArray rho
388
389 # add GridRange rapmin rapmax drap dphi
390 # rapmin - the minimum rapidity extent of the grid
391 # rapmax - the maximum rapidity extent of the grid
392 # drap - the grid spacing in rapidity
393 # dphi - the grid spacing in azimuth
394
395 add GridRange -5.0 -2.5 1.0 1.0
396 add GridRange -2.5 2.5 0.5 0.5
397 add GridRange 2.5 5.0 1.0 1.0
398
399}
400
401
402#####################
403# Neutrino Filter
404#####################
405
406module PdgCodeFilter NeutrinoFilter {
407
408 set InputArray Delphes/stableParticles
409 set OutputArray filteredParticles
410
411 set PTMin 0.0
412
413 add PdgCode {12}
414 add PdgCode {14}
415 add PdgCode {16}
416 add PdgCode {-12}
417 add PdgCode {-14}
418 add PdgCode {-16}
419
420}
421
422#####################
423# MC truth jet finder
424#####################
425
426module FastJetFinder GenJetFinder {
427 set InputArray NeutrinoFilter/filteredParticles
428
429 set OutputArray jets
430
431 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
432 set JetAlgorithm 6
433 set ParameterR 0.6
434
435 set JetPTMin 20.0
436}
437
438#########################
439# Gen Missing ET merger
440########################
441
442module Merger GenMissingET {
443# add InputArray InputArray
444 add InputArray NeutrinoFilter/filteredParticles
445 set MomentumOutputArray momentum
446}
447
448
449############
450# Jet finder
451############
452
453module FastJetFinder FastJetFinder {
454 set InputArray Calorimeter/towers
455
456 set OutputArray jets
457
458 # 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
459 set AreaAlgorithm 5
460
461 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
462 set JetAlgorithm 6
463 set ParameterR 0.6
464
465 set JetPTMin 20.0
466}
467
468###########################
469# Jet Pile-Up Subtraction
470###########################
471
472module JetPileUpSubtractor JetPileUpSubtractor {
473 set JetInputArray FastJetFinder/jets
474 set RhoInputArray Rho/rho
475
476 set OutputArray jets
477
478 set JetPTMin 20.0
479}
480
481##################
482# Jet Energy Scale
483##################
484
485module EnergyScale JetEnergyScale {
486 set InputArray JetPileUpSubtractor/jets
487 set OutputArray jets
488
489 # scale formula for jets
490 set ScaleFormula {1.0}
491}
492
493###################
494# Photon efficiency
495###################
496
497module Efficiency PhotonEfficiency {
498 set InputArray Calorimeter/eflowPhotons
499 set OutputArray photons
500
501 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
502
503 # efficiency formula for photons
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# Photon isolation
512##################
513
514module Isolation PhotonIsolation {
515 set CandidateInputArray PhotonEfficiency/photons
516 set IsolationInputArray EFlowMergerAllTracks/eflow
517 set RhoInputArray Rho/rho
518
519 set OutputArray photons
520
521 set DeltaRMax 0.5
522
523 set PTMin 0.5
524
525 set PTRatioMax 0.12
526}
527
528#####################
529# Electron efficiency
530#####################
531
532module Efficiency ElectronEfficiency {
533 set InputArray TrackPileUpSubtractor/electrons
534 set OutputArray electrons
535
536 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
537
538 # efficiency formula for electrons
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.5) * (pt > 10.0) * (0.85) +
542 (abs(eta) > 2.5) * (0.00)}
543}
544
545####################
546# Electron isolation
547####################
548
549module Isolation ElectronIsolation {
550 set CandidateInputArray ElectronEfficiency/electrons
551 set IsolationInputArray EFlowMergerAllTracks/eflow
552 set RhoInputArray Rho/rho
553
554 set OutputArray electrons
555
556 set DeltaRMax 0.5
557
558 set PTMin 0.5
559
560 set PTRatioMax 0.12
561}
562
563#################
564# Muon efficiency
565#################
566
567module Efficiency MuonEfficiency {
568 set InputArray TrackPileUpSubtractor/muons
569 set OutputArray muons
570
571 # set EfficiencyFormula {efficiency as a function of eta and pt}
572
573 # efficiency formula for muons
574 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
575 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
576 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) +
577 (abs(eta) > 2.7) * (0.00)}
578}
579
580################
581# Muon isolation
582################
583
584module Isolation MuonIsolation {
585 set CandidateInputArray MuonEfficiency/muons
586 set IsolationInputArray EFlowMergerAllTracks/eflow
587 set RhoInputArray Rho/rho
588
589 set OutputArray muons
590
591 set DeltaRMax 0.5
592
593 set PTMin 0.5
594
595 set PTRatioMax 0.25
596}
597
598###################
599# Missing ET merger
600###################
601
602module Merger MissingET {
603# add InputArray InputArray
604 add InputArray EFlowMergerAllTracks/eflow
605 set MomentumOutputArray momentum
606}
607
608
609##################
610# Scalar HT merger
611##################
612
613module Merger ScalarHT {
614# add InputArray InputArray
615 add InputArray UniqueObjectFinder/jets
616 add InputArray UniqueObjectFinder/electrons
617 add InputArray UniqueObjectFinder/photons
618 add InputArray UniqueObjectFinder/muons
619 set EnergyOutputArray energy
620}
621
622########################
623# Jet Flavor Association
624########################
625
626module JetFlavorAssociation JetFlavorAssociation {
627
628 set PartonInputArray Delphes/partons
629 set ParticleInputArray Delphes/allParticles
630 set ParticleLHEFInputArray Delphes/allParticlesLHEF
631 set JetInputArray JetEnergyScale/jets
632
633 set DeltaR 0.5
634 set PartonPTMin 1.0
635 set PartonEtaMax 2.5
636
637}
638
639###########
640# b-tagging
641###########
642
643module BTagging BTagging {
644 set JetInputArray JetEnergyScale/jets
645
646 set BitNumber 0
647
648 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
649 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
650 # gluon's PDG code has the lowest priority
651
652 # based on ATL-PHYS-PUB-2015-022
653
654 # default efficiency formula (misidentification rate)
655 add EfficiencyFormula {0} {0.002+7.3e-06*pt}
656
657 # efficiency formula for c-jets (misidentification rate)
658 add EfficiencyFormula {4} {0.20*tanh(0.02*pt)*(1/(1+0.0034*pt))}
659
660 # efficiency formula for b-jets
661 add EfficiencyFormula {5} {0.80*tanh(0.003*pt)*(30/(1+0.086*pt))}
662}
663
664#############
665# tau-tagging
666#############
667
668module TrackCountingTauTagging TauTagging {
669
670 set ParticleInputArray Delphes/allParticles
671 set PartonInputArray Delphes/partons
672 set TrackInputArray TrackMerger/tracks
673 set JetInputArray JetEnergyScale/jets
674
675 set DeltaR 0.2
676 set DeltaRTrack 0.2
677
678 set TrackPTMin 1.0
679
680 set TauPTMin 1.0
681 set TauEtaMax 2.5
682
683 # instructions: {n-prongs} {eff}
684
685 # 1 - one prong efficiency
686 # 2 - two or more efficiency
687 # -1 - one prong mistag rate
688 # -2 - two or more mistag rate
689
690 set BitNumber 0
691
692 # taken from ATL-PHYS-PUB-2015-045 (medium working point)
693 add EfficiencyFormula {1} {0.70}
694 add EfficiencyFormula {2} {0.60}
695 add EfficiencyFormula {-1} {0.02}
696 add EfficiencyFormula {-2} {0.01}
697
698}
699
700#####################################################
701# Find uniquely identified photons/electrons/tau/jets
702#####################################################
703
704module UniqueObjectFinder UniqueObjectFinder {
705# earlier arrays take precedence over later ones
706# add InputArray InputArray OutputArray
707 add InputArray PhotonIsolation/photons photons
708 add InputArray ElectronIsolation/electrons electrons
709 add InputArray MuonIsolation/muons muons
710 add InputArray JetEnergyScale/jets jets
711}
712
713##################
714# ROOT tree writer
715##################
716
717# tracks, towers and eflow objects are not stored by default in the output.
718# if needed (for jet constituent or other studies), uncomment the relevant
719# "add Branch ..." lines.
720
721module TreeWriter TreeWriter {
722# add Branch InputArray BranchName BranchClass
723 add Branch Delphes/allParticles Particle GenParticle
724
725# add Branch TrackMerger/tracks Track Track
726 add Branch Calorimeter/towers Tower Tower
727
728# add Branch Calorimeter/eflowTracks EFlowTrack Track
729# add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
730# add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
731
732 add Branch GenJetFinder/jets GenJet Jet
733 add Branch GenMissingET/momentum GenMissingET MissingET
734
735 add Branch UniqueObjectFinder/jets Jet Jet
736 add Branch UniqueObjectFinder/electrons Electron Electron
737 add Branch UniqueObjectFinder/photons Photon Photon
738 add Branch UniqueObjectFinder/muons Muon Muon
739 add Branch MissingET/momentum MissingET MissingET
740 add Branch ScalarHT/energy ScalarHT ScalarHT
741 add Branch Rho/rho Rho Rho
742 add Branch PileUpMerger/vertices Vertex Vertex
743
744}
Note: See TracBrowser for help on using the repository browser.