Fork me on GitHub

source: git/cards/delphes_card_CMS.tcl@ 96abea4

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 96abea4 was 998790a, checked in by Michele Selvaggi <michele.selvaggi@…>, 8 years ago

apply muon eff drop in tracking
  • Property mode set to 100644
File size: 20.9 KB
Line 
1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6 ParticlePropagator
7
8 ChargedHadronTrackingEfficiency
9 ElectronTrackingEfficiency
10 MuonTrackingEfficiency
11
12 ChargedHadronMomentumSmearing
13 ElectronMomentumSmearing
14 MuonMomentumSmearing
15
16 TrackMerger
17
18 ECal
19 HCal
20
21 Calorimeter
22 EFlowMerger
23
24 PhotonEfficiency
25 PhotonIsolation
26
27 ElectronFilter
28 ElectronEfficiency
29 ElectronIsolation
30
31 MuonEfficiency
32 MuonIsolation
33
34 MissingET
35
36 NeutrinoFilter
37 GenJetFinder
38 GenMissingET
39
40 FastJetFinder
41
42 JetEnergyScale
43
44 JetFlavorAssociation
45
46 BTagging
47 TauTagging
48
49 UniqueObjectFinder
50
51 ScalarHT
52
53 TreeWriter
54}
55
56#################################
57# Propagate particles in cylinder
58#################################
59
60module ParticlePropagator ParticlePropagator {
61 set InputArray Delphes/stableParticles
62
63 set OutputArray stableParticles
64 set ChargedHadronOutputArray chargedHadrons
65 set ElectronOutputArray electrons
66 set MuonOutputArray muons
67
68 # radius of the magnetic field coverage, in m
69 set Radius 1.29
70 # half-length of the magnetic field coverage, in m
71 set HalfLength 3.00
72
73 # magnetic field
74 set Bz 3.8
75}
76
77####################################
78# Charged hadron tracking efficiency
79####################################
80
81module Efficiency ChargedHadronTrackingEfficiency {
82 set InputArray ParticlePropagator/chargedHadrons
83 set OutputArray chargedHadrons
84
85 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
86
87 # tracking efficiency formula for charged hadrons
88 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
89 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
90 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
91 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
92 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
93 (abs(eta) > 2.5) * (0.00)}
94}
95
96##############################
97# Electron tracking efficiency
98##############################
99
100module Efficiency ElectronTrackingEfficiency {
101 set InputArray ParticlePropagator/electrons
102 set OutputArray electrons
103
104 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
105
106 # tracking efficiency formula for electrons
107 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
108 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
109 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
110 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
111 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
112 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
113 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
114 (abs(eta) > 2.5) * (0.00)}
115}
116
117##########################
118# Muon tracking efficiency
119##########################
120
121module Efficiency MuonTrackingEfficiency {
122 set InputArray ParticlePropagator/muons
123 set OutputArray muons
124
125 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
126
127 # tracking efficiency formula for muons
128 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
129 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
130 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e3) * (0.99) +
131 (abs(eta) <= 1.5) * (pt > 1.0e3 ) * (0.99 * exp(0.5 - pt*5.0e-4)) +
132
133 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
134 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e3) * (0.98) +
135 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3) * (0.98 * exp(0.5 - pt*5.0e-4)) +
136 (abs(eta) > 2.5) * (0.00)}
137}
138
139########################################
140# Momentum resolution for charged tracks
141########################################
142
143module MomentumSmearing ChargedHadronMomentumSmearing {
144 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
145 set OutputArray chargedHadrons
146
147 # set ResolutionFormula {resolution formula as a function of eta and pt}
148
149 # resolution formula for charged hadrons
150 # based on arXiv:1405.6569
151 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
152 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
153 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
154}
155
156###################################
157# Momentum resolution for electrons
158###################################
159
160module MomentumSmearing ElectronMomentumSmearing {
161 set InputArray ElectronTrackingEfficiency/electrons
162 set OutputArray electrons
163
164 # set ResolutionFormula {resolution formula as a function of eta and energy}
165
166 # resolution formula for electrons
167 # based on arXiv:1405.6569
168 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
169 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
170 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
171}
172
173###############################
174# Momentum resolution for muons
175###############################
176
177module MomentumSmearing MuonMomentumSmearing {
178 set InputArray MuonTrackingEfficiency/muons
179 set OutputArray muons
180
181 # set ResolutionFormula {resolution formula as a function of eta and pt}
182
183 # resolution formula for muons
184 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
185 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
186 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
187}
188
189##############
190# Track merger
191##############
192
193module Merger TrackMerger {
194# add InputArray InputArray
195 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
196 add InputArray ElectronMomentumSmearing/electrons
197 add InputArray MuonMomentumSmearing/muons
198 set OutputArray tracks
199}
200
201
202
203#############
204# ECAL
205#############
206
207module SimpleCalorimeter ECal {
208 set ParticleInputArray ParticlePropagator/stableParticles
209 set TrackInputArray TrackMerger/tracks
210
211 set TowerOutputArray ecalTowers
212 set EFlowTrackOutputArray eflowTracks
213 set EFlowTowerOutputArray eflowPhotons
214
215 set IsEcal true
216
217 set EnergyMin 0.5
218 set EnergySignificanceMin 2.0
219
220 set SmearTowerCenter true
221
222 set pi [expr {acos(-1)}]
223
224 # lists of the edges of each tower in eta and phi
225 # each list starts with the lower edge of the first tower
226 # the list ends with the higher edged of the last tower
227
228 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
229
230 set PhiBins {}
231 for {set i -180} {$i <= 180} {incr i} {
232 add PhiBins [expr {$i * $pi/180.0}]
233 }
234
235 # 0.02 unit in eta up to eta = 1.5 (barrel)
236 for {set i -85} {$i <= 86} {incr i} {
237 set eta [expr {$i * 0.0174}]
238 add EtaPhiBins $eta $PhiBins
239 }
240
241 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
242
243 set PhiBins {}
244 for {set i -180} {$i <= 180} {incr i} {
245 add PhiBins [expr {$i * $pi/180.0}]
246 }
247
248 # 0.02 unit in eta up to eta = 3
249 for {set i 1} {$i <= 84} {incr i} {
250 set eta [expr { -2.958 + $i * 0.0174}]
251 add EtaPhiBins $eta $PhiBins
252 }
253
254 for {set i 1} {$i <= 84} {incr i} {
255 set eta [expr { 1.4964 + $i * 0.0174}]
256 add EtaPhiBins $eta $PhiBins
257 }
258
259 # take present CMS granularity for HF
260
261 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
262 set PhiBins {}
263 for {set i -18} {$i <= 18} {incr i} {
264 add PhiBins [expr {$i * $pi/18.0}]
265 }
266
267 foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
268 add EtaPhiBins $eta $PhiBins
269 }
270
271
272 add EnergyFraction {0} {0.0}
273 # energy fractions for e, gamma and pi0
274 add EnergyFraction {11} {1.0}
275 add EnergyFraction {22} {1.0}
276 add EnergyFraction {111} {1.0}
277 # energy fractions for muon, neutrinos and neutralinos
278 add EnergyFraction {12} {0.0}
279 add EnergyFraction {13} {0.0}
280 add EnergyFraction {14} {0.0}
281 add EnergyFraction {16} {0.0}
282 add EnergyFraction {1000022} {0.0}
283 add EnergyFraction {1000023} {0.0}
284 add EnergyFraction {1000025} {0.0}
285 add EnergyFraction {1000035} {0.0}
286 add EnergyFraction {1000045} {0.0}
287 # energy fractions for K0short and Lambda
288 add EnergyFraction {310} {0.3}
289 add EnergyFraction {3122} {0.3}
290
291 # set ResolutionFormula {resolution formula as a function of eta and energy}
292
293 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
294
295 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
296 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
297 set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
298 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (2.16 + 5.6*(abs(eta)-2)^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
299 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
300
301}
302
303
304#############
305# HCAL
306#############
307
308module SimpleCalorimeter HCal {
309 set ParticleInputArray ParticlePropagator/stableParticles
310 set TrackInputArray ECal/eflowTracks
311
312 set TowerOutputArray hcalTowers
313 set EFlowTrackOutputArray eflowTracks
314 set EFlowTowerOutputArray eflowNeutralHadrons
315
316 set IsEcal false
317
318 set EnergyMin 1.0
319 set EnergySignificanceMin 2.0
320
321 set SmearTowerCenter true
322
323 set pi [expr {acos(-1)}]
324
325 # lists of the edges of each tower in eta and phi
326 # each list starts with the lower edge of the first tower
327 # the list ends with the higher edged of the last tower
328
329 # 5 degrees towers
330 set PhiBins {}
331 for {set i -36} {$i <= 36} {incr i} {
332 add PhiBins [expr {$i * $pi/36.0}]
333 }
334 foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
335 add EtaPhiBins $eta $PhiBins
336 }
337
338 # 10 degrees towers
339 set PhiBins {}
340 for {set i -18} {$i <= 18} {incr i} {
341 add PhiBins [expr {$i * $pi/18.0}]
342 }
343 foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
344 add EtaPhiBins $eta $PhiBins
345 }
346
347 # 20 degrees towers
348 set PhiBins {}
349 for {set i -9} {$i <= 9} {incr i} {
350 add PhiBins [expr {$i * $pi/9.0}]
351 }
352 foreach eta {-5 -4.7 -4.525 4.7 5} {
353 add EtaPhiBins $eta $PhiBins
354 }
355
356 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
357 add EnergyFraction {0} {1.0}
358 # energy fractions for e, gamma and pi0
359 add EnergyFraction {11} {0.0}
360 add EnergyFraction {22} {0.0}
361 add EnergyFraction {111} {0.0}
362 # energy fractions for muon, neutrinos and neutralinos
363 add EnergyFraction {12} {0.0}
364 add EnergyFraction {13} {0.0}
365 add EnergyFraction {14} {0.0}
366 add EnergyFraction {16} {0.0}
367 add EnergyFraction {1000022} {0.0}
368 add EnergyFraction {1000023} {0.0}
369 add EnergyFraction {1000025} {0.0}
370 add EnergyFraction {1000035} {0.0}
371 add EnergyFraction {1000045} {0.0}
372 # energy fractions for K0short and Lambda
373 add EnergyFraction {310} {0.7}
374 add EnergyFraction {3122} {0.7}
375
376 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
377 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
378 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
379
380}
381
382
383#################
384# Electron filter
385#################
386
387module PdgCodeFilter ElectronFilter {
388 set InputArray HCal/eflowTracks
389 set OutputArray electrons
390 set Invert true
391 add PdgCode {11}
392 add PdgCode {-11}
393}
394
395###################################################
396# Tower Merger (in case not using e-flow algorithm)
397###################################################
398
399module Merger Calorimeter {
400# add InputArray InputArray
401 add InputArray ECal/ecalTowers
402 add InputArray HCal/hcalTowers
403 set OutputArray towers
404}
405
406
407
408####################
409# Energy flow merger
410####################
411
412module Merger EFlowMerger {
413# add InputArray InputArray
414 add InputArray HCal/eflowTracks
415 add InputArray ECal/eflowPhotons
416 add InputArray HCal/eflowNeutralHadrons
417 set OutputArray eflow
418}
419
420###################
421# Photon efficiency
422###################
423
424module Efficiency PhotonEfficiency {
425 set InputArray ECal/eflowPhotons
426 set OutputArray photons
427
428 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
429
430 # efficiency formula for photons
431 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
432 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
433 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
434 (abs(eta) > 2.5) * (0.00)}
435}
436
437##################
438# Photon isolation
439##################
440
441module Isolation PhotonIsolation {
442 set CandidateInputArray PhotonEfficiency/photons
443 set IsolationInputArray EFlowMerger/eflow
444
445 set OutputArray photons
446
447 set DeltaRMax 0.5
448
449 set PTMin 0.5
450
451 set PTRatioMax 0.12
452}
453
454
455#####################
456# Electron efficiency
457#####################
458
459module Efficiency ElectronEfficiency {
460 set InputArray ElectronFilter/electrons
461 set OutputArray electrons
462
463 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
464
465 # efficiency formula for electrons
466 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
467 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
468 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
469 (abs(eta) > 2.5) * (0.00)}
470}
471
472####################
473# Electron isolation
474####################
475
476module Isolation ElectronIsolation {
477 set CandidateInputArray ElectronEfficiency/electrons
478 set IsolationInputArray EFlowMerger/eflow
479
480 set OutputArray electrons
481
482 set DeltaRMax 0.5
483
484 set PTMin 0.5
485
486 set PTRatioMax 0.12
487}
488
489#################
490# Muon efficiency
491#################
492
493module Efficiency MuonEfficiency {
494 set InputArray MuonMomentumSmearing/muons
495 set OutputArray muons
496
497 # set EfficiencyFormula {efficiency as a function of eta and pt}
498
499 # efficiency formula for muons
500 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
501 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
502 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0) * (0.95) +
503 (abs(eta) > 2.4) * (0.00)}
504}
505
506################
507# Muon isolation
508################
509
510module Isolation MuonIsolation {
511 set CandidateInputArray MuonEfficiency/muons
512 set IsolationInputArray EFlowMerger/eflow
513
514 set OutputArray muons
515
516 set DeltaRMax 0.5
517
518 set PTMin 0.5
519
520 set PTRatioMax 0.25
521}
522
523###################
524# Missing ET merger
525###################
526
527module Merger MissingET {
528# add InputArray InputArray
529 add InputArray EFlowMerger/eflow
530 set MomentumOutputArray momentum
531}
532
533##################
534# Scalar HT merger
535##################
536
537module Merger ScalarHT {
538# add InputArray InputArray
539 add InputArray UniqueObjectFinder/jets
540 add InputArray UniqueObjectFinder/electrons
541 add InputArray UniqueObjectFinder/photons
542 add InputArray UniqueObjectFinder/muons
543 set EnergyOutputArray energy
544}
545
546
547#####################
548# Neutrino Filter
549#####################
550
551module PdgCodeFilter NeutrinoFilter {
552
553 set InputArray Delphes/stableParticles
554 set OutputArray filteredParticles
555
556 set PTMin 0.0
557
558 add PdgCode {12}
559 add PdgCode {14}
560 add PdgCode {16}
561 add PdgCode {-12}
562 add PdgCode {-14}
563 add PdgCode {-16}
564
565}
566
567
568#####################
569# MC truth jet finder
570#####################
571
572module FastJetFinder GenJetFinder {
573 set InputArray NeutrinoFilter/filteredParticles
574
575 set OutputArray jets
576
577 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
578 set JetAlgorithm 6
579 set ParameterR 0.5
580
581 set JetPTMin 20.0
582}
583
584#########################
585# Gen Missing ET merger
586########################
587
588module Merger GenMissingET {
589# add InputArray InputArray
590 add InputArray NeutrinoFilter/filteredParticles
591 set MomentumOutputArray momentum
592}
593
594
595
596############
597# Jet finder
598############
599
600module FastJetFinder FastJetFinder {
601# set InputArray Calorimeter/towers
602 set InputArray EFlowMerger/eflow
603
604 set OutputArray jets
605
606 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
607 set JetAlgorithm 6
608 set ParameterR 0.5
609
610 set JetPTMin 20.0
611}
612
613##################
614# Jet Energy Scale
615##################
616
617module EnergyScale JetEnergyScale {
618 set InputArray FastJetFinder/jets
619 set OutputArray jets
620
621 # scale formula for jets
622 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
623}
624
625########################
626# Jet Flavor Association
627########################
628
629module JetFlavorAssociation JetFlavorAssociation {
630
631 set PartonInputArray Delphes/partons
632 set ParticleInputArray Delphes/allParticles
633 set ParticleLHEFInputArray Delphes/allParticlesLHEF
634 set JetInputArray JetEnergyScale/jets
635
636 set DeltaR 0.5
637 set PartonPTMin 1.0
638 set PartonEtaMax 2.5
639
640}
641
642###########
643# b-tagging
644###########
645
646module BTagging BTagging {
647 set JetInputArray JetEnergyScale/jets
648
649 set BitNumber 0
650
651 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
652 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
653 # gluon's PDG code has the lowest priority
654
655 # based on arXiv:1211.4462
656
657 # default efficiency formula (misidentification rate)
658 add EfficiencyFormula {0} {0.01+0.000038*pt}
659
660 # efficiency formula for c-jets (misidentification rate)
661 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
662
663 # efficiency formula for b-jets
664 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
665}
666
667#############
668# tau-tagging
669#############
670
671module TauTagging TauTagging {
672 set ParticleInputArray Delphes/allParticles
673 set PartonInputArray Delphes/partons
674 set JetInputArray JetEnergyScale/jets
675
676 set DeltaR 0.5
677
678 set TauPTMin 1.0
679
680 set TauEtaMax 2.5
681
682 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
683
684 # default efficiency formula (misidentification rate)
685 add EfficiencyFormula {0} {0.01}
686 # efficiency formula for tau-jets
687 add EfficiencyFormula {15} {0.6}
688}
689
690#####################################################
691# Find uniquely identified photons/electrons/tau/jets
692#####################################################
693
694module UniqueObjectFinder UniqueObjectFinder {
695# earlier arrays take precedence over later ones
696# add InputArray InputArray OutputArray
697 add InputArray PhotonIsolation/photons photons
698 add InputArray ElectronIsolation/electrons electrons
699 add InputArray MuonIsolation/muons muons
700 add InputArray JetEnergyScale/jets jets
701}
702
703##################
704# ROOT tree writer
705##################
706
707# tracks, towers and eflow objects are not stored by default in the output.
708# if needed (for jet constituent or other studies), uncomment the relevant
709# "add Branch ..." lines.
710
711module TreeWriter TreeWriter {
712# add Branch InputArray BranchName BranchClass
713 add Branch Delphes/allParticles Particle GenParticle
714
715 add Branch TrackMerger/tracks Track Track
716 add Branch Calorimeter/towers Tower Tower
717
718 add Branch HCal/eflowTracks EFlowTrack Track
719 add Branch ECal/eflowPhotons EFlowPhoton Tower
720 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
721
722 add Branch GenJetFinder/jets GenJet Jet
723 add Branch GenMissingET/momentum GenMissingET MissingET
724
725 add Branch UniqueObjectFinder/jets Jet Jet
726 add Branch UniqueObjectFinder/electrons Electron Electron
727 add Branch UniqueObjectFinder/photons Photon Photon
728 add Branch UniqueObjectFinder/muons Muon Muon
729 add Branch MissingET/momentum MissingET MissingET
730 add Branch ScalarHT/energy ScalarHT ScalarHT
731}
Note: See TracBrowser for help on using the repository browser.