Fork me on GitHub

source: git/cards/delphes_card_CMS_CSCCluster.tcl

Last change on this file was cc8716b, checked in by GitHub <noreply@…>, 3 years ago

Update to handle CMS endcap muon detector showers for long-lived particles (#103)

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