Fork me on GitHub

source: git/cards/CMS_PhaseII/CMS_PhaseII_0PU_v02.tcl@ 9e77ec8

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

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