Fork me on GitHub

source: git/cards/CMS_PhaseII_200PU.tcl@ b67a751

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

udated track reso. for tilted
  • Property mode set to 100644
File size: 106.8 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 MuonLooseIdEfficiency
52 MuonTightIdEfficiency
53
54 NeutrinoFilter
55
56 MissingET
57 GenMissingET
58 GenPileUpMissingET
59
60 GenJetFinder
61 FastJetFinder
62
63 ScalarHT
64
65 JetEnergyScale
66
67 JetFlavorAssociation
68
69 BTaggingLoose
70 BTaggingMedium
71 BTaggingTight
72
73 TauTagging
74
75 GenParticleFilter
76
77 TreeWriter
78}
79
80
81
82###############
83# PileUp Merger
84###############
85
86module PileUpMerger PileUpMerger {
87 set InputArray Delphes/stableParticles
88
89 set ParticleOutputArray stableParticles
90 set VertexOutputArray vertices
91
92 # pre-generated minbias input file
93 set PileUpFile ../eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup
94
95 # average expected pile up
96 set MeanPileUp 200
97
98 # maximum spread in the beam direction in m
99 set ZVertexSpread 0.15
100
101 # maximum spread in time in s
102 set TVertexSpread 1.5E-09
103
104 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
105
106 #set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
107
108
109 set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.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 (abs(t) > 1.0e-09) * (abs(z) > 0.15) * (0.00)}
113
114}
115
116
117#################################
118# Propagate particles in cylinder
119#################################
120
121module ParticlePropagator ParticlePropagator {
122 set InputArray PileUpMerger/stableParticles
123
124 set OutputArray stableParticles
125 set ChargedHadronOutputArray chargedHadrons
126 set ElectronOutputArray electrons
127 set MuonOutputArray muons
128
129 # radius of the magnetic field coverage, in m
130 set Radius 1.29
131 # half-length of the magnetic field coverage, in m
132 set HalfLength 3.0
133
134 # magnetic field
135 set Bz 3.8
136}
137
138
139####################################
140# Charged hadron tracking efficiency
141####################################
142
143module Efficiency ChargedHadronTrackingEfficiency {
144 ## particles after propagation
145 set InputArray ParticlePropagator/chargedHadrons
146 set OutputArray chargedHadrons
147 # tracking efficiency formula for charged hadrons
148 set EfficiencyFormula {
149 (pt <= 0.2) * (0.00) + \
150 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
151 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
152 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
153 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.87) + \
154 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
155 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.82) + \
156 (abs(eta) > 4.0) * (0.00)
157 }
158}
159
160
161#####################################
162# Electron tracking efficiency - ID
163####################################
164
165module Efficiency ElectronTrackingEfficiency {
166 set InputArray ParticlePropagator/electrons
167 set OutputArray electrons
168 # tracking efficiency formula for electrons
169 set EfficiencyFormula {
170 (pt <= 0.2) * (0.00) + \
171 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
172 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
173 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
174 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 10.0) * (0.82+pt*0.01) + \
175 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 10.0) * (0.90) + \
176 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
177 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.8+pt*0.01) + \
178 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.85) + \
179 (abs(eta) > 4.0) * (0.00)
180
181 }
182}
183
184##########################
185# Muon tracking efficiency
186##########################
187
188module Efficiency MuonTrackingEfficiency {
189 set InputArray ParticlePropagator/muons
190 set OutputArray muons
191 # tracking efficiency formula for muons
192 set EfficiencyFormula {
193 (pt <= 0.2) * (0.00) + \
194 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 1.00) + \
195 (abs(eta) <= 1.2) * (pt > 1.0) * (1.00) + \
196 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 0.2 && pt <= 1.0) * (pt*1.00) + \
197 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 1.0) * (1.00) + \
198 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
199 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
200 (abs(eta) > 4.0) * (0.00)
201
202 }
203}
204
205########################################
206# Momentum resolution for charged tracks
207########################################
208
209module MomentumSmearing ChargedHadronMomentumSmearing {
210 ## hadrons after having applied the tracking efficiency
211 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
212 set OutputArray chargedHadrons
213 # resolution formula for charged hadrons ,
214
215 # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
216
217 set ResolutionFormula { (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00467469) + \
218 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004675 + (pt-1.000000)* 0.000056) + \
219 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.005179 + (pt-10.000000)* 0.000064) + \
220 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.010955*pt/100.000000) + \
221 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00515885) + \
222 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005159 + (pt-1.000000)* 0.000049) + \
223 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005600 + (pt-10.000000)* 0.000062) + \
224 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.011149*pt/100.000000) + \
225 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00509775) + \
226 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005098 + (pt-1.000000)* 0.000055) + \
227 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005590 + (pt-10.000000)* 0.000063) + \
228 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.011251*pt/100.000000) + \
229 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00568785) + \
230 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005688 + (pt-1.000000)* 0.000061) + \
231 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.006237 + (pt-10.000000)* 0.000059) + \
232 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.011563*pt/100.000000) + \
233 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00668287) + \
234 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.006683 + (pt-1.000000)* 0.000065) + \
235 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007263 + (pt-10.000000)* 0.000067) + \
236 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.013251*pt/100.000000) + \
237 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01047734) + \
238 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010477 + (pt-1.000000)* 0.000059) + \
239 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.011005 + (pt-10.000000)* 0.000098) + \
240 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.019785*pt/100.000000) + \
241 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01430653) + \
242 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014307 + (pt-1.000000)* 0.000038) + \
243 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.014647 + (pt-10.000000)* 0.000087) + \
244 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.022499*pt/100.000000) + \
245 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01719240) + \
246 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017192 + (pt-1.000000)* 0.000020) + \
247 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.017372 + (pt-10.000000)* 0.000080) + \
248 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.024604*pt/100.000000) + \
249 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01783940) + \
250 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.017839 + (pt-1.000000)* 0.000045) + \
251 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.018244 + (pt-10.000000)* 0.000073) + \
252 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.024839*pt/100.000000) + \
253 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.01787758) + \
254 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.017878 + (pt-1.000000)* 0.000154) + \
255 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.019265 + (pt-10.000000)* 0.000141) + \
256 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.031979*pt/100.000000) + \
257 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01825549) + \
258 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.018255 + (pt-1.000000)* 0.000270) + \
259 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.020690 + (pt-10.000000)* 0.000212) + \
260 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.039811*pt/100.000000) + \
261 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01803308) + \
262 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.018033 + (pt-1.000000)* 0.000220) + \
263 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.020013 + (pt-10.000000)* 0.000247) + \
264 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.042262*pt/100.000000) + \
265 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.02156195) + \
266 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.021562 + (pt-1.000000)* 0.000225) + \
267 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.023591 + (pt-10.000000)* 0.000361) + \
268 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.056075*pt/100.000000) + \
269 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.02691276) + \
270 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.026913 + (pt-1.000000)* 0.000357) + \
271 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.030127 + (pt-10.000000)* 0.000483) + \
272 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.073578*pt/100.000000) + \
273 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.03253153) + \
274 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.032532 + (pt-1.000000)* 0.000382) + \
275 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.035969 + (pt-10.000000)* 0.000673) + \
276 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.096568*pt/100.000000) + \
277 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04551714) + \
278 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045517 + (pt-1.000000)* 0.002188) + \
279 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.065207 + (pt-10.000000)* 0.003254) + \
280 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.358089*pt/100.000000) + \
281 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.04724756) + \
282 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.047248 + (pt-1.000000)* 0.003322) + \
283 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.077142 + (pt-10.000000)* 0.004516) + \
284 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.483608*pt/100.000000) + \
285 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.04493020) + \
286 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.044930 + (pt-1.000000)* 0.002917) + \
287 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.071183 + (pt-10.000000)* 0.004186) + \
288 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.447915*pt/100.000000) + \
289 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.05464199) + \
290 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.054642 + (pt-1.000000)* 0.003456) + \
291 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.085748 + (pt-10.000000)* 0.004995) + \
292 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.535262*pt/100.000000) + \
293 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.06551185) + \
294 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.065512 + (pt-1.000000)* 0.006574) + \
295 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.124674 + (pt-10.000000)* 0.009560) + \
296 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.985040*pt/100.000000)
297}
298
299}
300
301########################################
302# Momentum resolution for electrons
303########################################
304
305module MomentumSmearing ElectronMomentumSmearing {
306 ## hadrons after having applied the tracking efficiency
307 set InputArray ElectronTrackingEfficiency/electrons
308 set OutputArray electrons
309 # resolution formula for electrons
310 # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
311
312 set ResolutionFormula { (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00467469) + \
313 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004675 + (pt-1.000000)* 0.000056) + \
314 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.005179 + (pt-10.000000)* 0.000064) + \
315 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.010955*pt/100.000000) + \
316 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00515885) + \
317 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005159 + (pt-1.000000)* 0.000049) + \
318 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005600 + (pt-10.000000)* 0.000062) + \
319 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.011149*pt/100.000000) + \
320 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00509775) + \
321 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005098 + (pt-1.000000)* 0.000055) + \
322 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005590 + (pt-10.000000)* 0.000063) + \
323 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.011251*pt/100.000000) + \
324 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00568785) + \
325 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005688 + (pt-1.000000)* 0.000061) + \
326 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.006237 + (pt-10.000000)* 0.000059) + \
327 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.011563*pt/100.000000) + \
328 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00668287) + \
329 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.006683 + (pt-1.000000)* 0.000065) + \
330 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007263 + (pt-10.000000)* 0.000067) + \
331 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.013251*pt/100.000000) + \
332 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01047734) + \
333 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010477 + (pt-1.000000)* 0.000059) + \
334 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.011005 + (pt-10.000000)* 0.000098) + \
335 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.019785*pt/100.000000) + \
336 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01430653) + \
337 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014307 + (pt-1.000000)* 0.000038) + \
338 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.014647 + (pt-10.000000)* 0.000087) + \
339 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.022499*pt/100.000000) + \
340 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01719240) + \
341 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017192 + (pt-1.000000)* 0.000020) + \
342 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.017372 + (pt-10.000000)* 0.000080) + \
343 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.024604*pt/100.000000) + \
344 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01783940) + \
345 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.017839 + (pt-1.000000)* 0.000045) + \
346 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.018244 + (pt-10.000000)* 0.000073) + \
347 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.024839*pt/100.000000) + \
348 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.01787758) + \
349 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.017878 + (pt-1.000000)* 0.000154) + \
350 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.019265 + (pt-10.000000)* 0.000141) + \
351 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.031979*pt/100.000000) + \
352 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01825549) + \
353 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.018255 + (pt-1.000000)* 0.000270) + \
354 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.020690 + (pt-10.000000)* 0.000212) + \
355 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.039811*pt/100.000000) + \
356 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01803308) + \
357 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.018033 + (pt-1.000000)* 0.000220) + \
358 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.020013 + (pt-10.000000)* 0.000247) + \
359 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.042262*pt/100.000000) + \
360 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.02156195) + \
361 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.021562 + (pt-1.000000)* 0.000225) + \
362 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.023591 + (pt-10.000000)* 0.000361) + \
363 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.056075*pt/100.000000) + \
364 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.02691276) + \
365 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.026913 + (pt-1.000000)* 0.000357) + \
366 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.030127 + (pt-10.000000)* 0.000483) + \
367 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.073578*pt/100.000000) + \
368 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.03253153) + \
369 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.032532 + (pt-1.000000)* 0.000382) + \
370 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.035969 + (pt-10.000000)* 0.000673) + \
371 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.096568*pt/100.000000) + \
372 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04551714) + \
373 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045517 + (pt-1.000000)* 0.002188) + \
374 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.065207 + (pt-10.000000)* 0.003254) + \
375 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.358089*pt/100.000000) + \
376 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.04724756) + \
377 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.047248 + (pt-1.000000)* 0.003322) + \
378 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.077142 + (pt-10.000000)* 0.004516) + \
379 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.483608*pt/100.000000) + \
380 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.04493020) + \
381 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.044930 + (pt-1.000000)* 0.002917) + \
382 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.071183 + (pt-10.000000)* 0.004186) + \
383 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.447915*pt/100.000000) + \
384 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.05464199) + \
385 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.054642 + (pt-1.000000)* 0.003456) + \
386 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.085748 + (pt-10.000000)* 0.004995) + \
387 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.535262*pt/100.000000) + \
388 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.06551185) + \
389 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.065512 + (pt-1.000000)* 0.006574) + \
390 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.124674 + (pt-10.000000)* 0.009560) + \
391 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.985040*pt/100.000000)
392}
393
394}
395
396
397###############################
398# Momentum resolution for muons
399###############################
400
401module MomentumSmearing MuonMomentumSmearing {
402 set InputArray MuonTrackingEfficiency/muons
403 set OutputArray muons
404 # resolution formula for muons
405
406 # up to |eta| < 2.8 take measurement from tracking + muon chambers
407 # for |eta| > 2.8 and pT < 5.0 take measurement from tracking alone taken from
408 # http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
409
410 set ResolutionFormula {
411
412 (pt > 5 && pt < 10) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.00519 + \
413 (pt > 5 && pt < 10) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.00575 + \
414 (pt > 5 && pt < 10) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.00629 + \
415 (pt > 5 && pt < 10) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.00626 + \
416 (pt > 5 && pt < 10) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.00722 + \
417 (pt > 5 && pt < 10) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.00968 + \
418 (pt > 5 && pt < 10) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.00983 + \
419 (pt > 5 && pt < 10) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.0102 + \
420 (pt > 5 && pt < 10) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.0126 + \
421 (pt > 5 && pt < 10) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.0178 + \
422 (pt > 5 && pt < 10) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.0248 + \
423 (pt > 5 && pt < 10) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.0347 + \
424 (pt > 5 && pt < 10) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.0359 + \
425 (pt > 5 && pt < 10) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.0412 + \
426 (pt > 5 && pt < 10) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.0422 + \
427 (pt > 5 && pt < 10) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.0496 + \
428
429 (pt > 10 && pt < 20) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.00533 + \
430 (pt > 10 && pt < 20) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.00592 + \
431 (pt > 10 && pt < 20) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.00602 + \
432 (pt > 10 && pt < 20) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.00614 + \
433 (pt > 10 && pt < 20) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.00736 + \
434 (pt > 10 && pt < 20) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.0104 + \
435 (pt > 10 && pt < 20) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.00968 + \
436 (pt > 10 && pt < 20) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.00991 + \
437 (pt > 10 && pt < 20) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.0106 + \
438 (pt > 10 && pt < 20) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.0158 + \
439 (pt > 10 && pt < 20) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.0207 + \
440 (pt > 10 && pt < 20) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.0327 + \
441 (pt > 10 && pt < 20) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.0369 + \
442 (pt > 10 && pt < 20) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.039 + \
443 (pt > 10 && pt < 20) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.0488 + \
444 (pt > 10 && pt < 20) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.0543 + \
445
446 (pt > 20 && pt < 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.00559 + \
447 (pt > 20 && pt < 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.00586 + \
448 (pt > 20 && pt < 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.00614 + \
449 (pt > 20 && pt < 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.00654 + \
450 (pt > 20 && pt < 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.00759 + \
451 (pt > 20 && pt < 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.0103 + \
452 (pt > 20 && pt < 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.00963 + \
453 (pt > 20 && pt < 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.0102 + \
454 (pt > 20 && pt < 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.011 + \
455 (pt > 20 && pt < 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.0142 + \
456 (pt > 20 && pt < 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.0197 + \
457 (pt > 20 && pt < 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.0315 + \
458 (pt > 20 && pt < 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.0409 + \
459 (pt > 20 && pt < 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.0471 + \
460 (pt > 20 && pt < 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.0579 + \
461 (pt > 20 && pt < 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.0679 + \
462
463 (pt > 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.00656 + \
464 (pt > 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.007 + \
465 (pt > 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.00733 + \
466 (pt > 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.00752 + \
467 (pt > 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.00851 + \
468 (pt > 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.0105 + \
469 (pt > 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.0105 + \
470 (pt > 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.0111 + \
471 (pt > 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.0119 + \
472 (pt > 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.0165 + \
473 (pt > 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.0242 + \
474 (pt > 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.0359 + \
475 (pt > 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.0474 + \
476 (pt > 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.052 + \
477 (pt > 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.0639 + \
478 (pt > 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.0763 + \
479
480 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.03253153) + \
481 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.032532 + (pt-1.000000)* 0.000382) + \
482 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.035969 + (pt-10.000000)* 0.000673) + \
483 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.096568*pt/100.000000) + \
484 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04551714) + \
485 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045517 + (pt-1.000000)* 0.002188) + \
486 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.065207 + (pt-10.000000)* 0.003254) + \
487 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.358089*pt/100.000000) + \
488 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.04724756) + \
489 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.047248 + (pt-1.000000)* 0.003322) + \
490 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.077142 + (pt-10.000000)* 0.004516) + \
491 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.483608*pt/100.000000) + \
492 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.04493020) + \
493 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.044930 + (pt-1.000000)* 0.002917) + \
494 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.071183 + (pt-10.000000)* 0.004186) + \
495 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.447915*pt/100.000000) + \
496 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.05464199) + \
497 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.054642 + (pt-1.000000)* 0.003456) + \
498 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.085748 + (pt-10.000000)* 0.004995) + \
499 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.535262*pt/100.000000) + \
500 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.06551185) + \
501 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.065512 + (pt-1.000000)* 0.006574) + \
502 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.124674 + (pt-10.000000)* 0.009560) + \
503 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.985040*pt/100.000000) + \
504
505 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00467469) + \
506 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 5.000) * (0.004675 + (pt-1.000000)* 0.000056) + \
507 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00515885) + \
508 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 5.000) * (0.005159 + (pt-1.000000)* 0.000049) + \
509 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00509775) + \
510 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 5.000) * (0.005098 + (pt-1.000000)* 0.000055) + \
511 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00568785) + \
512 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 5.000) * (0.005688 + (pt-1.000000)* 0.000061) + \
513 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00668287) + \
514 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 5.000) * (0.006683 + (pt-1.000000)* 0.000065) + \
515 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01047734) + \
516 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 5.000) * (0.010477 + (pt-1.000000)* 0.000059) + \
517 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01430653) + \
518 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 5.000) * (0.014307 + (pt-1.000000)* 0.000038) + \
519 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01719240) + \
520 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 5.000) * (0.017192 + (pt-1.000000)* 0.000020) + \
521 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01783940) + \
522 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 5.000) * (0.017839 + (pt-1.000000)* 0.000045) + \
523 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.01787758) + \
524 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 5.000) * (0.017878 + (pt-1.000000)* 0.000154) + \
525 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01825549) + \
526 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 5.000) * (0.018255 + (pt-1.000000)* 0.000270) + \
527 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01803308) + \
528 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 5.000) * (0.018033 + (pt-1.000000)* 0.000220) + \
529 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.02156195) + \
530 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 5.000) * (0.021562 + (pt-1.000000)* 0.000225) + \
531 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.02691276) + \
532 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 5.000) * (0.026913 + (pt-1.000000)* 0.000357) + \
533 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.03253153) + \
534 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 5.000) * (0.032532 + (pt-1.000000)* 0.000382) + \
535 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04551714) + \
536 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 5.000) * (0.045517 + (pt-1.000000)* 0.002188) + \
537 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.04724756) + \
538 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 5.000) * (0.047248 + (pt-1.000000)* 0.003322) + \
539 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.04493020) + \
540 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 5.000) * (0.044930 + (pt-1.000000)* 0.002917) + \
541 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.05464199) + \
542 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 5.000) * (0.054642 + (pt-1.000000)* 0.003456) + \
543 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.06551185) + \
544 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 5.000) * (0.065512 + (pt-1.000000)* 0.006574)
545
546 }
547}
548
549
550
551##############
552# Track merger
553##############
554
555module Merger TrackMerger {
556# add InputArray InputArray
557 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
558 add InputArray ElectronMomentumSmearing/electrons
559 add InputArray MuonMomentumSmearing/muons
560 set OutputArray tracks
561}
562
563#############
564# ECAL
565#############
566
567module SimpleCalorimeter ECal {
568 set ParticleInputArray ParticlePropagator/stableParticles
569 set TrackInputArray TrackMerger/tracks
570
571 set TowerOutputArray ecalTowers
572 set EFlowTrackOutputArray eflowTracks
573 set EFlowTowerOutputArray eflowPhotons
574
575 set IsEcal true
576
577 set EnergyMin 0.5
578 set EnergySignificanceMin 1.0
579
580 set SmearTowerCenter true
581
582 set pi [expr {acos(-1)}]
583
584 # lists of the edges of each tower in eta and phi
585 # each list starts with the lower edge of the first tower
586 # the list ends with the higher edged of the last tower
587
588 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
589
590 set PhiBins {}
591 for {set i -180} {$i <= 180} {incr i} {
592 add PhiBins [expr {$i * $pi/180.0}]
593 }
594
595 # 0.02 unit in eta up to eta = 1.5 (barrel)
596 for {set i -85} {$i <= 86} {incr i} {
597 set eta [expr {$i * 0.0174}]
598 add EtaPhiBins $eta $PhiBins
599 }
600
601 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
602
603 set PhiBins {}
604 for {set i -180} {$i <= 180} {incr i} {
605 add PhiBins [expr {$i * $pi/180.0}]
606 }
607
608 # 0.02 unit in eta up to eta = 3
609 for {set i 1} {$i <= 84} {incr i} {
610 set eta [expr { -2.958 + $i * 0.0174}]
611 add EtaPhiBins $eta $PhiBins
612 }
613
614 for {set i 1} {$i <= 84} {incr i} {
615 set eta [expr { 1.4964 + $i * 0.0174}]
616 add EtaPhiBins $eta $PhiBins
617 }
618
619 # take present CMS granularity for HF
620
621 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
622 set PhiBins {}
623 for {set i -18} {$i <= 18} {incr i} {
624 add PhiBins [expr {$i * $pi/18.0}]
625 }
626
627 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} {
628 add EtaPhiBins $eta $PhiBins
629 }
630
631
632 add EnergyFraction {0} {0.0}
633 # energy fractions for e, gamma and pi0
634 add EnergyFraction {11} {1.0}
635 add EnergyFraction {22} {1.0}
636 add EnergyFraction {111} {1.0}
637 # energy fractions for muon, neutrinos and neutralinos
638 add EnergyFraction {12} {0.0}
639 add EnergyFraction {13} {0.0}
640 add EnergyFraction {14} {0.0}
641 add EnergyFraction {16} {0.0}
642 add EnergyFraction {1000022} {0.0}
643 add EnergyFraction {1000023} {0.0}
644 add EnergyFraction {1000025} {0.0}
645 add EnergyFraction {1000035} {0.0}
646 add EnergyFraction {1000045} {0.0}
647 # energy fractions for K0short and Lambda
648 add EnergyFraction {310} {0.3}
649 add EnergyFraction {3122} {0.3}
650
651 # set ResolutionFormula {resolution formula as a function of eta and energy}
652
653 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016.
654 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117
655
656 set ResolutionFormula { (abs(eta) <= 1.50) * sqrt(energy^2*0.005^2 + energy*0.027^2 + 0.15^2) + \
657 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
658 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
659 (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
660 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
661
662}
663
664#############
665# HCAL
666#############
667
668module SimpleCalorimeter HCal {
669 set ParticleInputArray ParticlePropagator/stableParticles
670 set TrackInputArray ECal/eflowTracks
671
672 set TowerOutputArray hcalTowers
673 set EFlowTrackOutputArray eflowTracks
674 set EFlowTowerOutputArray eflowNeutralHadrons
675
676 set IsEcal false
677
678 set EnergyMin 1.0
679 set EnergySignificanceMin 1.0
680
681 set SmearTowerCenter true
682
683 set pi [expr {acos(-1)}]
684
685 # lists of the edges of each tower in eta and phi
686 # each list starts with the lower edge of the first tower
687 # the list ends with the higher edged of the last tower
688
689 # assume 0.087 x 0.087 resolution in eta,phi in the barrel |eta| < 1.5
690
691 set PhiBins {}
692 for {set i -36} {$i <= 36} {incr i} {
693 add PhiBins [expr {$i * $pi/36.0}]
694 }
695 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} {
696 add EtaPhiBins $eta $PhiBins
697 }
698
699 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- HCAL)
700
701 set PhiBins {}
702 for {set i -180} {$i <= 180} {incr i} {
703 add PhiBins [expr {$i * $pi/180.0}]
704 }
705
706 # 0.02 unit in eta up to eta = 3
707 for {set i 1} {$i <= 84} {incr i} {
708 set eta [expr { -2.958 + $i * 0.0174}]
709 add EtaPhiBins $eta $PhiBins
710 }
711
712 for {set i 1} {$i <= 84} {incr i} {
713 set eta [expr { 1.4964 + $i * 0.0174}]
714 add EtaPhiBins $eta $PhiBins
715 }
716
717 # take present CMS granularity for HF
718
719 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
720 set PhiBins {}
721 for {set i -18} {$i <= 18} {incr i} {
722 add PhiBins [expr {$i * $pi/18.0}]
723 }
724
725 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} {
726 add EtaPhiBins $eta $PhiBins
727 }
728
729
730 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
731 add EnergyFraction {0} {1.0}
732 # energy fractions for e, gamma and pi0
733 add EnergyFraction {11} {0.0}
734 add EnergyFraction {22} {0.0}
735 add EnergyFraction {111} {0.0}
736 # energy fractions for muon, neutrinos and neutralinos
737 add EnergyFraction {12} {0.0}
738 add EnergyFraction {13} {0.0}
739 add EnergyFraction {14} {0.0}
740 add EnergyFraction {16} {0.0}
741 add EnergyFraction {1000022} {0.0}
742 add EnergyFraction {1000023} {0.0}
743 add EnergyFraction {1000025} {0.0}
744 add EnergyFraction {1000035} {0.0}
745 add EnergyFraction {1000045} {0.0}
746 # energy fractions for K0short and Lambda
747 add EnergyFraction {310} {0.7}
748 add EnergyFraction {3122} {0.7}
749
750# set ResolutionFormula {resolution formula as a function of eta and energy}
751 set ResolutionFormula { (abs(eta) <= 1.5) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
752 (abs(eta) > 1.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
753 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.11^2 + energy*2.80^2)}
754
755}
756
757#################
758# Electron filter
759#################
760
761module PdgCodeFilter ElectronFilter {
762 set InputArray HCal/eflowTracks
763 set OutputArray electrons
764 set Invert true
765 add PdgCode {11}
766 add PdgCode {-11}
767}
768
769##########################
770# Track pile-up subtractor
771##########################
772
773module TrackPileUpSubtractor TrackPileUpSubtractor {
774# add InputArray InputArray OutputArray
775 add InputArray HCal/eflowTracks eflowTracks
776 add InputArray ElectronFilter/electrons electrons
777 add InputArray MuonMomentumSmearing/muons muons
778
779 set VertexInputArray PileUpMerger/vertices
780 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
781 # Z vertex resolution in m
782 set ZVertexResolution 0.0001
783}
784
785
786###################################################
787# Tower Merger (in case not using e-flow algorithm)
788###################################################
789
790module Merger TowerMerger {
791# add InputArray InputArray
792 add InputArray ECal/ecalTowers
793 add InputArray HCal/hcalTowers
794 set OutputArray towers
795}
796
797
798####################
799# Neutral eflow erger
800####################
801
802module Merger NeutralEFlowMerger {
803# add InputArray InputArray
804 add InputArray ECal/eflowPhotons
805 add InputArray HCal/eflowNeutralHadrons
806 set OutputArray eflowTowers
807}
808
809
810####################
811# Energy flow merger
812####################
813
814module Merger EFlowMerger {
815# add InputArray InputArray
816 add InputArray HCal/eflowTracks
817 add InputArray ECal/eflowPhotons
818 add InputArray HCal/eflowNeutralHadrons
819 set OutputArray eflow
820}
821
822##################################
823# Energy flow merger (all tracks)
824##################################
825
826module Merger EFlowMergerAllTracks {
827# add InputArray InputArray
828 add InputArray TrackMerger/tracks
829 add InputArray ECal/eflowPhotons
830 add InputArray HCal/eflowNeutralHadrons
831 set OutputArray eflow
832}
833
834#########################################
835### Run the puppi code (to be tuned) ###
836#########################################
837
838module RunPUPPI RunPUPPI {
839 ## input information
840 set TrackInputArray TrackMerger/tracks
841 set NeutralInputArray NeutralEFlowMerger/eflowTowers
842 set PVInputArray PileUpMerger/vertices
843 set MinPuppiWeight 0.05
844 set UseExp false
845
846 ## define puppi algorithm parameters (more than one for the same eta region is possible)
847 add EtaMinBin 0. 2.5 2.5 3.0 3.0
848 add EtaMaxBin 2.5 3.0 3.0 10.0 10.0
849 add PtMinBin 0. 0.5 0.5 0.5 0.5
850 add ConeSizeBin 0.25 0.25 0.25 0.25 0.25
851 add RMSPtMinBin 0.1 0.5 0.5 0.5 0.5
852 add RMSScaleFactorBin 1.0 1.0 1.0 1.0 1.0
853 add NeutralMinEBin 0.2 1.0 1.0 1.5 1.5
854 add NeutralPtSlope 0.02 0.02 0.02 0.02 0.02
855 add ApplyCHS true true true true true
856 add UseCharged true false false false false
857 add ApplyLowPUCorr true true true true true
858 add MetricId 5 5 1 5 1
859
860 ## output name
861 set OutputArray PuppiParticles
862 set OutputArrayTracks puppiTracks
863 set OutputArrayNeutrals puppiNeutrals
864}
865
866
867
868###################
869# Missing ET merger
870###################
871
872module Merger MissingET {
873# add InputArray InputArray
874# add InputArray RunPUPPI/PuppiParticles
875 add InputArray EFlowMerger/eflow
876 set MomentumOutputArray momentum
877}
878
879###################
880# Ger PileUp Missing ET
881###################
882
883module Merger GenPileUpMissingET {
884# add InputArray InputArray
885# add InputArray RunPUPPI/PuppiParticles
886 add InputArray ParticlePropagator/stableParticles
887 set MomentumOutputArray momentum
888}
889
890##################
891# Scalar HT merger
892##################
893
894module Merger ScalarHT {
895# add InputArray InputArray
896 add InputArray RunPUPPI/PuppiParticles
897 set EnergyOutputArray energy
898}
899
900#################
901# Neutrino Filter
902#################
903
904module PdgCodeFilter NeutrinoFilter {
905
906 set InputArray Delphes/stableParticles
907 set OutputArray filteredParticles
908
909 set PTMin 0.0
910
911 add PdgCode {12}
912 add PdgCode {14}
913 add PdgCode {16}
914 add PdgCode {-12}
915 add PdgCode {-14}
916 add PdgCode {-16}
917
918}
919
920
921#####################
922# MC truth jet finder
923#####################
924
925module FastJetFinder GenJetFinder {
926 set InputArray NeutrinoFilter/filteredParticles
927
928 set OutputArray jets
929
930 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
931 set JetAlgorithm 6
932 set ParameterR 0.4
933
934 set JetPTMin 15.0
935}
936
937#########################
938# Gen Missing ET merger
939########################
940
941module Merger GenMissingET {
942
943# add InputArray InputArray
944 add InputArray NeutrinoFilter/filteredParticles
945 set MomentumOutputArray momentum
946}
947
948
949
950############
951# Jet finder
952############
953
954module FastJetFinder FastJetFinder {
955# set InputArray TowerMerger/towers
956 set InputArray RunPUPPI/PuppiParticles
957
958 set OutputArray jets
959
960 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
961 set JetAlgorithm 6
962 set ParameterR 0.4
963
964 set JetPTMin 15.0
965}
966
967##################
968# Jet Energy Scale
969##################
970
971module EnergyScale JetEnergyScale {
972 set InputArray FastJetFinder/jets
973 set OutputArray jets
974
975 # scale formula for jets
976 set ScaleFormula {1.00}
977}
978
979
980####################
981# Photon isolation #
982####################
983
984module Isolation PhotonIsolation {
985
986 # particle for which calculate the isolation
987 set CandidateInputArray ECal/eflowPhotons
988
989 # isolation collection
990 set IsolationInputArray EFlowMerger/eflow
991
992 # output array
993 set OutputArray photons
994
995 # isolation cone
996 set DeltaRMax 0.3
997
998 # minimum pT
999 set PTMin 1.0
1000
1001 # iso ratio to cut
1002 set PTRatioMax 9999.
1003}
1004
1005
1006
1007#####################
1008# Photon efficiency #
1009#####################
1010
1011module Efficiency PhotonEfficiency {
1012
1013 ## input particles
1014 set InputArray PhotonIsolation/photons
1015 ## output particles
1016 set OutputArray photons
1017 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1018 # efficiency formula for photons
1019 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
1020 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
1021 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
1022 (abs(eta) > 4.0) * (0.00)}
1023
1024}
1025
1026
1027######################
1028# Electron isolation #
1029######################
1030
1031module Isolation ElectronIsolation {
1032
1033 set CandidateInputArray ElectronFilter/electrons
1034
1035 # isolation collection
1036 set IsolationInputArray EFlowMerger/eflow
1037
1038 set OutputArray electrons
1039
1040 set DeltaRMax 0.3
1041 set PTMin 1.0
1042 set PTRatioMax 9999.
1043
1044}
1045
1046
1047
1048#######################
1049# Electron efficiency #
1050#######################
1051
1052module Efficiency ElectronEfficiency {
1053
1054 set InputArray ElectronIsolation/electrons
1055 set OutputArray electrons
1056
1057 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1058 # efficiency formula for electrons
1059 set EfficiencyFormula {
1060 (pt <= 4.0) * (0.00) + \
1061 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
1062 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
1063 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
1064 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
1065 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
1066 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
1067 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
1068 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
1069 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
1070 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
1071 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
1072 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
1073 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
1074 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
1075 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
1076 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
1077 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1078 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1079 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1080 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1081 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
1082 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1083 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1084 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1085 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1086 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
1087 (abs(eta) > 4.0) * (0.00)
1088
1089 }
1090}
1091
1092##################
1093# Muon isolation #
1094##################
1095
1096module Isolation MuonIsolation {
1097 set CandidateInputArray MuonMomentumSmearing/muons
1098
1099 # isolation collection
1100 set IsolationInputArray EFlowMerger/eflow
1101
1102 set OutputArray muons
1103
1104 set DeltaRMax 0.3
1105 set PTMin 1.0
1106 set PTRatioMax 9999.
1107
1108}
1109
1110
1111##################
1112# Muon Loose Id #
1113##################
1114
1115module Efficiency MuonLooseIdEfficiency {
1116 set InputArray MuonIsolation/muons
1117 set OutputArray muons
1118 # tracking + LooseID efficiency formula for muons
1119 set EfficiencyFormula {
1120
1121 (pt <= 5.0) * (0.00) + \
1122 (pt > 5 && pt < 10) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.967 + \
1123 (pt > 5 && pt < 10) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.968 + \
1124 (pt > 5 && pt < 10) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.986 + \
1125 (pt > 5 && pt < 10) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.982 + \
1126 (pt > 5 && pt < 10) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.972 + \
1127 (pt > 5 && pt < 10) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.974 + \
1128 (pt > 5 && pt < 10) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.969 + \
1129 (pt > 5 && pt < 10) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.988 + \
1130 (pt > 5 && pt < 10) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.989 + \
1131 (pt > 5 && pt < 10) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.992 + \
1132 (pt > 5 && pt < 10) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.979 + \
1133 (pt > 5 && pt < 10) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.976 + \
1134 (pt > 5 && pt < 10) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.882 + \
1135 (pt > 5 && pt < 10) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.883 + \
1136 (pt > 5 && pt < 10) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.851 + \
1137 (pt > 5 && pt < 10) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.828 + \
1138
1139 (pt > 10 && pt < 20) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.989 + \
1140 (pt > 10 && pt < 20) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.989 + \
1141 (pt > 10 && pt < 20) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.995 + \
1142 (pt > 10 && pt < 20) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.995 + \
1143 (pt > 10 && pt < 20) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.987 + \
1144 (pt > 10 && pt < 20) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.989 + \
1145 (pt > 10 && pt < 20) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.981 + \
1146 (pt > 10 && pt < 20) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.99 + \
1147 (pt > 10 && pt < 20) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.997 + \
1148 (pt > 10 && pt < 20) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.997 + \
1149 (pt > 10 && pt < 20) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.989 + \
1150 (pt > 10 && pt < 20) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.995 + \
1151 (pt > 10 && pt < 20) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.888 + \
1152 (pt > 10 && pt < 20) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.894 + \
1153 (pt > 10 && pt < 20) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.886 + \
1154 (pt > 10 && pt < 20) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.883 + \
1155
1156 (pt > 20 && pt < 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.981 + \
1157 (pt > 20 && pt < 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.994 + \
1158 (pt > 20 && pt < 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.995 + \
1159 (pt > 20 && pt < 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.997 + \
1160 (pt > 20 && pt < 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.986 + \
1161 (pt > 20 && pt < 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.991 + \
1162 (pt > 20 && pt < 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.987 + \
1163 (pt > 20 && pt < 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.992 + \
1164 (pt > 20 && pt < 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.993 + \
1165 (pt > 20 && pt < 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.999 + \
1166 (pt > 20 && pt < 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.994 + \
1167 (pt > 20 && pt < 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.997 + \
1168 (pt > 20 && pt < 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.894 + \
1169 (pt > 20 && pt < 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.879 + \
1170 (pt > 20 && pt < 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.871 + \
1171 (pt > 20 && pt < 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.864 + \
1172
1173 (pt > 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.986 + \
1174 (pt > 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.989 + \
1175 (pt > 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.996 + \
1176 (pt > 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.994 + \
1177 (pt > 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.991 + \
1178 (pt > 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.989 + \
1179 (pt > 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.985 + \
1180 (pt > 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.993 + \
1181 (pt > 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.996 + \
1182 (pt > 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.997 + \
1183 (pt > 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.99 + \
1184 (pt > 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.993 + \
1185 (pt > 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.905 + \
1186 (pt > 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.86 + \
1187 (pt > 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.845 + \
1188 (pt > 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.855 + \
1189
1190 (abs(eta) > 2.8) * (0.00)
1191
1192 }
1193}
1194
1195
1196##################
1197# Muon Tight Id #
1198##################
1199
1200module Efficiency MuonTightIdEfficiency {
1201 set InputArray MuonIsolation/muons
1202 set OutputArray muons
1203 # tracking + TightID efficiency formula for muons
1204 set EfficiencyFormula {
1205
1206 (pt <= 5.0) * (0.00) + \
1207 (pt > 5 && pt < 10) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.949 + \
1208 (pt > 5 && pt < 10) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.894 + \
1209 (pt > 5 && pt < 10) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.966 + \
1210 (pt > 5 && pt < 10) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.964 + \
1211 (pt > 5 && pt < 10) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.923 + \
1212 (pt > 5 && pt < 10) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.94 + \
1213 (pt > 5 && pt < 10) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.959 + \
1214 (pt > 5 && pt < 10) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.972 + \
1215 (pt > 5 && pt < 10) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.962 + \
1216 (pt > 5 && pt < 10) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.972 + \
1217 (pt > 5 && pt < 10) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.947 + \
1218 (pt > 5 && pt < 10) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.911 + \
1219 (pt > 5 && pt < 10) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.825 + \
1220 (pt > 5 && pt < 10) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.821 + \
1221 (pt > 5 && pt < 10) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.783 + \
1222 (pt > 5 && pt < 10) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.746 + \
1223
1224 (pt > 10 && pt < 20) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.975 + \
1225 (pt > 10 && pt < 20) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.935 + \
1226 (pt > 10 && pt < 20) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.986 + \
1227 (pt > 10 && pt < 20) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.976 + \
1228 (pt > 10 && pt < 20) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.954 + \
1229 (pt > 10 && pt < 20) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.978 + \
1230 (pt > 10 && pt < 20) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.975 + \
1231 (pt > 10 && pt < 20) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.98 + \
1232 (pt > 10 && pt < 20) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.982 + \
1233 (pt > 10 && pt < 20) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.985 + \
1234 (pt > 10 && pt < 20) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.962 + \
1235 (pt > 10 && pt < 20) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.946 + \
1236 (pt > 10 && pt < 20) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.839 + \
1237 (pt > 10 && pt < 20) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.844 + \
1238 (pt > 10 && pt < 20) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.803 + \
1239 (pt > 10 && pt < 20) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.791 + \
1240
1241 (pt > 20 && pt < 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.965 + \
1242 (pt > 20 && pt < 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.938 + \
1243 (pt > 20 && pt < 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.99 + \
1244 (pt > 20 && pt < 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.977 + \
1245 (pt > 20 && pt < 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.959 + \
1246 (pt > 20 && pt < 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.976 + \
1247 (pt > 20 && pt < 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.984 + \
1248 (pt > 20 && pt < 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.981 + \
1249 (pt > 20 && pt < 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.978 + \
1250 (pt > 20 && pt < 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.985 + \
1251 (pt > 20 && pt < 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.969 + \
1252 (pt > 20 && pt < 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.964 + \
1253 (pt > 20 && pt < 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.863 + \
1254 (pt > 20 && pt < 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.828 + \
1255 (pt > 20 && pt < 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.807 + \
1256 (pt > 20 && pt < 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.764 + \
1257
1258 (pt > 40) * (abs(eta) > 0 && abs(eta) < 0.2) * 0.975 + \
1259 (pt > 40) * (abs(eta) > 0.2 && abs(eta) < 0.4) * 0.935 + \
1260 (pt > 40) * (abs(eta) > 0.4 && abs(eta) < 0.6) * 0.991 + \
1261 (pt > 40) * (abs(eta) > 0.6 && abs(eta) < 0.8) * 0.983 + \
1262 (pt > 40) * (abs(eta) > 0.8 && abs(eta) < 1) * 0.957 + \
1263 (pt > 40) * (abs(eta) > 1 && abs(eta) < 1.2) * 0.979 + \
1264 (pt > 40) * (abs(eta) > 1.2 && abs(eta) < 1.4) * 0.981 + \
1265 (pt > 40) * (abs(eta) > 1.4 && abs(eta) < 1.6) * 0.983 + \
1266 (pt > 40) * (abs(eta) > 1.6 && abs(eta) < 1.8) * 0.981 + \
1267 (pt > 40) * (abs(eta) > 1.8 && abs(eta) < 2) * 0.984 + \
1268 (pt > 40) * (abs(eta) > 2 && abs(eta) < 2.2) * 0.964 + \
1269 (pt > 40) * (abs(eta) > 2.2 && abs(eta) < 2.4) * 0.956 + \
1270 (pt > 40) * (abs(eta) > 2.4 && abs(eta) < 2.5) * 0.855 + \
1271 (pt > 40) * (abs(eta) > 2.5 && abs(eta) < 2.6) * 0.817 + \
1272 (pt > 40) * (abs(eta) > 2.6 && abs(eta) < 2.7) * 0.764 + \
1273 (pt > 40) * (abs(eta) > 2.7 && abs(eta) < 2.8) * 0.752 + \
1274
1275 (abs(eta) > 2.8) * (0.00)
1276
1277 }
1278}
1279
1280
1281########################
1282# Jet Flavor Association
1283########################
1284
1285module JetFlavorAssociation JetFlavorAssociation {
1286
1287 set PartonInputArray Delphes/partons
1288 set ParticleInputArray Delphes/allParticles
1289 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1290 set JetInputArray JetEnergyScale/jets
1291
1292 set DeltaR 0.5
1293 set PartonPTMin 10.0
1294 set PartonEtaMax 4.0
1295
1296}
1297
1298
1299#############
1300# b-tagging #
1301#############
1302module BTagging BTaggingLoose {
1303
1304 set JetInputArray JetEnergyScale/jets
1305
1306 set BitNumber 0
1307
1308 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
1309 (pt > 20.0) * (abs(eta) <= 3.4) * (0.1) + \
1310 (abs(eta) > 3.4) * (0.000)
1311 }
1312
1313
1314 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1315 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.341) + \
1316 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.389) + \
1317 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.416) + \
1318 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.432) + \
1319 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.447) + \
1320 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.451) + \
1321 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.451) + \
1322 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.447) + \
1323 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.447) + \
1324 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.444) + \
1325 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.430) + \
1326 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.424) + \
1327 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.419) + \
1328 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.407) + \
1329 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.384) + \
1330 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.364) + \
1331 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.353) + \
1332 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.303) + \
1333 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.273) + \
1334 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.249) + \
1335 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.233) + \
1336 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.209) + \
1337 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.194) + \
1338 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.227) + \
1339 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1340
1341 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1342 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.236) + \
1343 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.276) + \
1344 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.302) + \
1345 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.330) + \
1346 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.338) + \
1347 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.355) + \
1348 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.353) + \
1349 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.358) + \
1350 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.350) + \
1351 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.327) + \
1352 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.313) + \
1353 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.307) + \
1354 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.322) + \
1355 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.299) + \
1356 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.274) + \
1357 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.240) + \
1358 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.203) + \
1359 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.218) + \
1360 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.203) + \
1361 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.178) + \
1362 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1363
1364 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1365 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.179) + \
1366 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.198) + \
1367 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.219) + \
1368 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.225) + \
1369 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.232) + \
1370 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.251) + \
1371 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.267) + \
1372 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.252) + \
1373 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.243) + \
1374 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.244) + \
1375 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.242) + \
1376 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.246) + \
1377 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.217) + \
1378 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.223)+ \
1379 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.194) + \
1380 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.187) + \
1381 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1382
1383 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1384 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.09) + \
1385 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.10) + \
1386 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.11) + \
1387 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.11) + \
1388 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.12) + \
1389 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.13) + \
1390 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.13) + \
1391 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.13) + \
1392 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.12) + \
1393 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.12) + \
1394 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.12) + \
1395 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.12) + \
1396 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.11) + \
1397 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.11) + \
1398 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.10) + \
1399 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.09) + \
1400 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1401 (abs(eta) > 3.4) * (0.000)
1402 }
1403
1404 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1405 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.712) + \
1406 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.788) + \
1407 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.837) + \
1408 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.856) + \
1409 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.873) + \
1410 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.886) + \
1411 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.883) + \
1412 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.882) + \
1413 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.883) + \
1414 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.883) + \
1415 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.876) + \
1416 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.872) + \
1417 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.867) + \
1418 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.853) + \
1419 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.831) + \
1420 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.819) + \
1421 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.791) + \
1422 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.740) + \
1423 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.699) + \
1424 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.631) + \
1425 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.600) + \
1426 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.567) + \
1427 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.478) + \
1428 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.470) + \
1429 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1430
1431 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1432 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.543) + \
1433 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.627) + \
1434 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.701) + \
1435 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.741) + \
1436 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.754) + \
1437 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.796) + \
1438 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.778) + \
1439 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.798) + \
1440 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.788) + \
1441 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.775) + \
1442 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.765) + \
1443 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.752) + \
1444 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.740) + \
1445 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.716) + \
1446 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.671) + \
1447 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.637) + \
1448 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.577) + \
1449 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.572) + \
1450 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.514) + \
1451 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.516) + \
1452 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1453
1454 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1455 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.407) + \
1456 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.479) + \
1457 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.529) + \
1458 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.578) + \
1459 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.612) + \
1460 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.673) + \
1461 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.656) + \
1462 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.658) + \
1463 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.647) + \
1464 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.629) + \
1465 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.630) + \
1466 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.615) + \
1467 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.604) + \
1468 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.532)+ \
1469 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.532) + \
1470 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.503) + \
1471 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1472
1473 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1474 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.20) + \
1475 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.235) + \
1476 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.26) + \
1477 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.28) + \
1478 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.31) + \
1479 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.33) + \
1480 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.32) + \
1481 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.32) + \
1482 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.33) + \
1483 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.31) + \
1484 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.32) + \
1485 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.31) + \
1486 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.30) + \
1487 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.26) + \
1488 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.27) + \
1489 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.25) + \
1490 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1491 (abs(eta) > 3.4) * (0.000)
1492 }
1493}
1494
1495
1496#############
1497# b-tagging #
1498#############
1499module BTagging BTaggingMedium {
1500
1501 set JetInputArray JetEnergyScale/jets
1502
1503 set BitNumber 1
1504
1505 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
1506 (pt > 20.0) * (abs(eta) <= 3.4) * (0.01) + \
1507 (abs(eta) > 3.4) * (0.000)
1508 }
1509
1510 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1511 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.144) + \
1512 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.173) + \
1513 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.188) + \
1514 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.200) + \
1515 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.209) + \
1516 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.210) + \
1517 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.206) + \
1518 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.205) + \
1519 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.198) + \
1520 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.187) + \
1521 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.175) + \
1522 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.167) + \
1523 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.157) + \
1524 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.138) + \
1525 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.118) + \
1526 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.094) + \
1527 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.083) + \
1528 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.062) + \
1529 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.049) + \
1530 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.044) + \
1531 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.036) + \
1532 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.035) + \
1533 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.026) + \
1534 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.041) + \
1535 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1536
1537 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1538 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.070) + \
1539 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.089) + \
1540 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.104) + \
1541 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.121) + \
1542 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.121) + \
1543 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.127) + \
1544 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.127) + \
1545 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.117) + \
1546 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.117) + \
1547 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.096) + \
1548 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.099) + \
1549 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.086) + \
1550 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.075) + \
1551 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.066) + \
1552 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.056) + \
1553 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.048) + \
1554 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.031) + \
1555 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.042) + \
1556 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.030) + \
1557 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.030) + \
1558 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1559
1560 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1561 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.041) + \
1562 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.047) + \
1563 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.059) + \
1564 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.056) + \
1565 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.068) + \
1566 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.073) + \
1567 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.080) + \
1568 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.063) + \
1569 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.065) + \
1570 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.057) + \
1571 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.049) + \
1572 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.067) + \
1573 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.065) + \
1574 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.044)+ \
1575 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.031) + \
1576 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.026) + \
1577 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1578
1579 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1580 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.02) + \
1581 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.02) + \
1582 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.03) + \
1583 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.03) + \
1584 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.03) + \
1585 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.035) + \
1586 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.04) + \
1587 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.03) + \
1588 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.03) + \
1589 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.025) + \
1590 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.025) + \
1591 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.035) + \
1592 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.032) + \
1593 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.022) + \
1594 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.015) + \
1595 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.013) + \
1596 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1597 (abs(eta) > 3.4) * (0.000)
1598 }
1599
1600
1601 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1602 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.528) + \
1603 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.611) + \
1604 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.666) + \
1605 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.701) + \
1606 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.734) + \
1607 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.741) + \
1608 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.746) + \
1609 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.750) + \
1610 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.742) + \
1611 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.729) + \
1612 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.711) + \
1613 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.694) + \
1614 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.684) + \
1615 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.642) + \
1616 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.602) + \
1617 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.554) + \
1618 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.488) + \
1619 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.417) + \
1620 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.343) + \
1621 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.284) + \
1622 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.250) + \
1623 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.221) + \
1624 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.164) + \
1625 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.173) + \
1626 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1627
1628 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1629 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.322) + \
1630 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.414) + \
1631 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.488) + \
1632 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.521) + \
1633 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.554) + \
1634 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.590) + \
1635 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.571) + \
1636 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.583) + \
1637 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.575) + \
1638 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.557) + \
1639 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.536) + \
1640 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.516) + \
1641 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.483) + \
1642 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.432) + \
1643 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.366) + \
1644 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.323) + \
1645 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.241) + \
1646 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.252) + \
1647 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.237) + \
1648 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.243) + \
1649 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1650
1651 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1652 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.200) + \
1653 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.252) + \
1654 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.310) + \
1655 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.353) + \
1656 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.377) + \
1657 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.425) + \
1658 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.406) + \
1659 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.414) + \
1660 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.412) + \
1661 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.366) + \
1662 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.360) + \
1663 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.347) + \
1664 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.310) + \
1665 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.307)+ \
1666 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.224) + \
1667 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.194) + \
1668 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1669
1670 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1671 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.100) + \
1672 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.12) + \
1673 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.15) + \
1674 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.17) + \
1675 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.19) + \
1676 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.21) + \
1677 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.20) + \
1678 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.21) + \
1679 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.20) + \
1680 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.18) + \
1681 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.18) + \
1682 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.17) + \
1683 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.15) + \
1684 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.15) + \
1685 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.11) + \
1686 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.10) + \
1687 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1688 (abs(eta) > 3.4) * (0.000)
1689 }
1690}
1691
1692
1693
1694#############
1695# b-tagging #
1696#############
1697module BTagging BTaggingTight {
1698
1699 set JetInputArray JetEnergyScale/jets
1700
1701 set BitNumber 2
1702
1703 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
1704 (pt > 20.0) * (abs(eta) <= 3.4) * (0.001) + \
1705 (abs(eta) > 3.4) * (0.000)
1706 }
1707
1708
1709 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1710 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.038) + \
1711 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.047) + \
1712 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.046) + \
1713 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.044) + \
1714 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.041) + \
1715 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.044) + \
1716 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.041) + \
1717 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.038) + \
1718 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.034) + \
1719 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.027) + \
1720 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.022) + \
1721 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.021) + \
1722 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.016) + \
1723 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.012) + \
1724 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.009) + \
1725 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.005) + \
1726 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.004) + \
1727 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.003) + \
1728 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.002) + \
1729 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.002) + \
1730 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.001) + \
1731 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.002) + \
1732 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.001) + \
1733 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.006) + \
1734 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1735
1736 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1737 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.020) + \
1738 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.0242) + \
1739 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.0272) + \
1740 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.0349) + \
1741 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.0314) + \
1742 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.0220) + \
1743 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.028 ) + \
1744 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.0263) + \
1745 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.0245) + \
1746 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.0150) + \
1747 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.0088) + \
1748 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.0115) + \
1749 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.0063) + \
1750 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.0083) + \
1751 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.0027) + \
1752 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.0052) + \
1753 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.0007) + \
1754 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.0048) + \
1755 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.0013) + \
1756 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.0013) + \
1757 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1758
1759 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1760 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.0092) + \
1761 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.0094) + \
1762 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.0141) + \
1763 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.0108) + \
1764 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.0161) + \
1765 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.0137) + \
1766 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.0135) + \
1767 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.0176) + \
1768 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.0149) + \
1769 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.0137) + \
1770 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.0114) + \
1771 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.0079) + \
1772 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.0112) + \
1773 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.0081)+ \
1774 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.0033) + \
1775 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.0016) + \
1776 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1777
1778 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1779 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.005) + \
1780 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.005) + \
1781 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.005) + \
1782 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.007) + \
1783 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.005) + \
1784 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.008) + \
1785 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.006) + \
1786 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.007) + \
1787 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.008) + \
1788 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.007) + \
1789 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.006) + \
1790 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.004) + \
1791 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.005) + \
1792 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.004) + \
1793 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.001) + \
1794 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.0005) + \
1795 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1796 (abs(eta) > 3.4) * (0.000)
1797 }
1798
1799
1800
1801 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1802 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.349) + \
1803 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.435) + \
1804 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.485) + \
1805 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.504) + \
1806 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.519) + \
1807 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.537) + \
1808 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.541) + \
1809 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.528) + \
1810 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.514) + \
1811 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.473) + \
1812 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.447) + \
1813 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.408) + \
1814 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.373) + \
1815 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.315) + \
1816 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.267) + \
1817 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.181) + \
1818 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.138) + \
1819 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.101) + \
1820 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.067) + \
1821 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.048) + \
1822 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.032) + \
1823 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.029) + \
1824 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.022) + \
1825 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.038) + \
1826 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.000) + \
1827
1828 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1829 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.186) + \
1830 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.258) + \
1831 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.306) + \
1832 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.348) + \
1833 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.359) + \
1834 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.358) + \
1835 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.367) + \
1836 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.380) + \
1837 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.359) + \
1838 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.308) + \
1839 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.238) + \
1840 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.252) + \
1841 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.229) + \
1842 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.174) + \
1843 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.125) + \
1844 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.096) + \
1845 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.066) + \
1846 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.046) + \
1847 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.076) + \
1848 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 3000) * (0.072) + \
1849 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 3000.0) * (0.000) + \
1850
1851 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt <= 20.0) * (0.000) + \
1852 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 30) * (0.096) + \
1853 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 30.0 && pt <= 40) * (0.119) + \
1854 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 40.0 && pt <= 50) * (0.156) + \
1855 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 60) * (0.174) + \
1856 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 60.0 && pt <= 70) * (0.192) + \
1857 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 70.0 && pt <= 80) * (0.228) + \
1858 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 80.0 && pt <= 90) * (0.218) + \
1859 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 90.0 && pt <= 100) * (0.231) + \
1860 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 100.0 && pt <= 120) * (0.224) + \
1861 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 120.0 && pt <= 140) * (0.184) + \
1862 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 140.0 && pt <= 160) * (0.164) + \
1863 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 160.0 && pt <= 180) * (0.115) + \
1864 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 180.0 && pt <= 200) * (0.123) + \
1865 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 200.0 && pt <= 250) * (0.131)+ \
1866 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 250.0 && pt <= 350) * (0.063) + \
1867 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 350.0 && pt <= 3000) * (0.039) + \
1868 (abs(eta) > 2.4 && abs(eta) <= 3.0) * (pt > 3000.0) * (0.000) + \
1869
1870 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt <= 20.0) * (0.000) + \
1871 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 20.0 && pt <= 30) * (0.05) + \
1872 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 30.0 && pt <= 40) * (0.06) + \
1873 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 40.0 && pt <= 50) * (0.08) + \
1874 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 50.0 && pt <= 60) * (0.09) + \
1875 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 60.0 && pt <= 70) * (0.10) + \
1876 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 70.0 && pt <= 80) * (0.11) + \
1877 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 80.0 && pt <= 90) * (0.11) + \
1878 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 90.0 && pt <= 100) * (0.11) + \
1879 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 100.0 && pt <= 120) * (0.11) + \
1880 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 120.0 && pt <= 140) * (0.09) + \
1881 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 140.0 && pt <= 160) * (0.08) + \
1882 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 160.0 && pt <= 180) * (0.06) + \
1883 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 180.0 && pt <= 200) * (0.06) + \
1884 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 200.0 && pt <= 250) * (0.06) + \
1885 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 250.0 && pt <= 350) * (0.03) + \
1886 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 350.0 && pt <= 3000) * (0.02) + \
1887 (abs(eta) > 3.0 && abs(eta) <= 3.4) * (pt > 3000.0) * (0.000) + \
1888 (abs(eta) > 3.4) * (0.000)
1889 }
1890}
1891
1892
1893
1894#############
1895# tau-tagging
1896#############
1897
1898
1899module TauTagging TauTagging {
1900 set ParticleInputArray Delphes/allParticles
1901 set PartonInputArray Delphes/partons
1902 set JetInputArray JetEnergyScale/jets
1903
1904 set DeltaR 0.5
1905
1906 set TauPTMin 20.0
1907
1908 set TauEtaMax 2.3
1909
1910 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1911
1912 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)) + \
1913 (abs(eta) > 2.3) * (0.000)
1914 }
1915 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) ) + \
1916 (abs(eta) > 2.3) * (0.000)
1917 }
1918}
1919
1920
1921###############################################################################################################
1922# StatusPidFilter: this module removes all generated particles except electrons, muons, taus, and status == 3 #
1923###############################################################################################################
1924
1925module StatusPidFilter GenParticleFilter {
1926
1927 set InputArray Delphes/allParticles
1928 set OutputArray filteredParticles
1929 set PTMin 5.0
1930
1931}
1932
1933
1934##################
1935# ROOT tree writer
1936##################
1937
1938module TreeWriter TreeWriter {
1939# add Branch InputArray BranchName BranchClass
1940 add Branch GenParticleFilter/filteredParticles Particle GenParticle
1941 add Branch PileUpMerger/vertices Vertex Vertex
1942
1943 add Branch GenJetFinder/jets GenJet Jet
1944 add Branch GenMissingET/momentum GenMissingET MissingET
1945
1946# add Branch HCal/eflowTracks EFlowTrack Track
1947# add Branch ECal/eflowPhotons EFlowPhoton Tower
1948# add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1949
1950 add Branch PhotonEfficiency/photons Photon Photon
1951 add Branch ElectronEfficiency/electrons Electron Electron
1952 add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
1953 add Branch MuonTightIdEfficiency/muons MuonTight Muon
1954
1955 add Branch JetEnergyScale/jets Jet Jet
1956
1957 add Branch MissingET/momentum MissingET MissingET
1958 add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET
1959 add Branch ScalarHT/energy ScalarHT ScalarHT
1960}
1961
Note: See TracBrowser for help on using the repository browser.