Fork me on GitHub

source: git/cards/delphes_card_CMS_CSCCluster.tcl@ a5af1df

Last change on this file since a5af1df was a5af1df, checked in by christinaw97 <christina.wang@…>, 3 years ago

added CscCluster modules and classes
  • Property mode set to 100644
File size: 36.8 KB
Line 
1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6 ParticlePropagator
7
8
9
10 ChargedHadronTrackingEfficiency
11 ElectronTrackingEfficiency
12 MuonTrackingEfficiency
13
14 ChargedHadronMomentumSmearing
15 ElectronMomentumSmearing
16 MuonMomentumSmearing
17
18 TrackMerger
19
20 ECal
21 HCal
22
23 Calorimeter
24 EFlowMerger
25 EFlowFilter
26
27 PhotonEfficiency
28 PhotonIsolation
29
30 ElectronFilter
31 ElectronEfficiency
32 ElectronIsolation
33
34 ChargedHadronFilter
35
36 MuonEfficiency
37 MuonIsolation
38
39 MissingET
40
41 NeutrinoFilter
42 GenJetFinder
43 GenMissingET
44
45 FastJetFinder
46 FatJetFinder
47
48 JetEnergyScale
49
50 JetFlavorAssociation
51
52 BTagging
53 TauTagging
54
55 UniqueObjectFinder
56
57 ScalarHT
58
59 llpFilter
60 CSCFilter
61 CutBasedIDEfficiency
62 ClusterEfficiency
63
64
65 TreeWriter
66}
67
68#################################
69# Propagate particles in cylinder
70#################################
71
72module ParticlePropagator ParticlePropagator {
73 set InputArray Delphes/stableParticles
74
75 set OutputArray stableParticles
76 set ChargedHadronOutputArray chargedHadrons
77 set ElectronOutputArray electrons
78 set MuonOutputArray muons
79
80 # radius of the magnetic field coverage, in m
81 set Radius 1.29
82 # half-length of the magnetic field coverage, in m
83 set HalfLength 3.00
84
85 # magnetic field
86 set Bz 3.8
87}
88
89####################################
90# Charged hadron tracking efficiency
91####################################
92
93module Efficiency ChargedHadronTrackingEfficiency {
94 set InputArray ParticlePropagator/chargedHadrons
95 set OutputArray chargedHadrons
96
97 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
98
99 # tracking efficiency formula for charged hadrons
100 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
101 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
102 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
103 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
104 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
105 (abs(eta) > 2.5) * (0.00)}
106}
107
108##############################
109# Electron tracking efficiency
110##############################
111
112module Efficiency ElectronTrackingEfficiency {
113 set InputArray ParticlePropagator/electrons
114 set OutputArray electrons
115
116 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
117
118 # tracking efficiency formula for electrons
119 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
120 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
121 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
122 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
123 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
124 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
125 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
126 (abs(eta) > 2.5) * (0.00)}
127}
128
129##########################
130# Muon tracking efficiency
131##########################
132
133module Efficiency MuonTrackingEfficiency {
134 set InputArray ParticlePropagator/muons
135 set OutputArray muons
136
137 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
138
139 # tracking efficiency formula for muons
140 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
141 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
142 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e3) * (0.99) +
143 (abs(eta) <= 1.5) * (pt > 1.0e3 ) * (0.99 * exp(0.5 - pt*5.0e-4)) +
144
145 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
146 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e3) * (0.98) +
147 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3) * (0.98 * exp(0.5 - pt*5.0e-4)) +
148 (abs(eta) > 2.5) * (0.00)}
149}
150
151########################################
152# Momentum resolution for charged tracks
153########################################
154
155module MomentumSmearing ChargedHadronMomentumSmearing {
156 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
157 set OutputArray chargedHadrons
158
159 # set ResolutionFormula {resolution formula as a function of eta and pt}
160
161 # resolution formula for charged hadrons
162 # based on arXiv:1405.6569
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 electrons
170###################################
171
172module MomentumSmearing ElectronMomentumSmearing {
173 set InputArray ElectronTrackingEfficiency/electrons
174 set OutputArray electrons
175
176 # set ResolutionFormula {resolution formula as a function of eta and energy}
177
178 # resolution formula for electrons
179 # based on arXiv:1502.02701
180 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
181 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
182 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
183}
184
185###############################
186# Momentum resolution for muons
187###############################
188
189module MomentumSmearing MuonMomentumSmearing {
190 set InputArray MuonTrackingEfficiency/muons
191 set OutputArray muons
192
193 # set ResolutionFormula {resolution formula as a function of eta and pt}
194
195 # resolution formula for muons
196 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
197 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
198 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
199}
200
201##############
202# Track merger
203##############
204
205module Merger TrackMerger {
206# add InputArray InputArray
207 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
208 add InputArray ElectronMomentumSmearing/electrons
209 add InputArray MuonMomentumSmearing/muons
210 set OutputArray tracks
211}
212
213
214
215#############
216# ECAL
217#############
218
219module SimpleCalorimeter ECal {
220 set ParticleInputArray ParticlePropagator/stableParticles
221 set TrackInputArray TrackMerger/tracks
222
223 set TowerOutputArray ecalTowers
224 set EFlowTrackOutputArray eflowTracks
225 set EFlowTowerOutputArray eflowPhotons
226
227 set IsEcal true
228
229 set EnergyMin 0.5
230 set EnergySignificanceMin 2.0
231
232 set SmearTowerCenter true
233
234 set pi [expr {acos(-1)}]
235
236 # lists of the edges of each tower in eta and phi
237 # each list starts with the lower edge of the first tower
238 # the list ends with the higher edged of the last tower
239
240 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
241
242 set PhiBins {}
243 for {set i -180} {$i <= 180} {incr i} {
244 add PhiBins [expr {$i * $pi/180.0}]
245 }
246
247 # 0.02 unit in eta up to eta = 1.5 (barrel)
248 for {set i -85} {$i <= 86} {incr i} {
249 set eta [expr {$i * 0.0174}]
250 add EtaPhiBins $eta $PhiBins
251 }
252
253 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
254
255 set PhiBins {}
256 for {set i -180} {$i <= 180} {incr i} {
257 add PhiBins [expr {$i * $pi/180.0}]
258 }
259
260 # 0.02 unit in eta up to eta = 3
261 for {set i 1} {$i <= 84} {incr i} {
262 set eta [expr { -2.958 + $i * 0.0174}]
263 add EtaPhiBins $eta $PhiBins
264 }
265
266 for {set i 1} {$i <= 84} {incr i} {
267 set eta [expr { 1.4964 + $i * 0.0174}]
268 add EtaPhiBins $eta $PhiBins
269 }
270
271 # take present CMS granularity for HF
272
273 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
274 set PhiBins {}
275 for {set i -18} {$i <= 18} {incr i} {
276 add PhiBins [expr {$i * $pi/18.0}]
277 }
278
279 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} {
280 add EtaPhiBins $eta $PhiBins
281 }
282
283
284 add EnergyFraction {0} {0.0}
285 # energy fractions for e, gamma and pi0
286 add EnergyFraction {11} {1.0}
287 add EnergyFraction {22} {1.0}
288 add EnergyFraction {111} {1.0}
289 # energy fractions for muon, neutrinos and neutralinos
290 add EnergyFraction {12} {0.0}
291 add EnergyFraction {13} {0.0}
292 add EnergyFraction {14} {0.0}
293 add EnergyFraction {16} {0.0}
294 add EnergyFraction {1000022} {0.0}
295 add EnergyFraction {1000023} {0.0}
296 add EnergyFraction {1000025} {0.0}
297 add EnergyFraction {1000035} {0.0}
298 add EnergyFraction {1000045} {0.0}
299 # energy fractions for K0short and Lambda
300 add EnergyFraction {310} {0.3}
301 add EnergyFraction {3122} {0.3}
302
303 # set ResolutionFormula {resolution formula as a function of eta and energy}
304
305 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
306
307 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
308 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
309 set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
310 (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) +
311 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
312
313}
314
315
316#############
317# HCAL
318#############
319
320module SimpleCalorimeter HCal {
321 set ParticleInputArray ParticlePropagator/stableParticles
322 set TrackInputArray ECal/eflowTracks
323
324 set TowerOutputArray hcalTowers
325 set EFlowTrackOutputArray eflowTracks
326 set EFlowTowerOutputArray eflowNeutralHadrons
327
328 set IsEcal false
329
330 set EnergyMin 1.0
331 set EnergySignificanceMin 1.0
332
333 set SmearTowerCenter true
334
335 set pi [expr {acos(-1)}]
336
337 # lists of the edges of each tower in eta and phi
338 # each list starts with the lower edge of the first tower
339 # the list ends with the higher edged of the last tower
340
341 # 5 degrees towers
342 set PhiBins {}
343 for {set i -36} {$i <= 36} {incr i} {
344 add PhiBins [expr {$i * $pi/36.0}]
345 }
346 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} {
347 add EtaPhiBins $eta $PhiBins
348 }
349
350 # 10 degrees towers
351 set PhiBins {}
352 for {set i -18} {$i <= 18} {incr i} {
353 add PhiBins [expr {$i * $pi/18.0}]
354 }
355 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} {
356 add EtaPhiBins $eta $PhiBins
357 }
358
359 # 20 degrees towers
360 set PhiBins {}
361 for {set i -9} {$i <= 9} {incr i} {
362 add PhiBins [expr {$i * $pi/9.0}]
363 }
364 foreach eta {-5 -4.7 -4.525 4.7 5} {
365 add EtaPhiBins $eta $PhiBins
366 }
367
368 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
369 add EnergyFraction {0} {1.0}
370 # energy fractions for e, gamma and pi0
371 add EnergyFraction {11} {0.0}
372 add EnergyFraction {22} {0.0}
373 add EnergyFraction {111} {0.0}
374 # energy fractions for muon, neutrinos and neutralinos
375 add EnergyFraction {12} {0.0}
376 add EnergyFraction {13} {0.0}
377 add EnergyFraction {14} {0.0}
378 add EnergyFraction {16} {0.0}
379 add EnergyFraction {1000022} {0.0}
380 add EnergyFraction {1000023} {0.0}
381 add EnergyFraction {1000025} {0.0}
382 add EnergyFraction {1000035} {0.0}
383 add EnergyFraction {1000045} {0.0}
384 # energy fractions for K0short and Lambda
385 add EnergyFraction {310} {0.7}
386 add EnergyFraction {3122} {0.7}
387
388 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
389 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
390 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
391
392}
393
394
395#################
396# Electron filter
397#################
398
399module PdgCodeFilter ElectronFilter {
400 set InputArray HCal/eflowTracks
401 set OutputArray electrons
402 set Invert true
403 add PdgCode {11}
404 add PdgCode {-11}
405}
406#################
407# LLP filter
408#################
409
410module PdgCodeFilter HiggsFilter {
411set InputArray Delphes/allParticles
412 set OutputArray higgs
413 set Invert true
414 add PdgCode {25}
415}
416
417
418
419
420# filter out LLPs that decay in CSC
421module LLPFilter CSCFilter {
422 set InputArray Delphes/allParticles
423 set OutputArray LLP
424 set DecayRegion 1
425 # DecayRegion = 0: no cuts on decay region
426 # DecayRegion = 1: select LLP that decays in CSC volume
427 # DecayRegion = 2: select LLP that decays outside of calorimeters, for genMET calculation
428 set RequireStatus true
429 set Status 2
430 add PdgCode {9000005}
431
432}
433# filter out LLPs regardless of decay position
434module LLPFilter llpFilter {
435 set InputArray Delphes/allParticles
436 set OutputArray LLP
437 set DecayRegion 0
438 set RequireStatus true
439 set Status 2
440
441 # DecayRegion = 0: no cuts on decay region
442 # DecayRegion = 1: select LLP that decays in CSC volume
443 # DecayRegion = 2: select LLP that decays outside of calorimeters, for genMET calculation
444 add PdgCode {9000005}
445
446}
447
448
449
450module CscClusterEfficiency ClusterEfficiency {
451 set InputArray CutBasedIDEfficiency/cluster
452 set OutputArray cluster
453
454 # efficiency formula for Csc Cluster, as a function of LLP decay vertex in R, Z and hadronic and EM energy
455 set EfficiencyFormula {
456
457 (pt > 3900 && eta < 6710) * ((energy >= 0.0 && energy < 25.0 && 0.0 == phi)*(0.0049) +
458 (energy >= 0.0 && energy < 25.0&& phi > 0.0 && phi < 25.0)*(0.0130) +
459 (energy >= 0.0 && energy < 25.0&& phi >= 25.0 && phi < 50.0)*(0.0346) +
460 (energy >= 0.0 && energy < 25.0&& phi >= 50.0 && phi < 75.0)*(0.0623) +
461 (energy >= 0.0 && energy < 25.0&& phi >= 75.0 && phi < 100.0)*(0.0919) +
462 (energy >= 0.0 && energy < 25.0&& phi >= 100.0 && phi < 150.0)*(0.1086) +
463 (energy >= 0.0 && energy < 25.0&& phi >= 150.0 && phi < 200.0)*(0.1292) +
464 (energy >= 0.0 && energy < 25.0 && phi >= 200.0)*(0.1106) +
465 (energy >= 25.0 && energy < 50.0 && 0.0 == phi)*(0.0249) +
466 (energy >= 25.0 && energy < 50.0&& phi > 0.0 && phi < 25.0)*(0.0285) +
467 (energy >= 25.0 && energy < 50.0&& phi >= 25.0 && phi < 50.0)*(0.0501) +
468 (energy >= 25.0 && energy < 50.0&& phi >= 50.0 && phi < 75.0)*(0.0841) +
469 (energy >= 25.0 && energy < 50.0&& phi >= 75.0 && phi < 100.0)*(0.1021) +
470 (energy >= 25.0 && energy < 50.0&& phi >= 100.0 && phi < 150.0)*(0.1129) +
471 (energy >= 25.0 && energy < 50.0&& phi >= 150.0 && phi < 200.0)*(0.1141) +
472 (energy >= 25.0 && energy < 50.0 && phi >= 200.0)*(0.1370) +
473 (energy >= 50.0 && energy < 75.0 && 0.0 == phi)*(0.0282) +
474 (energy >= 50.0 && energy < 75.0&& phi > 0.0 && phi < 25.0)*(0.0445) +
475 (energy >= 50.0 && energy < 75.0&& phi >= 25.0 && phi < 50.0)*(0.0643) +
476 (energy >= 50.0 && energy < 75.0&& phi >= 50.0 && phi < 75.0)*(0.0903) +
477 (energy >= 50.0 && energy < 75.0&& phi >= 75.0 && phi < 100.0)*(0.0998) +
478 (energy >= 50.0 && energy < 75.0&& phi >= 100.0 && phi < 150.0)*(0.1420) +
479 (energy >= 50.0 && energy < 75.0&& phi >= 150.0 && phi < 200.0)*(0.1429) +
480 (energy >= 50.0 && energy < 75.0 && phi >= 200.0)*(0.0882) +
481 (energy >= 75.0 && energy < 100.0 && 0.0 == phi)*(0.0594) +
482 (energy >= 75.0 && energy < 100.0&& phi > 0.0 && phi < 25.0)*(0.0521) +
483 (energy >= 75.0 && energy < 100.0&& phi >= 25.0 && phi < 50.0)*(0.0605) +
484 (energy >= 75.0 && energy < 100.0&& phi >= 50.0 && phi < 75.0)*(0.0791) +
485 (energy >= 75.0 && energy < 100.0&& phi >= 75.0 && phi < 100.0)*(0.1117) +
486 (energy >= 75.0 && energy < 100.0&& phi >= 100.0 && phi < 150.0)*(0.0862) +
487 (energy >= 75.0 && energy < 100.0&& phi >= 150.0 && phi < 200.0)*(0.0698) +
488 (energy >= 75.0 && energy < 100.0 && phi >= 200.0)*(0.0500) +
489 (energy >= 100.0 && energy < 125.0 && 0.0 == phi)*(0.0758) +
490 (energy >= 100.0 && energy < 125.0&& phi > 0.0 && phi < 25.0)*(0.0414) +
491 (energy >= 100.0 && energy < 125.0&& phi >= 25.0 && phi < 50.0)*(0.0755) +
492 (energy >= 100.0 && energy < 125.0&& phi >= 50.0 && phi < 75.0)*(0.1027) +
493 (energy >= 100.0 && energy < 125.0&& phi >= 75.0 && phi < 100.0)*(0.0440) +
494 (energy >= 100.0 && energy < 125.0&& phi >= 100.0 && phi < 150.0)*(0.0811) +
495 (energy >= 100.0 && energy < 125.0&& phi >= 150.0 && phi < 200.0)*(0.1538) +
496 (energy >= 100.0 && energy < 125.0 && phi >= 200.0)*(0.0833) +
497 (energy >= 125.0 && energy < 150.0 && 0.0 == phi)*(0.0300) +
498 (energy >= 125.0 && energy < 150.0&& phi > 0.0 && phi < 25.0)*(0.0609) +
499 (energy >= 125.0 && energy < 150.0&& phi >= 25.0 && phi < 50.0)*(0.0745) +
500 (energy >= 125.0 && energy < 150.0&& phi >= 50.0 && phi < 75.0)*(0.0610) +
501 (energy >= 125.0 && energy < 150.0&& phi >= 75.0 && phi < 100.0)*(0.1224) +
502 (energy >= 125.0 && energy < 150.0&& phi >= 100.0 && phi < 150.0)*(0.1667) +
503 (energy >= 125.0 && energy < 150.0&& phi >= 150.0 && phi < 200.0)*(0.0000) +
504 (energy >= 125.0 && energy < 150.0 && phi >= 200.0)*(0.0000) +
505 (energy >= 150.0 && 0.0 == phi)*(0.0282) +
506 (energy >= 150.0&& phi > 0.0 && phi < 25.0)*(0.0809) +
507 (energy >= 150.0&& phi >= 25.0 && phi < 50.0)*(0.0352) +
508 (energy >= 150.0&& phi >= 50.0 && phi < 75.0)*(0.0984) +
509 (energy >= 150.0&& phi >= 75.0 && phi < 100.0)*(0.0968) +
510 (energy >= 150.0&& phi >= 100.0 && phi < 150.0)*(0.1282) +
511 (energy >= 150.0&& phi >= 150.0 && phi < 200.0)*(0.2105) +
512 (energy >= 150.0 && phi >= 200.0)*(0.0769)) +
513 (eta > 6710) * ((energy >= 0.0 && energy < 25.0 && 0.0 == phi)*(0.0184) +
514 (energy >= 0.0 && energy < 25.0&& phi > 0.0 && phi < 25.0)*(0.0772) +
515 (energy >= 0.0 && energy < 25.0&& phi >= 25.0 && phi < 50.0)*(0.2086) +
516 (energy >= 0.0 && energy < 25.0&& phi >= 50.0 && phi < 75.0)*(0.3091) +
517 (energy >= 0.0 && energy < 25.0&& phi >= 75.0 && phi < 100.0)*(0.3867) +
518 (energy >= 0.0 && energy < 25.0&& phi >= 100.0 && phi < 150.0)*(0.4500) +
519 (energy >= 0.0 && energy < 25.0&& phi >= 150.0 && phi < 200.0)*(0.4746) +
520 (energy >= 0.0 && energy < 25.0 && phi >= 200.0)*(0.4906) +
521 (energy >= 25.0 && energy < 50.0 && 0.0 == phi)*(0.0955) +
522 (energy >= 25.0 && energy < 50.0&& phi > 0.0 && phi < 25.0)*(0.1461) +
523 (energy >= 25.0 && energy < 50.0&& phi >= 25.0 && phi < 50.0)*(0.2594) +
524 (energy >= 25.0 && energy < 50.0&& phi >= 50.0 && phi < 75.0)*(0.3556) +
525 (energy >= 25.0 && energy < 50.0&& phi >= 75.0 && phi < 100.0)*(0.4165) +
526 (energy >= 25.0 && energy < 50.0&& phi >= 100.0 && phi < 150.0)*(0.4693) +
527 (energy >= 25.0 && energy < 50.0&& phi >= 150.0 && phi < 200.0)*(0.5054) +
528 (energy >= 25.0 && energy < 50.0 && phi >= 200.0)*(0.5219) +
529 (energy >= 50.0 && energy < 75.0 && 0.0 == phi)*(0.1472) +
530 (energy >= 50.0 && energy < 75.0&& phi > 0.0 && phi < 25.0)*(0.1970) +
531 (energy >= 50.0 && energy < 75.0&& phi >= 25.0 && phi < 50.0)*(0.2974) +
532 (energy >= 50.0 && energy < 75.0&& phi >= 50.0 && phi < 75.0)*(0.3783) +
533 (energy >= 50.0 && energy < 75.0&& phi >= 75.0 && phi < 100.0)*(0.4335) +
534 (energy >= 50.0 && energy < 75.0&& phi >= 100.0 && phi < 150.0)*(0.4736) +
535 (energy >= 50.0 && energy < 75.0&& phi >= 150.0 && phi < 200.0)*(0.4937) +
536 (energy >= 50.0 && energy < 75.0 && phi >= 200.0)*(0.5077) +
537 (energy >= 75.0 && energy < 100.0 && 0.0 == phi)*(0.2053) +
538 (energy >= 75.0 && energy < 100.0&& phi > 0.0 && phi < 25.0)*(0.2314) +
539 (energy >= 75.0 && energy < 100.0&& phi >= 25.0 && phi < 50.0)*(0.3114) +
540 (energy >= 75.0 && energy < 100.0&& phi >= 50.0 && phi < 75.0)*(0.3799) +
541 (energy >= 75.0 && energy < 100.0&& phi >= 75.0 && phi < 100.0)*(0.4420) +
542 (energy >= 75.0 && energy < 100.0&& phi >= 100.0 && phi < 150.0)*(0.4502) +
543 (energy >= 75.0 && energy < 100.0&& phi >= 150.0 && phi < 200.0)*(0.5348) +
544 (energy >= 75.0 && energy < 100.0 && phi >= 200.0)*(0.5115) +
545 (energy >= 100.0 && energy < 125.0 && 0.0 == phi)*(0.2198) +
546 (energy >= 100.0 && energy < 125.0&& phi > 0.0 && phi < 25.0)*(0.2404) +
547 (energy >= 100.0 && energy < 125.0&& phi >= 25.0 && phi < 50.0)*(0.3295) +
548 (energy >= 100.0 && energy < 125.0&& phi >= 50.0 && phi < 75.0)*(0.3932) +
549 (energy >= 100.0 && energy < 125.0&& phi >= 75.0 && phi < 100.0)*(0.4327) +
550 (energy >= 100.0 && energy < 125.0&& phi >= 100.0 && phi < 150.0)*(0.4377) +
551 (energy >= 100.0 && energy < 125.0&& phi >= 150.0 && phi < 200.0)*(0.5175) +
552 (energy >= 100.0 && energy < 125.0 && phi >= 200.0)*(0.6087) +
553 (energy >= 125.0 && energy < 150.0 && 0.0 == phi)*(0.2147) +
554 (energy >= 125.0 && energy < 150.0&& phi > 0.0 && phi < 25.0)*(0.2605) +
555 (energy >= 125.0 && energy < 150.0&& phi >= 25.0 && phi < 50.0)*(0.3442) +
556 (energy >= 125.0 && energy < 150.0&& phi >= 50.0 && phi < 75.0)*(0.3622) +
557 (energy >= 125.0 && energy < 150.0&& phi >= 75.0 && phi < 100.0)*(0.4407) +
558 (energy >= 125.0 && energy < 150.0&& phi >= 100.0 && phi < 150.0)*(0.5168) +
559 (energy >= 125.0 && energy < 150.0&& phi >= 150.0 && phi < 200.0)*(0.5056) +
560 (energy >= 125.0 && energy < 150.0 && phi >= 200.0)*(0.4559) +
561 (energy >= 150.0 && 0.0 == phi)*(0.2824) +
562 (energy >= 150.0&& phi > 0.0 && phi < 25.0)*(0.2447) +
563 (energy >= 150.0&& phi >= 25.0 && phi < 50.0)*(0.3519) +
564 (energy >= 150.0&& phi >= 50.0 && phi < 75.0)*(0.3772) +
565 (energy >= 150.0&& phi >= 75.0 && phi < 100.0)*(0.4447) +
566 (energy >= 150.0&& phi >= 100.0 && phi < 150.0)*(0.4703) +
567 (energy >= 150.0&& phi >= 150.0 && phi < 200.0)*(0.4460) +
568 (energy >= 150.0 && phi >= 200.0)*(0.4400)) +
569 (pt < 2700 && eta < 6710) * ((energy >= 0.0 && energy < 25.0 && 0.0 == phi)*(0.0002) +
570 (energy >= 0.0 && energy < 25.0&& phi > 0.0 && phi < 25.0)*(0.0001) +
571 (energy >= 0.0 && energy < 25.0&& phi >= 25.0 && phi < 50.0)*(0.0006) +
572 (energy >= 0.0 && energy < 25.0&& phi >= 50.0 && phi < 75.0)*(0.0014) +
573 (energy >= 0.0 && energy < 25.0&& phi >= 75.0 && phi < 100.0)*(0.0025) +
574 (energy >= 0.0 && energy < 25.0&& phi >= 100.0 && phi < 150.0)*(0.0046) +
575 (energy >= 0.0 && energy < 25.0&& phi >= 150.0 && phi < 200.0)*(0.0060) +
576 (energy >= 0.0 && energy < 25.0 && phi >= 200.0)*(0.0136) +
577 (energy >= 25.0 && energy < 50.0 && 0.0 == phi)*(0.0000) +
578 (energy >= 25.0 && energy < 50.0&& phi > 0.0 && phi < 25.0)*(0.0000) +
579 (energy >= 25.0 && energy < 50.0&& phi >= 25.0 && phi < 50.0)*(0.0006) +
580 (energy >= 25.0 && energy < 50.0&& phi >= 50.0 && phi < 75.0)*(0.0015) +
581 (energy >= 25.0 && energy < 50.0&& phi >= 75.0 && phi < 100.0)*(0.0033) +
582 (energy >= 25.0 && energy < 50.0&& phi >= 100.0 && phi < 150.0)*(0.0051) +
583 (energy >= 25.0 && energy < 50.0&& phi >= 150.0 && phi < 200.0)*(0.0098) +
584 (energy >= 25.0 && energy < 50.0 && phi >= 200.0)*(0.0146) +
585 (energy >= 50.0 && energy < 75.0 && 0.0 == phi)*(0.0000) +
586 (energy >= 50.0 && energy < 75.0&& phi > 0.0 && phi < 25.0)*(0.0001) +
587 (energy >= 50.0 && energy < 75.0&& phi >= 25.0 && phi < 50.0)*(0.0003) +
588 (energy >= 50.0 && energy < 75.0&& phi >= 50.0 && phi < 75.0)*(0.0015) +
589 (energy >= 50.0 && energy < 75.0&& phi >= 75.0 && phi < 100.0)*(0.0038) +
590 (energy >= 50.0 && energy < 75.0&& phi >= 100.0 && phi < 150.0)*(0.0052) +
591 (energy >= 50.0 && energy < 75.0&& phi >= 150.0 && phi < 200.0)*(0.0114) +
592 (energy >= 50.0 && energy < 75.0 && phi >= 200.0)*(0.0181) +
593 (energy >= 75.0 && energy < 100.0 && 0.0 == phi)*(0.0000) +
594 (energy >= 75.0 && energy < 100.0&& phi > 0.0 && phi < 25.0)*(0.0001) +
595 (energy >= 75.0 && energy < 100.0&& phi >= 25.0 && phi < 50.0)*(0.0005) +
596 (energy >= 75.0 && energy < 100.0&& phi >= 50.0 && phi < 75.0)*(0.0022) +
597 (energy >= 75.0 && energy < 100.0&& phi >= 75.0 && phi < 100.0)*(0.0067) +
598 (energy >= 75.0 && energy < 100.0&& phi >= 100.0 && phi < 150.0)*(0.0047) +
599 (energy >= 75.0 && energy < 100.0&& phi >= 150.0 && phi < 200.0)*(0.0113) +
600 (energy >= 75.0 && energy < 100.0 && phi >= 200.0)*(0.0145) +
601 (energy >= 100.0 && energy < 125.0 && 0.0 == phi)*(0.0000) +
602 (energy >= 100.0 && energy < 125.0&& phi > 0.0 && phi < 25.0)*(0.0001) +
603 (energy >= 100.0 && energy < 125.0&& phi >= 25.0 && phi < 50.0)*(0.0003) +
604 (energy >= 100.0 && energy < 125.0&& phi >= 50.0 && phi < 75.0)*(0.0016) +
605 (energy >= 100.0 && energy < 125.0&& phi >= 75.0 && phi < 100.0)*(0.0110) +
606 (energy >= 100.0 && energy < 125.0&& phi >= 100.0 && phi < 150.0)*(0.0029) +
607 (energy >= 100.0 && energy < 125.0&& phi >= 150.0 && phi < 200.0)*(0.0138) +
608 (energy >= 100.0 && energy < 125.0 && phi >= 200.0)*(0.0000) +
609 (energy >= 125.0 && energy < 150.0 && 0.0 == phi)*(0.0000) +
610 (energy >= 125.0 && energy < 150.0&& phi > 0.0 && phi < 25.0)*(0.0000) +
611 (energy >= 125.0 && energy < 150.0&& phi >= 25.0 && phi < 50.0)*(0.0000) +
612 (energy >= 125.0 && energy < 150.0&& phi >= 50.0 && phi < 75.0)*(0.0026) +
613 (energy >= 125.0 && energy < 150.0&& phi >= 75.0 && phi < 100.0)*(0.0047) +
614 (energy >= 125.0 && energy < 150.0&& phi >= 100.0 && phi < 150.0)*(0.0085) +
615 (energy >= 125.0 && energy < 150.0&& phi >= 150.0 && phi < 200.0)*(0.0152) +
616 (energy >= 125.0 && energy < 150.0 && phi >= 200.0)*(0.0164) +
617 (energy >= 150.0 && 0.0 == phi)*(0.0000) +
618 (energy >= 150.0&& phi > 0.0 && phi < 25.0)*(0.0000) +
619 (energy >= 150.0&& phi >= 25.0 && phi < 50.0)*(0.0000) +
620 (energy >= 150.0&& phi >= 50.0 && phi < 75.0)*(0.0000) +
621 (energy >= 150.0&& phi >= 75.0 && phi < 100.0)*(0.0000) +
622 (energy >= 150.0&& phi >= 100.0 && phi < 150.0)*(0.0000) +
623 (energy >= 150.0&& phi >= 150.0 && phi < 200.0)*(0.0080) +
624 (energy >= 150.0 && phi >= 200.0)*(0.0143))
625 }
626}
627
628module CscClusterId CutBasedIDEfficiency {
629set InputArray CSCFilter/LLP
630 set OutputArray cluster
631
632 # efficiency formula for Csc Cluster, as a function of LLP decay vertex in R, Z and hadronic and EM energy
633 set EfficiencyFormula {
634 (pt > 3900 && eta < 6710) * ((0.0 == phi)*(0.0656) +
635 (phi > 0.0 && phi < 25.0)*(0.0777) +
636 (phi >= 25.0 && phi < 50.0)*(0.1607) +
637 (phi >= 50.0 && phi < 75.0)*(0.2294) +
638 (phi >= 75.0 && phi < 100.0)*(0.3146) +
639 (phi >= 100.0 && phi < 150.0)*(0.3107) +
640 (phi >= 150.0 && phi < 200.0)*(0.3177) +
641 (phi >= 200.0)*(0.3229)) +
642 (eta > 6710) * ((0.0 == phi)*(0.2987) +
643 (phi > 0.0 && phi < 25.0)*(0.3100) +
644 (phi >= 25.0 && phi < 50.0)*(0.4476) +
645 (phi >= 50.0 && phi < 75.0)*(0.5335) +
646 (phi >= 75.0 && phi < 100.0)*(0.5961) +
647 (phi >= 100.0 && phi < 150.0)*(0.6368) +
648 (phi >= 150.0 && phi < 200.0)*(0.6814) +
649 (phi >= 200.0)*(0.6998)) +
650 (pt < 2700 && eta < 6710) * ((0.0 == phi)*(0.8604) +
651 (phi > 0.0 && phi < 25.0)*(0.3335) +
652 (phi >= 25.0 && phi < 50.0)*(0.2457) +
653 (phi >= 50.0 && phi < 75.0)*(0.1831) +
654 (phi >= 75.0 && phi < 100.0)*(0.2100) +
655 (phi >= 100.0 && phi < 150.0)*(0.2443) +
656 (phi >= 150.0 && phi < 200.0)*(0.2532) +
657 (phi >= 200.0)*(0.2404))
658 }
659}
660######################
661# ChargedHadronFilter
662######################
663
664module PdgCodeFilter ChargedHadronFilter {
665 set InputArray HCal/eflowTracks
666 set OutputArray chargedHadrons
667
668 add PdgCode {11}
669 add PdgCode {-11}
670 add PdgCode {13}
671 add PdgCode {-13}
672}
673
674
675###################################################
676# Tower Merger (in case not using e-flow algorithm)
677###################################################
678
679module Merger Calorimeter {
680# add InputArray InputArray
681 add InputArray ECal/ecalTowers
682 add InputArray HCal/hcalTowers
683 set OutputArray towers
684}
685
686
687
688####################
689# Energy flow merger
690####################
691
692module Merger EFlowMerger {
693# add InputArray InputArray
694 add InputArray HCal/eflowTracks
695 add InputArray ECal/eflowPhotons
696 add InputArray HCal/eflowNeutralHadrons
697 set OutputArray eflow
698}
699
700######################
701# EFlowFilter
702######################
703
704module PdgCodeFilter EFlowFilter {
705 set InputArray EFlowMerger/eflow
706 set OutputArray eflow
707
708 add PdgCode {11}
709 add PdgCode {-11}
710 add PdgCode {13}
711 add PdgCode {-13}
712}
713
714
715###################
716# Photon efficiency
717###################
718
719module Efficiency PhotonEfficiency {
720 set InputArray ECal/eflowPhotons
721 set OutputArray photons
722
723 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
724
725 # efficiency formula for photons
726 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
727 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
728 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
729 (abs(eta) > 2.5) * (0.00)}
730}
731
732##################
733# Photon isolation
734##################
735
736module Isolation PhotonIsolation {
737 set CandidateInputArray PhotonEfficiency/photons
738 set IsolationInputArray EFlowFilter/eflow
739
740 set OutputArray photons
741
742 set DeltaRMax 0.5
743
744 set PTMin 0.5
745
746 set PTRatioMax 0.12
747}
748
749
750#####################
751# Electron efficiency
752#####################
753
754module Efficiency ElectronEfficiency {
755 set InputArray ElectronFilter/electrons
756 set OutputArray electrons
757
758 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
759
760 # efficiency formula for electrons
761 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
762 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
763 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
764 (abs(eta) > 2.5) * (0.00)}
765}
766
767####################
768# Electron isolation
769####################
770
771module Isolation ElectronIsolation {
772 set CandidateInputArray ElectronEfficiency/electrons
773 set IsolationInputArray EFlowFilter/eflow
774
775 set OutputArray electrons
776
777 set DeltaRMax 0.5
778
779 set PTMin 0.5
780
781 set PTRatioMax 0.12
782}
783
784#################
785# Muon efficiency
786#################
787
788module Efficiency MuonEfficiency {
789 set InputArray MuonMomentumSmearing/muons
790 set OutputArray muons
791
792 # set EfficiencyFormula {efficiency as a function of eta and pt}
793
794 # efficiency formula for muons
795 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
796 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
797 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0) * (0.95) +
798 (abs(eta) > 2.4) * (0.00)}
799}
800
801################
802# Muon isolation
803################
804
805module Isolation MuonIsolation {
806 set CandidateInputArray MuonEfficiency/muons
807 set IsolationInputArray EFlowFilter/eflow
808
809 set OutputArray muons
810
811 set DeltaRMax 0.5
812
813 set PTMin 0.5
814
815 set PTRatioMax 0.25
816}
817
818###################
819# Missing ET merger
820###################
821
822module Merger MissingET {
823# add InputArray InputArray
824 add InputArray EFlowMerger/eflow
825 set MomentumOutputArray momentum
826}
827
828##################
829# Scalar HT merger
830##################
831
832module Merger ScalarHT {
833# add InputArray InputArray
834 add InputArray UniqueObjectFinder/jets
835 add InputArray UniqueObjectFinder/electrons
836 add InputArray UniqueObjectFinder/photons
837 add InputArray UniqueObjectFinder/muons
838 set EnergyOutputArray energy
839}
840
841
842#####################
843# Neutrino Filter
844#####################
845
846module PdgCodeFilter NeutrinoFilter {
847
848 set InputArray Delphes/stableParticles
849 set OutputArray filteredParticles
850
851 set PTMin 0.0
852
853 add PdgCode {12}
854 add PdgCode {14}
855 add PdgCode {16}
856 add PdgCode {-12}
857 add PdgCode {-14}
858 add PdgCode {-16}
859
860}
861
862
863#####################
864# MC truth jet finder
865#####################
866
867module FastJetFinder GenJetFinder {
868 set InputArray NeutrinoFilter/filteredParticles
869
870 set OutputArray jets
871
872 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
873 set JetAlgorithm 6
874 set ParameterR 0.5
875
876 set JetPTMin 20.0
877}
878
879#########################
880# Gen Missing ET merger
881########################
882
883module Merger GenMissingET {
884# add InputArray InputArray
885 add InputArray NeutrinoFilter/filteredParticles
886
887 set MomentumOutputArray momentum
888}
889
890
891
892############
893# Jet finder
894############
895
896module FastJetFinder FastJetFinder {
897# set InputArray Calorimeter/towers
898 set InputArray EFlowMerger/eflow
899
900 set OutputArray jets
901
902 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
903 set JetAlgorithm 6
904 set ParameterR 0.5
905
906 set JetPTMin 20.0
907}
908
909##################
910# Fat Jet finder
911##################
912
913module FastJetFinder FatJetFinder {
914 set InputArray EFlowMerger/eflow
915
916 set OutputArray jets
917
918 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
919 set JetAlgorithm 6
920 set ParameterR 0.8
921
922 set ComputeNsubjettiness 1
923 set Beta 1.0
924 set AxisMode 4
925
926 set ComputeTrimming 1
927 set RTrim 0.2
928 set PtFracTrim 0.05
929
930 set ComputePruning 1
931 set ZcutPrun 0.1
932 set RcutPrun 0.5
933 set RPrun 0.8
934
935 set ComputeSoftDrop 1
936 set BetaSoftDrop 0.0
937 set SymmetryCutSoftDrop 0.1
938 set R0SoftDrop 0.8
939
940 set JetPTMin 200.0
941}
942
943
944
945
946##################
947# Jet Energy Scale
948##################
949
950module EnergyScale JetEnergyScale {
951 set InputArray FastJetFinder/jets
952 set OutputArray jets
953
954 # scale formula for jets
955 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
956}
957
958########################
959# Jet Flavor Association
960########################
961
962module JetFlavorAssociation JetFlavorAssociation {
963
964 set PartonInputArray Delphes/partons
965 set ParticleInputArray Delphes/allParticles
966 set ParticleLHEFInputArray Delphes/allParticlesLHEF
967 set JetInputArray JetEnergyScale/jets
968
969 set DeltaR 0.5
970 set PartonPTMin 1.0
971 set PartonEtaMax 2.5
972
973}
974
975###########
976# b-tagging
977###########
978
979module BTagging BTagging {
980 set JetInputArray JetEnergyScale/jets
981
982 set BitNumber 0
983
984 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
985 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
986 # gluon's PDG code has the lowest priority
987
988 # based on arXiv:1211.4462
989
990 # default efficiency formula (misidentification rate)
991 add EfficiencyFormula {0} {0.01+0.000038*pt}
992
993 # efficiency formula for c-jets (misidentification rate)
994 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
995
996 # efficiency formula for b-jets
997 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
998}
999
1000#############
1001# tau-tagging
1002#############
1003
1004module TauTagging TauTagging {
1005 set ParticleInputArray Delphes/allParticles
1006 set PartonInputArray Delphes/partons
1007 set JetInputArray JetEnergyScale/jets
1008
1009 set DeltaR 0.5
1010
1011 set TauPTMin 1.0
1012
1013 set TauEtaMax 2.5
1014
1015 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1016
1017 # default efficiency formula (misidentification rate)
1018 add EfficiencyFormula {0} {0.01}
1019 # efficiency formula for tau-jets
1020 add EfficiencyFormula {15} {0.6}
1021}
1022
1023#####################################################
1024# Find uniquely identified photons/electrons/tau/jets
1025#####################################################
1026
1027module UniqueObjectFinder UniqueObjectFinder {
1028# earlier arrays take precedence over later ones
1029# add InputArray InputArray OutputArray
1030 add InputArray PhotonIsolation/photons photons
1031 add InputArray ElectronIsolation/electrons electrons
1032 add InputArray MuonIsolation/muons muons
1033 add InputArray JetEnergyScale/jets jets
1034}
1035
1036##################
1037# ROOT tree writer
1038##################
1039
1040# tracks, towers and eflow objects are not stored by default in the output.
1041# if needed (for jet constituent or other studies), uncomment the relevant
1042# "add Branch ..." lines.
1043
1044module TreeWriter TreeWriter {
1045# add Branch InputArray BranchName BranchClass
1046#add Branch Delphes/allParticles Particle GenParticle
1047 #add Branch Delphes/stableParticles Particle GenParticle
1048 #add Branch TrackMerger/tracks Track Track
1049 #add Branch Calorimeter/towers Tower Tower
1050
1051 #add Branch HCal/eflowTracks EFlowTrack Track
1052 #add Branch ECal/eflowPhotons EFlowPhoton Tower
1053 #add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1054
1055 #add Branch GenJetFinder/jets GenJet Jet
1056 #add Branch GenMissingET/momentum GenMissingET MissingET
1057 add Branch HiggsFilter/higgs higgs GenParticle
1058
1059
1060
1061 add Branch UniqueObjectFinder/jets Jet Jet
1062 add Branch UniqueObjectFinder/electrons Electron Electron
1063 #add Branch UniqueObjectFinder/photons Photon Photon
1064 add Branch UniqueObjectFinder/muons Muon Muon
1065
1066 #add Branch FatJetFinder/jets FatJet Jet
1067
1068 add Branch MissingET/momentum MissingET MissingET
1069 #add Branch ScalarHT/energy ScalarHT ScalarHT
1070 add Branch llpFilter/LLP llp CscCluster
1071 add Branch CSCFilter/LLP Cscllp CscCluster
1072
1073 add Branch ClusterEfficiency/cluster CscCluster130 CscCluster
1074}
Note: See TracBrowser for help on using the repository browser.