Fork me on GitHub

Ticket #725: test.tcl

File test.tcl, 21.0 KB (added by Michele Selvaggi, 9 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 ElectronEnergySmearing
14 MuonMomentumSmearing
15
16 TrackMerger
17 Calorimeter
18 EFlowMerger
19
20 PhotonEfficiency
21 PhotonIsolation
22
23 ElectronEfficiency
24 ElectronIsolation
25
26 MuonEfficiency
27 MuonIsolation
28
29 MissingET
30
31 NeutrinoFilter
32 GenJetFinder
33 FastJetFinder
34
35 JetFakeParticle
36 ElectronMerger
37 MuonMerger
38 PhotonMerger
39
40 JetEnergyScale
41
42 JetFlavorAssociation
43
44 BTagging
45 TauTagging
46
47 UniqueObjectFinder
48
49 ScalarHT
50
51 TreeWriter
52}
53
54#################################
55# Propagate particles in cylinder
56#################################
57
58module ParticlePropagator ParticlePropagator {
59 set InputArray Delphes/stableParticles
60
61 set OutputArray stableParticles
62 set ChargedHadronOutputArray chargedHadrons
63 set ElectronOutputArray electrons
64 set MuonOutputArray muons
65
66 # radius of the magnetic field coverage, in m
67 set Radius 1.29
68 # half-length of the magnetic field coverage, in m
69 set HalfLength 3.00
70
71 # magnetic field
72 set Bz 3.8
73}
74
75####################################
76# Charged hadron tracking efficiency
77####################################
78
79module Efficiency ChargedHadronTrackingEfficiency {
80 set InputArray ParticlePropagator/chargedHadrons
81 set OutputArray chargedHadrons
82
83 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
84
85 # tracking efficiency formula for charged hadrons
86 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
87 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
88 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
89 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
90 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
91 (abs(eta) > 2.5) * (0.00)}
92}
93
94##############################
95# Electron tracking efficiency
96##############################
97
98module Efficiency ElectronTrackingEfficiency {
99 set InputArray ParticlePropagator/electrons
100 set OutputArray electrons
101
102 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
103
104 # tracking efficiency formula for electrons
105 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
106 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
107 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
108 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
109 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
110 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
111 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
112 (abs(eta) > 2.5) * (0.00)}
113}
114
115##########################
116# Muon tracking efficiency
117##########################
118
119module Efficiency MuonTrackingEfficiency {
120 set InputArray ParticlePropagator/muons
121 set OutputArray muons
122
123 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
124
125 # tracking efficiency formula for muons
126 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
127 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
128 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) +
129 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
130 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) +
131 (abs(eta) > 2.5) * (0.00)}
132}
133
134########################################
135# Momentum resolution for charged tracks
136########################################
137
138module MomentumSmearing ChargedHadronMomentumSmearing {
139 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
140 set OutputArray chargedHadrons
141
142 # set ResolutionFormula {resolution formula as a function of eta and pt}
143
144 # resolution formula for charged hadrons
145 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) +
146 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) +
147 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) +
148 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) +
149 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) +
150 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) +
151 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) +
152 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
153}
154
155#################################
156# Energy resolution for electrons
157#################################
158
159module EnergySmearing ElectronEnergySmearing {
160 set InputArray ElectronTrackingEfficiency/electrons
161 set OutputArray electrons
162
163 # set ResolutionFormula {resolution formula as a function of eta and energy}
164
165 # resolution formula for electrons
166 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.0e1) * (energy*0.0225) +
167 (abs(eta) <= 2.5) * (energy > 2.0e1) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
168 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
169 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
170
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 && pt <= 5.0) * (0.02) +
185 (abs(eta) <= 0.5) * (pt > 5.0 && pt <= 1.0e2) * (0.015) +
186 (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) +
187 (abs(eta) <= 0.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) +
188 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1 && pt <= 5.0) * (0.03) +
189 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0 && pt <= 1.0e2) * (0.02) +
190 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) +
191 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4) +
192 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 5.0) * (0.04) +
193 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0 && pt <= 1.0e2) * (0.035) +
194 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) +
195 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05 + pt*1.e-4)}
196}
197
198##############
199# Track merger
200##############
201
202module Merger TrackMerger {
203# add InputArray InputArray
204 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
205 add InputArray ElectronEnergySmearing/electrons
206 add InputArray MuonMomentumSmearing/muons
207 set OutputArray tracks
208}
209
210#############
211# Calorimeter
212#############
213
214module Calorimeter Calorimeter {
215 set ParticleInputArray ParticlePropagator/stableParticles
216 set TrackInputArray TrackMerger/tracks
217
218 set TowerOutputArray towers
219 set PhotonOutputArray photons
220
221 set EFlowTrackOutputArray eflowTracks
222 set EFlowPhotonOutputArray eflowPhotons
223 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
224
225 set ECalEnergyMin 0.5
226 set HCalEnergyMin 1.0
227
228 set ECalEnergySignificanceMin 1.0
229 set HCalEnergySignificanceMin 1.0
230
231 set SmearTowerCenter true
232
233 set pi [expr {acos(-1)}]
234
235 # lists of the edges of each tower in eta and phi
236 # each list starts with the lower edge of the first tower
237 # the list ends with the higher edged of the last tower
238
239 # 5 degrees towers
240 set PhiBins {}
241 for {set i -36} {$i <= 36} {incr i} {
242 add PhiBins [expr {$i * $pi/36.0}]
243 }
244 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} {
245 add EtaPhiBins $eta $PhiBins
246 }
247
248 # 10 degrees towers
249 set PhiBins {}
250 for {set i -18} {$i <= 18} {incr i} {
251 add PhiBins [expr {$i * $pi/18.0}]
252 }
253 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} {
254 add EtaPhiBins $eta $PhiBins
255 }
256
257 # 20 degrees towers
258 set PhiBins {}
259 for {set i -9} {$i <= 9} {incr i} {
260 add PhiBins [expr {$i * $pi/9.0}]
261 }
262 foreach eta {-5 -4.7 -4.525 4.7 5} {
263 add EtaPhiBins $eta $PhiBins
264 }
265
266 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
267 add EnergyFraction {0} {0.0 1.0}
268 # energy fractions for e, gamma and pi0
269 add EnergyFraction {11} {1.0 0.0}
270 add EnergyFraction {22} {1.0 0.0}
271 add EnergyFraction {111} {1.0 0.0}
272 # energy fractions for muon, neutrinos and neutralinos
273 add EnergyFraction {12} {0.0 0.0}
274 add EnergyFraction {13} {0.0 0.0}
275 add EnergyFraction {14} {0.0 0.0}
276 add EnergyFraction {16} {0.0 0.0}
277 add EnergyFraction {1000022} {0.0 0.0}
278 add EnergyFraction {1000023} {0.0 0.0}
279 add EnergyFraction {1000025} {0.0 0.0}
280 add EnergyFraction {1000035} {0.0 0.0}
281 add EnergyFraction {1000045} {0.0 0.0}
282 # energy fractions for K0short and Lambda
283 add EnergyFraction {310} {0.3 0.7}
284 add EnergyFraction {3122} {0.3 0.7}
285
286 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
287 set ECalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
288 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
289
290 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
291 set HCalResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
292 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
293}
294
295####################
296# Energy flow merger
297####################
298
299module Merger EFlowMerger {
300# add InputArray InputArray
301 add InputArray Calorimeter/eflowTracks
302 add InputArray Calorimeter/eflowPhotons
303 add InputArray Calorimeter/eflowNeutralHadrons
304 set OutputArray eflow
305}
306
307###################
308# Photon efficiency
309###################
310
311module Efficiency PhotonEfficiency {
312 set InputArray Calorimeter/eflowPhotons
313 set OutputArray photons
314
315 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
316
317 # efficiency formula for photons
318 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
319 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
320 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
321 (abs(eta) > 2.5) * (0.00)}
322}
323
324##################
325# Photon isolation
326##################
327
328module Isolation PhotonIsolation {
329 set CandidateInputArray PhotonEfficiency/photons
330 set IsolationInputArray EFlowMerger/eflow
331
332 set OutputArray photons
333
334 set DeltaRMax 0.5
335
336 set PTMin 0.5
337
338 set PTRatioMax 0.1
339}
340
341#####################
342# Electron efficiency
343#####################
344
345module Efficiency ElectronEfficiency {
346 set InputArray ElectronEnergySmearing/electrons
347 set OutputArray electrons
348
349 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
350
351 # efficiency formula for electrons
352 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
353 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
354 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
355 (abs(eta) > 2.5) * (0.00)}
356}
357
358####################
359# Electron isolation
360####################
361
362module Isolation ElectronIsolation {
363 set CandidateInputArray ElectronEfficiency/electrons
364 set IsolationInputArray EFlowMerger/eflow
365
366 set OutputArray electrons
367
368 set DeltaRMax 0.5
369
370 set PTMin 0.5
371
372 set PTRatioMax 0.1
373}
374
375#################
376# Muon efficiency
377#################
378
379module Efficiency MuonEfficiency {
380 set InputArray MuonMomentumSmearing/muons
381 set OutputArray muons
382
383 # set EfficiencyFormula {efficiency as a function of eta and pt}
384
385 # efficiency formula for muons
386 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
387 (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
388 (abs(eta) <= 1.5) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
389 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
390 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3) * (0.95 * exp(0.5 - pt*5.0e-4)) +
391 (abs(eta) > 2.4) * (0.00)}
392}
393
394################
395# Muon isolation
396################
397
398module Isolation MuonIsolation {
399 set CandidateInputArray MuonEfficiency/muons
400 set IsolationInputArray EFlowMerger/eflow
401
402 set OutputArray muons
403
404 set DeltaRMax 0.5
405
406 set PTMin 0.5
407
408 set PTRatioMax 0.1
409}
410
411###################
412# Missing ET merger
413###################
414
415module Merger MissingET {
416# add InputArray InputArray
417 add InputArray EFlowMerger/eflow
418 set MomentumOutputArray momentum
419}
420
421##################
422# Scalar HT merger
423##################
424
425module Merger ScalarHT {
426# add InputArray InputArray
427 add InputArray UniqueObjectFinder/jets
428 add InputArray UniqueObjectFinder/electrons
429 add InputArray UniqueObjectFinder/photons
430 add InputArray UniqueObjectFinder/muons
431 set EnergyOutputArray energy
432}
433
434
435#####################
436# Neutrino Filter
437#####################
438
439module PdgCodeFilter NeutrinoFilter {
440
441 set InputArray Delphes/stableParticles
442 set OutputArray filteredParticles
443
444 set PTMin 0.0
445
446 add PdgCode {12}
447 add PdgCode {14}
448 add PdgCode {16}
449 add PdgCode {-12}
450 add PdgCode {-14}
451 add PdgCode {-16}
452
453}
454
455
456#####################
457# MC truth jet finder
458#####################
459
460module FastJetFinder GenJetFinder {
461 set InputArray NeutrinoFilter/filteredParticles
462
463 set OutputArray jets
464
465 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
466 set JetAlgorithm 6
467 set ParameterR 0.5
468
469 set JetPTMin 20.0
470}
471
472
473############
474# Jet finder
475############
476
477module FastJetFinder FastJetFinder {
478# set InputArray Calorimeter/towers
479 set InputArray EFlowMerger/eflow
480
481 set OutputArray jets
482
483 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
484 set JetAlgorithm 6
485 set ParameterR 0.5
486
487 set JetPTMin 20.0
488}
489
490###############
491# Jet Fakes
492###############
493
494
495module JetFakeParticle JetFakeParticle {
496 set InputArray FastJetFinder/jets
497
498
499 set PhotonOutputArray photons
500 set ElectronOutputArray electrons
501 set MuonOutputArray muons
502
503 set OutputArray jets
504
505 # {PID out} {formula}
506 # make sure output particles are only electrons, muons or photons
507
508 add EfficiencyFormula {11} { (eta <= 2.5) * (pt > 20.0) * (0.30)}
509 add EfficiencyFormula {13} { (eta <= 2.5) * (pt > 20.0) * (0.30)}
510 add EfficiencyFormula {22} { (eta <= 2.5) * (pt > 20.0) * (0.30)}
511
512}
513
514##################
515# Electron merger
516##################
517
518module Merger ElectronMerger {
519# add InputArray InputArray
520 add InputArray ElectronIsolation/electrons
521 add InputArray JetFakeParticle/electrons
522 set OutputArray electrons
523}
524
525##################
526# Muon merger
527##################
528
529module Merger MuonMerger {
530# add InputArray InputArray
531 add InputArray MuonIsolation/muons
532 add InputArray JetFakeParticle/muons
533 set OutputArray muons
534}
535
536##################
537# Photon merger
538##################
539
540module Merger PhotonMerger {
541# add InputArray InputArray
542 add InputArray PhotonIsolation/photons
543 add InputArray JetFakeParticle/photons
544 set OutputArray photons
545}
546
547
548
549##################
550# Jet Energy Scale
551##################
552
553module EnergyScale JetEnergyScale {
554 set InputArray JetFakeParticle/jets
555 set OutputArray jets
556
557 # scale formula for jets
558 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
559}
560
561########################
562# Jet Flavor Association
563########################
564
565module JetFlavorAssociation JetFlavorAssociation {
566
567 set PartonInputArray Delphes/partons
568 set ParticleInputArray Delphes/allParticles
569 set ParticleLHEFInputArray Delphes/allParticlesLHEF
570 set JetInputArray JetEnergyScale/jets
571
572 set DeltaR 0.5
573 set PartonPTMin 1.0
574 set PartonEtaMax 2.5
575
576}
577
578###########
579# b-tagging
580###########
581
582module BTagging BTagging {
583 set JetInputArray JetEnergyScale/jets
584
585 set BitNumber 0
586
587 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
588 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
589 # gluon's PDG code has the lowest priority
590
591 # https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsBTV
592 # default efficiency formula (misidentification rate)
593 add EfficiencyFormula {0} {0.001}
594
595 # efficiency formula for c-jets (misidentification rate)
596 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) +
597 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) +
598 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) +
599 (abs(eta) > 2.5) * (0.000)}
600
601 # efficiency formula for b-jets
602 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) +
603 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) +
604 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) +
605 (abs(eta) > 2.5) * (0.000)}
606}
607
608#############
609# tau-tagging
610#############
611
612module TauTagging TauTagging {
613 set ParticleInputArray Delphes/allParticles
614 set PartonInputArray Delphes/partons
615 set JetInputArray JetEnergyScale/jets
616
617 set DeltaR 0.5
618
619 set TauPTMin 1.0
620
621 set TauEtaMax 2.5
622
623 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
624
625 # default efficiency formula (misidentification rate)
626 add EfficiencyFormula {0} {0.001}
627 # efficiency formula for tau-jets
628 add EfficiencyFormula {15} {0.4}
629}
630
631#####################################################
632# Find uniquely identified photons/electrons/tau/jets
633#####################################################
634
635module UniqueObjectFinder UniqueObjectFinder {
636# earlier arrays take precedence over later ones
637# add InputArray InputArray OutputArray
638 add InputArray PhotonMerger/photons photons
639 add InputArray ElectronMerger/electrons electrons
640 add InputArray MuonMerger/muons muons
641 add InputArray JetFakeParticle/jets jets
642}
643
644##################
645# ROOT tree writer
646##################
647
648# tracks, towers and eflow objects are not stored by default in the output.
649# if needed (for jet constituent or other studies), uncomment the relevant
650# "add Branch ..." lines.
651
652module TreeWriter TreeWriter {
653# add Branch InputArray BranchName BranchClass
654 add Branch Delphes/allParticles Particle GenParticle
655
656 add Branch TrackMerger/tracks Track Track
657 add Branch Calorimeter/towers Tower Tower
658
659 add Branch Calorimeter/eflowTracks EFlowTrack Track
660 add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
661 add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
662
663 add Branch GenJetFinder/jets GenJet Jet
664 add Branch UniqueObjectFinder/jets Jet Jet
665 add Branch UniqueObjectFinder/electrons Electron Electron
666 add Branch UniqueObjectFinder/photons Photon Photon
667 add Branch UniqueObjectFinder/muons Muon Muon
668 add Branch MissingET/momentum MissingET MissingET
669 add Branch ScalarHT/energy ScalarHT ScalarHT
670}