Fork me on GitHub

source: git/cards/CMS_PhaseII/CMS_PhaseII_200PU_v03.tcl@ cdcceb7

ImprovedOutputFile Timing dual_readout llp
Last change on this file since cdcceb7 was cdcceb7, checked in by Basil Schneider <basil.schneider@…>, 6 years ago

Add TauTagging for AK8, PUPPI and PUPPIAK8 jets
  • Property mode set to 100644
File size: 133.8 KB
Line 
1#
2# Main authors: Michele Selvaggi (CERN)
3#
4# Released on: Jun 26 - 2017
5#
6# Version: v02
7#
8#
9#######################################
10# Order of execution of various modules
11#######################################
12
13set ExecutionPath {
14
15 PileUpMerger
16 ParticlePropagator
17
18 ChargedHadronTrackingEfficiency
19 ElectronTrackingEfficiency
20 MuonTrackingEfficiency
21
22 ChargedHadronMomentumSmearing
23 ElectronEnergySmearing
24 MuonMomentumSmearing
25
26 DenseProp
27 DenseMergeTracks
28 TrackMergerProp
29 TrackMerger
30
31 ECal
32 HCal
33
34 PhotonEnergySmearing
35 ElectronFilter
36
37 TrackPileUpSubtractor
38 RecoPuFilter
39
40 TowerMerger
41 NeutralEFlowMerger
42
43 EFlowMerger
44 EFlowMergerCHS
45 Rho
46
47 LeptonFilterNoLep
48 LeptonFilterLep
49 RunPUPPIBase
50 RunPUPPIMerger
51 RunPUPPI
52
53 EFlowFilterPuppi
54 EFlowFilterCHS
55
56 PhotonFilter
57
58 PhotonCloner
59 PhotonIsolation
60 PhotonIsolationCHS
61 PhotonEfficiency
62 PhotonEfficiencyCHS
63
64 ElectronCloner
65 ElectronIsolation
66 ElectronIsolationCHS
67
68 ElectronEfficiency
69 ElectronEfficiencyCHS
70
71 MuonCloner
72 MuonIsolation
73 MuonIsolationCHS
74
75 MuonLooseIdEfficiency
76 MuonTightIdEfficiency
77
78 MuonLooseIdEfficiencyCHS
79 MuonTightIdEfficiencyCHS
80
81 NeutrinoFilter
82
83 MissingET
84 PuppiMissingET
85 GenMissingET
86 GenPileUpMissingET
87
88 GenJetFinder
89 GenJetFinderAK8
90 FastJetFinder
91 FastJetFinderAK8
92 JetPileUpSubtractor
93 JetPileUpSubtractorAK8
94 FastJetFinderPUPPI
95 FastJetFinderPUPPIAK8
96
97 ScalarHT
98
99 JetEnergyScale
100 JetEnergyScaleAK8
101 JetEnergyScalePUPPI
102 JetEnergyScalePUPPIAK8
103
104 JetFlavorAssociation
105 JetFlavorAssociationAK8
106 JetFlavorAssociationPUPPI
107 JetFlavorAssociationPUPPIAK8
108
109 BTagging
110 BTaggingAK8
111 BTaggingPUPPI
112 BTaggingPUPPIAK8
113
114 TauTagging
115 TauTaggingAK8
116 TauTaggingPUPPI
117 TauTaggingPUPPIAK8
118
119 GenParticleFilter
120
121 TreeWriter
122}
123
124
125###############
126# PileUp Merger
127###############
128
129module PileUpMerger PileUpMerger {
130 set InputArray Delphes/stableParticles
131
132 set ParticleOutputArray stableParticles
133 set VertexOutputArray vertices
134
135 # pre-generated minbias input file
136 set PileUpFile /eos/cms/store/group/upgrade/delphes/PhaseII/MinBias_100k.pileup
137 #set PileUpFile MinBias.pileup
138
139 # average expected pile up
140 set MeanPileUp 200
141
142 # maximum spread in the beam direction in m
143 set ZVertexSpread 0.25
144
145 # maximum spread in time in s
146 set TVertexSpread 800E-12
147
148 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s) - {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
149 set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
150
151}
152
153
154
155#################################
156# Dense Track propagation calo
157#################################
158
159module ParticlePropagator ParticlePropagator {
160 set InputArray PileUpMerger/stableParticles
161
162 set OutputArray stableParticles
163 set NeutralOutputArray neutralParticles
164 set ChargedHadronOutputArray chargedHadrons
165 set ElectronOutputArray electrons
166 set MuonOutputArray muons
167
168 # radius of the magnetic field coverage, in m
169 set Radius 1.29
170 # half-length of the magnetic field coverage, in m
171 set HalfLength 3.0
172
173 # magnetic field
174 set Bz 3.8
175}
176
177
178
179
180####################################
181# Charged hadron tracking efficiency
182####################################
183
184module Efficiency ChargedHadronTrackingEfficiency {
185 ## particles after propagation
186 set InputArray ParticlePropagator/chargedHadrons
187 set OutputArray chargedHadrons
188 # tracking efficiency formula for charged hadrons
189 set EfficiencyFormula {
190 (pt <= 0.2) * (0.00) + \
191 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
192 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
193 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
194 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.87) + \
195 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
196 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.82) + \
197 (abs(eta) > 4.0) * (0.00)
198 }
199}
200
201
202#####################################
203# Electron tracking efficiency - ID
204####################################
205
206module Efficiency ElectronTrackingEfficiency {
207 set InputArray ParticlePropagator/electrons
208 set OutputArray electrons
209 # tracking efficiency formula for electrons
210 set EfficiencyFormula {
211 (pt <= 0.2) * (0.00) + \
212 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
213 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
214 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
215 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 10.0) * (0.82+pt*0.01) + \
216 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 10.0) * (0.90) + \
217 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
218 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.8+pt*0.01) + \
219 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.85) + \
220 (abs(eta) > 4.0) * (0.00)
221
222 }
223}
224
225##########################
226# Muon tracking efficiency
227##########################
228
229module Efficiency MuonTrackingEfficiency {
230 set InputArray ParticlePropagator/muons
231 set OutputArray muons
232 # tracking efficiency formula for muons
233 set EfficiencyFormula {
234 (pt <= 0.2) * (0.00) + \
235 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 1.00) + \
236 (abs(eta) <= 1.2) * (pt > 1.0) * (1.00) + \
237 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 0.2 && pt <= 1.0) * (pt*1.00) + \
238 (abs(eta) > 1.2 && abs(eta) <= 2.8) * (pt > 1.0) * (1.00) + \
239 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
240 (abs(eta) > 2.8 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
241 (abs(eta) > 4.0) * (0.00)
242
243 }
244}
245
246
247########################################
248# Momentum resolution for charged tracks
249########################################
250
251module MomentumSmearing ChargedHadronMomentumSmearing {
252 ## hadrons after having applied the tracking efficiency
253 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
254 set OutputArray chargedHadrons
255 # resolution formula for charged hadrons ,
256
257 #
258 # Automatically generated tracker resolution formula for layout: OT612IT4025
259 #
260 # By Unknown author on: 2017-06-30.17:03:00
261 #
262 set ResolutionFormula { (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.00457888) + \
263 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.004579 + (pt-1.000000)* 0.000045) + \
264 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.004983 + (pt-10.000000)* 0.000047) + \
265 (abs(eta) >= 0.0000 && abs(eta) < 0.2000) * (pt >= 100.0000) * (0.009244*pt/100.000000) + \
266 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.00505011) + \
267 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.005050 + (pt-1.000000)* 0.000033) + \
268 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.005343 + (pt-10.000000)* 0.000043) + \
269 (abs(eta) >= 0.2000 && abs(eta) < 0.4000) * (pt >= 100.0000) * (0.009172*pt/100.000000) + \
270 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.00510573) + \
271 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.005106 + (pt-1.000000)* 0.000023) + \
272 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.005317 + (pt-10.000000)* 0.000042) + \
273 (abs(eta) >= 0.4000 && abs(eta) < 0.6000) * (pt >= 100.0000) * (0.009077*pt/100.000000) + \
274 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.00578020) + \
275 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.005780 + (pt-1.000000)* -0.000000) + \
276 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.005779 + (pt-10.000000)* 0.000038) + \
277 (abs(eta) >= 0.6000 && abs(eta) < 0.8000) * (pt >= 100.0000) * (0.009177*pt/100.000000) + \
278 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.00728723) + \
279 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.007287 + (pt-1.000000)* -0.000031) + \
280 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.007011 + (pt-10.000000)* 0.000038) + \
281 (abs(eta) >= 0.8000 && abs(eta) < 1.0000) * (pt >= 100.0000) * (0.010429*pt/100.000000) + \
282 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.01045117) + \
283 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.010451 + (pt-1.000000)* -0.000051) + \
284 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.009989 + (pt-10.000000)* 0.000043) + \
285 (abs(eta) >= 1.0000 && abs(eta) < 1.2000) * (pt >= 100.0000) * (0.013867*pt/100.000000) + \
286 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.01477199) + \
287 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.014772 + (pt-1.000000)* -0.000128) + \
288 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.013616 + (pt-10.000000)* 0.000035) + \
289 (abs(eta) >= 1.2000 && abs(eta) < 1.4000) * (pt >= 100.0000) * (0.016800*pt/100.000000) + \
290 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.01731474) + \
291 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.017315 + (pt-1.000000)* -0.000208) + \
292 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.015439 + (pt-10.000000)* 0.000030) + \
293 (abs(eta) >= 1.4000 && abs(eta) < 1.6000) * (pt >= 100.0000) * (0.018161*pt/100.000000) + \
294 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.01942025) + \
295 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.019420 + (pt-1.000000)* -0.000417) + \
296 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.015669 + (pt-10.000000)* 0.000026) + \
297 (abs(eta) >= 1.6000 && abs(eta) < 1.8000) * (pt >= 100.0000) * (0.018039*pt/100.000000) + \
298 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.02201432) + \
299 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.022014 + (pt-1.000000)* -0.000667) + \
300 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.016012 + (pt-10.000000)* 0.000045) + \
301 (abs(eta) >= 1.8000 && abs(eta) < 2.0000) * (pt >= 100.0000) * (0.020098*pt/100.000000) + \
302 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.02574300) + \
303 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.025743 + (pt-1.000000)* -0.001118) + \
304 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.015681 + (pt-10.000000)* 0.000051) + \
305 (abs(eta) >= 2.0000 && abs(eta) < 2.2000) * (pt >= 100.0000) * (0.020289*pt/100.000000) + \
306 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.02885821) + \
307 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.028858 + (pt-1.000000)* -0.001345) + \
308 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.016753 + (pt-10.000000)* 0.000053) + \
309 (abs(eta) >= 2.2000 && abs(eta) < 2.4000) * (pt >= 100.0000) * (0.021524*pt/100.000000) + \
310 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.03204812) + \
311 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.032048 + (pt-1.000000)* -0.001212) + \
312 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.021138 + (pt-10.000000)* 0.000037) + \
313 (abs(eta) >= 2.4000 && abs(eta) < 2.6000) * (pt >= 100.0000) * (0.024477*pt/100.000000) + \
314 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.03950405) + \
315 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.039504 + (pt-1.000000)* -0.001386) + \
316 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.027026 + (pt-10.000000)* 0.000037) + \
317 (abs(eta) >= 2.6000 && abs(eta) < 2.8000) * (pt >= 100.0000) * (0.030392*pt/100.000000) + \
318 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.04084751) + \
319 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.040848 + (pt-1.000000)* -0.001780) + \
320 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.024824 + (pt-10.000000)* 0.000029) + \
321 (abs(eta) >= 2.8000 && abs(eta) < 3.0000) * (pt >= 100.0000) * (0.027445*pt/100.000000) + \
322 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 0.0000 && pt < 1.0000) * (0.04532425) + \
323 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 1.0000 && pt < 10.0000) * (0.045324 + (pt-1.000000)* -0.002497) + \
324 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 10.0000 && pt < 100.0000) * (0.022851 + (pt-10.000000)* 0.000024) + \
325 (abs(eta) >= 3.0000 && abs(eta) < 3.2000) * (pt >= 100.0000) * (0.025053*pt/100.000000) + \
326 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 0.0000 && pt < 1.0000) * (0.06418925) + \
327 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 1.0000 && pt < 10.0000) * (0.064189 + (pt-1.000000)* -0.004055) + \
328 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 10.0000 && pt < 100.0000) * (0.027691 + (pt-10.000000)* 0.000034) + \
329 (abs(eta) >= 3.2000 && abs(eta) < 3.4000) * (pt >= 100.0000) * (0.030710*pt/100.000000) + \
330 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 0.0000 && pt < 1.0000) * (0.07682500) + \
331 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 1.0000 && pt < 10.0000) * (0.076825 + (pt-1.000000)* -0.004510) + \
332 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 10.0000 && pt < 100.0000) * (0.036234 + (pt-10.000000)* 0.000049) + \
333 (abs(eta) >= 3.4000 && abs(eta) < 3.6000) * (pt >= 100.0000) * (0.040629*pt/100.000000) + \
334 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 0.0000 && pt < 1.0000) * (0.09796358) + \
335 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 1.0000 && pt < 10.0000) * (0.097964 + (pt-1.000000)* -0.005758) + \
336 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 10.0000 && pt < 100.0000) * (0.046145 + (pt-10.000000)* 0.000069) + \
337 (abs(eta) >= 3.6000 && abs(eta) < 3.8000) * (pt >= 100.0000) * (0.052345*pt/100.000000) + \
338 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 0.0000 && pt < 1.0000) * (0.13415929) + \
339 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 1.0000 && pt < 10.0000) * (0.134159 + (pt-1.000000)* -0.008283) + \
340 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 10.0000 && pt < 100.0000) * (0.059612 + (pt-10.000000)* 0.000111) + \
341 (abs(eta) >= 3.8000 && abs(eta) < 4.0000) * (pt >= 100.0000) * (0.069617*pt/100.000000)
342 }
343
344
345}
346
347#################################
348# Energy resolution for electrons
349#################################
350
351module EnergySmearing ElectronEnergySmearing {
352 set InputArray ElectronTrackingEfficiency/electrons
353 set OutputArray electrons
354
355 # set ResolutionFormula {resolution formula as a function of eta and energy}
356
357 # resolution formula for electrons
358
359 # taking something flat in energy for now, ECAL will take over at high energy anyway.
360 # inferred from hep-ex/1306.2016 and 1502.02701
361 set ResolutionFormula {
362
363 (abs(eta) <= 1.5) * (energy*0.028) +
364 (abs(eta) > 1.5 && abs(eta) <= 1.75) * (energy*0.037) +
365 (abs(eta) > 1.75 && abs(eta) <= 2.15) * (energy*0.038) +
366 (abs(eta) > 2.15 && abs(eta) <= 3.00) * (energy*0.044) +
367 (abs(eta) > 3.00 && abs(eta) <= 4.00) * (energy*0.10)}
368
369}
370
371###############################
372# Momentum resolution for muons
373###############################
374
375module MomentumSmearing MuonMomentumSmearing {
376 set InputArray MuonTrackingEfficiency/muons
377 set OutputArray muons
378 # resolution formula for muons
379
380 # up to |eta| < 2.8 take measurement from tracking + muon chambers
381 # for |eta| > 2.8 and pT < 5.0 take measurement from tracking alone taken from
382 # http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
383 source muonMomentumResolution.tcl
384}
385
386##############
387# Track merger
388##############
389
390module Merger TrackMergerProp {
391# add InputArray InputArray
392 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
393 add InputArray ElectronEnergySmearing/electrons
394 add InputArray MuonMomentumSmearing/muons
395 set OutputArray tracks
396}
397
398
399#################################
400# Dense Track propagation
401#################################
402
403module ParticlePropagator DenseProp {
404 set InputArray PileUpMerger/stableParticles
405
406 set OutputArray stableParticles
407 set NeutralOutputArray neutralParticles
408 set ChargedHadronOutputArray chargedHadrons
409 set ElectronOutputArray electrons
410 set MuonOutputArray muons
411
412 # radius of the first pixel layer
413 set Radius 0.3
414 set RadiusMax 1.29
415 # half-length of the magnetic field coverage, in m
416 set HalfLength 0.7
417 set HalfLengthMax 3.0
418
419 # magnetic field
420 set Bz 3.8
421}
422
423####################
424# Dense Track merger
425###################
426
427module Merger DenseMergeTracks {
428# add InputArray InputArray
429 add InputArray DenseProp/chargedHadrons
430 add InputArray DenseProp/electrons
431 add InputArray DenseProp/muons
432 set OutputArray tracks
433}
434
435######################
436# Dense Track Filter
437######################
438
439module DenseTrackFilter TrackMerger {
440 set TrackInputArray TrackMergerProp/tracks
441 set DenseChargedInputArray DenseMergeTracks/tracks
442
443 set TrackOutputArray tracks
444
445 set EtaPhiRes 0.003
446 set EtaMax 4.0
447
448 set pi [expr {acos(-1)}]
449
450 set nbins_phi [expr {$pi/$EtaPhiRes} ]
451 set nbins_phi [expr {int($nbins_phi)} ]
452
453 set PhiBins {}
454 for {set i -$nbins_phi} {$i <= $nbins_phi} {incr i} {
455 add PhiBins [expr {$i * $pi/$nbins_phi}]
456 }
457
458 set nbins_eta [expr {$EtaMax/$EtaPhiRes} ]
459 set nbins_eta [expr {int($nbins_eta)} ]
460
461 for {set i -$nbins_eta} {$i <= $nbins_eta} {incr i} {
462 set eta [expr {$i * $EtaPhiRes}]
463 add EtaPhiBins $eta $PhiBins
464 }
465
466}
467
468
469#############
470# ECAL
471#############
472
473module SimpleCalorimeter ECal {
474 set ParticleInputArray ParticlePropagator/stableParticles
475 set TrackInputArray TrackMerger/tracks
476
477 set TowerOutputArray ecalTowers
478 set EFlowTrackOutputArray eflowTracks
479 set EFlowTowerOutputArray eflowPhotons
480
481 set IsEcal true
482
483 set EnergyMin 0.5
484 set EnergySignificanceMin 1.0
485
486 set SmearTowerCenter true
487
488 set pi [expr {acos(-1)}]
489
490 # lists of the edges of each tower in eta and phi
491 # each list starts with the lower edge of the first tower
492 # the list ends with the higher edged of the last tower
493
494 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
495
496 set PhiBins {}
497 for {set i -180} {$i <= 180} {incr i} {
498 add PhiBins [expr {$i * $pi/180.0}]
499 }
500
501 # 0.02 unit in eta up to eta = 1.5 (barrel)
502 for {set i -85} {$i <= 86} {incr i} {
503 set eta [expr {$i * 0.0174}]
504 add EtaPhiBins $eta $PhiBins
505 }
506
507 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
508
509 set PhiBins {}
510 for {set i -180} {$i <= 180} {incr i} {
511 add PhiBins [expr {$i * $pi/180.0}]
512 }
513
514 # 0.02 unit in eta up to eta = 3
515 for {set i 1} {$i <= 84} {incr i} {
516 set eta [expr { -2.958 + $i * 0.0174}]
517 add EtaPhiBins $eta $PhiBins
518 }
519
520 for {set i 1} {$i <= 84} {incr i} {
521 set eta [expr { 1.4964 + $i * 0.0174}]
522 add EtaPhiBins $eta $PhiBins
523 }
524
525 # take present CMS granularity for HF
526
527 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
528 set PhiBins {}
529 for {set i -18} {$i <= 18} {incr i} {
530 add PhiBins [expr {$i * $pi/18.0}]
531 }
532
533 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} {
534 add EtaPhiBins $eta $PhiBins
535 }
536
537
538 add EnergyFraction {0} {0.0}
539 # energy fractions for e, gamma and pi0
540 add EnergyFraction {11} {1.0}
541 add EnergyFraction {22} {1.0}
542 add EnergyFraction {111} {1.0}
543 # energy fractions for muon, neutrinos and neutralinos
544 add EnergyFraction {12} {0.0}
545 add EnergyFraction {13} {0.0}
546 add EnergyFraction {14} {0.0}
547 add EnergyFraction {16} {0.0}
548 add EnergyFraction {1000022} {0.0}
549 add EnergyFraction {1000023} {0.0}
550 add EnergyFraction {1000025} {0.0}
551 add EnergyFraction {1000035} {0.0}
552 add EnergyFraction {1000045} {0.0}
553 # energy fractions for K0short and Lambda
554 add EnergyFraction {310} {0.3}
555 add EnergyFraction {3122} {0.3}
556
557 # set ResolutionFormula {resolution formula as a function of eta and energy}
558
559 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
560 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117
561
562 set ResolutionFormula { (abs(eta) <= 1.50) * sqrt(energy^2*0.009^2 + energy*0.12^2 + 0.45^2) +
563 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
564 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
565 (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
566 (abs(eta) >= 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
567
568}
569
570#############
571# HCAL
572#############
573
574module SimpleCalorimeter HCal {
575 set ParticleInputArray ParticlePropagator/stableParticles
576 set TrackInputArray ECal/eflowTracks
577
578 set TowerOutputArray hcalTowers
579 set EFlowTrackOutputArray eflowTracks
580 set EFlowTowerOutputArray eflowNeutralHadrons
581
582 set IsEcal false
583
584 set EnergyMin 1.0
585 set EnergySignificanceMin 1.0
586
587 set SmearTowerCenter true
588
589 set pi [expr {acos(-1)}]
590
591 # lists of the edges of each tower in eta and phi
592 # each list starts with the lower edge of the first tower
593 # the list ends with the higher edged of the last tower
594
595 # assume 0.087 x 0.087 resolution in eta,phi in the barrel |eta| < 1.5
596
597 set PhiBins {}
598 for {set i -36} {$i <= 36} {incr i} {
599 add PhiBins [expr {$i * $pi/36.0}]
600 }
601 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} {
602 add EtaPhiBins $eta $PhiBins
603 }
604
605 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- HCAL)
606
607 set PhiBins {}
608 for {set i -45} {$i <= 45} {incr i} {
609 add PhiBins [expr {$i * $pi/45.0}]
610 }
611
612 # 0.02 unit in eta up to eta = 3
613 for {set i 1} {$i <= 21} {incr i} {
614 set eta [expr { -2.958 + $i * 0.0696}]
615 add EtaPhiBins $eta $PhiBins
616 }
617
618 for {set i 1} {$i <= 21} {incr i} {
619 set eta [expr { 1.4964 + $i * 0.0696}]
620 add EtaPhiBins $eta $PhiBins
621 }
622
623 # take present CMS granularity for HF
624
625 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
626 set PhiBins {}
627 for {set i -18} {$i <= 18} {incr i} {
628 add PhiBins [expr {$i * $pi/18.0}]
629 }
630
631 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} {
632 add EtaPhiBins $eta $PhiBins
633 }
634
635
636 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
637 add EnergyFraction {0} {1.0}
638 # energy fractions for e, gamma and pi0
639 add EnergyFraction {11} {0.0}
640 add EnergyFraction {22} {0.0}
641 add EnergyFraction {111} {0.0}
642 # energy fractions for muon, neutrinos and neutralinos
643 add EnergyFraction {12} {0.0}
644 add EnergyFraction {13} {0.0}
645 add EnergyFraction {14} {0.0}
646 add EnergyFraction {16} {0.0}
647 add EnergyFraction {1000022} {0.0}
648 add EnergyFraction {1000023} {0.0}
649 add EnergyFraction {1000025} {0.0}
650 add EnergyFraction {1000035} {0.0}
651 add EnergyFraction {1000045} {0.0}
652 # energy fractions for K0short and Lambda
653 add EnergyFraction {310} {0.7}
654 add EnergyFraction {3122} {0.7}
655
656# set ResolutionFormula {resolution formula as a function of eta and energy}
657 set ResolutionFormula { (abs(eta) <= 1.5) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
658 (abs(eta) > 1.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
659 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.11^2 + energy*2.80^2)}
660
661}
662
663#################################
664# Energy resolution for electrons
665#################################
666
667module EnergySmearing PhotonEnergySmearing {
668 set InputArray ECal/eflowPhotons
669 set OutputArray eflowPhotons
670
671 # adding 1% extra photon smearing
672 set ResolutionFormula {energy*0.01}
673
674}
675
676
677
678#################
679# Electron filter
680#################
681
682module PdgCodeFilter ElectronFilter {
683 set InputArray HCal/eflowTracks
684 set OutputArray electrons
685 set Invert true
686 add PdgCode {11}
687 add PdgCode {-11}
688}
689
690
691##########################
692# Track pile-up subtractor
693##########################
694
695module TrackPileUpSubtractor TrackPileUpSubtractor {
696# add InputArray InputArray OutputArray
697 add InputArray HCal/eflowTracks eflowTracks
698 add InputArray ElectronFilter/electrons electrons
699 add InputArray MuonMomentumSmearing/muons muons
700
701 set VertexInputArray PileUpMerger/vertices
702 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
703 # Z vertex resolution in m
704 set ZVertexResolution {0.0001}
705}
706
707########################
708# Reco PU filter
709########################
710
711module RecoPuFilter RecoPuFilter {
712 set InputArray HCal/eflowTracks
713 set OutputArray eflowTracks
714}
715
716###################################################
717# Tower Merger (in case not using e-flow algorithm)
718###################################################
719
720module Merger TowerMerger {
721# add InputArray InputArray
722 add InputArray ECal/ecalTowers
723 add InputArray HCal/hcalTowers
724 set OutputArray towers
725}
726
727####################
728# Neutral eflow erger
729####################
730
731module Merger NeutralEFlowMerger {
732# add InputArray InputArray
733 add InputArray PhotonEnergySmearing/eflowPhotons
734 add InputArray HCal/eflowNeutralHadrons
735 set OutputArray eflowTowers
736}
737
738#####################
739# Energy flow merger
740#####################
741
742module Merger EFlowMerger {
743# add InputArray InputArray
744 add InputArray HCal/eflowTracks
745 add InputArray PhotonEnergySmearing/eflowPhotons
746 add InputArray HCal/eflowNeutralHadrons
747 set OutputArray eflow
748}
749
750############################
751# Energy flow merger no PU
752############################
753
754module Merger EFlowMergerCHS {
755# add InputArray InputArray
756 add InputArray RecoPuFilter/eflowTracks
757 add InputArray PhotonEnergySmearing/eflowPhotons
758 add InputArray HCal/eflowNeutralHadrons
759 set OutputArray eflow
760}
761
762#########################################
763### Run the puppi code (to be tuned) ###
764#########################################
765
766module PdgCodeFilter LeptonFilterNoLep {
767 set InputArray HCal/eflowTracks
768 set OutputArray eflowTracksNoLeptons
769 set Invert false
770 add PdgCode {13}
771 add PdgCode {-13}
772 add PdgCode {11}
773 add PdgCode {-11}
774}
775
776module PdgCodeFilter LeptonFilterLep {
777 set InputArray HCal/eflowTracks
778 set OutputArray eflowTracksLeptons
779 set Invert true
780 add PdgCode {11}
781 add PdgCode {-11}
782 add PdgCode {13}
783 add PdgCode {-13}
784}
785
786module RunPUPPI RunPUPPIBase {
787 ## input information
788 set TrackInputArray LeptonFilterNoLep/eflowTracksNoLeptons
789 set NeutralInputArray NeutralEFlowMerger/eflowTowers
790 set PVInputArray PileUpMerger/vertices
791 set MinPuppiWeight 0.05
792 set UseExp false
793 set UseNoLep false
794
795 ## define puppi algorithm parameters (more than one for the same eta region is possible)
796 add EtaMinBin 0.0 1.5 4.0
797 add EtaMaxBin 1.5 4.0 10.0
798 add PtMinBin 0.0 0.0 0.0
799 add ConeSizeBin 0.2 0.2 0.2
800 add RMSPtMinBin 0.1 0.5 0.5
801 add RMSScaleFactorBin 1.0 1.0 1.0
802 add NeutralMinEBin 0.2 0.2 0.5
803 add NeutralPtSlope 0.006 0.013 0.067
804 add ApplyCHS true true true
805 add UseCharged true true false
806 add ApplyLowPUCorr true true true
807 add MetricId 5 5 5
808 add CombId 0 0 0
809
810 ## output name
811 set OutputArray PuppiParticles
812 set OutputArrayTracks puppiTracks
813 set OutputArrayNeutrals puppiNeutrals
814}
815
816module Merger RunPUPPIMerger {
817 add InputArray RunPUPPIBase/PuppiParticles
818 add InputArray LeptonFilterLep/eflowTracksLeptons
819 set OutputArray PuppiParticles
820}
821
822# need this because of leptons that were added back
823module RecoPuFilter RunPUPPI {
824 set InputArray RunPUPPIMerger/PuppiParticles
825 set OutputArray PuppiParticles
826}
827
828######################
829# EFlowFilterPuppi
830######################
831
832module PdgCodeFilter EFlowFilterPuppi {
833 set InputArray RunPUPPI/PuppiParticles
834 set OutputArray eflow
835
836 add PdgCode {11}
837 add PdgCode {-11}
838 add PdgCode {13}
839 add PdgCode {-13}
840}
841
842######################
843# EFlowFilterCHS
844######################
845
846module PdgCodeFilter EFlowFilterCHS {
847 set InputArray EFlowMergerCHS/eflow
848 set OutputArray eflow
849
850 add PdgCode {11}
851 add PdgCode {-11}
852 add PdgCode {13}
853 add PdgCode {-13}
854}
855
856
857###################
858# Missing ET merger
859###################
860
861module Merger MissingET {
862# add InputArray InputArray
863# add InputArray RunPUPPI/PuppiParticles
864 add InputArray EFlowMerger/eflow
865 set MomentumOutputArray momentum
866}
867
868module Merger PuppiMissingET {
869 #add InputArray InputArray
870 add InputArray RunPUPPI/PuppiParticles
871 #add InputArray EFlowMerger/eflow
872 set MomentumOutputArray momentum
873}
874
875###################
876# Ger PileUp Missing ET
877###################
878
879module Merger GenPileUpMissingET {
880# add InputArray InputArray
881# add InputArray RunPUPPI/PuppiParticles
882 add InputArray ParticlePropagator/stableParticles
883 set MomentumOutputArray momentum
884}
885
886##################
887# Scalar HT merger
888##################
889
890module Merger ScalarHT {
891# add InputArray InputArray
892 add InputArray RunPUPPI/PuppiParticles
893 set EnergyOutputArray energy
894}
895
896#################
897# Neutrino Filter
898#################
899
900module PdgCodeFilter NeutrinoFilter {
901
902 set InputArray Delphes/stableParticles
903 set OutputArray filteredParticles
904
905 set PTMin 0.0
906
907 add PdgCode {12}
908 add PdgCode {14}
909 add PdgCode {16}
910 add PdgCode {-12}
911 add PdgCode {-14}
912 add PdgCode {-16}
913
914}
915
916#####################
917# MC truth jet finder
918#####################
919
920module FastJetFinder GenJetFinder {
921 set InputArray NeutrinoFilter/filteredParticles
922
923 set OutputArray jets
924
925 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
926 set JetAlgorithm 6
927 set ParameterR 0.4
928
929 set JetPTMin 15.0
930}
931
932module FastJetFinder GenJetFinderAK8 {
933 set InputArray NeutrinoFilter/filteredParticles
934
935 set OutputArray jetsAK8
936
937 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
938 set JetAlgorithm 6
939 set ParameterR 0.8
940
941 set JetPTMin 200.0
942}
943
944#########################
945# Gen Missing ET merger
946########################
947
948module Merger GenMissingET {
949
950# add InputArray InputArray
951 add InputArray NeutrinoFilter/filteredParticles
952 set MomentumOutputArray momentum
953}
954
955
956#############
957# Rho pile-up
958#############
959
960module FastJetGridMedianEstimator Rho {
961
962 set InputArray EFlowMergerCHS/eflow
963 set RhoOutputArray rho
964
965 # add GridRange rapmin rapmax drap dphi
966 # rapmin - the minimum rapidity extent of the grid
967 # rapmax - the maximum rapidity extent of the grid
968 # drap - the grid spacing in rapidity
969 # dphi - the grid spacing in azimuth
970
971 add GridRange -5.0 -4.0 1.0 1.0
972 add GridRange -4.0 -1.5 1.0 1.0
973 add GridRange -1.5 1.5 1.0 1.0
974 add GridRange 1.5 4.0 1.0 1.0
975 add GridRange 4.0 5.0 1.0 1.0
976
977}
978
979
980##############
981# Jet finder
982##############
983
984module FastJetFinder FastJetFinder {
985# set InputArray TowerMerger/towers
986 set InputArray EFlowMergerCHS/eflow
987
988 set OutputArray jets
989
990 set AreaAlgorithm 5
991
992 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
993 set JetAlgorithm 6
994 set ParameterR 0.4
995
996 set JetPTMin 15.0
997}
998
999#module Class Name
1000module FastJetFinder FastJetFinderAK8 {
1001# set InputArray TowerMerger/towers
1002 set InputArray EFlowMergerCHS/eflow
1003
1004 set OutputArray jets
1005
1006 set AreaAlgorithm 5
1007
1008 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1009 set JetAlgorithm 6
1010 set ParameterR 0.8
1011
1012 set ComputeNsubjettiness 1
1013 set Beta 1.0
1014 set AxisMode 4
1015
1016 set ComputeTrimming 1
1017 set RTrim 0.2
1018 set PtFracTrim 0.05
1019
1020 set ComputePruning 1
1021 set ZcutPrun 0.1
1022 set RcutPrun 0.5
1023 set RPrun 0.8
1024
1025 set ComputeSoftDrop 1
1026 set BetaSoftDrop 0.0
1027 set SymmetryCutSoftDrop 0.1
1028 set R0SoftDrop 0.8
1029
1030 set JetPTMin 200.0
1031}
1032
1033###########################
1034# Jet Pile-Up Subtraction
1035###########################
1036
1037module JetPileUpSubtractor JetPileUpSubtractor {
1038 set JetInputArray FastJetFinder/jets
1039 set RhoInputArray Rho/rho
1040
1041 set OutputArray jets
1042
1043 set JetPTMin 15.0
1044}
1045
1046##############################
1047# Jet Pile-Up Subtraction AK8
1048##############################
1049
1050module JetPileUpSubtractor JetPileUpSubtractorAK8 {
1051 set JetInputArray FastJetFinderAK8/jets
1052 set RhoInputArray Rho/rho
1053
1054 set OutputArray jets
1055
1056 set JetPTMin 15.0
1057}
1058
1059module FastJetFinder FastJetFinderPUPPI {
1060# set InputArray TowerMerger/towers
1061 set InputArray RunPUPPI/PuppiParticles
1062
1063 set OutputArray jets
1064
1065 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
1066 set JetAlgorithm 6
1067 set ParameterR 0.4
1068
1069 set JetPTMin 15.0
1070}
1071
1072
1073module FastJetFinder FastJetFinderPUPPIAK8 {
1074# set InputArray TowerMerger/towers
1075 set InputArray RunPUPPI/PuppiParticles
1076
1077 set OutputArray jets
1078
1079 set JetAlgorithm 6
1080 set ParameterR 0.8
1081
1082 set ComputeNsubjettiness 1
1083 set Beta 1.0
1084 set AxisMode 4
1085
1086 set ComputeTrimming 1
1087 set RTrim 0.2
1088 set PtFracTrim 0.05
1089
1090 set ComputePruning 1
1091 set ZcutPrun 0.1
1092 set RcutPrun 0.5
1093 set RPrun 0.8
1094
1095 set ComputeSoftDrop 1
1096 set BetaSoftDrop 0.0
1097 set SymmetryCutSoftDrop 0.1
1098 set R0SoftDrop 0.8
1099
1100 set JetPTMin 200.0
1101}
1102
1103##################
1104# Jet Energy Scale
1105##################
1106
1107module EnergyScale JetEnergyScale {
1108 set InputArray JetPileUpSubtractor/jets
1109 set OutputArray jets
1110
1111 # scale formula for jets
1112 set ScaleFormula {1.00}
1113}
1114
1115module EnergyScale JetEnergyScaleAK8 {
1116 set InputArray JetPileUpSubtractorAK8/jets
1117 set OutputArray jets
1118
1119 # scale formula for jets
1120 set ScaleFormula {1.00}
1121}
1122
1123module EnergyScale JetEnergyScalePUPPI {
1124 set InputArray FastJetFinderPUPPI/jets
1125 set OutputArray jets
1126
1127 # scale formula for jets
1128 set ScaleFormula {1.00}
1129}
1130
1131module EnergyScale JetEnergyScalePUPPIAK8 {
1132 set InputArray FastJetFinderPUPPIAK8/jets
1133 set OutputArray jets
1134
1135 # scale formula for jets
1136 set ScaleFormula {1.00}
1137}
1138
1139#################
1140# Photon filter
1141#################
1142
1143module PdgCodeFilter PhotonFilter {
1144 set InputArray PhotonEnergySmearing/eflowPhotons
1145 set OutputArray photons
1146 set Invert true
1147 set PTMin 5.0
1148 add PdgCode {22}
1149}
1150
1151
1152##################
1153# Muon cloner #
1154##################
1155
1156module Cloner MuonCloner {
1157 set InputArray MuonMomentumSmearing/muons
1158 set OutputArray muons
1159}
1160
1161####################
1162# Electron cloner #
1163####################
1164
1165module Cloner ElectronCloner {
1166 set InputArray ElectronFilter/electrons
1167 set OutputArray electrons
1168}
1169
1170##################
1171# Photon cloner #
1172##################
1173
1174module Cloner PhotonCloner {
1175 set InputArray PhotonFilter/photons
1176 set OutputArray photons
1177}
1178
1179
1180####################
1181# Photon isolation #
1182####################
1183
1184module Isolation PhotonIsolation {
1185
1186 # particle for which calculate the isolation
1187 set CandidateInputArray PhotonFilter/photons
1188
1189 # isolation collection
1190 set IsolationInputArray EFlowFilterPuppi/eflow
1191
1192 # output array
1193 set OutputArray photons
1194
1195 # veto isolation cand. based on proximity to input cand.
1196 set DeltaRMin 0.01
1197 set UseMiniCone true
1198
1199 # isolation cone
1200 set DeltaRMax 0.3
1201
1202 # minimum pT
1203 set PTMin 0.0
1204
1205 # iso ratio to cut
1206 set PTRatioMax 9999.
1207
1208}
1209
1210
1211########################
1212# Photon isolation CHS #
1213########################
1214
1215module Isolation PhotonIsolationCHS {
1216
1217 # particle for which calculate the isolation
1218 set CandidateInputArray PhotonCloner/photons
1219
1220 # isolation collection
1221 set IsolationInputArray EFlowFilterCHS/eflow
1222 set RhoInputArray Rho/rho
1223
1224 # output array
1225 set OutputArray photons
1226
1227 # isolation cone
1228 set DeltaRMax 0.3
1229
1230 # minimum pT
1231 set PTMin 0.0
1232
1233 # iso ratio to cut
1234 set PTRatioMax 9999.
1235
1236}
1237
1238
1239#####################
1240# Photon efficiency #
1241#####################
1242
1243module Efficiency PhotonEfficiency {
1244
1245 ## input particles
1246 set InputArray PhotonIsolation/photons
1247 ## output particles
1248 set OutputArray photons
1249 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1250 # efficiency formula for photons
1251 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
1252 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
1253 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
1254 (abs(eta) > 4.0) * (0.00)}
1255
1256}
1257
1258
1259#####################
1260# Photon efficiency #
1261#####################
1262
1263module Efficiency PhotonEfficiencyCHS {
1264
1265 ## input particles
1266 set InputArray PhotonIsolationCHS/photons
1267 ## output particles
1268 set OutputArray photons
1269 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1270 # efficiency formula for photons
1271 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
1272 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
1273 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
1274 (abs(eta) > 4.0) * (0.00)}
1275
1276}
1277
1278######################
1279# Electron isolation #
1280######################
1281
1282module Isolation ElectronIsolation {
1283
1284 set CandidateInputArray ElectronFilter/electrons
1285
1286 # isolation collection
1287 set IsolationInputArray EFlowFilterPuppi/eflow
1288
1289 set OutputArray electrons
1290
1291 set DeltaRMax 0.3
1292 set PTMin 0.0
1293 set PTRatioMax 9999.
1294
1295}
1296
1297
1298##########################
1299# Electron isolation CHS #
1300##########################
1301
1302module Isolation ElectronIsolationCHS {
1303
1304 set CandidateInputArray ElectronCloner/electrons
1305
1306 # isolation collection
1307 set IsolationInputArray EFlowFilterCHS/eflow
1308 set RhoInputArray Rho/rho
1309
1310 set OutputArray electrons
1311
1312 # veto isolation cand. based on proximity to input cand.
1313 set DeltaRMin 0.01
1314 set UseMiniCone true
1315
1316 set DeltaRMax 0.3
1317 set PTMin 0.0
1318 set PTRatioMax 9999.
1319
1320}
1321
1322
1323#######################
1324# Electron efficiency #
1325#######################
1326
1327module Efficiency ElectronEfficiency {
1328
1329 set InputArray ElectronIsolation/electrons
1330 set OutputArray electrons
1331
1332 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1333 # efficiency formula for electrons
1334 # efficiency for low pT leptons is set artifically to 100%. Analyzers should rescale according to proper lepton Id
1335 set EfficiencyFormula {
1336 (pt <= 4.0) * (0.00) +
1337 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 4.0 && pt <= 6.0) * (0.018) +
1338 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 6.0 && pt <= 8.0) * (0.252) +
1339 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 8.0 && pt <= 10.0) * (0.480) +
1340 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 10.0 && pt <= 20.0) * (0.681) +
1341 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 20.0 && pt <= 35.0) * (0.792) +
1342 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 35.0 && pt <= 50.0) * (0.862) +
1343 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 50.0 && pt <= 14000.0) * (0.859) +
1344 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 4.0 && pt <= 6.0) * (0.016) +
1345 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 6.0 && pt <= 8.0) * (0.198) +
1346 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 8.0 && pt <= 10.0) * (0.446) +
1347 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 10.0 && pt <= 20.0) * (0.598) +
1348 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 20.0 && pt <= 35.0) * (0.759) +
1349 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 35.0 && pt <= 50.0) * (0.847) +
1350 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 50.0 && pt <= 14000.0) * (0.872) +
1351 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 4.0 && pt <= 6.0) * (0.005) +
1352 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 6.0 && pt <= 8.0) * (0.029) +
1353 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 8.0 && pt <= 10.0) * (0.108) +
1354 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 10.0 && pt <= 20.0) * (0.289) +
1355 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 20.0 && pt <= 35.0) * (0.570) +
1356 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 35.0 && pt <= 50.0) * (0.743) +
1357 (abs(eta) > 1.0 && abs(eta) <= 1.45) * (pt > 50.0 && pt <= 14000.0) * (0.828) +
1358 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 6.0) * (0.026) +
1359 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 6.0 && pt <= 8.0) * (0.045) +
1360 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 8.0 && pt <= 10.0) * (0.133) +
1361 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 20.0) * (0.411) +
1362 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 20.0 && pt <= 35.0) * (0.629) +
1363 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 35.0 && pt <= 50.0) * (0.761) +
1364 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 50.0 && pt <= 14000.0) * (0.752) +
1365 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 4.0 && pt <= 6.0) * (0.061) +
1366 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 6.0 && pt <= 8.0) * (0.191) +
1367 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 8.0 && pt <= 10.0) * (0.337) +
1368 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.475) +
1369 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 35.0) * (0.605) +
1370 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 35.0 && pt <= 50.0) * (0.713) +
1371 (abs(eta) > 1.55 && abs(eta) <= 2.0) * (pt > 50.0 && pt <= 14000.0) * (0.794) +
1372 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 4.0 && pt <= 6.0) * (0.100) +
1373 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 6.0 && pt <= 8.0) * (0.223) +
1374 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 8.0 && pt <= 10.0) * (0.427) +
1375 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.590) +
1376 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 35.0) * (0.720) +
1377 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 35.0 && pt <= 50.0) * (0.800) +
1378 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 50.0 && pt <= 14000.0) * (0.840) +
1379 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 4.0 && pt <= 6.0) * (0.049) +
1380 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 6.0 && pt <= 8.0) * (0.152) +
1381 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 8.0 && pt <= 10.0) * (0.436) +
1382 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 20.0) * (0.679) +
1383 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 20.0 && pt <= 35.0) * (0.778) +
1384 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 35.0 && pt <= 50.0) * (0.830) +
1385 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 50.0 && pt <= 14000.0) * (0.919) +
1386 (abs(eta) > 3.0) * (0.00)
1387 }
1388}
1389
1390###########################
1391# Electron efficiency CHS #
1392###########################
1393
1394module Efficiency ElectronEfficiencyCHS {
1395
1396 set InputArray ElectronIsolationCHS/electrons
1397 set OutputArray electrons
1398
1399 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1400 # efficiency formula for electrons
1401 # efficiency for low pT leptons is set artifically to 100%. Analyzers should rescale according to proper lepton Id
1402 set EfficiencyFormula {
1403 (pt <= 4.0) * (1.00) + \
1404 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
1405 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
1406 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
1407 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
1408 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
1409 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
1410 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
1411 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
1412 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
1413 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
1414 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
1415 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
1416 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
1417 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
1418 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
1419 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
1420 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1421 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1422 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1423 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1424 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
1425 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1426 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1427 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1428 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1429 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
1430 (abs(eta) > 4.0) * (0.00)
1431
1432 }
1433}
1434
1435
1436
1437##################
1438# Muon isolation #
1439##################
1440
1441module Isolation MuonIsolation {
1442 set CandidateInputArray MuonMomentumSmearing/muons
1443
1444 # isolation collection
1445 set IsolationInputArray EFlowFilterPuppi/eflow
1446
1447 set OutputArray muons
1448
1449 set DeltaRMax 0.3
1450 set PTMin 0.0
1451 set PTRatioMax 9999.
1452
1453}
1454
1455######################
1456# Muon isolation CHS #
1457######################
1458
1459module Isolation MuonIsolationCHS {
1460 set CandidateInputArray MuonCloner/muons
1461
1462 # isolation collection
1463 set IsolationInputArray EFlowFilterCHS/eflow
1464 set RhoInputArray Rho/rho
1465
1466 set OutputArray muons
1467
1468
1469 # veto isolation cand. based on proximity to input cand.
1470 set DeltaRMin 0.01
1471 set UseMiniCone true
1472
1473 set DeltaRMax 0.3
1474 set PTMin 0.0
1475 set PTRatioMax 9999.
1476
1477}
1478
1479
1480#####################
1481# Muon Loose Id #
1482#####################
1483
1484module Efficiency MuonLooseIdEfficiency {
1485 set InputArray MuonIsolation/muons
1486 set OutputArray muons
1487 # tracking + LooseID efficiency formula for muons
1488 source muonLooseId.tcl
1489
1490}
1491
1492
1493##################
1494# Muon Tight Id #
1495##################
1496
1497module Efficiency MuonTightIdEfficiency {
1498 set InputArray MuonIsolation/muons
1499 set OutputArray muons
1500 # tracking + TightID efficiency formula for muons
1501 #source muonTightId.tcl
1502 set EfficiencyFormula {
1503 (pt <= 2.0) * (0.00) +
1504 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 2.0 && pt <= 4.0) * (0.04) +
1505 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 4.0 && pt <= 6.0) * (0.43) +
1506 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 6.0 && pt <= 8.0) * (0.53) +
1507 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 8.0 && pt <= 10.0) * (0.67) +
1508 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 10.0 && pt <= 20.0) * (0.81) +
1509 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 20.0 && pt <= 35.0) * (0.90) +
1510 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 35.0 && pt <= 50.0) * (0.92) +
1511 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 50.0 && pt <= 14000.0) * (0.90) +
1512 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 2.0 && pt <= 4.0) * (0.05) +
1513 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 4.0 && pt <= 6.0) * (0.46) +
1514 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 6.0 && pt <= 8.0) * (0.56) +
1515 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 8.0 && pt <= 10.0) * (0.65) +
1516 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 10.0 && pt <= 20.0) * (0.79) +
1517 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 20.0 && pt <= 35.0) * (0.91) +
1518 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 35.0 && pt <= 50.0) * (0.93) +
1519 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 50.0 && pt <= 14000.0) * (0.92) +
1520 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 2.0 && pt <= 4.0) * (0.15) +
1521 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 4.0 && pt <= 6.0) * (0.47) +
1522 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 6.0 && pt <= 8.0) * (0.55) +
1523 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 8.0 && pt <= 10.0) * (0.64) +
1524 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 10.0 && pt <= 20.0) * (0.78) +
1525 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 20.0 && pt <= 35.0) * (0.89) +
1526 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 35.0 && pt <= 50.0) * (0.96) +
1527 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1528 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 2.0 && pt <= 4.0) * (0.23) +
1529 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 4.0 && pt <= 6.0) * (0.44) +
1530 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 6.0 && pt <= 8.0) * (0.53) +
1531 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 8.0 && pt <= 10.0) * (0.68) +
1532 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.78) +
1533 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 35.0) * (0.89) +
1534 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 35.0 && pt <= 50.0) * (0.95) +
1535 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 50.0 && pt <= 14000.0) * (0.88) +
1536 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 2.0 && pt <= 4.0) * (0.22) +
1537 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 4.0 && pt <= 6.0) * (0.36) +
1538 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 6.0 && pt <= 8.0) * (0.44) +
1539 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 8.0 && pt <= 10.0) * (0.57) +
1540 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.63) +
1541 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 35.0) * (0.71) +
1542 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 35.0 && pt <= 50.0) * (0.76) +
1543 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 50.0 && pt <= 14000.0) * (0.82) +
1544 (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 2.0) * (0.80) +
1545 (abs(eta) > 2.8) * (0.00)
1546 }
1547}
1548
1549
1550#####################
1551# Muon Loose Id CHS #
1552#####################
1553
1554module Efficiency MuonLooseIdEfficiencyCHS {
1555 set InputArray MuonIsolationCHS/muons
1556 set OutputArray muons
1557 # tracking + LooseID efficiency formula for muons
1558 source muonLooseId.tcl
1559
1560}
1561
1562
1563######################
1564# Muon Tight Id CHS #
1565######################
1566
1567module Efficiency MuonTightIdEfficiencyCHS {
1568 set InputArray MuonIsolationCHS/muons
1569 set OutputArray muons
1570 # tracking + TightID efficiency formula for muons
1571 #source muonTightId.tcl
1572 set EfficiencyFormula {
1573 (pt <= 2.0) * (0.00) +
1574 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 2.0 && pt <= 4.0) * (0.04) +
1575 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 4.0 && pt <= 6.0) * (0.43) +
1576 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 6.0 && pt <= 8.0) * (0.53) +
1577 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 8.0 && pt <= 10.0) * (0.67) +
1578 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 10.0 && pt <= 20.0) * (0.81) +
1579 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 20.0 && pt <= 35.0) * (0.90) +
1580 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 35.0 && pt <= 50.0) * (0.92) +
1581 (abs(eta) > 0.0 && abs(eta) <= 0.5) * (pt > 50.0 && pt <= 14000.0) * (0.90) +
1582 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 2.0 && pt <= 4.0) * (0.05) +
1583 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 4.0 && pt <= 6.0) * (0.46) +
1584 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 6.0 && pt <= 8.0) * (0.56) +
1585 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 8.0 && pt <= 10.0) * (0.65) +
1586 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 10.0 && pt <= 20.0) * (0.79) +
1587 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 20.0 && pt <= 35.0) * (0.91) +
1588 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 35.0 && pt <= 50.0) * (0.93) +
1589 (abs(eta) > 0.5 && abs(eta) <= 1.0) * (pt > 50.0 && pt <= 14000.0) * (0.92) +
1590 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 2.0 && pt <= 4.0) * (0.15) +
1591 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 4.0 && pt <= 6.0) * (0.47) +
1592 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 6.0 && pt <= 8.0) * (0.55) +
1593 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 8.0 && pt <= 10.0) * (0.64) +
1594 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 10.0 && pt <= 20.0) * (0.78) +
1595 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 20.0 && pt <= 35.0) * (0.89) +
1596 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 35.0 && pt <= 50.0) * (0.96) +
1597 (abs(eta) > 1.0 && abs(eta) <= 1.5) * (pt > 50.0 && pt <= 14000.0) * (0.91) +
1598 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 2.0 && pt <= 4.0) * (0.23) +
1599 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 4.0 && pt <= 6.0) * (0.44) +
1600 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 6.0 && pt <= 8.0) * (0.53) +
1601 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 8.0 && pt <= 10.0) * (0.68) +
1602 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 10.0 && pt <= 20.0) * (0.78) +
1603 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 20.0 && pt <= 35.0) * (0.89) +
1604 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 35.0 && pt <= 50.0) * (0.95) +
1605 (abs(eta) > 1.5 && abs(eta) <= 2.0) * (pt > 50.0 && pt <= 14000.0) * (0.88) +
1606 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 2.0 && pt <= 4.0) * (0.22) +
1607 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 4.0 && pt <= 6.0) * (0.36) +
1608 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 6.0 && pt <= 8.0) * (0.44) +
1609 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 8.0 && pt <= 10.0) * (0.57) +
1610 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 10.0 && pt <= 20.0) * (0.63) +
1611 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 20.0 && pt <= 35.0) * (0.71) +
1612 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 35.0 && pt <= 50.0) * (0.76) +
1613 (abs(eta) > 2.0 && abs(eta) <= 2.5) * (pt > 50.0 && pt <= 14000.0) * (0.82) +
1614 (abs(eta) > 2.5 && abs(eta) <= 2.8) * (pt > 2.0) * (0.80) +
1615 (abs(eta) > 2.8) * (0.00)
1616 }
1617
1618
1619
1620
1621}
1622
1623
1624########################
1625# Jet Flavor Association
1626########################
1627
1628module JetFlavorAssociation JetFlavorAssociation {
1629
1630 set PartonInputArray Delphes/partons
1631 set ParticleInputArray Delphes/allParticles
1632 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1633 set JetInputArray JetEnergyScale/jets
1634
1635 set DeltaR 0.5
1636 set PartonPTMin 10.0
1637 set PartonEtaMax 4.0
1638
1639}
1640
1641module JetFlavorAssociation JetFlavorAssociationAK8 {
1642
1643 set PartonInputArray Delphes/partons
1644 set ParticleInputArray Delphes/allParticles
1645 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1646 set JetInputArray JetEnergyScaleAK8/jets
1647
1648 set DeltaR 0.8
1649 set PartonPTMin 100.0
1650 set PartonEtaMax 4.0
1651
1652}
1653
1654module JetFlavorAssociation JetFlavorAssociationPUPPI {
1655
1656 set PartonInputArray Delphes/partons
1657 set ParticleInputArray Delphes/allParticles
1658 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1659 set JetInputArray JetEnergyScalePUPPI/jets
1660
1661 set DeltaR 0.5
1662 set PartonPTMin 10.0
1663 set PartonEtaMax 4.0
1664
1665}
1666
1667module JetFlavorAssociation JetFlavorAssociationPUPPIAK8 {
1668
1669 set PartonInputArray Delphes/partons
1670 set ParticleInputArray Delphes/allParticles
1671 set ParticleLHEFInputArray Delphes/allParticlesLHEF
1672 set JetInputArray JetEnergyScalePUPPIAK8/jets
1673
1674 set DeltaR 0.8
1675 set PartonPTMin 100.0
1676 set PartonEtaMax 4.0
1677
1678}
1679
1680
1681#############
1682# b-tagging #
1683#############
1684module BTagging BTagging {
1685
1686 set JetInputArray JetEnergyScale/jets
1687 set BitNumber 0
1688
1689 add EfficiencyFormula {0} {0.001}
1690
1691 add EfficiencyFormula {5} {
1692 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 20.00 && pt <= 30.00) * (0.63) +
1693 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 30.00 && pt <= 40.00) * (0.70) +
1694 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 40.00 && pt <= 50.00) * (0.74) +
1695 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 50.00 && pt <= 60.00) * (0.74) +
1696 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 60.00 && pt <= 70.00) * (0.75) +
1697 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 70.00 && pt <= 80.00) * (0.75) +
1698 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 80.00 && pt <= 90.00) * (0.75) +
1699 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 90.00 && pt <= 100.00) * (0.74) +
1700 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 100.00 && pt <= 120.00) * (0.73) +
1701 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 120.00 && pt <= 140.00) * (0.73) +
1702 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 140.00 && pt <= 160.00) * (0.72) +
1703 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 160.00 && pt <= 180.00) * (0.69) +
1704 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 180.00 && pt <= 200.00) * (0.68) +
1705 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 200.00 && pt <= 250.00) * (0.66) +
1706 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 250.00 && pt <= 300.00) * (0.64) +
1707 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 300.00 && pt <= 350.00) * (0.59) +
1708 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 350.00 && pt <= 400.00) * (0.56) +
1709 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 400.00 && pt <= 500.00) * (0.50) +
1710 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 500.00 && pt <= 600.00) * (0.44) +
1711 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 600.00 && pt <= 700.00) * (0.40) +
1712 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 700.00 && pt <= 800.00) * (0.32) +
1713 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 800.00 && pt <= 1000.00) * (0.26) +
1714 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1000.00 && pt <= 1400.00) * (0.21) +
1715 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1400.00 && pt <= 2000.00) * (0.11) +
1716 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 2000.00 && pt <= 3000.00) * (0.17) +
1717 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 3000.00) * (0.00) +
1718 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 20.00 && pt <= 30.00) * (0.43) +
1719 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 30.00 && pt <= 40.00) * (0.53) +
1720 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 40.00 && pt <= 50.00) * (0.56) +
1721 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 50.00 && pt <= 60.00) * (0.59) +
1722 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 60.00 && pt <= 70.00) * (0.60) +
1723 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 70.00 && pt <= 80.00) * (0.61) +
1724 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 80.00 && pt <= 90.00) * (0.60) +
1725 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 90.00 && pt <= 100.00) * (0.59) +
1726 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 100.00 && pt <= 120.00) * (0.59) +
1727 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 120.00 && pt <= 140.00) * (0.60) +
1728 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 140.00 && pt <= 160.00) * (0.58) +
1729 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 160.00 && pt <= 180.00) * (0.56) +
1730 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 180.00 && pt <= 200.00) * (0.55) +
1731 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 200.00 && pt <= 250.00) * (0.53) +
1732 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 250.00 && pt <= 300.00) * (0.49) +
1733 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 300.00 && pt <= 350.00) * (0.45) +
1734 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 350.00 && pt <= 400.00) * (0.42) +
1735 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 400.00 && pt <= 500.00) * (0.38) +
1736 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 500.00 && pt <= 600.00) * (0.32) +
1737 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 600.00 && pt <= 700.00) * (0.36) +
1738 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 700.00 && pt <= 800.00) * (0.34) +
1739 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 800.00 && pt <= 1000.00) * (0.29) +
1740 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1000.00 && pt <= 1400.00) * (0.20) +
1741 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1400.00 && pt <= 2000.00) * (0.00) +
1742 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 2000.00 && pt <= 3000.00) * (0.00) +
1743 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 3000.00) * (0.00) +
1744 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 20.00 && pt <= 30.00) * (0.25) +
1745 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 30.00 && pt <= 40.00) * (0.33) +
1746 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 40.00 && pt <= 50.00) * (0.37) +
1747 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 50.00 && pt <= 60.00) * (0.39) +
1748 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 60.00 && pt <= 70.00) * (0.41) +
1749 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 70.00 && pt <= 80.00) * (0.43) +
1750 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 80.00 && pt <= 90.00) * (0.41) +
1751 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 90.00 && pt <= 100.00) * (0.44) +
1752 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 100.00 && pt <= 120.00) * (0.44) +
1753 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 120.00 && pt <= 140.00) * (0.42) +
1754 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 140.00 && pt <= 160.00) * (0.40) +
1755 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 160.00 && pt <= 180.00) * (0.41) +
1756 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 180.00 && pt <= 200.00) * (0.39) +
1757 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 200.00 && pt <= 250.00) * (0.34) +
1758 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 250.00 && pt <= 300.00) * (0.30) +
1759 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 300.00 && pt <= 350.00) * (0.23) +
1760 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 350.00 && pt <= 400.00) * (0.33) +
1761 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 400.00 && pt <= 500.00) * (0.14) +
1762 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 500.00 && pt <= 600.00) * (0.00) +
1763 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 600.00 && pt <= 700.00) * (0.00) +
1764 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 700.00 && pt <= 800.00) * (0.00) +
1765 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 800.00 && pt <= 1000.00) * (0.00) +
1766 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1000.00 && pt <= 1400.00) * (0.00) +
1767 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1400.00 && pt <= 2000.00) * (0.00) +
1768 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 2000.00 && pt <= 3000.00) * (0.00) +
1769 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 3000.00) * (0.00)
1770 }
1771
1772 add EfficiencyFormula {4} {
1773 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 20.00 && pt <= 30.00) * (0.17) +
1774 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 30.00 && pt <= 40.00) * (0.19) +
1775 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 40.00 && pt <= 50.00) * (0.20) +
1776 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 50.00 && pt <= 60.00) * (0.21) +
1777 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 60.00 && pt <= 70.00) * (0.20) +
1778 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 70.00 && pt <= 80.00) * (0.20) +
1779 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 80.00 && pt <= 90.00) * (0.20) +
1780 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 90.00 && pt <= 100.00) * (0.19) +
1781 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 100.00 && pt <= 120.00) * (0.18) +
1782 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 120.00 && pt <= 140.00) * (0.18) +
1783 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 140.00 && pt <= 160.00) * (0.17) +
1784 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 160.00 && pt <= 180.00) * (0.16) +
1785 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 180.00 && pt <= 200.00) * (0.16) +
1786 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 200.00 && pt <= 250.00) * (0.15) +
1787 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 250.00 && pt <= 300.00) * (0.12) +
1788 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 300.00 && pt <= 350.00) * (0.12) +
1789 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 350.00 && pt <= 400.00) * (0.10) +
1790 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 400.00 && pt <= 500.00) * (0.08) +
1791 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 500.00 && pt <= 600.00) * (0.07) +
1792 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 600.00 && pt <= 700.00) * (0.06) +
1793 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 700.00 && pt <= 800.00) * (0.05) +
1794 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 800.00 && pt <= 1000.00) * (0.04) +
1795 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1000.00 && pt <= 1400.00) * (0.03) +
1796 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1400.00 && pt <= 2000.00) * (0.02) +
1797 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 2000.00 && pt <= 3000.00) * (0.02) +
1798 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 3000.00) * (0.00) +
1799 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 20.00 && pt <= 30.00) * (0.08) +
1800 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 30.00 && pt <= 40.00) * (0.10) +
1801 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 40.00 && pt <= 50.00) * (0.10) +
1802 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 50.00 && pt <= 60.00) * (0.12) +
1803 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 60.00 && pt <= 70.00) * (0.12) +
1804 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 70.00 && pt <= 80.00) * (0.13) +
1805 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 80.00 && pt <= 90.00) * (0.13) +
1806 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 90.00 && pt <= 100.00) * (0.14) +
1807 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 100.00 && pt <= 120.00) * (0.14) +
1808 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 120.00 && pt <= 140.00) * (0.14) +
1809 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 140.00 && pt <= 160.00) * (0.15) +
1810 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 160.00 && pt <= 180.00) * (0.15) +
1811 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 180.00 && pt <= 200.00) * (0.15) +
1812 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 200.00 && pt <= 250.00) * (0.15) +
1813 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 250.00 && pt <= 300.00) * (0.13) +
1814 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 300.00 && pt <= 350.00) * (0.12) +
1815 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 350.00 && pt <= 400.00) * (0.11) +
1816 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 400.00 && pt <= 500.00) * (0.09) +
1817 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 500.00 && pt <= 600.00) * (0.07) +
1818 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 600.00 && pt <= 700.00) * (0.07) +
1819 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 700.00 && pt <= 800.00) * (0.06) +
1820 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 800.00 && pt <= 1000.00) * (0.03) +
1821 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1000.00 && pt <= 1400.00) * (0.03) +
1822 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1400.00 && pt <= 2000.00) * (0.00) +
1823 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 2000.00 && pt <= 3000.00) * (0.00) +
1824 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 3000.00) * (0.00) +
1825 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 20.00 && pt <= 30.00) * (0.06) +
1826 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 30.00 && pt <= 40.00) * (0.07) +
1827 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 40.00 && pt <= 50.00) * (0.08) +
1828 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 50.00 && pt <= 60.00) * (0.08) +
1829 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 60.00 && pt <= 70.00) * (0.09) +
1830 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 70.00 && pt <= 80.00) * (0.09) +
1831 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 80.00 && pt <= 90.00) * (0.10) +
1832 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 90.00 && pt <= 100.00) * (0.09) +
1833 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 100.00 && pt <= 120.00) * (0.10) +
1834 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 120.00 && pt <= 140.00) * (0.10) +
1835 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 140.00 && pt <= 160.00) * (0.10) +
1836 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 160.00 && pt <= 180.00) * (0.09) +
1837 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 180.00 && pt <= 200.00) * (0.09) +
1838 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 200.00 && pt <= 250.00) * (0.10) +
1839 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 250.00 && pt <= 300.00) * (0.08) +
1840 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 300.00 && pt <= 350.00) * (0.07) +
1841 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 350.00 && pt <= 400.00) * (0.07) +
1842 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 400.00 && pt <= 500.00) * (0.07) +
1843 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 500.00 && pt <= 600.00) * (0.00) +
1844 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 600.00 && pt <= 700.00) * (0.00) +
1845 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 700.00 && pt <= 800.00) * (0.00) +
1846 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 800.00 && pt <= 1000.00) * (0.00) +
1847 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1000.00 && pt <= 1400.00) * (0.00) +
1848 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1400.00 && pt <= 2000.00) * (0.00) +
1849 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 2000.00 && pt <= 3000.00) * (0.00) +
1850 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 3000.00) * (0.00)
1851 }
1852
1853
1854
1855}
1856
1857module BTagging BTaggingAK8 {
1858
1859 set JetInputArray JetEnergyScaleAK8/jets
1860 set BitNumber 0
1861 add EfficiencyFormula {0} {0.001}
1862
1863 add EfficiencyFormula {5} {
1864 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 20.00 && pt <= 30.00) * (0.63) +
1865 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 30.00 && pt <= 40.00) * (0.70) +
1866 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 40.00 && pt <= 50.00) * (0.74) +
1867 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 50.00 && pt <= 60.00) * (0.74) +
1868 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 60.00 && pt <= 70.00) * (0.75) +
1869 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 70.00 && pt <= 80.00) * (0.75) +
1870 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 80.00 && pt <= 90.00) * (0.75) +
1871 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 90.00 && pt <= 100.00) * (0.74) +
1872 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 100.00 && pt <= 120.00) * (0.73) +
1873 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 120.00 && pt <= 140.00) * (0.73) +
1874 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 140.00 && pt <= 160.00) * (0.72) +
1875 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 160.00 && pt <= 180.00) * (0.69) +
1876 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 180.00 && pt <= 200.00) * (0.68) +
1877 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 200.00 && pt <= 250.00) * (0.66) +
1878 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 250.00 && pt <= 300.00) * (0.64) +
1879 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 300.00 && pt <= 350.00) * (0.59) +
1880 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 350.00 && pt <= 400.00) * (0.56) +
1881 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 400.00 && pt <= 500.00) * (0.50) +
1882 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 500.00 && pt <= 600.00) * (0.44) +
1883 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 600.00 && pt <= 700.00) * (0.40) +
1884 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 700.00 && pt <= 800.00) * (0.32) +
1885 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 800.00 && pt <= 1000.00) * (0.26) +
1886 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1000.00 && pt <= 1400.00) * (0.21) +
1887 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1400.00 && pt <= 2000.00) * (0.11) +
1888 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 2000.00 && pt <= 3000.00) * (0.17) +
1889 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 3000.00) * (0.00) +
1890 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 20.00 && pt <= 30.00) * (0.43) +
1891 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 30.00 && pt <= 40.00) * (0.53) +
1892 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 40.00 && pt <= 50.00) * (0.56) +
1893 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 50.00 && pt <= 60.00) * (0.59) +
1894 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 60.00 && pt <= 70.00) * (0.60) +
1895 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 70.00 && pt <= 80.00) * (0.61) +
1896 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 80.00 && pt <= 90.00) * (0.60) +
1897 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 90.00 && pt <= 100.00) * (0.59) +
1898 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 100.00 && pt <= 120.00) * (0.59) +
1899 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 120.00 && pt <= 140.00) * (0.60) +
1900 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 140.00 && pt <= 160.00) * (0.58) +
1901 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 160.00 && pt <= 180.00) * (0.56) +
1902 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 180.00 && pt <= 200.00) * (0.55) +
1903 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 200.00 && pt <= 250.00) * (0.53) +
1904 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 250.00 && pt <= 300.00) * (0.49) +
1905 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 300.00 && pt <= 350.00) * (0.45) +
1906 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 350.00 && pt <= 400.00) * (0.42) +
1907 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 400.00 && pt <= 500.00) * (0.38) +
1908 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 500.00 && pt <= 600.00) * (0.32) +
1909 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 600.00 && pt <= 700.00) * (0.36) +
1910 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 700.00 && pt <= 800.00) * (0.34) +
1911 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 800.00 && pt <= 1000.00) * (0.29) +
1912 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1000.00 && pt <= 1400.00) * (0.20) +
1913 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1400.00 && pt <= 2000.00) * (0.00) +
1914 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 2000.00 && pt <= 3000.00) * (0.00) +
1915 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 3000.00) * (0.00) +
1916 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 20.00 && pt <= 30.00) * (0.25) +
1917 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 30.00 && pt <= 40.00) * (0.33) +
1918 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 40.00 && pt <= 50.00) * (0.37) +
1919 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 50.00 && pt <= 60.00) * (0.39) +
1920 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 60.00 && pt <= 70.00) * (0.41) +
1921 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 70.00 && pt <= 80.00) * (0.43) +
1922 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 80.00 && pt <= 90.00) * (0.41) +
1923 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 90.00 && pt <= 100.00) * (0.44) +
1924 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 100.00 && pt <= 120.00) * (0.44) +
1925 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 120.00 && pt <= 140.00) * (0.42) +
1926 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 140.00 && pt <= 160.00) * (0.40) +
1927 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 160.00 && pt <= 180.00) * (0.41) +
1928 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 180.00 && pt <= 200.00) * (0.39) +
1929 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 200.00 && pt <= 250.00) * (0.34) +
1930 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 250.00 && pt <= 300.00) * (0.30) +
1931 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 300.00 && pt <= 350.00) * (0.23) +
1932 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 350.00 && pt <= 400.00) * (0.33) +
1933 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 400.00 && pt <= 500.00) * (0.14) +
1934 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 500.00 && pt <= 600.00) * (0.00) +
1935 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 600.00 && pt <= 700.00) * (0.00) +
1936 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 700.00 && pt <= 800.00) * (0.00) +
1937 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 800.00 && pt <= 1000.00) * (0.00) +
1938 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1000.00 && pt <= 1400.00) * (0.00) +
1939 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1400.00 && pt <= 2000.00) * (0.00) +
1940 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 2000.00 && pt <= 3000.00) * (0.00) +
1941 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 3000.00) * (0.00)
1942 }
1943
1944 add EfficiencyFormula {4} {
1945 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 20.00 && pt <= 30.00) * (0.17) +
1946 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 30.00 && pt <= 40.00) * (0.19) +
1947 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 40.00 && pt <= 50.00) * (0.20) +
1948 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 50.00 && pt <= 60.00) * (0.21) +
1949 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 60.00 && pt <= 70.00) * (0.20) +
1950 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 70.00 && pt <= 80.00) * (0.20) +
1951 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 80.00 && pt <= 90.00) * (0.20) +
1952 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 90.00 && pt <= 100.00) * (0.19) +
1953 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 100.00 && pt <= 120.00) * (0.18) +
1954 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 120.00 && pt <= 140.00) * (0.18) +
1955 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 140.00 && pt <= 160.00) * (0.17) +
1956 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 160.00 && pt <= 180.00) * (0.16) +
1957 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 180.00 && pt <= 200.00) * (0.16) +
1958 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 200.00 && pt <= 250.00) * (0.15) +
1959 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 250.00 && pt <= 300.00) * (0.12) +
1960 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 300.00 && pt <= 350.00) * (0.12) +
1961 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 350.00 && pt <= 400.00) * (0.10) +
1962 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 400.00 && pt <= 500.00) * (0.08) +
1963 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 500.00 && pt <= 600.00) * (0.07) +
1964 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 600.00 && pt <= 700.00) * (0.06) +
1965 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 700.00 && pt <= 800.00) * (0.05) +
1966 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 800.00 && pt <= 1000.00) * (0.04) +
1967 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1000.00 && pt <= 1400.00) * (0.03) +
1968 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1400.00 && pt <= 2000.00) * (0.02) +
1969 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 2000.00 && pt <= 3000.00) * (0.02) +
1970 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 3000.00) * (0.00) +
1971 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 20.00 && pt <= 30.00) * (0.08) +
1972 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 30.00 && pt <= 40.00) * (0.10) +
1973 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 40.00 && pt <= 50.00) * (0.10) +
1974 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 50.00 && pt <= 60.00) * (0.12) +
1975 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 60.00 && pt <= 70.00) * (0.12) +
1976 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 70.00 && pt <= 80.00) * (0.13) +
1977 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 80.00 && pt <= 90.00) * (0.13) +
1978 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 90.00 && pt <= 100.00) * (0.14) +
1979 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 100.00 && pt <= 120.00) * (0.14) +
1980 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 120.00 && pt <= 140.00) * (0.14) +
1981 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 140.00 && pt <= 160.00) * (0.15) +
1982 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 160.00 && pt <= 180.00) * (0.15) +
1983 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 180.00 && pt <= 200.00) * (0.15) +
1984 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 200.00 && pt <= 250.00) * (0.15) +
1985 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 250.00 && pt <= 300.00) * (0.13) +
1986 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 300.00 && pt <= 350.00) * (0.12) +
1987 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 350.00 && pt <= 400.00) * (0.11) +
1988 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 400.00 && pt <= 500.00) * (0.09) +
1989 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 500.00 && pt <= 600.00) * (0.07) +
1990 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 600.00 && pt <= 700.00) * (0.07) +
1991 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 700.00 && pt <= 800.00) * (0.06) +
1992 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 800.00 && pt <= 1000.00) * (0.03) +
1993 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1000.00 && pt <= 1400.00) * (0.03) +
1994 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1400.00 && pt <= 2000.00) * (0.00) +
1995 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 2000.00 && pt <= 3000.00) * (0.00) +
1996 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 3000.00) * (0.00) +
1997 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 20.00 && pt <= 30.00) * (0.06) +
1998 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 30.00 && pt <= 40.00) * (0.07) +
1999 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 40.00 && pt <= 50.00) * (0.08) +
2000 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 50.00 && pt <= 60.00) * (0.08) +
2001 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 60.00 && pt <= 70.00) * (0.09) +
2002 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 70.00 && pt <= 80.00) * (0.09) +
2003 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 80.00 && pt <= 90.00) * (0.10) +
2004 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 90.00 && pt <= 100.00) * (0.09) +
2005 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 100.00 && pt <= 120.00) * (0.10) +
2006 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 120.00 && pt <= 140.00) * (0.10) +
2007 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 140.00 && pt <= 160.00) * (0.10) +
2008 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 160.00 && pt <= 180.00) * (0.09) +
2009 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 180.00 && pt <= 200.00) * (0.09) +
2010 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 200.00 && pt <= 250.00) * (0.10) +
2011 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 250.00 && pt <= 300.00) * (0.08) +
2012 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 300.00 && pt <= 350.00) * (0.07) +
2013 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 350.00 && pt <= 400.00) * (0.07) +
2014 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 400.00 && pt <= 500.00) * (0.07) +
2015 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 500.00 && pt <= 600.00) * (0.00) +
2016 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 600.00 && pt <= 700.00) * (0.00) +
2017 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 700.00 && pt <= 800.00) * (0.00) +
2018 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 800.00 && pt <= 1000.00) * (0.00) +
2019 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1000.00 && pt <= 1400.00) * (0.00) +
2020 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1400.00 && pt <= 2000.00) * (0.00) +
2021 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 2000.00 && pt <= 3000.00) * (0.00) +
2022 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 3000.00) * (0.00)
2023 }
2024
2025
2026}
2027
2028module BTagging BTaggingPUPPI {
2029
2030 set JetInputArray JetEnergyScalePUPPI/jets
2031 set BitNumber 0
2032 add EfficiencyFormula {0} {0.001}
2033
2034 add EfficiencyFormula {5} {
2035 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 20.00 && pt <= 30.00) * (0.63) +
2036 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 30.00 && pt <= 40.00) * (0.70) +
2037 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 40.00 && pt <= 50.00) * (0.74) +
2038 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 50.00 && pt <= 60.00) * (0.74) +
2039 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 60.00 && pt <= 70.00) * (0.75) +
2040 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 70.00 && pt <= 80.00) * (0.75) +
2041 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 80.00 && pt <= 90.00) * (0.75) +
2042 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 90.00 && pt <= 100.00) * (0.74) +
2043 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 100.00 && pt <= 120.00) * (0.73) +
2044 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 120.00 && pt <= 140.00) * (0.73) +
2045 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 140.00 && pt <= 160.00) * (0.72) +
2046 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 160.00 && pt <= 180.00) * (0.69) +
2047 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 180.00 && pt <= 200.00) * (0.68) +
2048 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 200.00 && pt <= 250.00) * (0.66) +
2049 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 250.00 && pt <= 300.00) * (0.64) +
2050 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 300.00 && pt <= 350.00) * (0.59) +
2051 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 350.00 && pt <= 400.00) * (0.56) +
2052 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 400.00 && pt <= 500.00) * (0.50) +
2053 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 500.00 && pt <= 600.00) * (0.44) +
2054 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 600.00 && pt <= 700.00) * (0.40) +
2055 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 700.00 && pt <= 800.00) * (0.32) +
2056 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 800.00 && pt <= 1000.00) * (0.26) +
2057 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1000.00 && pt <= 1400.00) * (0.21) +
2058 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1400.00 && pt <= 2000.00) * (0.11) +
2059 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 2000.00 && pt <= 3000.00) * (0.17) +
2060 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 3000.00) * (0.00) +
2061 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 20.00 && pt <= 30.00) * (0.43) +
2062 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 30.00 && pt <= 40.00) * (0.53) +
2063 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 40.00 && pt <= 50.00) * (0.56) +
2064 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 50.00 && pt <= 60.00) * (0.59) +
2065 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 60.00 && pt <= 70.00) * (0.60) +
2066 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 70.00 && pt <= 80.00) * (0.61) +
2067 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 80.00 && pt <= 90.00) * (0.60) +
2068 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 90.00 && pt <= 100.00) * (0.59) +
2069 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 100.00 && pt <= 120.00) * (0.59) +
2070 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 120.00 && pt <= 140.00) * (0.60) +
2071 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 140.00 && pt <= 160.00) * (0.58) +
2072 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 160.00 && pt <= 180.00) * (0.56) +
2073 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 180.00 && pt <= 200.00) * (0.55) +
2074 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 200.00 && pt <= 250.00) * (0.53) +
2075 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 250.00 && pt <= 300.00) * (0.49) +
2076 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 300.00 && pt <= 350.00) * (0.45) +
2077 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 350.00 && pt <= 400.00) * (0.42) +
2078 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 400.00 && pt <= 500.00) * (0.38) +
2079 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 500.00 && pt <= 600.00) * (0.32) +
2080 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 600.00 && pt <= 700.00) * (0.36) +
2081 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 700.00 && pt <= 800.00) * (0.34) +
2082 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 800.00 && pt <= 1000.00) * (0.29) +
2083 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1000.00 && pt <= 1400.00) * (0.20) +
2084 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1400.00 && pt <= 2000.00) * (0.00) +
2085 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 2000.00 && pt <= 3000.00) * (0.00) +
2086 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 3000.00) * (0.00) +
2087 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 20.00 && pt <= 30.00) * (0.25) +
2088 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 30.00 && pt <= 40.00) * (0.33) +
2089 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 40.00 && pt <= 50.00) * (0.37) +
2090 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 50.00 && pt <= 60.00) * (0.39) +
2091 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 60.00 && pt <= 70.00) * (0.41) +
2092 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 70.00 && pt <= 80.00) * (0.43) +
2093 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 80.00 && pt <= 90.00) * (0.41) +
2094 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 90.00 && pt <= 100.00) * (0.44) +
2095 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 100.00 && pt <= 120.00) * (0.44) +
2096 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 120.00 && pt <= 140.00) * (0.42) +
2097 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 140.00 && pt <= 160.00) * (0.40) +
2098 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 160.00 && pt <= 180.00) * (0.41) +
2099 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 180.00 && pt <= 200.00) * (0.39) +
2100 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 200.00 && pt <= 250.00) * (0.34) +
2101 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 250.00 && pt <= 300.00) * (0.30) +
2102 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 300.00 && pt <= 350.00) * (0.23) +
2103 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 350.00 && pt <= 400.00) * (0.33) +
2104 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 400.00 && pt <= 500.00) * (0.14) +
2105 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 500.00 && pt <= 600.00) * (0.00) +
2106 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 600.00 && pt <= 700.00) * (0.00) +
2107 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 700.00 && pt <= 800.00) * (0.00) +
2108 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 800.00 && pt <= 1000.00) * (0.00) +
2109 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1000.00 && pt <= 1400.00) * (0.00) +
2110 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1400.00 && pt <= 2000.00) * (0.00) +
2111 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 2000.00 && pt <= 3000.00) * (0.00) +
2112 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 3000.00) * (0.00)
2113 }
2114
2115 add EfficiencyFormula {4} {
2116 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 20.00 && pt <= 30.00) * (0.17) +
2117 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 30.00 && pt <= 40.00) * (0.19) +
2118 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 40.00 && pt <= 50.00) * (0.20) +
2119 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 50.00 && pt <= 60.00) * (0.21) +
2120 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 60.00 && pt <= 70.00) * (0.20) +
2121 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 70.00 && pt <= 80.00) * (0.20) +
2122 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 80.00 && pt <= 90.00) * (0.20) +
2123 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 90.00 && pt <= 100.00) * (0.19) +
2124 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 100.00 && pt <= 120.00) * (0.18) +
2125 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 120.00 && pt <= 140.00) * (0.18) +
2126 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 140.00 && pt <= 160.00) * (0.17) +
2127 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 160.00 && pt <= 180.00) * (0.16) +
2128 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 180.00 && pt <= 200.00) * (0.16) +
2129 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 200.00 && pt <= 250.00) * (0.15) +
2130 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 250.00 && pt <= 300.00) * (0.12) +
2131 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 300.00 && pt <= 350.00) * (0.12) +
2132 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 350.00 && pt <= 400.00) * (0.10) +
2133 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 400.00 && pt <= 500.00) * (0.08) +
2134 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 500.00 && pt <= 600.00) * (0.07) +
2135 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 600.00 && pt <= 700.00) * (0.06) +
2136 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 700.00 && pt <= 800.00) * (0.05) +
2137 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 800.00 && pt <= 1000.00) * (0.04) +
2138 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1000.00 && pt <= 1400.00) * (0.03) +
2139 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1400.00 && pt <= 2000.00) * (0.02) +
2140 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 2000.00 && pt <= 3000.00) * (0.02) +
2141 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 3000.00) * (0.00) +
2142 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 20.00 && pt <= 30.00) * (0.08) +
2143 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 30.00 && pt <= 40.00) * (0.10) +
2144 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 40.00 && pt <= 50.00) * (0.10) +
2145 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 50.00 && pt <= 60.00) * (0.12) +
2146 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 60.00 && pt <= 70.00) * (0.12) +
2147 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 70.00 && pt <= 80.00) * (0.13) +
2148 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 80.00 && pt <= 90.00) * (0.13) +
2149 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 90.00 && pt <= 100.00) * (0.14) +
2150 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 100.00 && pt <= 120.00) * (0.14) +
2151 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 120.00 && pt <= 140.00) * (0.14) +
2152 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 140.00 && pt <= 160.00) * (0.15) +
2153 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 160.00 && pt <= 180.00) * (0.15) +
2154 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 180.00 && pt <= 200.00) * (0.15) +
2155 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 200.00 && pt <= 250.00) * (0.15) +
2156 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 250.00 && pt <= 300.00) * (0.13) +
2157 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 300.00 && pt <= 350.00) * (0.12) +
2158 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 350.00 && pt <= 400.00) * (0.11) +
2159 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 400.00 && pt <= 500.00) * (0.09) +
2160 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 500.00 && pt <= 600.00) * (0.07) +
2161 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 600.00 && pt <= 700.00) * (0.07) +
2162 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 700.00 && pt <= 800.00) * (0.06) +
2163 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 800.00 && pt <= 1000.00) * (0.03) +
2164 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1000.00 && pt <= 1400.00) * (0.03) +
2165 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1400.00 && pt <= 2000.00) * (0.00) +
2166 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 2000.00 && pt <= 3000.00) * (0.00) +
2167 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 3000.00) * (0.00) +
2168 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 20.00 && pt <= 30.00) * (0.06) +
2169 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 30.00 && pt <= 40.00) * (0.07) +
2170 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 40.00 && pt <= 50.00) * (0.08) +
2171 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 50.00 && pt <= 60.00) * (0.08) +
2172 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 60.00 && pt <= 70.00) * (0.09) +
2173 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 70.00 && pt <= 80.00) * (0.09) +
2174 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 80.00 && pt <= 90.00) * (0.10) +
2175 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 90.00 && pt <= 100.00) * (0.09) +
2176 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 100.00 && pt <= 120.00) * (0.10) +
2177 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 120.00 && pt <= 140.00) * (0.10) +
2178 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 140.00 && pt <= 160.00) * (0.10) +
2179 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 160.00 && pt <= 180.00) * (0.09) +
2180 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 180.00 && pt <= 200.00) * (0.09) +
2181 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 200.00 && pt <= 250.00) * (0.10) +
2182 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 250.00 && pt <= 300.00) * (0.08) +
2183 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 300.00 && pt <= 350.00) * (0.07) +
2184 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 350.00 && pt <= 400.00) * (0.07) +
2185 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 400.00 && pt <= 500.00) * (0.07) +
2186 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 500.00 && pt <= 600.00) * (0.00) +
2187 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 600.00 && pt <= 700.00) * (0.00) +
2188 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 700.00 && pt <= 800.00) * (0.00) +
2189 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 800.00 && pt <= 1000.00) * (0.00) +
2190 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1000.00 && pt <= 1400.00) * (0.00) +
2191 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1400.00 && pt <= 2000.00) * (0.00) +
2192 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 2000.00 && pt <= 3000.00) * (0.00) +
2193 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 3000.00) * (0.00)
2194 }
2195
2196
2197}
2198
2199module BTagging BTaggingPUPPIAK8 {
2200
2201 set JetInputArray JetEnergyScalePUPPIAK8/jets
2202 set BitNumber 0
2203 add EfficiencyFormula {0} {0.001}
2204
2205 add EfficiencyFormula {5} {
2206 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 20.00 && pt <= 30.00) * (0.63) +
2207 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 30.00 && pt <= 40.00) * (0.70) +
2208 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 40.00 && pt <= 50.00) * (0.74) +
2209 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 50.00 && pt <= 60.00) * (0.74) +
2210 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 60.00 && pt <= 70.00) * (0.75) +
2211 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 70.00 && pt <= 80.00) * (0.75) +
2212 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 80.00 && pt <= 90.00) * (0.75) +
2213 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 90.00 && pt <= 100.00) * (0.74) +
2214 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 100.00 && pt <= 120.00) * (0.73) +
2215 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 120.00 && pt <= 140.00) * (0.73) +
2216 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 140.00 && pt <= 160.00) * (0.72) +
2217 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 160.00 && pt <= 180.00) * (0.69) +
2218 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 180.00 && pt <= 200.00) * (0.68) +
2219 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 200.00 && pt <= 250.00) * (0.66) +
2220 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 250.00 && pt <= 300.00) * (0.64) +
2221 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 300.00 && pt <= 350.00) * (0.59) +
2222 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 350.00 && pt <= 400.00) * (0.56) +
2223 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 400.00 && pt <= 500.00) * (0.50) +
2224 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 500.00 && pt <= 600.00) * (0.44) +
2225 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 600.00 && pt <= 700.00) * (0.40) +
2226 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 700.00 && pt <= 800.00) * (0.32) +
2227 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 800.00 && pt <= 1000.00) * (0.26) +
2228 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1000.00 && pt <= 1400.00) * (0.21) +
2229 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1400.00 && pt <= 2000.00) * (0.11) +
2230 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 2000.00 && pt <= 3000.00) * (0.17) +
2231 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 3000.00) * (0.00) +
2232 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 20.00 && pt <= 30.00) * (0.43) +
2233 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 30.00 && pt <= 40.00) * (0.53) +
2234 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 40.00 && pt <= 50.00) * (0.56) +
2235 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 50.00 && pt <= 60.00) * (0.59) +
2236 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 60.00 && pt <= 70.00) * (0.60) +
2237 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 70.00 && pt <= 80.00) * (0.61) +
2238 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 80.00 && pt <= 90.00) * (0.60) +
2239 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 90.00 && pt <= 100.00) * (0.59) +
2240 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 100.00 && pt <= 120.00) * (0.59) +
2241 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 120.00 && pt <= 140.00) * (0.60) +
2242 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 140.00 && pt <= 160.00) * (0.58) +
2243 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 160.00 && pt <= 180.00) * (0.56) +
2244 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 180.00 && pt <= 200.00) * (0.55) +
2245 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 200.00 && pt <= 250.00) * (0.53) +
2246 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 250.00 && pt <= 300.00) * (0.49) +
2247 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 300.00 && pt <= 350.00) * (0.45) +
2248 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 350.00 && pt <= 400.00) * (0.42) +
2249 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 400.00 && pt <= 500.00) * (0.38) +
2250 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 500.00 && pt <= 600.00) * (0.32) +
2251 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 600.00 && pt <= 700.00) * (0.36) +
2252 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 700.00 && pt <= 800.00) * (0.34) +
2253 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 800.00 && pt <= 1000.00) * (0.29) +
2254 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1000.00 && pt <= 1400.00) * (0.20) +
2255 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1400.00 && pt <= 2000.00) * (0.00) +
2256 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 2000.00 && pt <= 3000.00) * (0.00) +
2257 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 3000.00) * (0.00) +
2258 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 20.00 && pt <= 30.00) * (0.25) +
2259 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 30.00 && pt <= 40.00) * (0.33) +
2260 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 40.00 && pt <= 50.00) * (0.37) +
2261 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 50.00 && pt <= 60.00) * (0.39) +
2262 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 60.00 && pt <= 70.00) * (0.41) +
2263 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 70.00 && pt <= 80.00) * (0.43) +
2264 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 80.00 && pt <= 90.00) * (0.41) +
2265 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 90.00 && pt <= 100.00) * (0.44) +
2266 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 100.00 && pt <= 120.00) * (0.44) +
2267 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 120.00 && pt <= 140.00) * (0.42) +
2268 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 140.00 && pt <= 160.00) * (0.40) +
2269 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 160.00 && pt <= 180.00) * (0.41) +
2270 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 180.00 && pt <= 200.00) * (0.39) +
2271 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 200.00 && pt <= 250.00) * (0.34) +
2272 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 250.00 && pt <= 300.00) * (0.30) +
2273 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 300.00 && pt <= 350.00) * (0.23) +
2274 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 350.00 && pt <= 400.00) * (0.33) +
2275 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 400.00 && pt <= 500.00) * (0.14) +
2276 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 500.00 && pt <= 600.00) * (0.00) +
2277 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 600.00 && pt <= 700.00) * (0.00) +
2278 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 700.00 && pt <= 800.00) * (0.00) +
2279 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 800.00 && pt <= 1000.00) * (0.00) +
2280 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1000.00 && pt <= 1400.00) * (0.00) +
2281 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1400.00 && pt <= 2000.00) * (0.00) +
2282 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 2000.00 && pt <= 3000.00) * (0.00) +
2283 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 3000.00) * (0.00)
2284 }
2285
2286 add EfficiencyFormula {4} {
2287 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 20.00 && pt <= 30.00) * (0.17) +
2288 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 30.00 && pt <= 40.00) * (0.19) +
2289 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 40.00 && pt <= 50.00) * (0.20) +
2290 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 50.00 && pt <= 60.00) * (0.21) +
2291 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 60.00 && pt <= 70.00) * (0.20) +
2292 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 70.00 && pt <= 80.00) * (0.20) +
2293 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 80.00 && pt <= 90.00) * (0.20) +
2294 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 90.00 && pt <= 100.00) * (0.19) +
2295 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 100.00 && pt <= 120.00) * (0.18) +
2296 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 120.00 && pt <= 140.00) * (0.18) +
2297 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 140.00 && pt <= 160.00) * (0.17) +
2298 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 160.00 && pt <= 180.00) * (0.16) +
2299 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 180.00 && pt <= 200.00) * (0.16) +
2300 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 200.00 && pt <= 250.00) * (0.15) +
2301 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 250.00 && pt <= 300.00) * (0.12) +
2302 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 300.00 && pt <= 350.00) * (0.12) +
2303 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 350.00 && pt <= 400.00) * (0.10) +
2304 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 400.00 && pt <= 500.00) * (0.08) +
2305 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 500.00 && pt <= 600.00) * (0.07) +
2306 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 600.00 && pt <= 700.00) * (0.06) +
2307 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 700.00 && pt <= 800.00) * (0.05) +
2308 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 800.00 && pt <= 1000.00) * (0.04) +
2309 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1000.00 && pt <= 1400.00) * (0.03) +
2310 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 1400.00 && pt <= 2000.00) * (0.02) +
2311 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 2000.00 && pt <= 3000.00) * (0.02) +
2312 (abs(eta) > 0.00 && abs(eta) <= 1.50) * (pt > 3000.00) * (0.00) +
2313 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 20.00 && pt <= 30.00) * (0.08) +
2314 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 30.00 && pt <= 40.00) * (0.10) +
2315 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 40.00 && pt <= 50.00) * (0.10) +
2316 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 50.00 && pt <= 60.00) * (0.12) +
2317 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 60.00 && pt <= 70.00) * (0.12) +
2318 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 70.00 && pt <= 80.00) * (0.13) +
2319 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 80.00 && pt <= 90.00) * (0.13) +
2320 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 90.00 && pt <= 100.00) * (0.14) +
2321 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 100.00 && pt <= 120.00) * (0.14) +
2322 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 120.00 && pt <= 140.00) * (0.14) +
2323 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 140.00 && pt <= 160.00) * (0.15) +
2324 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 160.00 && pt <= 180.00) * (0.15) +
2325 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 180.00 && pt <= 200.00) * (0.15) +
2326 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 200.00 && pt <= 250.00) * (0.15) +
2327 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 250.00 && pt <= 300.00) * (0.13) +
2328 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 300.00 && pt <= 350.00) * (0.12) +
2329 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 350.00 && pt <= 400.00) * (0.11) +
2330 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 400.00 && pt <= 500.00) * (0.09) +
2331 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 500.00 && pt <= 600.00) * (0.07) +
2332 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 600.00 && pt <= 700.00) * (0.07) +
2333 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 700.00 && pt <= 800.00) * (0.06) +
2334 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 800.00 && pt <= 1000.00) * (0.03) +
2335 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1000.00 && pt <= 1400.00) * (0.03) +
2336 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 1400.00 && pt <= 2000.00) * (0.00) +
2337 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 2000.00 && pt <= 3000.00) * (0.00) +
2338 (abs(eta) > 1.50 && abs(eta) <= 2.50) * (pt > 3000.00) * (0.00) +
2339 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 20.00 && pt <= 30.00) * (0.06) +
2340 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 30.00 && pt <= 40.00) * (0.07) +
2341 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 40.00 && pt <= 50.00) * (0.08) +
2342 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 50.00 && pt <= 60.00) * (0.08) +
2343 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 60.00 && pt <= 70.00) * (0.09) +
2344 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 70.00 && pt <= 80.00) * (0.09) +
2345 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 80.00 && pt <= 90.00) * (0.10) +
2346 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 90.00 && pt <= 100.00) * (0.09) +
2347 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 100.00 && pt <= 120.00) * (0.10) +
2348 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 120.00 && pt <= 140.00) * (0.10) +
2349 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 140.00 && pt <= 160.00) * (0.10) +
2350 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 160.00 && pt <= 180.00) * (0.09) +
2351 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 180.00 && pt <= 200.00) * (0.09) +
2352 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 200.00 && pt <= 250.00) * (0.10) +
2353 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 250.00 && pt <= 300.00) * (0.08) +
2354 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 300.00 && pt <= 350.00) * (0.07) +
2355 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 350.00 && pt <= 400.00) * (0.07) +
2356 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 400.00 && pt <= 500.00) * (0.07) +
2357 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 500.00 && pt <= 600.00) * (0.00) +
2358 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 600.00 && pt <= 700.00) * (0.00) +
2359 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 700.00 && pt <= 800.00) * (0.00) +
2360 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 800.00 && pt <= 1000.00) * (0.00) +
2361 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1000.00 && pt <= 1400.00) * (0.00) +
2362 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 1400.00 && pt <= 2000.00) * (0.00) +
2363 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 2000.00 && pt <= 3000.00) * (0.00) +
2364 (abs(eta) > 2.50 && abs(eta) <= 3.50) * (pt > 3000.00) * (0.00)
2365 }
2366
2367
2368}
2369
2370#############
2371# tau-tagging
2372#############
2373
2374
2375module TauTagging TauTagging {
2376 set ParticleInputArray Delphes/allParticles
2377 set PartonInputArray Delphes/partons
2378 set JetInputArray JetEnergyScale/jets
2379
2380 set DeltaR 0.5
2381
2382 set TauPTMin 20.0
2383
2384 set TauEtaMax 2.3
2385
2386 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
2387
2388 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)) + \
2389 (abs(eta) > 2.3) * (0.000)
2390 }
2391 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) ) + \
2392 (abs(eta) > 2.3) * (0.000)
2393 }
2394}
2395
2396module TauTagging TauTaggingAK8 {
2397 set ParticleInputArray Delphes/allParticles
2398 set PartonInputArray Delphes/partons
2399 set JetInputArray JetEnergyScaleAK8/jets
2400
2401 set DeltaR 0.5
2402
2403 set TauPTMin 20.0
2404
2405 set TauEtaMax 2.3
2406
2407 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
2408
2409 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)) + \
2410 (abs(eta) > 2.3) * (0.000)
2411 }
2412 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) ) + \
2413 (abs(eta) > 2.3) * (0.000)
2414 }
2415}
2416
2417module TauTagging TauTaggingPUPPI {
2418 set ParticleInputArray Delphes/allParticles
2419 set PartonInputArray Delphes/partons
2420 set JetInputArray JetEnergyScalePUPPI/jets
2421
2422 set DeltaR 0.5
2423
2424 set TauPTMin 20.0
2425
2426 set TauEtaMax 2.3
2427
2428 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
2429
2430 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)) + \
2431 (abs(eta) > 2.3) * (0.000)
2432 }
2433 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) ) + \
2434 (abs(eta) > 2.3) * (0.000)
2435 }
2436}
2437
2438module TauTagging TauTaggingPUPPIAK8 {
2439 set ParticleInputArray Delphes/allParticles
2440 set PartonInputArray Delphes/partons
2441 set JetInputArray JetEnergyScalePUPPIAK8/jets
2442
2443 set DeltaR 0.5
2444
2445 set TauPTMin 20.0
2446
2447 set TauEtaMax 2.3
2448
2449 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
2450
2451 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)) + \
2452 (abs(eta) > 2.3) * (0.000)
2453 }
2454 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) ) + \
2455 (abs(eta) > 2.3) * (0.000)
2456 }
2457}
2458
2459
2460###############################################################################################################
2461# StatusPidFilter: this module removes all generated particles except electrons, muons, taus, and status == 3 #
2462###############################################################################################################
2463
2464module StatusPidFilter GenParticleFilter {
2465
2466 set InputArray Delphes/allParticles
2467 set OutputArray filteredParticles
2468 set PTMin 5.0
2469
2470}
2471
2472
2473##################
2474# ROOT tree writer
2475##################
2476
2477module TreeWriter TreeWriter {
2478
2479
2480
2481# add Branch InputArray BranchName BranchClass
2482 add Branch GenParticleFilter/filteredParticles Particle GenParticle
2483 #add Branch Delphes/allParticles Particle GenParticle
2484 add Branch PileUpMerger/vertices Vertex Vertex
2485
2486 add Branch GenJetFinder/jets GenJet Jet
2487 add Branch GenJetFinderAK8/jetsAK8 GenJetAK8 Jet
2488 add Branch GenMissingET/momentum GenMissingET MissingET
2489
2490# add Branch HCal/eflowTracks EFlowTrack Track
2491# add Branch ECal/eflowPhotons EFlowPhoton Tower
2492# add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
2493
2494 add Branch PhotonEfficiency/photons Photon Photon
2495 add Branch ElectronEfficiency/electrons Electron Electron
2496 add Branch MuonLooseIdEfficiency/muons MuonLoose Muon
2497 add Branch MuonTightIdEfficiency/muons MuonTight Muon
2498
2499 add Branch PhotonEfficiencyCHS/photons PhotonCHS Photon
2500 add Branch ElectronEfficiencyCHS/electrons ElectronCHS Electron
2501 add Branch MuonLooseIdEfficiencyCHS/muons MuonLooseCHS Muon
2502 add Branch MuonTightIdEfficiencyCHS/muons MuonTightCHS Muon
2503
2504 add Branch JetEnergyScale/jets Jet Jet
2505 add Branch JetEnergyScalePUPPI/jets JetPUPPI Jet
2506 add Branch JetEnergyScaleAK8/jets JetAK8 Jet
2507 add Branch JetEnergyScalePUPPIAK8/jets JetPUPPIAK8 Jet
2508
2509 add Branch Rho/rho Rho Rho
2510
2511 add Branch MissingET/momentum MissingET MissingET
2512 add Branch PuppiMissingET/momentum PuppiMissingET MissingET
2513 add Branch GenPileUpMissingET/momentum GenPileUpMissingET MissingET
2514 add Branch ScalarHT/energy ScalarHT ScalarHT
2515
2516}
Note: See TracBrowser for help on using the repository browser.