Fork me on GitHub

Ticket #870: test.tcl

File test.tcl, 20.0 KB (added by Li Huang, 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 ElectronMomentumSmearing
14 MuonMomentumSmearing
15
16 TrackMerger
17 Calorimeter
18 EFlowMerger
19
20 PhotonEfficiency
21 PhotonIsolation
22
23 ElectronFilter
24 ElectronEfficiency
25 ElectronIsolation
26
27 MuonEfficiency
28 MuonIsolation
29
30 MissingET
31
32 NeutrinoFilter
33 GenJetFinder
34 GenMissingET
35
36 FastJetFinder
37 FastJetFinder1
38
39 JetEnergyScale
40 JetEnergyScale1
41
42 JetFlavorAssociation
43
44 BTagging
45 TauTagging
46 BTagging1
47 TauTagging1
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.15
70 # half-length of the magnetic field coverage, in m
71 set HalfLength 3.51
72
73 # magnetic field
74 set Bz 2.0
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) * (0.99) +
131 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
132 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) +
133 (abs(eta) > 2.5) * (0.00)}
134}
135
136########################################
137# Momentum resolution for charged tracks
138########################################
139
140module MomentumSmearing ChargedHadronMomentumSmearing {
141 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
142 set OutputArray chargedHadrons
143
144 # set ResolutionFormula {resolution formula as a function of eta and pt}
145
146 # resolution formula for charged hadrons
147 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
148 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
149 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
150}
151
152###################################
153# Momentum resolution for electrons
154###################################
155
156module MomentumSmearing ElectronMomentumSmearing {
157 set InputArray ElectronTrackingEfficiency/electrons
158 set OutputArray electrons
159
160 # set ResolutionFormula {resolution formula as a function of eta and energy}
161
162 # resolution formula for electrons
163 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
164 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
165 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
166}
167
168###############################
169# Momentum resolution for muons
170###############################
171
172module MomentumSmearing MuonMomentumSmearing {
173 set InputArray MuonTrackingEfficiency/muons
174 set OutputArray muons
175
176 # set ResolutionFormula {resolution formula as a function of eta and pt}
177
178 # resolution formula for muons
179 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.02^2 + pt^2*2.0e-4^2) +
180 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*3.0e-4^2) +
181 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*6.0e-4^2)}
182}
183
184##############
185# Track merger
186##############
187
188module Merger TrackMerger {
189# add InputArray InputArray
190 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
191 add InputArray ElectronMomentumSmearing/electrons
192 add InputArray MuonMomentumSmearing/muons
193 set OutputArray tracks
194}
195
196#############
197# Calorimeter
198#############
199
200module Calorimeter Calorimeter {
201 set ParticleInputArray ParticlePropagator/stableParticles
202 set TrackInputArray TrackMerger/tracks
203
204 set TowerOutputArray towers
205 set PhotonOutputArray photons
206
207 set EFlowTrackOutputArray eflowTracks
208 set EFlowPhotonOutputArray eflowPhotons
209 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
210
211 set ECalEnergyMin 0.5
212 set HCalEnergyMin 1.0
213
214 set ECalEnergySignificanceMin 1.0
215 set HCalEnergySignificanceMin 1.0
216
217 set SmearTowerCenter true
218
219 set pi [expr {acos(-1)}]
220
221 # lists of the edges of each tower in eta and phi
222 # each list starts with the lower edge of the first tower
223 # the list ends with the higher edged of the last tower
224
225 # 10 degrees towers
226 set PhiBins {}
227 for {set i -18} {$i <= 18} {incr i} {
228 add PhiBins [expr {$i * $pi/18.0}]
229 }
230 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} {
231 add EtaPhiBins $eta $PhiBins
232 }
233
234 # 20 degrees towers
235 set PhiBins {}
236 for {set i -9} {$i <= 9} {incr i} {
237 add PhiBins [expr {$i * $pi/9.0}]
238 }
239 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} {
240 add EtaPhiBins $eta $PhiBins
241 }
242
243 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
244 add EnergyFraction {0} {0.0 1.0}
245 # energy fractions for e, gamma and pi0
246 add EnergyFraction {11} {1.0 0.0}
247 add EnergyFraction {22} {1.0 0.0}
248 add EnergyFraction {111} {1.0 0.0}
249 # energy fractions for muon, neutrinos and neutralinos
250 add EnergyFraction {12} {0.0 0.0}
251 add EnergyFraction {13} {0.0 0.0}
252 add EnergyFraction {14} {0.0 0.0}
253 add EnergyFraction {16} {0.0 0.0}
254 add EnergyFraction {1000022} {0.0 0.0}
255 add EnergyFraction {1000023} {0.0 0.0}
256 add EnergyFraction {1000025} {0.0 0.0}
257 add EnergyFraction {1000035} {0.0 0.0}
258 add EnergyFraction {1000045} {0.0 0.0}
259 # energy fractions for K0short and Lambda
260 add EnergyFraction {310} {0.3 0.7}
261 add EnergyFraction {3122} {0.3 0.7}
262
263 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
264 # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
265 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
266 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
267 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) +
268 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
269
270 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
271 # http://arxiv.org/pdf/hep-ex/0004009v1
272 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
273 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
274 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
275 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
276}
277
278####################
279# Energy flow merger
280####################
281
282module Merger EFlowMerger {
283# add InputArray InputArray
284 add InputArray Calorimeter/eflowTracks
285 add InputArray Calorimeter/eflowPhotons
286 add InputArray Calorimeter/eflowNeutralHadrons
287 set OutputArray eflow
288}
289
290###################
291# Photon efficiency
292###################
293
294module Efficiency PhotonEfficiency {
295 set InputArray Calorimeter/eflowPhotons
296 set OutputArray photons
297
298 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
299
300 # efficiency formula for photons
301 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
302 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
303 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
304 (abs(eta) > 2.5) * (0.00)}
305}
306
307##################
308# Photon isolation
309##################
310
311module Isolation PhotonIsolation {
312 set CandidateInputArray PhotonEfficiency/photons
313 set IsolationInputArray EFlowMerger/eflow
314
315 set OutputArray photons
316
317 set DeltaRMax 0.5
318
319 set PTMin 0.5
320
321 set PTRatioMax 0.12
322}
323
324#################
325# Electron filter
326#################
327
328module PdgCodeFilter ElectronFilter {
329 set InputArray Calorimeter/eflowTracks
330 set OutputArray electrons
331 set Invert true
332 add PdgCode {11}
333 add PdgCode {-11}
334}
335
336#####################
337# Electron efficiency
338#####################
339
340module Efficiency ElectronEfficiency {
341 set InputArray ElectronFilter/electrons
342 set OutputArray electrons
343
344 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
345
346 # efficiency formula for electrons
347 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
348 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
349 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
350 (abs(eta) > 2.5) * (0.00)}
351}
352
353####################
354# Electron isolation
355####################
356
357module Isolation ElectronIsolation {
358 set CandidateInputArray ElectronEfficiency/electrons
359 set IsolationInputArray EFlowMerger/eflow
360
361 set OutputArray electrons
362
363 set DeltaRMax 0.5
364
365 set PTMin 0.5
366
367 set PTRatioMax 0.12
368}
369
370#################
371# Muon efficiency
372#################
373
374module Efficiency MuonEfficiency {
375 set InputArray MuonMomentumSmearing/muons
376 set OutputArray muons
377
378 # set EfficiencyFormula {efficiency as a function of eta and pt}
379
380 # efficiency formula for muons
381 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
382 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
383 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) +
384 (abs(eta) > 2.7) * (0.00)}
385}
386
387################
388# Muon isolation
389################
390
391module Isolation MuonIsolation {
392 set CandidateInputArray MuonEfficiency/muons
393 set IsolationInputArray EFlowMerger/eflow
394
395 set OutputArray muons
396
397 set DeltaRMax 0.5
398
399 set PTMin 0.5
400
401 set PTRatioMax 0.25
402}
403
404###################
405# Missing ET merger
406###################
407
408module Merger MissingET {
409# add InputArray InputArray
410 add InputArray EFlowMerger/eflow
411 set MomentumOutputArray momentum
412}
413
414##################
415# Scalar HT merger
416##################
417
418module Merger ScalarHT {
419# add InputArray InputArray
420 add InputArray UniqueObjectFinder/jets
421 add InputArray UniqueObjectFinder/electrons
422 add InputArray UniqueObjectFinder/photons
423 add InputArray UniqueObjectFinder/muons
424 set EnergyOutputArray energy
425}
426
427
428#####################
429# Neutrino Filter
430#####################
431
432module PdgCodeFilter NeutrinoFilter {
433
434 set InputArray Delphes/stableParticles
435 set OutputArray filteredParticles
436
437 set PTMin 0.0
438
439 add PdgCode {12}
440 add PdgCode {14}
441 add PdgCode {16}
442 add PdgCode {-12}
443 add PdgCode {-14}
444 add PdgCode {-16}
445
446}
447
448#####################
449# MC truth jet finder
450#####################
451
452module FastJetFinder GenJetFinder {
453 set InputArray NeutrinoFilter/filteredParticles
454
455 set OutputArray jets
456
457 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
458 set JetAlgorithm 6
459 set ParameterR 0.6
460
461 set JetPTMin 20.0
462}
463
464
465#########################
466# Gen Missing ET merger
467########################
468
469module Merger GenMissingET {
470# add InputArray InputArray
471 add InputArray NeutrinoFilter/filteredParticles
472 set MomentumOutputArray momentum
473}
474
475
476
477############
478# Jet finder
479############
480
481module FastJetFinder FastJetFinder {
482 set InputArray Calorimeter/towers
483
484 set OutputArray jets
485
486 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
487 set JetAlgorithm 6
488 set ParameterR 0.6
489
490 set JetPTMin 20.0
491}
492
493module FastJetFinder FastJetFinder1 {
494 set InputArray Calorimeter/towers
495
496 set OutputArray jets
497
498 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
499 set JetAlgorithm 6
500 set ParameterR 0.3
501
502 set JetPTMin 20.0
503}
504
505##################
506# Jet Energy Scale
507##################
508
509module EnergyScale JetEnergyScale {
510 set InputArray FastJetFinder/jets
511 set OutputArray jets
512
513 # scale formula for jets
514 set ScaleFormula { sqrt( (3.0 - 0.2*(abs(eta)))^2 / pt + 1.0 ) }
515}
516
517module EnergyScale JetEnergyScale1 {
518 set InputArray FastJetFinder1/jets
519 set OutputArray jets
520
521 # scale formula for jets
522 set ScaleFormula { sqrt( (3.0 - 0.2*(abs(eta)))^2 / pt + 1.0 ) }
523}
524
525########################
526# Jet Flavor Association
527########################
528
529module JetFlavorAssociation JetFlavorAssociation {
530
531 set PartonInputArray Delphes/partons
532 set ParticleInputArray Delphes/allParticles
533 set ParticleLHEFInputArray Delphes/allParticlesLHEF
534 set JetInputArray JetEnergyScale/jets
535 set JetInputArray JetEnergyScale1/jets
536
537 set DeltaR 0.5
538 set PartonPTMin 1.0
539 set PartonEtaMax 2.5
540
541}
542
543###########
544# b-tagging
545###########
546
547module BTagging BTagging {
548 set JetInputArray JetEnergyScale/jets
549
550 set BitNumber 0
551
552 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
553 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
554 # gluon's PDG code has the lowest priority
555
556 # based on ATL-PHYS-PUB-2015-022
557
558 # default efficiency formula (misidentification rate)
559 add EfficiencyFormula {0} {0.002+7.3e-06*pt}
560
561 # efficiency formula for c-jets (misidentification rate)
562 add EfficiencyFormula {4} {0.20*tanh(0.02*pt)*(1/(1+0.0034*pt))}
563
564 # efficiency formula for b-jets
565 add EfficiencyFormula {5} {0.80*tanh(0.003*pt)*(30/(1+0.086*pt))}
566}
567
568module BTagging BTagging1 {
569 set JetInputArray JetEnergyScale1/jets
570
571 set BitNumber 0
572
573 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
574 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
575 # gluon's PDG code has the lowest priority
576
577 # based on ATL-PHYS-PUB-2015-022
578
579 # default efficiency formula (misidentification rate)
580 add EfficiencyFormula {0} {0.002+7.3e-06*pt}
581
582 # efficiency formula for c-jets (misidentification rate)
583 add EfficiencyFormula {4} {0.20*tanh(0.02*pt)*(1/(1+0.0034*pt))}
584
585 # efficiency formula for b-jets
586 add EfficiencyFormula {5} {0.80*tanh(0.003*pt)*(30/(1+0.086*pt))}
587}
588
589#############
590# tau-tagging
591#############
592
593module TrackCountingTauTagging TauTagging {
594
595 set ParticleInputArray Delphes/allParticles
596 set PartonInputArray Delphes/partons
597 set TrackInputArray TrackMerger/tracks
598 set JetInputArray JetEnergyScale/jets
599
600 set DeltaR 0.2
601 set DeltaRTrack 0.2
602
603 set TrackPTMin 1.0
604
605 set TauPTMin 1.0
606 set TauEtaMax 2.5
607
608 # instructions: {n-prongs} {eff}
609
610 # 1 - one prong efficiency
611 # 2 - two or more efficiency
612 # -1 - one prong mistag rate
613 # -2 - two or more mistag rate
614
615 set BitNumber 0
616
617 # taken from ATL-PHYS-PUB-2015-045 (medium working point)
618 add EfficiencyFormula {1} {0.70}
619 add EfficiencyFormula {2} {0.60}
620 add EfficiencyFormula {-1} {0.02}
621 add EfficiencyFormula {-2} {0.01}
622
623}
624
625module TrackCountingTauTagging TauTagging1 {
626
627 set ParticleInputArray Delphes/allParticles
628 set PartonInputArray Delphes/partons
629 set TrackInputArray TrackMerger/tracks
630 set JetInputArray JetEnergyScale1/jets
631
632 set DeltaR 0.2
633 set DeltaRTrack 0.2
634
635 set TrackPTMin 1.0
636
637 set TauPTMin 1.0
638 set TauEtaMax 2.5
639
640 # instructions: {n-prongs} {eff}
641
642 # 1 - one prong efficiency
643 # 2 - two or more efficiency
644 # -1 - one prong mistag rate
645 # -2 - two or more mistag rate
646
647 set BitNumber 0
648
649 # taken from ATL-PHYS-PUB-2015-045 (medium working point)
650 add EfficiencyFormula {1} {0.70}
651 add EfficiencyFormula {2} {0.60}
652 add EfficiencyFormula {-1} {0.02}
653 add EfficiencyFormula {-2} {0.01}
654
655}
656
657#####################################################
658# Find uniquely identified photons/electrons/tau/jets
659#####################################################
660
661module UniqueObjectFinder UniqueObjectFinder {
662# earlier arrays take precedence over later ones
663# add InputArray InputArray OutputArray
664 add InputArray PhotonIsolation/photons photons
665 add InputArray ElectronIsolation/electrons electrons
666 add InputArray MuonIsolation/muons muons
667 add InputArray JetEnergyScale/jets jets
668 add InputArray JetEnergyScale1/jets jetsN
669}
670
671##################
672# ROOT tree writer
673##################
674
675# tracks, towers and eflow objects are not stored by default in the output.
676# if needed (for jet constituent or other studies), uncomment the relevant
677# "add Branch ..." lines.
678
679module TreeWriter TreeWriter {
680# add Branch InputArray BranchName BranchClass
681 add Branch Delphes/allParticles Particle GenParticle
682
683 add Branch TrackMerger/tracks Track Track
684 add Branch Calorimeter/towers Tower Tower
685
686 add Branch Calorimeter/eflowTracks EFlowTrack Track
687 add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
688 add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
689
690 add Branch GenJetFinder/jets GenJet Jet
691 add Branch GenMissingET/momentum GenMissingET MissingET
692
693 add Branch UniqueObjectFinder/jets Jet Jet
694 add Branch UniqueObjectFinder/jetsN JetN Jet
695 add Branch UniqueObjectFinder/electrons Electron Electron
696 add Branch UniqueObjectFinder/photons Photon Photon
697 add Branch UniqueObjectFinder/muons Muon Muon
698 add Branch MissingET/momentum MissingET MissingET
699 add Branch ScalarHT/energy ScalarHT ScalarHT
700}
701