Fork me on GitHub

source: git/cards/CMS_PhaseII_200PU.tcl@ 61c0925

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 61c0925 was 61c0925, checked in by Michele Selvaggi <michele.selvaggi@…>, 9 years ago

set HGCAL/HCAL granularity same as in HGCAL/ECAL
  • Property mode set to 100644
File size: 75.0 KB
Line 
1#
2# Phase II - Pile-Up
3#
4# Main authors: Michele Selvaggi (UCL)
5#
6# Released on:
7#
8# Version: v01
9#
10#
11#######################################
12# Order of execution of various modules
13#######################################
14
15set ExecutionPath {
16
17 PileUpMerger
18
19 ParticlePropagator
20
21 ChargedHadronTrackingEfficiency
22 ElectronTrackingEfficiency
23 MuonTrackingEfficiency
24
25 ChargedHadronMomentumSmearing
26 ElectronMomentumSmearing
27 MuonMomentumSmearing
28
29 TrackMerger
30
31 ECal
32 HCal
33
34 ElectronFilter
35 TrackPileUpSubtractor
36
37 TowerMerger
38 NeutralEFlowMerger
39 EFlowMergerAllTracks
40 EFlowMerger
41
42 RunPUPPI
43
44 PhotonEfficiency
45 PhotonIsolation
46
47 ElectronEfficiency
48 ElectronIsolation
49
50 MuonEfficiency
51 MuonIsolation
52
53 NeutrinoFilter
54
55 MissingET
56 GenMissingET
57 GenPileUpMissingET
58
59 GenJetFinder
60 FastJetFinder
61
62 ScalarHT
63
64 JetEnergyScale
65
66 JetFlavorAssociation
67
68 BTaggingLoose
69 BTaggingMedium
70 BTaggingTight
71
72 TauTagging
73
74 GenParticleFilter
75
76 TreeWriter
77}
78
79
80
81###############
82# PileUp Merger
83###############
84
85module PileUpMerger PileUpMerger {
86 set InputArray Delphes/stableParticles
87
88 set ParticleOutputArray stableParticles
89 set VertexOutputArray vertices
90
91 # pre-generated minbias input file
92 set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup
93
94 # average expected pile up
95 set MeanPileUp 200
96
97 # maximum spread in the beam direction in m
98 set ZVertexSpread 0.15
99
100 # maximum spread in time in s
101 set TVertexSpread 1.5E-09
102
103 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
104
105 #set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
106
107
108 set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) +
109 (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) +
110 (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) +
111 (abs(t) > 1.0e-09) * (abs(z) > 0.15) * (0.00)}
112
113}
114
115
116#################################
117# Propagate particles in cylinder
118#################################
119
120module ParticlePropagator ParticlePropagator {
121 set InputArray PileUpMerger/stableParticles
122
123 set OutputArray stableParticles
124 set ChargedHadronOutputArray chargedHadrons
125 set ElectronOutputArray electrons
126 set MuonOutputArray muons
127
128 # radius of the magnetic field coverage, in m
129 set Radius 1.29
130 # half-length of the magnetic field coverage, in m
131 set HalfLength 3.0
132
133 # magnetic field
134 set Bz 3.8
135}
136
137
138####################################
139# Charged hadron tracking efficiency
140####################################
141
142module Efficiency ChargedHadronTrackingEfficiency {
143 ## particles after propagation
144 set InputArray ParticlePropagator/chargedHadrons
145 set OutputArray chargedHadrons
146 # tracking efficiency formula for charged hadrons
147 set EfficiencyFormula {
148 (pt <= 0.2) * (0.00) + \
149 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
150 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
151 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
152 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.87) + \
153 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
154 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.82) + \
155 (abs(eta) > 4.0) * (0.00)
156 }
157}
158
159
160#####################################
161# Electron tracking efficiency - ID
162####################################
163
164module Efficiency ElectronTrackingEfficiency {
165 set InputArray ParticlePropagator/electrons
166 set OutputArray electrons
167 # tracking efficiency formula for electrons
168 set EfficiencyFormula {
169 (pt <= 0.2) * (0.00) + \
170 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
171 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
172 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
173 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 10.0) * (0.82+pt*0.01) + \
174 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 10.0) * (0.90) + \
175 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
176 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.8+pt*0.01) + \
177 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.85) + \
178 (abs(eta) > 4.0) * (0.00)
179
180 }
181}
182
183##########################
184# Muon tracking efficiency
185##########################
186
187module Efficiency MuonTrackingEfficiency {
188 set InputArray ParticlePropagator/muons
189 set OutputArray muons
190 # tracking efficiency formula for muons
191 set EfficiencyFormula {
192 (pt <= 0.2) * (0.00) + \
193 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.998) + \
194 (abs(eta) <= 1.2) * (pt > 1.0) * (0.998) + \
195 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.99) + \
196 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.99) + \
197 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
198 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
199 (abs(eta) > 4.0) * (0.00)
200
201 }
202}
203
204########################################
205# Momentum resolution for charged tracks
206########################################
207
208module MomentumSmearing ChargedHadronMomentumSmearing {
209 ## hadrons after having applied the tracking efficiency
210 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
211 set OutputArray chargedHadrons
212 # resolution formula for charged hadrons
213 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
214 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.013) + \
215 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.02) + \
216 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
217 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
218 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.015) + \
219 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
220 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05) + \
221 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 0.1 && pt <= 1.0) * (0.017) + \
222 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 1.0 && pt <= 10.0) * (0.03) + \
223 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 100.0) * (0.05) + \
224 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 100.0) * (0.30) + \
225 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
226 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 1.0 && pt <= 10.0) * (0.04) + \
227 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 10.0 && pt <= 100.0) * (0.07) + \
228 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 100.0) * (0.30) + \
229 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.025) + \
230 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.05) + \
231 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 100.0) * (0.20) + \
232 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 100.0) * (0.80)
233 }
234}
235
236########################################
237# Momentum resolution for electrons
238########################################
239
240module MomentumSmearing ElectronMomentumSmearing {
241 ## hadrons after having applied the tracking efficiency
242 set InputArray ElectronTrackingEfficiency/electrons
243 set OutputArray electrons
244 # resolution formula for electrons
245 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
246 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.013) + \
247 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.02) + \
248 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
249 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
250 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.015) + \
251 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
252 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05) + \
253 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 0.1 && pt <= 1.0) * (0.017) + \
254 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 1.0 && pt <= 10.0) * (0.03) + \
255 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 100.0) * (0.05) + \
256 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 100.0) * (0.30) + \
257 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
258 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 1.0 && pt <= 10.0) * (0.04) + \
259 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 10.0 && pt <= 100.0) * (0.07) + \
260 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 100.0) * (0.30) + \
261 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.025) + \
262 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.05) + \
263 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 100.0) * (0.20) + \
264 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 100.0) * (0.80)
265 }
266}
267
268
269###############################
270# Momentum resolution for muons
271###############################
272
273module MomentumSmearing MuonMomentumSmearing {
274 set InputArray MuonTrackingEfficiency/muons
275 set OutputArray muons
276 # resolution formula for muons
277 set ResolutionFormula {
278 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.012) + \
279 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
280 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.012) + \
281 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.030) + \
282 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.011) + \
283 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.011) + \
284 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.022) + \
285 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.030) + \
286 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 0.1 && pt <= 1.0) * (0.017) + \
287 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 1.0 && pt <= 10.0) * (0.03) + \
288 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 100.0) * (0.05) + \
289 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 100.0) * (0.30) + \
290 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
291 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 1.0 && pt <= 10.0) * (0.04) + \
292 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 10.0 && pt <= 100.0) * (0.07) + \
293 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 100.0) * (0.30) + \
294 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.025) + \
295 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.05) + \
296 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 100.0) * (0.20) + \
297 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 100.0) * (0.80)
298
299 }
300}
301
302
303##############
304# Track merger
305##############
306
307module Merger TrackMerger {
308# add InputArray InputArray
309 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
310 add InputArray ElectronMomentumSmearing/electrons
311 add InputArray MuonMomentumSmearing/muons
312 set OutputArray tracks
313}
314
315#############
316# ECAL
317#############
318
319module SimpleCalorimeter ECal {
320 set ParticleInputArray ParticlePropagator/stableParticles
321 set TrackInputArray TrackMerger/tracks
322
323 set TowerOutputArray ecalTowers
324 set EFlowTrackOutputArray eflowTracks
325 set EFlowTowerOutputArray eflowPhotons
326
327 set IsEcal true
328
329 set EnergyMin 0.5
330 set EnergySignificanceMin 1.0
331
332 set SmearTowerCenter true
333
334 set pi [expr {acos(-1)}]
335
336 # lists of the edges of each tower in eta and phi
337 # each list starts with the lower edge of the first tower
338 # the list ends with the higher edged of the last tower
339
340 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
341
342 set PhiBins {}
343 for {set i -180} {$i <= 180} {incr i} {
344 add PhiBins [expr {$i * $pi/180.0}]
345 }
346
347 # 0.02 unit in eta up to eta = 1.5 (barrel)
348 for {set i -85} {$i <= 86} {incr i} {
349 set eta [expr {$i * 0.0174}]
350 add EtaPhiBins $eta $PhiBins
351 }
352
353 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
354
355 set PhiBins {}
356 for {set i -180} {$i <= 180} {incr i} {
357 add PhiBins [expr {$i * $pi/180.0}]
358 }
359
360 # 0.02 unit in eta up to eta = 3
361 for {set i 1} {$i <= 84} {incr i} {
362 set eta [expr { -2.958 + $i * 0.0174}]
363 add EtaPhiBins $eta $PhiBins
364 }
365
366 for {set i 1} {$i <= 84} {incr i} {
367 set eta [expr { 1.4964 + $i * 0.0174}]
368 add EtaPhiBins $eta $PhiBins
369 }
370
371 # take present CMS granularity for HF
372
373 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
374 set PhiBins {}
375 for {set i -18} {$i <= 18} {incr i} {
376 add PhiBins [expr {$i * $pi/18.0}]
377 }
378
379 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} {
380 add EtaPhiBins $eta $PhiBins
381 }
382
383
384 add EnergyFraction {0} {0.0}
385 # energy fractions for e, gamma and pi0
386 add EnergyFraction {11} {1.0}
387 add EnergyFraction {22} {1.0}
388 add EnergyFraction {111} {1.0}
389 # energy fractions for muon, neutrinos and neutralinos
390 add EnergyFraction {12} {0.0}
391 add EnergyFraction {13} {0.0}
392 add EnergyFraction {14} {0.0}
393 add EnergyFraction {16} {0.0}
394 add EnergyFraction {1000022} {0.0}
395 add EnergyFraction {1000023} {0.0}
396 add EnergyFraction {1000025} {0.0}
397 add EnergyFraction {1000035} {0.0}
398 add EnergyFraction {1000045} {0.0}
399 # energy fractions for K0short and Lambda
400 add EnergyFraction {310} {0.3}
401 add EnergyFraction {3122} {0.3}
402
403 # set ResolutionFormula {resolution formula as a function of eta and energy}
404
405 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016.
406 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117
407
408 set ResolutionFormula { (abs(eta) <= 1.50) * sqrt(energy^2*0.005^2 + energy*0.027^2 + 0.15^2) + \
409 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
410 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
411 (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
412 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
413
414}
415
416#############
417# HCAL
418#############
419
420module SimpleCalorimeter HCal {
421 set ParticleInputArray ParticlePropagator/stableParticles
422 set TrackInputArray ECal/eflowTracks
423
424 set TowerOutputArray hcalTowers
425 set EFlowTrackOutputArray eflowTracks
426 set EFlowTowerOutputArray eflowNeutralHadrons
427
428 set IsEcal false
429
430 set EnergyMin 1.0
431 set EnergySignificanceMin 1.0
432
433 set SmearTowerCenter true
434
435 set pi [expr {acos(-1)}]
436
437 # lists of the edges of each tower in eta and phi
438 # each list starts with the lower edge of the first tower
439 # the list ends with the higher edged of the last tower
440
441 # assume 0.087 x 0.087 resolution in eta,phi in the barrel |eta| < 1.5
442
443 set PhiBins {}
444 for {set i -36} {$i <= 36} {incr i} {
445 add PhiBins [expr {$i * $pi/36.0}]
446 }
447 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.65} {
448 add EtaPhiBins $eta $PhiBins
449 }
450
451 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- HCAL)
452
453 set PhiBins {}
454 for {set i -180} {$i <= 180} {incr i} {
455 add PhiBins [expr {$i * $pi/180.0}]
456 }
457
458 # 0.02 unit in eta up to eta = 3
459 for {set i 1} {$i <= 84} {incr i} {
460 set eta [expr { -2.958 + $i * 0.0174}]
461 add EtaPhiBins $eta $PhiBins
462 }
463
464 for {set i 1} {$i <= 84} {incr i} {
465 set eta [expr { 1.4964 + $i * 0.0174}]
466 add EtaPhiBins $eta $PhiBins
467 }
468
469 # take present CMS granularity for HF
470
471 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
472 set PhiBins {}
473 for {set i -18} {$i <= 18} {incr i} {
474 add PhiBins [expr {$i * $pi/18.0}]
475 }
476
477 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} {
478 add EtaPhiBins $eta $PhiBins
479 }
480
481
482 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
483 add EnergyFraction {0} {1.0}
484 # energy fractions for e, gamma and pi0
485 add EnergyFraction {11} {0.0}
486 add EnergyFraction {22} {0.0}
487 add EnergyFraction {111} {0.0}
488 # energy fractions for muon, neutrinos and neutralinos
489 add EnergyFraction {12} {0.0}
490 add EnergyFraction {13} {0.0}
491 add EnergyFraction {14} {0.0}
492 add EnergyFraction {16} {0.0}
493 add EnergyFraction {1000022} {0.0}
494 add EnergyFraction {1000023} {0.0}
495 add EnergyFraction {1000025} {0.0}
496 add EnergyFraction {1000035} {0.0}
497 add EnergyFraction {1000045} {0.0}
498 # energy fractions for K0short and Lambda
499 add EnergyFraction {310} {0.7}
500 add EnergyFraction {3122} {0.7}
501
502# set ResolutionFormula {resolution formula as a function of eta and energy}
503 set ResolutionFormula { (abs(eta) <= 1.5) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
504 (abs(eta) > 1.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
505 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.11^2 + energy*2.80^2)}
506
507}
508
509#################
510# Electron filter
511#################
512
513module PdgCodeFilter ElectronFilter {
514 set InputArray HCal/eflowTracks
515 set OutputArray electrons
516 set Invert true
517 add PdgCode {11}
518 add PdgCode {-11}
519}
520
521##########################
522# Track pile-up subtractor
523##########################
524
525module TrackPileUpSubtractor TrackPileUpSubtractor {
526# add InputArray InputArray OutputArray
527 add InputArray HCal/eflowTracks eflowTracks
528 add InputArray ElectronFilter/electrons electrons
529 add InputArray MuonMomentumSmearing/muons muons
530
531 set VertexInputArray PileUpMerger/vertices
532 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
533 # Z vertex resolution in m
534 set ZVertexResolution 0.0001
535}
536
537
538###################################################
539# Tower Merger (in case not using e-flow algorithm)
540###################################################
541
542module Merger TowerMerger {
543# add InputArray InputArray
544 add InputArray ECal/ecalTowers
545 add InputArray HCal/hcalTowers
546 set OutputArray towers
547}
548
549
550####################
551# Neutral eflow erger
552####################
553
554module Merger NeutralEFlowMerger {
555# add InputArray InputArray
556 add InputArray ECal/eflowPhotons
557 add InputArray HCal/eflowNeutralHadrons
558 set OutputArray eflowTowers
559}
560
561
562####################
563# Energy flow merger
564####################
565
566module Merger EFlowMerger {
567# add InputArray InputArray
568 add InputArray HCal/eflowTracks
569 add InputArray ECal/eflowPhotons
570 add InputArray HCal/eflowNeutralHadrons
571 set OutputArray eflow
572}
573
574##################################
575# Energy flow merger (all tracks)
576##################################
577
578module Merger EFlowMergerAllTracks {
579# add InputArray InputArray
580 add InputArray TrackMerger/tracks
581 add InputArray ECal/eflowPhotons
582 add InputArray HCal/eflowNeutralHadrons
583 set OutputArray eflow
584}
585
586#########################################
587### Run the puppi code (to be tuned) ###
588#########################################
589
590module RunPUPPI RunPUPPI {
591 ## input information
592 set TrackInputArray TrackMerger/tracks
593 set NeutralInputArray NeutralEFlowMerger/eflowTowers
594 set PVInputArray PileUpMerger/vertices
595 set MinPuppiWeight 0.05
596 set UseExp false
597
598 ## define puppi algorithm parameters (more than one for the same eta region is possible)
599 add EtaMinBin 0. 2.5 2.5 3.0 3.0
600 add EtaMaxBin 2.5 3.0 3.0 10.0 10.0
601 add PtMinBin 0. 0.5 0.5 0.5 0.5
602 add ConeSizeBin 0.25 0.25 0.25 0.25 0.25
603 add RMSPtMinBin 0.1 0.5 0.5 0.5 0.5
604 add RMSScaleFactorBin 1.0 1.0 1.0 1.0 1.0
605 add NeutralMinEBin 0.2 1.0 1.0 1.5 1.5
606 add NeutralPtSlope 0.02 0.02 0.02 0.02 0.02
607 add ApplyCHS true true true true true
608 add UseCharged true false false false false
609 add ApplyLowPUCorr true true true true true
610 add MetricId 5 5 1 5 1
611
612 ## output name
613 set OutputArray PuppiParticles
614 set OutputArrayTracks puppiTracks
615 set OutputArrayNeutrals puppiNeutrals
616}
617
618
619
620###################
621# Missing ET merger
622###################
623
624module Merger MissingET {
625# add InputArray InputArray
626# add InputArray RunPUPPI/PuppiParticles
627 add InputArray EFlowMerger/eflow
628 set MomentumOutputArray momentum
629}
630
631###################
632# Ger PileUp Missing ET
633###################
634
635module Merger GenPileUpMissingET {
636# add InputArray InputArray
637# add InputArray RunPUPPI/PuppiParticles
638 add InputArray ParticlePropagator/stableParticles
639 set MomentumOutputArray momentum
640}
641
642##################
643# Scalar HT merger
644##################
645
646module Merger ScalarHT {
647# add InputArray InputArray
648 add InputArray RunPUPPI/PuppiParticles
649 set EnergyOutputArray energy
650}
651
652#################
653# Neutrino Filter
654#################
655
656module PdgCodeFilter NeutrinoFilter {
657
658 set InputArray Delphes/stableParticles
659 set OutputArray filteredParticles
660
661 set PTMin 0.0
662
663 add PdgCode {12}
664 add PdgCode {14}
665 add PdgCode {16}
666 add PdgCode {-12}
667 add PdgCode {-14}
668 add PdgCode {-16}
669
670}
671
672
673#####################
674# MC truth jet finder
675#####################
676
677module FastJetFinder GenJetFinder {
678 set InputArray NeutrinoFilter/filteredParticles
679
680 set OutputArray jets
681
682 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
683 set JetAlgorithm 6
684 set ParameterR 0.4
685
686 set JetPTMin 15.0
687}
688
689#########################
690# Gen Missing ET merger
691########################
692
693module Merger GenMissingET {
694
695# add InputArray InputArray
696 add InputArray NeutrinoFilter/filteredParticles
697 set MomentumOutputArray momentum
698}
699
700
701
702############
703# Jet finder
704############
705
706module FastJetFinder FastJetFinder {
707# set InputArray TowerMerger/towers
708 set InputArray RunPUPPI/PuppiParticles
709
710 set OutputArray jets
711
712 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
713 set JetAlgorithm 6
714 set ParameterR 0.4
715
716 set JetPTMin 15.0
717}
718
719##################
720# Jet Energy Scale
721##################
722
723module EnergyScale JetEnergyScale {
724 set InputArray FastJetFinder/jets
725 set OutputArray jets
726
727 # scale formula for jets
728 set ScaleFormula {1.00}
729}
730
731
732#####################
733# Photon efficiency #
734#####################
735
736module Efficiency PhotonEfficiency {
737
738 ## input particles
739 set InputArray ECal/eflowPhotons
740 ## output particles
741 set OutputArray photons
742 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
743 # efficiency formula for photons
744 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
745 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
746 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
747 (abs(eta) > 4.0) * (0.00)}
748
749}
750
751####################
752# Photon isolation #
753####################
754
755module Isolation PhotonIsolation {
756
757 # particle for which calculate the isolation
758 set CandidateInputArray PhotonEfficiency/photons
759
760 # isolation collection
761 set IsolationInputArray RunPUPPI/PuppiParticles
762
763 # output array
764 set OutputArray photons
765
766 # isolation cone
767 set DeltaRMax 0.3
768
769 # minimum pT
770 set PTMin 1.0
771
772 # iso ratio to cut
773 set PTRatioMax 9999.
774}
775
776
777#######################
778# Electron efficiency #
779#######################
780
781module Efficiency ElectronEfficiency {
782
783 set InputArray ElectronFilter/electrons
784 set OutputArray electrons
785
786 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
787 # efficiency formula for electrons
788 set EfficiencyFormula {
789 (pt <= 4.0) * (0.00) + \
790 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
791 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
792 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
793 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
794 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
795 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
796 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
797 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
798 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
799 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
800 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
801 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
802 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
803 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
804 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
805 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
806 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
807 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
808 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
809 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
810 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
811 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
812 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
813 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
814 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
815 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
816 (abs(eta) > 4.0) * (0.00)
817
818 }
819}
820
821######################
822# Electron isolation #
823######################
824
825module Isolation ElectronIsolation {
826
827 set CandidateInputArray ElectronEfficiency/electrons
828
829 # isolation collection
830 set IsolationInputArray RunPUPPI/PuppiParticles
831
832 set OutputArray electrons
833
834 set DeltaRMax 0.3
835 set PTMin 1.0
836 set PTRatioMax 9999.
837
838}
839
840
841
842###################
843# Muon efficiency #
844###################
845
846module Efficiency MuonEfficiency {
847
848 set InputArray MuonMomentumSmearing/muons
849
850 set OutputArray muons
851 # set EfficiencyFormula {efficiency as a function of eta and pt}
852 # efficiency formula for muons
853 set EfficiencyFormula { (pt <= 2.0) * (0.00) + \
854 (abs(eta) <= 4.00) * (pt > 2.0 && pt <= 3.0) * (0.51) + \
855 (abs(eta) <= 4.00) * (pt > 3.0 && pt <= 4.0) * (0.85) + \
856 (abs(eta) <= 4.00) * (pt > 4.0 && pt <= 11.0) * (0.93) + \
857 (abs(eta) <= 4.00) * (pt > 11. && pt <= 50.) * (0.96) + \
858 (abs(eta) <= 4.00) * (pt > 50. && pt <= 70.) * (0.98) + \
859 (abs(eta) <= 4.00) * (pt > 70.0 ) * (1.00) + \
860 (abs(eta) > 4.00) * (0.00)
861 }
862
863}
864
865
866
867##################
868# Muon isolation #
869##################
870
871module Isolation MuonIsolation {
872 set CandidateInputArray MuonEfficiency/muons
873
874 # isolation collection
875 set IsolationInputArray RunPUPPI/PuppiParticles
876
877 set OutputArray muons
878
879 set DeltaRMax 0.3
880 set PTMin 1.0
881 set PTRatioMax 9999.
882
883}
884
885
886
887########################
888# Jet Flavor Association
889########################
890
891module JetFlavorAssociation JetFlavorAssociation {
892
893 set PartonInputArray Delphes/partons
894 set ParticleInputArray Delphes/allParticles
895 set ParticleLHEFInputArray Delphes/allParticlesLHEF
896 set JetInputArray JetEnergyScale/jets
897
898 set DeltaR 0.5
899 set PartonPTMin 10.0
900 set PartonEtaMax 4.0
901
902}
903
904
905#############
906# b-tagging #
907#############
908module BTagging BTaggingLoose {
909
910 set JetInputArray JetEnergyScale/jets
911
912 set BitNumber 0
913
914 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
915 (pt > 20.0) * (abs(eta) <= 3.4) * (0.1) + \
916 (abs(eta) > 3.4) * (0.000)
917 }
918
919
920 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
921 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.341) + \
922 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.389) + \
923 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.416) + \
924 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.432) + \
925 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.447) + \
926 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.451) + \
927 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.451) + \
928 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.447) + \
929 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.447) + \
930 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.444) + \
931 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.430) + \
932 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.424) + \
933 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.419) + \
934 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.407) + \
935 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.384) + \
936 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.364) + \
937 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.353) + \
938 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.303) + \
939 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.273) + \
940 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.249) + \
941 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.233) + \
942 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.209) + \
943 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.194) + \
944 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.227) + \
945 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
946
947 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
948 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.236) + \
949 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.276) + \
950 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.302) + \
951 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.330) + \
952 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.338) + \
953 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.355) + \
954 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.353) + \
955 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.358) + \
956 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.350) + \
957 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.327) + \
958 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.313) + \
959 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.307) + \
960 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.322) + \
961 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.299) + \
962 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.274) + \
963 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.240) + \
964 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.203) + \
965 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.218) + \
966 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.203) + \
967 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.178) + \
968 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
969
970 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
971 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.179) + \
972 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.198) + \
973 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.219) + \
974 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.225) + \
975 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.232) + \
976 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.251) + \
977 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.267) + \
978 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.252) + \
979 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.243) + \
980 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.244) + \
981 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.242) + \
982 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.246) + \
983 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.217) + \
984 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.223)+ \
985 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.194) + \
986 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.187) + \
987 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
988
989 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
990 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.09) + \
991 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.10) + \
992 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.11) + \
993 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.11) + \
994 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.12) + \
995 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.13) + \
996 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.13) + \
997 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.13) + \
998 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.12) + \
999 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.12) + \
1000 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.12) + \
1001 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.12) + \
1002 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.11) + \
1003 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.11) + \
1004 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.10) + \
1005 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.09) + \
1006 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1007 (abs(eta) > 3.4) * (0.000)
1008 }
1009
1010 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1011 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.712) + \
1012 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.788) + \
1013 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.837) + \
1014 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.856) + \
1015 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.873) + \
1016 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.886) + \
1017 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.883) + \
1018 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.882) + \
1019 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.883) + \
1020 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.883) + \
1021 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.876) + \
1022 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.872) + \
1023 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.867) + \
1024 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.853) + \
1025 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.831) + \
1026 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.819) + \
1027 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.791) + \
1028 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.740) + \
1029 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.699) + \
1030 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.631) + \
1031 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.600) + \
1032 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.567) + \
1033 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.478) + \
1034 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.470) + \
1035 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1036
1037 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1038 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.543) + \
1039 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.627) + \
1040 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.701) + \
1041 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.741) + \
1042 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.754) + \
1043 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.796) + \
1044 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.778) + \
1045 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.798) + \
1046 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.788) + \
1047 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.775) + \
1048 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.765) + \
1049 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.752) + \
1050 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.740) + \
1051 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.716) + \
1052 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.671) + \
1053 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.637) + \
1054 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.577) + \
1055 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.572) + \
1056 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.514) + \
1057 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.516) + \
1058 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1059
1060 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1061 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.407) + \
1062 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.479) + \
1063 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.529) + \
1064 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.578) + \
1065 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.612) + \
1066 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.673) + \
1067 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.656) + \
1068 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.658) + \
1069 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.647) + \
1070 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.629) + \
1071 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.630) + \
1072 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.615) + \
1073 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.604) + \
1074 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.532)+ \
1075 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.532) + \
1076 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.503) + \
1077 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1078
1079 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1080 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.20) + \
1081 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.235) + \
1082 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.26) + \
1083 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.28) + \
1084 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.31) + \
1085 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.33) + \
1086 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.32) + \
1087 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.32) + \
1088 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.33) + \
1089 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.31) + \
1090 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.32) + \
1091 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.31) + \
1092 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.30) + \
1093 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.26) + \
1094 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.27) + \
1095 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.25) + \
1096 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1097 (abs(eta) > 3.4) * (0.000)
1098 }
1099}
1100
1101
1102#############
1103# b-tagging #
1104#############
1105module BTagging BTaggingMedium {
1106
1107 set JetInputArray JetEnergyScale/jets
1108
1109 set BitNumber 1
1110
1111 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
1112 (pt > 20.0) * (abs(eta) <= 3.4) * (0.01) + \
1113 (abs(eta) > 3.4) * (0.000)
1114 }
1115
1116 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1117 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.144) + \
1118 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.173) + \
1119 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.188) + \
1120 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.200) + \
1121 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.209) + \
1122 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.210) + \
1123 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.206) + \
1124 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.205) + \
1125 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.198) + \
1126 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.187) + \
1127 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.175) + \
1128 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.167) + \
1129 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.157) + \
1130 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.138) + \
1131 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.118) + \
1132 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.094) + \
1133 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.083) + \
1134 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.062) + \
1135 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.049) + \
1136 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.044) + \
1137 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.036) + \
1138 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.035) + \
1139 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.026) + \
1140 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.041) + \
1141 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1142
1143 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1144 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.070) + \
1145 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.089) + \
1146 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.104) + \
1147 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.121) + \
1148 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.121) + \
1149 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.127) + \
1150 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.127) + \
1151 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.117) + \
1152 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.117) + \
1153 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.096) + \
1154 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.099) + \
1155 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.086) + \
1156 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.075) + \
1157 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.066) + \
1158 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.056) + \
1159 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.048) + \
1160 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.031) + \
1161 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.042) + \
1162 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.030) + \
1163 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.030) + \
1164 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1165
1166 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1167 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.041) + \
1168 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.047) + \
1169 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.059) + \
1170 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.056) + \
1171 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.068) + \
1172 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.073) + \
1173 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.080) + \
1174 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.063) + \
1175 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.065) + \
1176 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.057) + \
1177 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.049) + \
1178 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.067) + \
1179 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.065) + \
1180 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.044)+ \
1181 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.031) + \
1182 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.026) + \
1183 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1184
1185 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1186 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.02) + \
1187 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.02) + \
1188 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.03) + \
1189 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.03) + \
1190 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.03) + \
1191 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.035) + \
1192 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.04) + \
1193 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.03) + \
1194 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.03) + \
1195 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.025) + \
1196 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.025) + \
1197 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.035) + \
1198 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.032) + \
1199 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.022) + \
1200 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.015) + \
1201 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.013) + \
1202 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1203 (abs(eta) > 3.4) * (0.000)
1204 }
1205
1206
1207 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1208 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.528) + \
1209 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.611) + \
1210 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.666) + \
1211 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.701) + \
1212 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.734) + \
1213 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.741) + \
1214 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.746) + \
1215 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.750) + \
1216 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.742) + \
1217 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.729) + \
1218 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.711) + \
1219 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.694) + \
1220 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.684) + \
1221 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.642) + \
1222 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.602) + \
1223 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.554) + \
1224 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.488) + \
1225 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.417) + \
1226 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.343) + \
1227 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.284) + \
1228 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.250) + \
1229 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.221) + \
1230 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.164) + \
1231 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.173) + \
1232 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1233
1234 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1235 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.322) + \
1236 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.414) + \
1237 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.488) + \
1238 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.521) + \
1239 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.554) + \
1240 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.590) + \
1241 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.571) + \
1242 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.583) + \
1243 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.575) + \
1244 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.557) + \
1245 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.536) + \
1246 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.516) + \
1247 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.483) + \
1248 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.432) + \
1249 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.366) + \
1250 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.323) + \
1251 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.241) + \
1252 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.252) + \
1253 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.237) + \
1254 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.243) + \
1255 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1256
1257 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1258 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.200) + \
1259 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.252) + \
1260 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.310) + \
1261 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.353) + \
1262 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.377) + \
1263 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.425) + \
1264 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.406) + \
1265 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.414) + \
1266 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.412) + \
1267 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.366) + \
1268 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.360) + \
1269 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.347) + \
1270 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.310) + \
1271 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.307)+ \
1272 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.224) + \
1273 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.194) + \
1274 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1275
1276 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1277 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.100) + \
1278 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.12) + \
1279 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.15) + \
1280 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.17) + \
1281 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.19) + \
1282 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.21) + \
1283 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.20) + \
1284 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.21) + \
1285 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.20) + \
1286 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.18) + \
1287 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.18) + \
1288 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.17) + \
1289 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.15) + \
1290 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.15) + \
1291 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.11) + \
1292 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.10) + \
1293 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1294 (abs(eta) > 3.4) * (0.000)
1295 }
1296}
1297
1298
1299
1300#############
1301# b-tagging #
1302#############
1303module BTagging BTaggingTight {
1304
1305 set JetInputArray JetEnergyScale/jets
1306
1307 set BitNumber 2
1308
1309 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
1310 (pt > 20.0) * (abs(eta) <= 3.4) * (0.001) + \
1311 (abs(eta) > 3.4) * (0.000)
1312 }
1313
1314
1315 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1316 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.038) + \
1317 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.047) + \
1318 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.046) + \
1319 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.044) + \
1320 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.041) + \
1321 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.044) + \
1322 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.041) + \
1323 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.038) + \
1324 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.034) + \
1325 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.027) + \
1326 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.022) + \
1327 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.021) + \
1328 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.016) + \
1329 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.012) + \
1330 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.009) + \
1331 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.005) + \
1332 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.004) + \
1333 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.003) + \
1334 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.002) + \
1335 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.002) + \
1336 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.001) + \
1337 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.002) + \
1338 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.001) + \
1339 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.006) + \
1340 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1341
1342 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1343 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.020) + \
1344 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.0242) + \
1345 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.0272) + \
1346 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.0349) + \
1347 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.0314) + \
1348 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.0220) + \
1349 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.028 ) + \
1350 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.0263) + \
1351 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.0245) + \
1352 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.0150) + \
1353 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.0088) + \
1354 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.0115) + \
1355 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.0063) + \
1356 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.0083) + \
1357 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.0027) + \
1358 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.0052) + \
1359 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.0007) + \
1360 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.0048) + \
1361 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.0013) + \
1362 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.0013) + \
1363 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1364
1365 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1366 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.0092) + \
1367 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.0094) + \
1368 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.0141) + \
1369 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.0108) + \
1370 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.0161) + \
1371 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.0137) + \
1372 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.0135) + \
1373 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.0176) + \
1374 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.0149) + \
1375 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.0137) + \
1376 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.0114) + \
1377 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.0079) + \
1378 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.0112) + \
1379 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.0081)+ \
1380 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.0033) + \
1381 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.0016) + \
1382 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1383
1384 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1385 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.005) + \
1386 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.005) + \
1387 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.005) + \
1388 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.007) + \
1389 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.005) + \
1390 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.008) + \
1391 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.006) + \
1392 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.007) + \
1393 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.008) + \
1394 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.007) + \
1395 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.006) + \
1396 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.004) + \
1397 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.005) + \
1398 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.004) + \
1399 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.001) + \
1400 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.0005) + \
1401 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1402 (abs(eta) > 3.4) * (0.000)
1403 }
1404
1405
1406
1407 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1408 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.349) + \
1409 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.435) + \
1410 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.485) + \
1411 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.504) + \
1412 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.519) + \
1413 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.537) + \
1414 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.541) + \
1415 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.528) + \
1416 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.514) + \
1417 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.473) + \
1418 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.447) + \
1419 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.408) + \
1420 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.373) + \
1421 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.315) + \
1422 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.267) + \
1423 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.181) + \
1424 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.138) + \
1425 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.101) + \
1426 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.067) + \
1427 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.048) + \
1428 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.032) + \
1429 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.029) + \
1430 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.022) + \
1431 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.038) + \
1432 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1433
1434 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1435 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.186) + \
1436 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.258) + \
1437 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.306) + \
1438 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.348) + \
1439 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.359) + \
1440 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.358) + \
1441 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.367) + \
1442 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.380) + \
1443 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.359) + \
1444 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.308) + \
1445 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.238) + \
1446 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.252) + \
1447 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.229) + \
1448 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.174) + \
1449 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.125) + \
1450 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.096) + \
1451 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.066) + \
1452 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.046) + \
1453 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.076) + \
1454 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.072) + \
1455 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1456
1457 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1458 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.096) + \
1459 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.119) + \
1460 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.156) + \
1461 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.174) + \
1462 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.192) + \
1463 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.228) + \
1464 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.218) + \
1465 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.231) + \
1466 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.224) + \
1467 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.184) + \
1468 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.164) + \
1469 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.115) + \
1470 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.123) + \
1471 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.131)+ \
1472 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.063) + \
1473 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.039) + \
1474 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1475
1476 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1477 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.05) + \
1478 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.06) + \
1479 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.08) + \
1480 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.09) + \
1481 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.10) + \
1482 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.11) + \
1483 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.11) + \
1484 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.11) + \
1485 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.11) + \
1486 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.09) + \
1487 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.08) + \
1488 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.06) + \
1489 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.06) + \
1490 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.06) + \
1491 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.03) + \
1492 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.02) + \
1493 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1494 (abs(eta) > 3.4) * (0.000)
1495 }
1496}
1497
1498
1499
1500#############
1501# tau-tagging
1502#############
1503
1504
1505module TauTagging TauTagging {
1506 set ParticleInputArray Delphes/allParticles
1507 set PartonInputArray Delphes/partons
1508 set JetInputArray JetEnergyScale/jets
1509
1510 set DeltaR 0.5
1511
1512 set TauPTMin 20.0
1513
1514 set TauEtaMax 2.3
1515
1516 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1517
1518 add EfficiencyFormula {0} { (abs(eta) < 2.3) * ((( -0.00621816+0.00130097*pt-2.19642e-5*pt^2+1.49393e-7*pt^3-4.58972e-10*pt^4+5.27983e-13*pt^5 )) * (pt<250) + 0.0032*(pt>250)) + \
1519 (abs(eta) > 2.3) * (0.000)
1520 }
1521 add EfficiencyFormula {15} { (abs(eta) < 2.3) * 0.97*0.77*( (0.32 + 0.01*pt - 0.000054*pt*pt )*(pt<100)+0.78*(pt>100) ) + \
1522 (abs(eta) > 2.3) * (0.000)
1523 }
1524}
1525
1526
1527###############################################################################################################
1528# StatusPidFilter: this module removes all generated particles except electrons, muons, taus, and status == 3 #
1529###############################################################################################################
1530
1531module StatusPidFilter GenParticleFilter {
1532
1533 set InputArray Delphes/allParticles
1534 set OutputArray filteredParticles
1535 set PTMin 5.0
1536
1537}
1538
1539
1540##################
1541# ROOT tree writer
1542##################
1543
1544module TreeWriter TreeWriter {
1545# add Branch InputArray BranchName BranchClass
1546 add Branch GenParticleFilter/filteredParticles Particle GenParticle
1547 add Branch PileUpMerger/vertices Vertex Vertex
1548
1549 add Branch GenJetFinder/jets GenJet Jet
1550 add Branch GenMissingET/momentum GenMissingET MissingET
1551
1552# add Branch HCal/eflowTracks EFlowTrack Track
1553# add Branch ECal/eflowPhotons EFlowPhoton Tower
1554# add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1555
1556 add Branch PhotonIsolation/photons Photon Photon
1557 add Branch ElectronIsolation/electrons Electron Electron
1558 add Branch MuonIsolation/muons Muon Muon
1559
1560 add Branch JetEnergyScale/jets Jet Jet
1561
1562 add Branch MissingET/momentum MissingET MissingET
1563 add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET
1564 add Branch ScalarHT/energy ScalarHT ScalarHT
1565}
1566
Note: See TracBrowser for help on using the repository browser.