Fork me on GitHub

source: git/cards/delphes_card_CMS_CSCCluster.tcl@ a09b75f

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

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