Fork me on GitHub

source: git/cards/CMS_PhaseII_0PU.tcl@ 576428a

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

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