Fork me on GitHub

source: git/cards/delphes_card_CMS_CSCCluster.tcl@ c61b5ce

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

cleaning
  • Property mode set to 100644
File size: 35.8 KB
Line 
1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6 ParticlePropagator
7
8 ChargedHadronTrackingEfficiency
9 ElectronTrackingEfficiency
10 MuonTrackingEfficiency
11
12 ChargedHadronMomentumSmearing
13 ElectronMomentumSmearing
14 MuonMomentumSmearing
15
16 TrackMerger
17
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 # DecayRegion = 0: no cuts on decay region
414 # DecayRegion = 1: select LLP that decays in CSC volume
415 # DecayRegion = 2: select LLP that decays outside of calorimeters, for genMET calculation
416 set DecayRegion 1
417 set RequireStatus false
418 add PdgCode {1500001}
419
420
421}
422# filter out LLPs regardless of decay position
423module LLPFilter llpFilter {
424 set InputArray Delphes/allParticles
425 set OutputArray LLP
426 # DecayRegion = 0: no cuts on decay region
427 # DecayRegion = 1: select LLP that decays in CSC volume
428 # DecayRegion = 2: select LLP that decays outside of calorimeters, for genMET calculation
429 set DecayRegion 0
430 set RequireStatus false
431 add PdgCode {1500001}
432
433}
434
435
436
437module CscClusterEfficiency ClusterEfficiency {
438 set InputArray CutBasedIDEfficiency/cluster
439 set OutputArray cluster
440
441 # efficiency formula for Csc Cluster, as a function of LLP decay vertex in R, Z and hadronic and EM energy
442 set EfficiencyFormula {
443
444 (decayR > 3900 && decayZ < 6710) * ((Eem >= 0.0 && Eem < 25.0 && 0.0 == Ehad)*(0.0049) +
445 (Eem >= 0.0 && Eem < 25.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0130) +
446 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0346) +
447 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0623) +
448 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0919) +
449 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.1086) +
450 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.1292) +
451 (Eem >= 0.0 && Eem < 25.0 && Ehad >= 200.0)*(0.1106) +
452 (Eem >= 25.0 && Eem < 50.0 && 0.0 == Ehad)*(0.0249) +
453 (Eem >= 25.0 && Eem < 50.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0285) +
454 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0501) +
455 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0841) +
456 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.1021) +
457 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.1129) +
458 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.1141) +
459 (Eem >= 25.0 && Eem < 50.0 && Ehad >= 200.0)*(0.1370) +
460 (Eem >= 50.0 && Eem < 75.0 && 0.0 == Ehad)*(0.0282) +
461 (Eem >= 50.0 && Eem < 75.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0445) +
462 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0643) +
463 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0903) +
464 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0998) +
465 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.1420) +
466 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.1429) +
467 (Eem >= 50.0 && Eem < 75.0 && Ehad >= 200.0)*(0.0882) +
468 (Eem >= 75.0 && Eem < 100.0 && 0.0 == Ehad)*(0.0594) +
469 (Eem >= 75.0 && Eem < 100.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0521) +
470 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0605) +
471 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0791) +
472 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.1117) +
473 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0862) +
474 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0698) +
475 (Eem >= 75.0 && Eem < 100.0 && Ehad >= 200.0)*(0.0500) +
476 (Eem >= 100.0 && Eem < 125.0 && 0.0 == Ehad)*(0.0758) +
477 (Eem >= 100.0 && Eem < 125.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0414) +
478 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0755) +
479 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.1027) +
480 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0440) +
481 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0811) +
482 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.1538) +
483 (Eem >= 100.0 && Eem < 125.0 && Ehad >= 200.0)*(0.0833) +
484 (Eem >= 125.0 && Eem < 150.0 && 0.0 == Ehad)*(0.0300) +
485 (Eem >= 125.0 && Eem < 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0609) +
486 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0745) +
487 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0610) +
488 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.1224) +
489 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.1667) +
490 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0000) +
491 (Eem >= 125.0 && Eem < 150.0 && Ehad >= 200.0)*(0.0000) +
492 (Eem >= 150.0 && 0.0 == Ehad)*(0.0282) +
493 (Eem >= 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0809) +
494 (Eem >= 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0352) +
495 (Eem >= 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0984) +
496 (Eem >= 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0968) +
497 (Eem >= 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.1282) +
498 (Eem >= 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.2105) +
499 (Eem >= 150.0 && Ehad >= 200.0)*(0.0769)) +
500 (decayZ > 6710) * ((Eem >= 0.0 && Eem < 25.0 && 0.0 == Ehad)*(0.0184) +
501 (Eem >= 0.0 && Eem < 25.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0772) +
502 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.2086) +
503 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3091) +
504 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.3867) +
505 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4500) +
506 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.4746) +
507 (Eem >= 0.0 && Eem < 25.0 && Ehad >= 200.0)*(0.4906) +
508 (Eem >= 25.0 && Eem < 50.0 && 0.0 == Ehad)*(0.0955) +
509 (Eem >= 25.0 && Eem < 50.0&& Ehad > 0.0 && Ehad < 25.0)*(0.1461) +
510 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.2594) +
511 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3556) +
512 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4165) +
513 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4693) +
514 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.5054) +
515 (Eem >= 25.0 && Eem < 50.0 && Ehad >= 200.0)*(0.5219) +
516 (Eem >= 50.0 && Eem < 75.0 && 0.0 == Ehad)*(0.1472) +
517 (Eem >= 50.0 && Eem < 75.0&& Ehad > 0.0 && Ehad < 25.0)*(0.1970) +
518 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.2974) +
519 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3783) +
520 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4335) +
521 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4736) +
522 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.4937) +
523 (Eem >= 50.0 && Eem < 75.0 && Ehad >= 200.0)*(0.5077) +
524 (Eem >= 75.0 && Eem < 100.0 && 0.0 == Ehad)*(0.2053) +
525 (Eem >= 75.0 && Eem < 100.0&& Ehad > 0.0 && Ehad < 25.0)*(0.2314) +
526 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.3114) +
527 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3799) +
528 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4420) +
529 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4502) +
530 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.5348) +
531 (Eem >= 75.0 && Eem < 100.0 && Ehad >= 200.0)*(0.5115) +
532 (Eem >= 100.0 && Eem < 125.0 && 0.0 == Ehad)*(0.2198) +
533 (Eem >= 100.0 && Eem < 125.0&& Ehad > 0.0 && Ehad < 25.0)*(0.2404) +
534 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.3295) +
535 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3932) +
536 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4327) +
537 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4377) +
538 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.5175) +
539 (Eem >= 100.0 && Eem < 125.0 && Ehad >= 200.0)*(0.6087) +
540 (Eem >= 125.0 && Eem < 150.0 && 0.0 == Ehad)*(0.2147) +
541 (Eem >= 125.0 && Eem < 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.2605) +
542 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.3442) +
543 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3622) +
544 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4407) +
545 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.5168) +
546 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.5056) +
547 (Eem >= 125.0 && Eem < 150.0 && Ehad >= 200.0)*(0.4559) +
548 (Eem >= 150.0 && 0.0 == Ehad)*(0.2824) +
549 (Eem >= 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.2447) +
550 (Eem >= 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.3519) +
551 (Eem >= 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.3772) +
552 (Eem >= 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.4447) +
553 (Eem >= 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.4703) +
554 (Eem >= 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.4460) +
555 (Eem >= 150.0 && Ehad >= 200.0)*(0.4400)) +
556 (decayR < 2700 && decayZ < 6710) * ((Eem >= 0.0 && Eem < 25.0 && 0.0 == Ehad)*(0.0002) +
557 (Eem >= 0.0 && Eem < 25.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0001) +
558 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0006) +
559 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0014) +
560 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0025) +
561 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0046) +
562 (Eem >= 0.0 && Eem < 25.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0060) +
563 (Eem >= 0.0 && Eem < 25.0 && Ehad >= 200.0)*(0.0136) +
564 (Eem >= 25.0 && Eem < 50.0 && 0.0 == Ehad)*(0.0000) +
565 (Eem >= 25.0 && Eem < 50.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0000) +
566 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0006) +
567 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0015) +
568 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0033) +
569 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0051) +
570 (Eem >= 25.0 && Eem < 50.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0098) +
571 (Eem >= 25.0 && Eem < 50.0 && Ehad >= 200.0)*(0.0146) +
572 (Eem >= 50.0 && Eem < 75.0 && 0.0 == Ehad)*(0.0000) +
573 (Eem >= 50.0 && Eem < 75.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0001) +
574 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0003) +
575 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0015) +
576 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0038) +
577 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0052) +
578 (Eem >= 50.0 && Eem < 75.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0114) +
579 (Eem >= 50.0 && Eem < 75.0 && Ehad >= 200.0)*(0.0181) +
580 (Eem >= 75.0 && Eem < 100.0 && 0.0 == Ehad)*(0.0000) +
581 (Eem >= 75.0 && Eem < 100.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0001) +
582 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0005) +
583 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0022) +
584 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0067) +
585 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0047) +
586 (Eem >= 75.0 && Eem < 100.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0113) +
587 (Eem >= 75.0 && Eem < 100.0 && Ehad >= 200.0)*(0.0145) +
588 (Eem >= 100.0 && Eem < 125.0 && 0.0 == Ehad)*(0.0000) +
589 (Eem >= 100.0 && Eem < 125.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0001) +
590 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0003) +
591 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0016) +
592 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0110) +
593 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0029) +
594 (Eem >= 100.0 && Eem < 125.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0138) +
595 (Eem >= 100.0 && Eem < 125.0 && Ehad >= 200.0)*(0.0000) +
596 (Eem >= 125.0 && Eem < 150.0 && 0.0 == Ehad)*(0.0000) +
597 (Eem >= 125.0 && Eem < 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0000) +
598 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0000) +
599 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0026) +
600 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0047) +
601 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0085) +
602 (Eem >= 125.0 && Eem < 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0152) +
603 (Eem >= 125.0 && Eem < 150.0 && Ehad >= 200.0)*(0.0164) +
604 (Eem >= 150.0 && 0.0 == Ehad)*(0.0000) +
605 (Eem >= 150.0&& Ehad > 0.0 && Ehad < 25.0)*(0.0000) +
606 (Eem >= 150.0&& Ehad >= 25.0 && Ehad < 50.0)*(0.0000) +
607 (Eem >= 150.0&& Ehad >= 50.0 && Ehad < 75.0)*(0.0000) +
608 (Eem >= 150.0&& Ehad >= 75.0 && Ehad < 100.0)*(0.0000) +
609 (Eem >= 150.0&& Ehad >= 100.0 && Ehad < 150.0)*(0.0000) +
610 (Eem >= 150.0&& Ehad >= 150.0 && Ehad < 200.0)*(0.0080) +
611 (Eem >= 150.0 && Ehad >= 200.0)*(0.0143)) }
612}
613
614module CscClusterId CutBasedIDEfficiency {
615set InputArray CSCFilter/LLP
616 set OutputArray cluster
617
618 # efficiency formula for Csc Cluster, as a function of LLP decay vertex in R, Z and hadronic and EM energy
619 set EfficiencyFormula {
620 (decayR > 3900 && decayZ < 6710) * ((0.0 == Ehad)*(0.0656) +
621 (Ehad > 0.0 && Ehad < 25.0)*(0.0777) +
622 (Ehad >= 25.0 && Ehad < 50.0)*(0.1607) +
623 (Ehad >= 50.0 && Ehad < 75.0)*(0.2294) +
624 (Ehad >= 75.0 && Ehad < 100.0)*(0.3146) +
625 (Ehad >= 100.0 && Ehad < 150.0)*(0.3107) +
626 (Ehad >= 150.0 && Ehad < 200.0)*(0.3177) +
627 (Ehad >= 200.0)*(0.3229)) +
628 (decayZ > 6710) * ((0.0 == Ehad)*(0.2987) +
629 (Ehad > 0.0 && Ehad < 25.0)*(0.3100) +
630 (Ehad >= 25.0 && Ehad < 50.0)*(0.4476) +
631 (Ehad >= 50.0 && Ehad < 75.0)*(0.5335) +
632 (Ehad >= 75.0 && Ehad < 100.0)*(0.5961) +
633 (Ehad >= 100.0 && Ehad < 150.0)*(0.6368) +
634 (Ehad >= 150.0 && Ehad < 200.0)*(0.6814) +
635 (Ehad >= 200.0)*(0.6998)) +
636 (decayR < 2700 && decayZ < 6710) * ((0.0 == Ehad)*(0.8604) +
637 (Ehad > 0.0 && Ehad < 25.0)*(0.3335) +
638 (Ehad >= 25.0 && Ehad < 50.0)*(0.2457) +
639 (Ehad >= 50.0 && Ehad < 75.0)*(0.1831) +
640 (Ehad >= 75.0 && Ehad < 100.0)*(0.2100) +
641 (Ehad >= 100.0 && Ehad < 150.0)*(0.2443) +
642 (Ehad >= 150.0 && Ehad < 200.0)*(0.2532) +
643 (Ehad >= 200.0)*(0.2404))
644 }
645}
646######################
647# ChargedHadronFilter
648######################
649
650module PdgCodeFilter ChargedHadronFilter {
651 set InputArray HCal/eflowTracks
652 set OutputArray chargedHadrons
653
654 add PdgCode {11}
655 add PdgCode {-11}
656 add PdgCode {13}
657 add PdgCode {-13}
658}
659
660
661###################################################
662# Tower Merger (in case not using e-flow algorithm)
663###################################################
664
665module Merger Calorimeter {
666# add InputArray InputArray
667 add InputArray ECal/ecalTowers
668 add InputArray HCal/hcalTowers
669 set OutputArray towers
670}
671
672
673
674####################
675# Energy flow merger
676####################
677
678module Merger EFlowMerger {
679# add InputArray InputArray
680 add InputArray HCal/eflowTracks
681 add InputArray ECal/eflowPhotons
682 add InputArray HCal/eflowNeutralHadrons
683 set OutputArray eflow
684}
685
686######################
687# EFlowFilter
688######################
689
690module PdgCodeFilter EFlowFilter {
691 set InputArray EFlowMerger/eflow
692 set OutputArray eflow
693
694 add PdgCode {11}
695 add PdgCode {-11}
696 add PdgCode {13}
697 add PdgCode {-13}
698}
699
700
701###################
702# Photon efficiency
703###################
704
705module Efficiency PhotonEfficiency {
706 set InputArray ECal/eflowPhotons
707 set OutputArray photons
708
709 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
710
711 # efficiency formula for photons
712 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
713 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
714 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
715 (abs(eta) > 2.5) * (0.00)}
716}
717
718##################
719# Photon isolation
720##################
721
722module Isolation PhotonIsolation {
723 set CandidateInputArray PhotonEfficiency/photons
724 set IsolationInputArray EFlowFilter/eflow
725
726 set OutputArray photons
727
728 set DeltaRMax 0.5
729
730 set PTMin 0.5
731
732 set PTRatioMax 0.12
733}
734
735
736#####################
737# Electron efficiency
738#####################
739
740module Efficiency ElectronEfficiency {
741 set InputArray ElectronFilter/electrons
742 set OutputArray electrons
743
744 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
745
746 # efficiency formula for electrons
747 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
748 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
749 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
750 (abs(eta) > 2.5) * (0.00)}
751}
752
753####################
754# Electron isolation
755####################
756
757module Isolation ElectronIsolation {
758 set CandidateInputArray ElectronEfficiency/electrons
759 set IsolationInputArray EFlowFilter/eflow
760
761 set OutputArray electrons
762
763 set DeltaRMax 0.5
764
765 set PTMin 0.5
766
767 set PTRatioMax 0.12
768}
769
770#################
771# Muon efficiency
772#################
773
774module Efficiency MuonEfficiency {
775 set InputArray MuonMomentumSmearing/muons
776 set OutputArray muons
777
778 # set EfficiencyFormula {efficiency as a function of eta and pt}
779
780 # efficiency formula for muons
781 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
782 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
783 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0) * (0.95) +
784 (abs(eta) > 2.4) * (0.00)}
785}
786
787################
788# Muon isolation
789################
790
791module Isolation MuonIsolation {
792 set CandidateInputArray MuonEfficiency/muons
793 set IsolationInputArray EFlowFilter/eflow
794
795 set OutputArray muons
796
797 set DeltaRMax 0.5
798
799 set PTMin 0.5
800
801 set PTRatioMax 0.25
802}
803
804###################
805# Missing ET merger
806###################
807
808module Merger MissingET {
809# add InputArray InputArray
810 add InputArray EFlowMerger/eflow
811 set MomentumOutputArray momentum
812}
813
814##################
815# Scalar HT merger
816##################
817
818module Merger ScalarHT {
819# add InputArray InputArray
820 add InputArray UniqueObjectFinder/jets
821 add InputArray UniqueObjectFinder/electrons
822 add InputArray UniqueObjectFinder/photons
823 add InputArray UniqueObjectFinder/muons
824 set EnergyOutputArray energy
825}
826
827
828#####################
829# Neutrino Filter
830#####################
831
832module PdgCodeFilter NeutrinoFilter {
833
834 set InputArray Delphes/stableParticles
835 set OutputArray filteredParticles
836
837 set PTMin 0.0
838
839 add PdgCode {12}
840 add PdgCode {14}
841 add PdgCode {16}
842 add PdgCode {-12}
843 add PdgCode {-14}
844 add PdgCode {-16}
845
846}
847
848
849#####################
850# MC truth jet finder
851#####################
852
853module FastJetFinder GenJetFinder {
854 set InputArray NeutrinoFilter/filteredParticles
855
856 set OutputArray jets
857
858 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
859 set JetAlgorithm 6
860 set ParameterR 0.5
861
862 set JetPTMin 20.0
863}
864
865#########################
866# Gen Missing ET merger
867########################
868
869module Merger GenMissingET {
870# add InputArray InputArray
871 add InputArray NeutrinoFilter/filteredParticles
872
873 set MomentumOutputArray momentum
874}
875
876
877
878############
879# Jet finder
880############
881
882module FastJetFinder FastJetFinder {
883# set InputArray Calorimeter/towers
884 set InputArray EFlowMerger/eflow
885
886 set OutputArray jets
887
888 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
889 set JetAlgorithm 6
890 set ParameterR 0.5
891
892 set JetPTMin 20.0
893}
894
895##################
896# Fat Jet finder
897##################
898
899module FastJetFinder FatJetFinder {
900 set InputArray EFlowMerger/eflow
901
902 set OutputArray jets
903
904 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
905 set JetAlgorithm 6
906 set ParameterR 0.8
907
908 set ComputeNsubjettiness 1
909 set Beta 1.0
910 set AxisMode 4
911
912 set ComputeTrimming 1
913 set RTrim 0.2
914 set PtFracTrim 0.05
915
916 set ComputePruning 1
917 set ZcutPrun 0.1
918 set RcutPrun 0.5
919 set RPrun 0.8
920
921 set ComputeSoftDrop 1
922 set BetaSoftDrop 0.0
923 set SymmetryCutSoftDrop 0.1
924 set R0SoftDrop 0.8
925
926 set JetPTMin 200.0
927}
928
929
930
931
932##################
933# Jet Energy Scale
934##################
935
936module EnergyScale JetEnergyScale {
937 set InputArray FastJetFinder/jets
938 set OutputArray jets
939
940 # scale formula for jets
941 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
942}
943
944########################
945# Jet Flavor Association
946########################
947
948module JetFlavorAssociation JetFlavorAssociation {
949
950 set PartonInputArray Delphes/partons
951 set ParticleInputArray Delphes/allParticles
952 set ParticleLHEFInputArray Delphes/allParticlesLHEF
953 set JetInputArray JetEnergyScale/jets
954
955 set DeltaR 0.5
956 set PartonPTMin 1.0
957 set PartonEtaMax 2.5
958
959}
960
961###########
962# b-tagging
963###########
964
965module BTagging BTagging {
966 set JetInputArray JetEnergyScale/jets
967
968 set BitNumber 0
969
970 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
971 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
972 # gluon's PDG code has the lowest priority
973
974 # based on arXiv:1211.4462
975
976 # default efficiency formula (misidentification rate)
977 add EfficiencyFormula {0} {0.01+0.000038*pt}
978
979 # efficiency formula for c-jets (misidentification rate)
980 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
981
982 # efficiency formula for b-jets
983 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
984}
985
986#############
987# tau-tagging
988#############
989
990module TauTagging TauTagging {
991 set ParticleInputArray Delphes/allParticles
992 set PartonInputArray Delphes/partons
993 set JetInputArray JetEnergyScale/jets
994
995 set DeltaR 0.5
996
997 set TauPTMin 1.0
998
999 set TauEtaMax 2.5
1000
1001 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1002
1003 # default efficiency formula (misidentification rate)
1004 add EfficiencyFormula {0} {0.01}
1005 # efficiency formula for tau-jets
1006 add EfficiencyFormula {15} {0.6}
1007}
1008
1009#####################################################
1010# Find uniquely identified photons/electrons/tau/jets
1011#####################################################
1012
1013module UniqueObjectFinder UniqueObjectFinder {
1014# earlier arrays take precedence over later ones
1015# add InputArray InputArray OutputArray
1016 add InputArray PhotonIsolation/photons photons
1017 add InputArray ElectronIsolation/electrons electrons
1018 add InputArray MuonIsolation/muons muons
1019 add InputArray JetEnergyScale/jets jets
1020}
1021
1022##################
1023# ROOT tree writer
1024##################
1025
1026# tracks, towers and eflow objects are not stored by default in the output.
1027# if needed (for jet constituent or other studies), uncomment the relevant
1028# "add Branch ..." lines.
1029
1030module TreeWriter TreeWriter {
1031# add Branch InputArray BranchName BranchClass
1032# add Branch Delphes/allParticles Particle GenParticle
1033
1034 add Branch UniqueObjectFinder/jets Jet Jet
1035 add Branch UniqueObjectFinder/electrons Electron Electron
1036 add Branch UniqueObjectFinder/muons Muon Muon
1037
1038 add Branch MissingET/momentum MissingET MissingET
1039 add Branch llpFilter/LLP llp CscCluster
1040 add Branch CSCFilter/LLP Cscllp CscCluster
1041 add Branch ClusterEfficiency/cluster CscCluster CscCluster
1042}
Note: See TracBrowser for help on using the repository browser.