Fork me on GitHub

source: git/cards/CMS_PhaseII_200PU.tcl@ 2aaf039

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

added muon id efficiencies
  • Property mode set to 100644
File size: 88.2 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 PhotonIsolation
45 PhotonEfficiency
46
47 ElectronIsolation
48 ElectronEfficiency
49
50 MuonIsolation
51 MuonEfficiency
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 * 1.00) + \
194 (abs(eta) <= 1.2) * (pt > 1.0) * (1.00) + \
195 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 0.2 && pt <= 1.0) * (pt*1.00) + \
196 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 1.0) * (1.00) + \
197 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
198 (abs(eta) > 2.8 && 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 # up to |eta| < 2.8 take measurement from tracking + muon chambers
277 # for |eta| > 2.8 and pT < 5.0 take measurement from tracking alone
278
279 set ResolutionFormula {
280
281 (pt > 5 && pt < 10) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.00519 + \
282 (pt > 5 && pt < 10) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.00575 + \
283 (pt > 5 && pt < 10) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.00629 + \
284 (pt > 5 && pt < 10) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.00626 + \
285 (pt > 5 && pt < 10) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.00722 + \
286 (pt > 5 && pt < 10) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.00968 + \
287 (pt > 5 && pt < 10) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.00983 + \
288 (pt > 5 && pt < 10) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.0102 + \
289 (pt > 5 && pt < 10) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.0126 + \
290 (pt > 5 && pt < 10) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.0178 + \
291 (pt > 5 && pt < 10) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.0248 + \
292 (pt > 5 && pt < 10) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.0347 + \
293 (pt > 5 && pt < 10) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.0359 + \
294 (pt > 5 && pt < 10) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.0412 + \
295 (pt > 5 && pt < 10) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.0422 + \
296 (pt > 5 && pt < 10) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.0496 + \
297
298 (pt > 10 && pt < 20) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.00533 + \
299 (pt > 10 && pt < 20) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.00592 + \
300 (pt > 10 && pt < 20) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.00602 + \
301 (pt > 10 && pt < 20) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.00614 + \
302 (pt > 10 && pt < 20) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.00736 + \
303 (pt > 10 && pt < 20) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.0104 + \
304 (pt > 10 && pt < 20) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.00968 + \
305 (pt > 10 && pt < 20) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.00991 + \
306 (pt > 10 && pt < 20) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.0106 + \
307 (pt > 10 && pt < 20) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.0158 + \
308 (pt > 10 && pt < 20) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.0207 + \
309 (pt > 10 && pt < 20) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.0327 + \
310 (pt > 10 && pt < 20) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.0369 + \
311 (pt > 10 && pt < 20) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.039 + \
312 (pt > 10 && pt < 20) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.0488 + \
313 (pt > 10 && pt < 20) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.0543 + \
314
315 (pt > 20 && pt < 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.00559 + \
316 (pt > 20 && pt < 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.00586 + \
317 (pt > 20 && pt < 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.00614 + \
318 (pt > 20 && pt < 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.00654 + \
319 (pt > 20 && pt < 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.00759 + \
320 (pt > 20 && pt < 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.0103 + \
321 (pt > 20 && pt < 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.00963 + \
322 (pt > 20 && pt < 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.0102 + \
323 (pt > 20 && pt < 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.011 + \
324 (pt > 20 && pt < 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.0142 + \
325 (pt > 20 && pt < 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.0197 + \
326 (pt > 20 && pt < 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.0315 + \
327 (pt > 20 && pt < 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.0409 + \
328 (pt > 20 && pt < 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.0471 + \
329 (pt > 20 && pt < 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.0579 + \
330 (pt > 20 && pt < 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.0679 + \
331
332 (pt > 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.00656 + \
333 (pt > 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.007 + \
334 (pt > 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.00733 + \
335 (pt > 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.00752 + \
336 (pt > 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.00851 + \
337 (pt > 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.0105 + \
338 (pt > 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.0105 + \
339 (pt > 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.0111 + \
340 (pt > 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.0119 + \
341 (pt > 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.0165 + \
342 (pt > 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.0242 + \
343 (pt > 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.0359 + \
344 (pt > 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.0474 + \
345 (pt > 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.052 + \
346 (pt > 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.0639 + \
347 (pt > 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.0763 + \
348
349 (abs(eta) > 2.8 && abs(eta) <= 3.0) * (pt > 0.1 && pt <= 1.0) * (0.017) + \
350 (abs(eta) > 2.8 && abs(eta) <= 3.0) * (pt > 1.0 && pt <= 10.0) * (0.03) + \
351 (abs(eta) > 2.8 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 100.0) * (0.05) + \
352 (abs(eta) > 2.8 && abs(eta) <= 3.0) * (pt > 100.0) * (0.30) + \
353 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
354 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 1.0 && pt <= 10.0) * (0.04) + \
355 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 10.0 && pt <= 100.0) * (0.07) + \
356 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 100.0) * (0.30) + \
357 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.025) + \
358 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.05) + \
359 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 100.0) * (0.20) + \
360 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 100.0) * (0.80) + \
361
362 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
363 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 5.0) * (0.013) + \
364 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
365 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 5.0) * (0.015) + \
366 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 0.1 && pt <= 1.0) * (0.017) + \
367 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 1.0 && pt <= 5.0) * (0.03) + \
368 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
369 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 1.0 && pt <= 5.0) * (0.04) + \
370 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.025) + \
371 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 5.0) * (0.05)
372
373 }
374}
375
376
377##############
378# Track merger
379##############
380
381module Merger TrackMerger {
382# add InputArray InputArray
383 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
384 add InputArray ElectronMomentumSmearing/electrons
385 add InputArray MuonMomentumSmearing/muons
386 set OutputArray tracks
387}
388
389#############
390# ECAL
391#############
392
393module SimpleCalorimeter ECal {
394 set ParticleInputArray ParticlePropagator/stableParticles
395 set TrackInputArray TrackMerger/tracks
396
397 set TowerOutputArray ecalTowers
398 set EFlowTrackOutputArray eflowTracks
399 set EFlowTowerOutputArray eflowPhotons
400
401 set IsEcal true
402
403 set EnergyMin 0.5
404 set EnergySignificanceMin 1.0
405
406 set SmearTowerCenter true
407
408 set pi [expr {acos(-1)}]
409
410 # lists of the edges of each tower in eta and phi
411 # each list starts with the lower edge of the first tower
412 # the list ends with the higher edged of the last tower
413
414 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
415
416 set PhiBins {}
417 for {set i -180} {$i <= 180} {incr i} {
418 add PhiBins [expr {$i * $pi/180.0}]
419 }
420
421 # 0.02 unit in eta up to eta = 1.5 (barrel)
422 for {set i -85} {$i <= 86} {incr i} {
423 set eta [expr {$i * 0.0174}]
424 add EtaPhiBins $eta $PhiBins
425 }
426
427 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
428
429 set PhiBins {}
430 for {set i -180} {$i <= 180} {incr i} {
431 add PhiBins [expr {$i * $pi/180.0}]
432 }
433
434 # 0.02 unit in eta up to eta = 3
435 for {set i 1} {$i <= 84} {incr i} {
436 set eta [expr { -2.958 + $i * 0.0174}]
437 add EtaPhiBins $eta $PhiBins
438 }
439
440 for {set i 1} {$i <= 84} {incr i} {
441 set eta [expr { 1.4964 + $i * 0.0174}]
442 add EtaPhiBins $eta $PhiBins
443 }
444
445 # take present CMS granularity for HF
446
447 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
448 set PhiBins {}
449 for {set i -18} {$i <= 18} {incr i} {
450 add PhiBins [expr {$i * $pi/18.0}]
451 }
452
453 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} {
454 add EtaPhiBins $eta $PhiBins
455 }
456
457
458 add EnergyFraction {0} {0.0}
459 # energy fractions for e, gamma and pi0
460 add EnergyFraction {11} {1.0}
461 add EnergyFraction {22} {1.0}
462 add EnergyFraction {111} {1.0}
463 # energy fractions for muon, neutrinos and neutralinos
464 add EnergyFraction {12} {0.0}
465 add EnergyFraction {13} {0.0}
466 add EnergyFraction {14} {0.0}
467 add EnergyFraction {16} {0.0}
468 add EnergyFraction {1000022} {0.0}
469 add EnergyFraction {1000023} {0.0}
470 add EnergyFraction {1000025} {0.0}
471 add EnergyFraction {1000035} {0.0}
472 add EnergyFraction {1000045} {0.0}
473 # energy fractions for K0short and Lambda
474 add EnergyFraction {310} {0.3}
475 add EnergyFraction {3122} {0.3}
476
477 # set ResolutionFormula {resolution formula as a function of eta and energy}
478
479 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016.
480 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117
481
482 set ResolutionFormula { (abs(eta) <= 1.50) * sqrt(energy^2*0.005^2 + energy*0.027^2 + 0.15^2) + \
483 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
484 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
485 (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
486 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
487
488}
489
490#############
491# HCAL
492#############
493
494module SimpleCalorimeter HCal {
495 set ParticleInputArray ParticlePropagator/stableParticles
496 set TrackInputArray ECal/eflowTracks
497
498 set TowerOutputArray hcalTowers
499 set EFlowTrackOutputArray eflowTracks
500 set EFlowTowerOutputArray eflowNeutralHadrons
501
502 set IsEcal false
503
504 set EnergyMin 1.0
505 set EnergySignificanceMin 1.0
506
507 set SmearTowerCenter true
508
509 set pi [expr {acos(-1)}]
510
511 # lists of the edges of each tower in eta and phi
512 # each list starts with the lower edge of the first tower
513 # the list ends with the higher edged of the last tower
514
515 # assume 0.087 x 0.087 resolution in eta,phi in the barrel |eta| < 1.5
516
517 set PhiBins {}
518 for {set i -36} {$i <= 36} {incr i} {
519 add PhiBins [expr {$i * $pi/36.0}]
520 }
521 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} {
522 add EtaPhiBins $eta $PhiBins
523 }
524
525 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- HCAL)
526
527 set PhiBins {}
528 for {set i -180} {$i <= 180} {incr i} {
529 add PhiBins [expr {$i * $pi/180.0}]
530 }
531
532 # 0.02 unit in eta up to eta = 3
533 for {set i 1} {$i <= 84} {incr i} {
534 set eta [expr { -2.958 + $i * 0.0174}]
535 add EtaPhiBins $eta $PhiBins
536 }
537
538 for {set i 1} {$i <= 84} {incr i} {
539 set eta [expr { 1.4964 + $i * 0.0174}]
540 add EtaPhiBins $eta $PhiBins
541 }
542
543 # take present CMS granularity for HF
544
545 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
546 set PhiBins {}
547 for {set i -18} {$i <= 18} {incr i} {
548 add PhiBins [expr {$i * $pi/18.0}]
549 }
550
551 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} {
552 add EtaPhiBins $eta $PhiBins
553 }
554
555
556 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
557 add EnergyFraction {0} {1.0}
558 # energy fractions for e, gamma and pi0
559 add EnergyFraction {11} {0.0}
560 add EnergyFraction {22} {0.0}
561 add EnergyFraction {111} {0.0}
562 # energy fractions for muon, neutrinos and neutralinos
563 add EnergyFraction {12} {0.0}
564 add EnergyFraction {13} {0.0}
565 add EnergyFraction {14} {0.0}
566 add EnergyFraction {16} {0.0}
567 add EnergyFraction {1000022} {0.0}
568 add EnergyFraction {1000023} {0.0}
569 add EnergyFraction {1000025} {0.0}
570 add EnergyFraction {1000035} {0.0}
571 add EnergyFraction {1000045} {0.0}
572 # energy fractions for K0short and Lambda
573 add EnergyFraction {310} {0.7}
574 add EnergyFraction {3122} {0.7}
575
576# set ResolutionFormula {resolution formula as a function of eta and energy}
577 set ResolutionFormula { (abs(eta) <= 1.5) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
578 (abs(eta) > 1.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
579 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.11^2 + energy*2.80^2)}
580
581}
582
583#################
584# Electron filter
585#################
586
587module PdgCodeFilter ElectronFilter {
588 set InputArray HCal/eflowTracks
589 set OutputArray electrons
590 set Invert true
591 add PdgCode {11}
592 add PdgCode {-11}
593}
594
595##########################
596# Track pile-up subtractor
597##########################
598
599module TrackPileUpSubtractor TrackPileUpSubtractor {
600# add InputArray InputArray OutputArray
601 add InputArray HCal/eflowTracks eflowTracks
602 add InputArray ElectronFilter/electrons electrons
603 add InputArray MuonMomentumSmearing/muons muons
604
605 set VertexInputArray PileUpMerger/vertices
606 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
607 # Z vertex resolution in m
608 set ZVertexResolution 0.0001
609}
610
611
612###################################################
613# Tower Merger (in case not using e-flow algorithm)
614###################################################
615
616module Merger TowerMerger {
617# add InputArray InputArray
618 add InputArray ECal/ecalTowers
619 add InputArray HCal/hcalTowers
620 set OutputArray towers
621}
622
623
624####################
625# Neutral eflow erger
626####################
627
628module Merger NeutralEFlowMerger {
629# add InputArray InputArray
630 add InputArray ECal/eflowPhotons
631 add InputArray HCal/eflowNeutralHadrons
632 set OutputArray eflowTowers
633}
634
635
636####################
637# Energy flow merger
638####################
639
640module Merger EFlowMerger {
641# add InputArray InputArray
642 add InputArray HCal/eflowTracks
643 add InputArray ECal/eflowPhotons
644 add InputArray HCal/eflowNeutralHadrons
645 set OutputArray eflow
646}
647
648##################################
649# Energy flow merger (all tracks)
650##################################
651
652module Merger EFlowMergerAllTracks {
653# add InputArray InputArray
654 add InputArray TrackMerger/tracks
655 add InputArray ECal/eflowPhotons
656 add InputArray HCal/eflowNeutralHadrons
657 set OutputArray eflow
658}
659
660#########################################
661### Run the puppi code (to be tuned) ###
662#########################################
663
664module RunPUPPI RunPUPPI {
665 ## input information
666 set TrackInputArray TrackMerger/tracks
667 set NeutralInputArray NeutralEFlowMerger/eflowTowers
668 set PVInputArray PileUpMerger/vertices
669 set MinPuppiWeight 0.05
670 set UseExp false
671
672 ## define puppi algorithm parameters (more than one for the same eta region is possible)
673 add EtaMinBin 0. 2.5 2.5 3.0 3.0
674 add EtaMaxBin 2.5 3.0 3.0 10.0 10.0
675 add PtMinBin 0. 0.5 0.5 0.5 0.5
676 add ConeSizeBin 0.25 0.25 0.25 0.25 0.25
677 add RMSPtMinBin 0.1 0.5 0.5 0.5 0.5
678 add RMSScaleFactorBin 1.0 1.0 1.0 1.0 1.0
679 add NeutralMinEBin 0.2 1.0 1.0 1.5 1.5
680 add NeutralPtSlope 0.02 0.02 0.02 0.02 0.02
681 add ApplyCHS true true true true true
682 add UseCharged true false false false false
683 add ApplyLowPUCorr true true true true true
684 add MetricId 5 5 1 5 1
685
686 ## output name
687 set OutputArray PuppiParticles
688 set OutputArrayTracks puppiTracks
689 set OutputArrayNeutrals puppiNeutrals
690}
691
692
693
694###################
695# Missing ET merger
696###################
697
698module Merger MissingET {
699# add InputArray InputArray
700# add InputArray RunPUPPI/PuppiParticles
701 add InputArray EFlowMerger/eflow
702 set MomentumOutputArray momentum
703}
704
705###################
706# Ger PileUp Missing ET
707###################
708
709module Merger GenPileUpMissingET {
710# add InputArray InputArray
711# add InputArray RunPUPPI/PuppiParticles
712 add InputArray ParticlePropagator/stableParticles
713 set MomentumOutputArray momentum
714}
715
716##################
717# Scalar HT merger
718##################
719
720module Merger ScalarHT {
721# add InputArray InputArray
722 add InputArray RunPUPPI/PuppiParticles
723 set EnergyOutputArray energy
724}
725
726#################
727# Neutrino Filter
728#################
729
730module PdgCodeFilter NeutrinoFilter {
731
732 set InputArray Delphes/stableParticles
733 set OutputArray filteredParticles
734
735 set PTMin 0.0
736
737 add PdgCode {12}
738 add PdgCode {14}
739 add PdgCode {16}
740 add PdgCode {-12}
741 add PdgCode {-14}
742 add PdgCode {-16}
743
744}
745
746
747#####################
748# MC truth jet finder
749#####################
750
751module FastJetFinder GenJetFinder {
752 set InputArray NeutrinoFilter/filteredParticles
753
754 set OutputArray jets
755
756 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
757 set JetAlgorithm 6
758 set ParameterR 0.4
759
760 set JetPTMin 15.0
761}
762
763#########################
764# Gen Missing ET merger
765########################
766
767module Merger GenMissingET {
768
769# add InputArray InputArray
770 add InputArray NeutrinoFilter/filteredParticles
771 set MomentumOutputArray momentum
772}
773
774
775
776############
777# Jet finder
778############
779
780module FastJetFinder FastJetFinder {
781# set InputArray TowerMerger/towers
782 set InputArray RunPUPPI/PuppiParticles
783
784 set OutputArray jets
785
786 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
787 set JetAlgorithm 6
788 set ParameterR 0.4
789
790 set JetPTMin 15.0
791}
792
793##################
794# Jet Energy Scale
795##################
796
797module EnergyScale JetEnergyScale {
798 set InputArray FastJetFinder/jets
799 set OutputArray jets
800
801 # scale formula for jets
802 set ScaleFormula {1.00}
803}
804
805
806####################
807# Photon isolation #
808####################
809
810module Isolation PhotonIsolation {
811
812 # particle for which calculate the isolation
813 set CandidateInputArray ECal/eflowPhotons
814
815 # isolation collection
816 set IsolationInputArray EFlowMerger/eflow
817
818 # output array
819 set OutputArray photons
820
821 # isolation cone
822 set DeltaRMax 0.3
823
824 # minimum pT
825 set PTMin 1.0
826
827 # iso ratio to cut
828 set PTRatioMax 9999.
829}
830
831
832
833#####################
834# Photon efficiency #
835#####################
836
837module Efficiency PhotonEfficiency {
838
839 ## input particles
840 set InputArray PhotonIsolation/photons
841 ## output particles
842 set OutputArray photons
843 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
844 # efficiency formula for photons
845 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
846 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
847 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
848 (abs(eta) > 4.0) * (0.00)}
849
850}
851
852
853######################
854# Electron isolation #
855######################
856
857module Isolation ElectronIsolation {
858
859 set CandidateInputArray ElectronFilter/electrons
860
861 # isolation collection
862 set IsolationInputArray EFlowMerger/eflow
863
864 set OutputArray electrons
865
866 set DeltaRMax 0.3
867 set PTMin 1.0
868 set PTRatioMax 9999.
869
870}
871
872
873
874#######################
875# Electron efficiency #
876#######################
877
878module Efficiency ElectronEfficiency {
879
880 set InputArray ElectronIsolation/electrons
881 set OutputArray electrons
882
883 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
884 # efficiency formula for electrons
885 set EfficiencyFormula {
886 (pt <= 4.0) * (0.00) + \
887 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
888 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
889 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
890 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
891 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
892 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
893 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
894 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
895 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
896 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
897 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
898 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
899 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
900 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
901 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
902 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
903 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
904 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
905 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
906 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
907 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
908 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
909 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
910 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
911 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
912 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
913 (abs(eta) > 4.0) * (0.00)
914
915 }
916}
917
918##################
919# Muon isolation #
920##################
921
922module Isolation MuonIsolation {
923 set CandidateInputArray MuonMomentumSmearing/muons
924
925 # isolation collection
926 set IsolationInputArray EFlowMerger/eflow
927
928 set OutputArray muons
929
930 set DeltaRMax 0.3
931 set PTMin 1.0
932 set PTRatioMax 9999.
933
934}
935
936
937##################
938# Muon Loose Id #
939##################
940
941module Efficiency MuonLooseIdEfficiency {
942 set InputArray MuonIsolation/muons
943 set OutputArray muons
944 # tracking + LooseID efficiency formula for muons
945 set EfficiencyFormula {
946
947 (pt <= 5.0) * (0.00) + \
948 (pt > 5 && pt < 10) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.967 + \
949 (pt > 5 && pt < 10) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.968 + \
950 (pt > 5 && pt < 10) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.986 + \
951 (pt > 5 && pt < 10) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.982 + \
952 (pt > 5 && pt < 10) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.972 + \
953 (pt > 5 && pt < 10) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.974 + \
954 (pt > 5 && pt < 10) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.969 + \
955 (pt > 5 && pt < 10) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.988 + \
956 (pt > 5 && pt < 10) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.989 + \
957 (pt > 5 && pt < 10) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.992 + \
958 (pt > 5 && pt < 10) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.979 + \
959 (pt > 5 && pt < 10) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.976 + \
960 (pt > 5 && pt < 10) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.882 + \
961 (pt > 5 && pt < 10) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.883 + \
962 (pt > 5 && pt < 10) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.851 + \
963 (pt > 5 && pt < 10) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.828 + \
964
965 (pt > 10 && pt < 20) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.989 + \
966 (pt > 10 && pt < 20) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.989 + \
967 (pt > 10 && pt < 20) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.995 + \
968 (pt > 10 && pt < 20) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.995 + \
969 (pt > 10 && pt < 20) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.987 + \
970 (pt > 10 && pt < 20) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.989 + \
971 (pt > 10 && pt < 20) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.981 + \
972 (pt > 10 && pt < 20) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.99 + \
973 (pt > 10 && pt < 20) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.997 + \
974 (pt > 10 && pt < 20) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.997 + \
975 (pt > 10 && pt < 20) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.989 + \
976 (pt > 10 && pt < 20) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.995 + \
977 (pt > 10 && pt < 20) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.888 + \
978 (pt > 10 && pt < 20) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.894 + \
979 (pt > 10 && pt < 20) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.886 + \
980 (pt > 10 && pt < 20) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.883 + \
981
982 (pt > 20 && pt < 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.981 + \
983 (pt > 20 && pt < 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.994 + \
984 (pt > 20 && pt < 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.995 + \
985 (pt > 20 && pt < 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.997 + \
986 (pt > 20 && pt < 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.986 + \
987 (pt > 20 && pt < 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.991 + \
988 (pt > 20 && pt < 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.987 + \
989 (pt > 20 && pt < 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.992 + \
990 (pt > 20 && pt < 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.993 + \
991 (pt > 20 && pt < 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.999 + \
992 (pt > 20 && pt < 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.994 + \
993 (pt > 20 && pt < 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.997 + \
994 (pt > 20 && pt < 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.894 + \
995 (pt > 20 && pt < 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.879 + \
996 (pt > 20 && pt < 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.871 + \
997 (pt > 20 && pt < 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.864 + \
998
999 (pt > 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.986 + \
1000 (pt > 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.989 + \
1001 (pt > 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.996 + \
1002 (pt > 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.994 + \
1003 (pt > 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.991 + \
1004 (pt > 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.989 + \
1005 (pt > 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.985 + \
1006 (pt > 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.993 + \
1007 (pt > 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.996 + \
1008 (pt > 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.997 + \
1009 (pt > 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.99 + \
1010 (pt > 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.993 + \
1011 (pt > 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.905 + \
1012 (pt > 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.86 + \
1013 (pt > 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.845 + \
1014 (pt > 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.855 + \
1015
1016 (abs(eta) > 2.8) * (0.00)
1017
1018 }
1019}
1020
1021
1022##################
1023# Muon Tight Id #
1024##################
1025
1026module Efficiency MuonTightIdEfficiency {
1027 set InputArray MuonIsolation/muons
1028 set OutputArray muons
1029 # tracking + TightID efficiency formula for muons
1030 set EfficiencyFormula {
1031
1032 (pt <= 5.0) * (0.00) + \
1033 (pt > 5 && pt < 10) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.949 + \
1034 (pt > 5 && pt < 10) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.894 + \
1035 (pt > 5 && pt < 10) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.966 + \
1036 (pt > 5 && pt < 10) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.964 + \
1037 (pt > 5 && pt < 10) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.923 + \
1038 (pt > 5 && pt < 10) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.94 + \
1039 (pt > 5 && pt < 10) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.959 + \
1040 (pt > 5 && pt < 10) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.972 + \
1041 (pt > 5 && pt < 10) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.962 + \
1042 (pt > 5 && pt < 10) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.972 + \
1043 (pt > 5 && pt < 10) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.947 + \
1044 (pt > 5 && pt < 10) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.911 + \
1045 (pt > 5 && pt < 10) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.825 + \
1046 (pt > 5 && pt < 10) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.821 + \
1047 (pt > 5 && pt < 10) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.783 + \
1048 (pt > 5 && pt < 10) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.746 + \
1049
1050 (pt > 10 && pt < 20) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.975 + \
1051 (pt > 10 && pt < 20) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.935 + \
1052 (pt > 10 && pt < 20) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.986 + \
1053 (pt > 10 && pt < 20) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.976 + \
1054 (pt > 10 && pt < 20) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.954 + \
1055 (pt > 10 && pt < 20) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.978 + \
1056 (pt > 10 && pt < 20) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.975 + \
1057 (pt > 10 && pt < 20) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.98 + \
1058 (pt > 10 && pt < 20) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.982 + \
1059 (pt > 10 && pt < 20) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.985 + \
1060 (pt > 10 && pt < 20) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.962 + \
1061 (pt > 10 && pt < 20) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.946 + \
1062 (pt > 10 && pt < 20) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.839 + \
1063 (pt > 10 && pt < 20) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.844 + \
1064 (pt > 10 && pt < 20) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.803 + \
1065 (pt > 10 && pt < 20) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.791 + \
1066
1067 (pt > 20 && pt < 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.965 + \
1068 (pt > 20 && pt < 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.938 + \
1069 (pt > 20 && pt < 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.99 + \
1070 (pt > 20 && pt < 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.977 + \
1071 (pt > 20 && pt < 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.959 + \
1072 (pt > 20 && pt < 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.976 + \
1073 (pt > 20 && pt < 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.984 + \
1074 (pt > 20 && pt < 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.981 + \
1075 (pt > 20 && pt < 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.978 + \
1076 (pt > 20 && pt < 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.985 + \
1077 (pt > 20 && pt < 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.969 + \
1078 (pt > 20 && pt < 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.964 + \
1079 (pt > 20 && pt < 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.863 + \
1080 (pt > 20 && pt < 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.828 + \
1081 (pt > 20 && pt < 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.807 + \
1082 (pt > 20 && pt < 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.764 + \
1083
1084 (pt > 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.975 + \
1085 (pt > 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.935 + \
1086 (pt > 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.991 + \
1087 (pt > 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.983 + \
1088 (pt > 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.957 + \
1089 (pt > 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.979 + \
1090 (pt > 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.981 + \
1091 (pt > 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.983 + \
1092 (pt > 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.981 + \
1093 (pt > 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.984 + \
1094 (pt > 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.964 + \
1095 (pt > 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.956 + \
1096 (pt > 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.855 + \
1097 (pt > 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.817 + \
1098 (pt > 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.764 + \
1099 (pt > 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.752 + \
1100
1101 (abs(eta) > 2.8) * (0.00)
1102
1103 }
1104}
1105
1106
1107########################
1108# Jet Flavor Association
1109########################
1110
1111module JetFlavorAssociation JetFlavorAssociation {
1112
1113 set PartonInputArray Delphes/partons
1114 set ParticleInputArray Delphes/allParticles
1115 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1116 set JetInputArray JetEnergyScale/jets
1117
1118 set DeltaR 0.5
1119 set PartonPTMin 10.0
1120 set PartonEtaMax 4.0
1121
1122}
1123
1124
1125#############
1126# b-tagging #
1127#############
1128module BTagging BTaggingLoose {
1129
1130 set JetInputArray JetEnergyScale/jets
1131
1132 set BitNumber 0
1133
1134 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
1135 (pt > 20.0) * (abs(eta) <= 3.4) * (0.1) + \
1136 (abs(eta) > 3.4) * (0.000)
1137 }
1138
1139
1140 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1141 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.341) + \
1142 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.389) + \
1143 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.416) + \
1144 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.432) + \
1145 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.447) + \
1146 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.451) + \
1147 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.451) + \
1148 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.447) + \
1149 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.447) + \
1150 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.444) + \
1151 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.430) + \
1152 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.424) + \
1153 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.419) + \
1154 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.407) + \
1155 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.384) + \
1156 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.364) + \
1157 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.353) + \
1158 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.303) + \
1159 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.273) + \
1160 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.249) + \
1161 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.233) + \
1162 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.209) + \
1163 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.194) + \
1164 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.227) + \
1165 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1166
1167 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1168 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.236) + \
1169 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.276) + \
1170 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.302) + \
1171 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.330) + \
1172 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.338) + \
1173 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.355) + \
1174 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.353) + \
1175 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.358) + \
1176 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.350) + \
1177 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.327) + \
1178 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.313) + \
1179 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.307) + \
1180 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.322) + \
1181 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.299) + \
1182 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.274) + \
1183 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.240) + \
1184 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.203) + \
1185 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.218) + \
1186 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.203) + \
1187 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.178) + \
1188 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1189
1190 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1191 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.179) + \
1192 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.198) + \
1193 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.219) + \
1194 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.225) + \
1195 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.232) + \
1196 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.251) + \
1197 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.267) + \
1198 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.252) + \
1199 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.243) + \
1200 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.244) + \
1201 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.242) + \
1202 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.246) + \
1203 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.217) + \
1204 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.223)+ \
1205 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.194) + \
1206 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.187) + \
1207 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1208
1209 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1210 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.09) + \
1211 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.10) + \
1212 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.11) + \
1213 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.11) + \
1214 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.12) + \
1215 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.13) + \
1216 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.13) + \
1217 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.13) + \
1218 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.12) + \
1219 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.12) + \
1220 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.12) + \
1221 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.12) + \
1222 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.11) + \
1223 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.11) + \
1224 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.10) + \
1225 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.09) + \
1226 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1227 (abs(eta) > 3.4) * (0.000)
1228 }
1229
1230 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1231 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.712) + \
1232 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.788) + \
1233 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.837) + \
1234 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.856) + \
1235 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.873) + \
1236 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.886) + \
1237 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.883) + \
1238 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.882) + \
1239 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.883) + \
1240 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.883) + \
1241 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.876) + \
1242 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.872) + \
1243 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.867) + \
1244 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.853) + \
1245 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.831) + \
1246 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.819) + \
1247 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.791) + \
1248 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.740) + \
1249 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.699) + \
1250 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.631) + \
1251 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.600) + \
1252 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.567) + \
1253 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.478) + \
1254 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.470) + \
1255 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1256
1257 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1258 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.543) + \
1259 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.627) + \
1260 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.701) + \
1261 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.741) + \
1262 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.754) + \
1263 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.796) + \
1264 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.778) + \
1265 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.798) + \
1266 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.788) + \
1267 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.775) + \
1268 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.765) + \
1269 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.752) + \
1270 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.740) + \
1271 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.716) + \
1272 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.671) + \
1273 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.637) + \
1274 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.577) + \
1275 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.572) + \
1276 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.514) + \
1277 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.516) + \
1278 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1279
1280 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1281 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.407) + \
1282 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.479) + \
1283 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.529) + \
1284 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.578) + \
1285 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.612) + \
1286 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.673) + \
1287 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.656) + \
1288 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.658) + \
1289 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.647) + \
1290 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.629) + \
1291 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.630) + \
1292 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.615) + \
1293 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.604) + \
1294 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.532)+ \
1295 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.532) + \
1296 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.503) + \
1297 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1298
1299 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1300 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.20) + \
1301 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.235) + \
1302 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.26) + \
1303 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.28) + \
1304 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.31) + \
1305 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.33) + \
1306 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.32) + \
1307 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.32) + \
1308 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.33) + \
1309 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.31) + \
1310 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.32) + \
1311 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.31) + \
1312 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.30) + \
1313 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.26) + \
1314 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.27) + \
1315 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.25) + \
1316 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1317 (abs(eta) > 3.4) * (0.000)
1318 }
1319}
1320
1321
1322#############
1323# b-tagging #
1324#############
1325module BTagging BTaggingMedium {
1326
1327 set JetInputArray JetEnergyScale/jets
1328
1329 set BitNumber 1
1330
1331 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
1332 (pt > 20.0) * (abs(eta) <= 3.4) * (0.01) + \
1333 (abs(eta) > 3.4) * (0.000)
1334 }
1335
1336 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1337 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.144) + \
1338 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.173) + \
1339 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.188) + \
1340 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.200) + \
1341 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.209) + \
1342 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.210) + \
1343 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.206) + \
1344 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.205) + \
1345 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.198) + \
1346 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.187) + \
1347 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.175) + \
1348 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.167) + \
1349 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.157) + \
1350 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.138) + \
1351 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.118) + \
1352 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.094) + \
1353 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.083) + \
1354 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.062) + \
1355 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.049) + \
1356 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.044) + \
1357 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.036) + \
1358 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.035) + \
1359 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.026) + \
1360 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.041) + \
1361 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1362
1363 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1364 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.070) + \
1365 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.089) + \
1366 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.104) + \
1367 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.121) + \
1368 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.121) + \
1369 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.127) + \
1370 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.127) + \
1371 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.117) + \
1372 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.117) + \
1373 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.096) + \
1374 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.099) + \
1375 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.086) + \
1376 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.075) + \
1377 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.066) + \
1378 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.056) + \
1379 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.048) + \
1380 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.031) + \
1381 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.042) + \
1382 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.030) + \
1383 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.030) + \
1384 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1385
1386 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1387 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.041) + \
1388 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.047) + \
1389 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.059) + \
1390 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.056) + \
1391 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.068) + \
1392 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.073) + \
1393 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.080) + \
1394 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.063) + \
1395 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.065) + \
1396 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.057) + \
1397 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.049) + \
1398 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.067) + \
1399 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.065) + \
1400 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.044)+ \
1401 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.031) + \
1402 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.026) + \
1403 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1404
1405 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1406 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.02) + \
1407 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.02) + \
1408 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.03) + \
1409 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.03) + \
1410 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.03) + \
1411 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.035) + \
1412 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.04) + \
1413 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.03) + \
1414 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.03) + \
1415 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.025) + \
1416 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.025) + \
1417 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.035) + \
1418 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.032) + \
1419 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.022) + \
1420 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.015) + \
1421 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.013) + \
1422 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1423 (abs(eta) > 3.4) * (0.000)
1424 }
1425
1426
1427 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1428 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.528) + \
1429 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.611) + \
1430 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.666) + \
1431 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.701) + \
1432 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.734) + \
1433 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.741) + \
1434 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.746) + \
1435 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.750) + \
1436 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.742) + \
1437 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.729) + \
1438 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.711) + \
1439 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.694) + \
1440 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.684) + \
1441 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.642) + \
1442 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.602) + \
1443 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.554) + \
1444 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.488) + \
1445 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.417) + \
1446 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.343) + \
1447 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.284) + \
1448 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.250) + \
1449 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.221) + \
1450 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.164) + \
1451 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.173) + \
1452 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1453
1454 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1455 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.322) + \
1456 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.414) + \
1457 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.488) + \
1458 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.521) + \
1459 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.554) + \
1460 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.590) + \
1461 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.571) + \
1462 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.583) + \
1463 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.575) + \
1464 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.557) + \
1465 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.536) + \
1466 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.516) + \
1467 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.483) + \
1468 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.432) + \
1469 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.366) + \
1470 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.323) + \
1471 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.241) + \
1472 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.252) + \
1473 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.237) + \
1474 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.243) + \
1475 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1476
1477 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1478 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.200) + \
1479 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.252) + \
1480 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.310) + \
1481 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.353) + \
1482 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.377) + \
1483 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.425) + \
1484 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.406) + \
1485 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.414) + \
1486 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.412) + \
1487 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.366) + \
1488 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.360) + \
1489 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.347) + \
1490 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.310) + \
1491 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.307)+ \
1492 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.224) + \
1493 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.194) + \
1494 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1495
1496 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1497 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.100) + \
1498 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.12) + \
1499 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.15) + \
1500 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.17) + \
1501 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.19) + \
1502 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.21) + \
1503 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.20) + \
1504 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.21) + \
1505 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.20) + \
1506 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.18) + \
1507 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.18) + \
1508 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.17) + \
1509 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.15) + \
1510 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.15) + \
1511 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.11) + \
1512 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.10) + \
1513 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1514 (abs(eta) > 3.4) * (0.000)
1515 }
1516}
1517
1518
1519
1520#############
1521# b-tagging #
1522#############
1523module BTagging BTaggingTight {
1524
1525 set JetInputArray JetEnergyScale/jets
1526
1527 set BitNumber 2
1528
1529 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
1530 (pt > 20.0) * (abs(eta) <= 3.4) * (0.001) + \
1531 (abs(eta) > 3.4) * (0.000)
1532 }
1533
1534
1535 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1536 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.038) + \
1537 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.047) + \
1538 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.046) + \
1539 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.044) + \
1540 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.041) + \
1541 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.044) + \
1542 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.041) + \
1543 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.038) + \
1544 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.034) + \
1545 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.027) + \
1546 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.022) + \
1547 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.021) + \
1548 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.016) + \
1549 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.012) + \
1550 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.009) + \
1551 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.005) + \
1552 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.004) + \
1553 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.003) + \
1554 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.002) + \
1555 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.002) + \
1556 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.001) + \
1557 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.002) + \
1558 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.001) + \
1559 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.006) + \
1560 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1561
1562 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1563 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.020) + \
1564 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.0242) + \
1565 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.0272) + \
1566 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.0349) + \
1567 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.0314) + \
1568 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.0220) + \
1569 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.028 ) + \
1570 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.0263) + \
1571 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.0245) + \
1572 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.0150) + \
1573 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.0088) + \
1574 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.0115) + \
1575 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.0063) + \
1576 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.0083) + \
1577 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.0027) + \
1578 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.0052) + \
1579 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.0007) + \
1580 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.0048) + \
1581 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.0013) + \
1582 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.0013) + \
1583 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1584
1585 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1586 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.0092) + \
1587 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.0094) + \
1588 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.0141) + \
1589 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.0108) + \
1590 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.0161) + \
1591 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.0137) + \
1592 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.0135) + \
1593 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.0176) + \
1594 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.0149) + \
1595 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.0137) + \
1596 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.0114) + \
1597 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.0079) + \
1598 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.0112) + \
1599 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.0081)+ \
1600 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.0033) + \
1601 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.0016) + \
1602 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1603
1604 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1605 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.005) + \
1606 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.005) + \
1607 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.005) + \
1608 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.007) + \
1609 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.005) + \
1610 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.008) + \
1611 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.006) + \
1612 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.007) + \
1613 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.008) + \
1614 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.007) + \
1615 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.006) + \
1616 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.004) + \
1617 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.005) + \
1618 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.004) + \
1619 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.001) + \
1620 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.0005) + \
1621 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1622 (abs(eta) > 3.4) * (0.000)
1623 }
1624
1625
1626
1627 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1628 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.349) + \
1629 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.435) + \
1630 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.485) + \
1631 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.504) + \
1632 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.519) + \
1633 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.537) + \
1634 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.541) + \
1635 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.528) + \
1636 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.514) + \
1637 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.473) + \
1638 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.447) + \
1639 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.408) + \
1640 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.373) + \
1641 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.315) + \
1642 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.267) + \
1643 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.181) + \
1644 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.138) + \
1645 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.101) + \
1646 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.067) + \
1647 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.048) + \
1648 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.032) + \
1649 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.029) + \
1650 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.022) + \
1651 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.038) + \
1652 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1653
1654 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1655 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.186) + \
1656 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.258) + \
1657 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.306) + \
1658 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.348) + \
1659 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.359) + \
1660 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.358) + \
1661 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.367) + \
1662 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.380) + \
1663 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.359) + \
1664 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.308) + \
1665 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.238) + \
1666 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.252) + \
1667 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.229) + \
1668 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.174) + \
1669 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.125) + \
1670 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.096) + \
1671 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.066) + \
1672 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.046) + \
1673 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.076) + \
1674 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.072) + \
1675 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1676
1677 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1678 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.096) + \
1679 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.119) + \
1680 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.156) + \
1681 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.174) + \
1682 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.192) + \
1683 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.228) + \
1684 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.218) + \
1685 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.231) + \
1686 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.224) + \
1687 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.184) + \
1688 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.164) + \
1689 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.115) + \
1690 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.123) + \
1691 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.131)+ \
1692 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.063) + \
1693 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.039) + \
1694 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1695
1696 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1697 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.05) + \
1698 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.06) + \
1699 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.08) + \
1700 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.09) + \
1701 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.10) + \
1702 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.11) + \
1703 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.11) + \
1704 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.11) + \
1705 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.11) + \
1706 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.09) + \
1707 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.08) + \
1708 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.06) + \
1709 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.06) + \
1710 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.06) + \
1711 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.03) + \
1712 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.02) + \
1713 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1714 (abs(eta) > 3.4) * (0.000)
1715 }
1716}
1717
1718
1719
1720#############
1721# tau-tagging
1722#############
1723
1724
1725module TauTagging TauTagging {
1726 set ParticleInputArray Delphes/allParticles
1727 set PartonInputArray Delphes/partons
1728 set JetInputArray JetEnergyScale/jets
1729
1730 set DeltaR 0.5
1731
1732 set TauPTMin 20.0
1733
1734 set TauEtaMax 2.3
1735
1736 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1737
1738 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)) + \
1739 (abs(eta) > 2.3) * (0.000)
1740 }
1741 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) ) + \
1742 (abs(eta) > 2.3) * (0.000)
1743 }
1744}
1745
1746
1747###############################################################################################################
1748# StatusPidFilter: this module removes all generated particles except electrons, muons, taus, and status == 3 #
1749###############################################################################################################
1750
1751module StatusPidFilter GenParticleFilter {
1752
1753 set InputArray Delphes/allParticles
1754 set OutputArray filteredParticles
1755 set PTMin 5.0
1756
1757}
1758
1759
1760##################
1761# ROOT tree writer
1762##################
1763
1764module TreeWriter TreeWriter {
1765# add Branch InputArray BranchName BranchClass
1766 add Branch GenParticleFilter/filteredParticles Particle GenParticle
1767 add Branch PileUpMerger/vertices Vertex Vertex
1768
1769 add Branch GenJetFinder/jets GenJet Jet
1770 add Branch GenMissingET/momentum GenMissingET MissingET
1771
1772# add Branch HCal/eflowTracks EFlowTrack Track
1773# add Branch ECal/eflowPhotons EFlowPhoton Tower
1774# add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1775
1776 add Branch PhotonEfficiency/photons Photon Photon
1777 add Branch ElectronEfficiency/electrons Electron Electron
1778 add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
1779 add Branch MuonTightIdEfficiency/muons MuonTight Muon
1780
1781 add Branch JetEnergyScale/jets Jet Jet
1782
1783 add Branch MissingET/momentum MissingET MissingET
1784 add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET
1785 add Branch ScalarHT/energy ScalarHT ScalarHT
1786}
1787
Note: See TracBrowser for help on using the repository browser.