Fork me on GitHub

Ticket #318: delphes_card_Snowmass_NoPileUp_updateMay8_fixed_V2.tcl

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