Fork me on GitHub

source: git/cards/delphes_card_CMS_CSCCluster.tcl@ 157fe13

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

LLPFilter hard-coded decay region requirements are moved to configuration cards

  • Property mode set to 100644
File size: 35.7 KB
RevLine 
[a5af1df]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
[29b722a]65
[a5af1df]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
[157fe13]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
[29b722a]425 set RequireStatus false
426 add PdgCode {1500001}
427
[a5af1df]428
429}
430# filter out LLPs regardless of decay position
431module LLPFilter llpFilter {
432 set InputArray Delphes/allParticles
433 set OutputArray LLP
[157fe13]434 set RequireDecayRegion false
435
436
[c61b5ce]437 set RequireStatus false
[29b722a]438 add PdgCode {1500001}
[a5af1df]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
[29b722a]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)) }
[a5af1df]619}
620
621module CscClusterId CutBasedIDEfficiency {
622set InputArray CSCFilter/LLP
623 set OutputArray cluster
624
625 # efficiency formula for Csc Cluster, as a function of LLP decay vertex in R, Z and hadronic and EM energy
626 set EfficiencyFormula {
[29b722a]627 (decayR > 3900 && decayZ < 6710) * ((0.0 == Ehad)*(0.0656) +
628 (Ehad > 0.0 && Ehad < 25.0)*(0.0777) +
629 (Ehad >= 25.0 && Ehad < 50.0)*(0.1607) +
630 (Ehad >= 50.0 && Ehad < 75.0)*(0.2294) +
631 (Ehad >= 75.0 && Ehad < 100.0)*(0.3146) +
632 (Ehad >= 100.0 && Ehad < 150.0)*(0.3107) +
633 (Ehad >= 150.0 && Ehad < 200.0)*(0.3177) +
634 (Ehad >= 200.0)*(0.3229)) +
635 (decayZ > 6710) * ((0.0 == Ehad)*(0.2987) +
636 (Ehad > 0.0 && Ehad < 25.0)*(0.3100) +
637 (Ehad >= 25.0 && Ehad < 50.0)*(0.4476) +
638 (Ehad >= 50.0 && Ehad < 75.0)*(0.5335) +
639 (Ehad >= 75.0 && Ehad < 100.0)*(0.5961) +
640 (Ehad >= 100.0 && Ehad < 150.0)*(0.6368) +
641 (Ehad >= 150.0 && Ehad < 200.0)*(0.6814) +
642 (Ehad >= 200.0)*(0.6998)) +
643 (decayR < 2700 && decayZ < 6710) * ((0.0 == Ehad)*(0.8604) +
644 (Ehad > 0.0 && Ehad < 25.0)*(0.3335) +
645 (Ehad >= 25.0 && Ehad < 50.0)*(0.2457) +
646 (Ehad >= 50.0 && Ehad < 75.0)*(0.1831) +
647 (Ehad >= 75.0 && Ehad < 100.0)*(0.2100) +
648 (Ehad >= 100.0 && Ehad < 150.0)*(0.2443) +
649 (Ehad >= 150.0 && Ehad < 200.0)*(0.2532) +
650 (Ehad >= 200.0)*(0.2404))
[a5af1df]651 }
652}
653######################
654# ChargedHadronFilter
655######################
656
657module PdgCodeFilter ChargedHadronFilter {
658 set InputArray HCal/eflowTracks
659 set OutputArray chargedHadrons
660
661 add PdgCode {11}
662 add PdgCode {-11}
663 add PdgCode {13}
664 add PdgCode {-13}
665}
666
667
668###################################################
669# Tower Merger (in case not using e-flow algorithm)
670###################################################
671
672module Merger Calorimeter {
673# add InputArray InputArray
674 add InputArray ECal/ecalTowers
675 add InputArray HCal/hcalTowers
676 set OutputArray towers
677}
678
679
680
681####################
682# Energy flow merger
683####################
684
685module Merger EFlowMerger {
686# add InputArray InputArray
687 add InputArray HCal/eflowTracks
688 add InputArray ECal/eflowPhotons
689 add InputArray HCal/eflowNeutralHadrons
690 set OutputArray eflow
691}
692
693######################
694# EFlowFilter
695######################
696
697module PdgCodeFilter EFlowFilter {
698 set InputArray EFlowMerger/eflow
699 set OutputArray eflow
700
701 add PdgCode {11}
702 add PdgCode {-11}
703 add PdgCode {13}
704 add PdgCode {-13}
705}
706
707
708###################
709# Photon efficiency
710###################
711
712module Efficiency PhotonEfficiency {
713 set InputArray ECal/eflowPhotons
714 set OutputArray photons
715
716 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
717
718 # efficiency formula for photons
719 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
720 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
721 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
722 (abs(eta) > 2.5) * (0.00)}
723}
724
725##################
726# Photon isolation
727##################
728
729module Isolation PhotonIsolation {
730 set CandidateInputArray PhotonEfficiency/photons
731 set IsolationInputArray EFlowFilter/eflow
732
733 set OutputArray photons
734
735 set DeltaRMax 0.5
736
737 set PTMin 0.5
738
739 set PTRatioMax 0.12
740}
741
742
743#####################
744# Electron efficiency
745#####################
746
747module Efficiency ElectronEfficiency {
748 set InputArray ElectronFilter/electrons
749 set OutputArray electrons
750
751 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
752
753 # efficiency formula for electrons
754 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
755 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
756 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
757 (abs(eta) > 2.5) * (0.00)}
758}
759
760####################
761# Electron isolation
762####################
763
764module Isolation ElectronIsolation {
765 set CandidateInputArray ElectronEfficiency/electrons
766 set IsolationInputArray EFlowFilter/eflow
767
768 set OutputArray electrons
769
770 set DeltaRMax 0.5
771
772 set PTMin 0.5
773
774 set PTRatioMax 0.12
775}
776
777#################
778# Muon efficiency
779#################
780
781module Efficiency MuonEfficiency {
782 set InputArray MuonMomentumSmearing/muons
783 set OutputArray muons
784
785 # set EfficiencyFormula {efficiency as a function of eta and pt}
786
787 # efficiency formula for muons
788 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
789 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
790 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0) * (0.95) +
791 (abs(eta) > 2.4) * (0.00)}
792}
793
794################
795# Muon isolation
796################
797
798module Isolation MuonIsolation {
799 set CandidateInputArray MuonEfficiency/muons
800 set IsolationInputArray EFlowFilter/eflow
801
802 set OutputArray muons
803
804 set DeltaRMax 0.5
805
806 set PTMin 0.5
807
808 set PTRatioMax 0.25
809}
810
811###################
812# Missing ET merger
813###################
814
815module Merger MissingET {
816# add InputArray InputArray
817 add InputArray EFlowMerger/eflow
818 set MomentumOutputArray momentum
819}
820
821##################
822# Scalar HT merger
823##################
824
825module Merger ScalarHT {
826# add InputArray InputArray
827 add InputArray UniqueObjectFinder/jets
828 add InputArray UniqueObjectFinder/electrons
829 add InputArray UniqueObjectFinder/photons
830 add InputArray UniqueObjectFinder/muons
831 set EnergyOutputArray energy
832}
833
834
835#####################
836# Neutrino Filter
837#####################
838
839module PdgCodeFilter NeutrinoFilter {
840
841 set InputArray Delphes/stableParticles
842 set OutputArray filteredParticles
843
844 set PTMin 0.0
845
846 add PdgCode {12}
847 add PdgCode {14}
848 add PdgCode {16}
849 add PdgCode {-12}
850 add PdgCode {-14}
851 add PdgCode {-16}
852
853}
854
855
856#####################
857# MC truth jet finder
858#####################
859
860module FastJetFinder GenJetFinder {
861 set InputArray NeutrinoFilter/filteredParticles
862
863 set OutputArray jets
864
865 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
866 set JetAlgorithm 6
867 set ParameterR 0.5
868
869 set JetPTMin 20.0
870}
871
872#########################
873# Gen Missing ET merger
874########################
875
876module Merger GenMissingET {
877# add InputArray InputArray
878 add InputArray NeutrinoFilter/filteredParticles
879
880 set MomentumOutputArray momentum
881}
882
883
884
885############
886# Jet finder
887############
888
889module FastJetFinder FastJetFinder {
890# set InputArray Calorimeter/towers
891 set InputArray EFlowMerger/eflow
892
893 set OutputArray jets
894
895 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
896 set JetAlgorithm 6
897 set ParameterR 0.5
898
899 set JetPTMin 20.0
900}
901
902##################
903# Fat Jet finder
904##################
905
906module FastJetFinder FatJetFinder {
907 set InputArray EFlowMerger/eflow
908
909 set OutputArray jets
910
911 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
912 set JetAlgorithm 6
913 set ParameterR 0.8
914
915 set ComputeNsubjettiness 1
916 set Beta 1.0
917 set AxisMode 4
918
919 set ComputeTrimming 1
920 set RTrim 0.2
921 set PtFracTrim 0.05
922
923 set ComputePruning 1
924 set ZcutPrun 0.1
925 set RcutPrun 0.5
926 set RPrun 0.8
927
928 set ComputeSoftDrop 1
929 set BetaSoftDrop 0.0
930 set SymmetryCutSoftDrop 0.1
931 set R0SoftDrop 0.8
932
933 set JetPTMin 200.0
934}
935
936
937
938
939##################
940# Jet Energy Scale
941##################
942
943module EnergyScale JetEnergyScale {
944 set InputArray FastJetFinder/jets
945 set OutputArray jets
946
947 # scale formula for jets
948 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
949}
950
951########################
952# Jet Flavor Association
953########################
954
955module JetFlavorAssociation JetFlavorAssociation {
956
957 set PartonInputArray Delphes/partons
958 set ParticleInputArray Delphes/allParticles
959 set ParticleLHEFInputArray Delphes/allParticlesLHEF
960 set JetInputArray JetEnergyScale/jets
961
962 set DeltaR 0.5
963 set PartonPTMin 1.0
964 set PartonEtaMax 2.5
965
966}
967
968###########
969# b-tagging
970###########
971
972module BTagging BTagging {
973 set JetInputArray JetEnergyScale/jets
974
975 set BitNumber 0
976
977 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
978 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
979 # gluon's PDG code has the lowest priority
980
981 # based on arXiv:1211.4462
982
983 # default efficiency formula (misidentification rate)
984 add EfficiencyFormula {0} {0.01+0.000038*pt}
985
986 # efficiency formula for c-jets (misidentification rate)
987 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
988
989 # efficiency formula for b-jets
990 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
991}
992
993#############
994# tau-tagging
995#############
996
997module TauTagging TauTagging {
998 set ParticleInputArray Delphes/allParticles
999 set PartonInputArray Delphes/partons
1000 set JetInputArray JetEnergyScale/jets
1001
1002 set DeltaR 0.5
1003
1004 set TauPTMin 1.0
1005
1006 set TauEtaMax 2.5
1007
1008 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1009
1010 # default efficiency formula (misidentification rate)
1011 add EfficiencyFormula {0} {0.01}
1012 # efficiency formula for tau-jets
1013 add EfficiencyFormula {15} {0.6}
1014}
1015
1016#####################################################
1017# Find uniquely identified photons/electrons/tau/jets
1018#####################################################
1019
1020module UniqueObjectFinder UniqueObjectFinder {
1021# earlier arrays take precedence over later ones
1022# add InputArray InputArray OutputArray
1023 add InputArray PhotonIsolation/photons photons
1024 add InputArray ElectronIsolation/electrons electrons
1025 add InputArray MuonIsolation/muons muons
1026 add InputArray JetEnergyScale/jets jets
1027}
1028
1029##################
1030# ROOT tree writer
1031##################
1032
1033# tracks, towers and eflow objects are not stored by default in the output.
1034# if needed (for jet constituent or other studies), uncomment the relevant
1035# "add Branch ..." lines.
1036
1037module TreeWriter TreeWriter {
1038# add Branch InputArray BranchName BranchClass
[c61b5ce]1039# add Branch Delphes/allParticles Particle GenParticle
[a5af1df]1040
1041 add Branch UniqueObjectFinder/jets Jet Jet
1042 add Branch UniqueObjectFinder/electrons Electron Electron
1043 add Branch UniqueObjectFinder/muons Muon Muon
1044
1045 add Branch MissingET/momentum MissingET MissingET
1046 add Branch llpFilter/LLP llp CscCluster
1047 add Branch CSCFilter/LLP Cscllp CscCluster
[29b722a]1048 add Branch ClusterEfficiency/cluster CscCluster CscCluster
[a5af1df]1049}
Note: See TracBrowser for help on using the repository browser.