MA5SandBox: delphes_card_atlas_susy_2018_04.tcl

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