Fork me on GitHub

source: git/cards/CMS_PhaseII_140PU_conf4.tcl@ d67b86b

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

Update and rename delphes_card_CMS_Phase_II_140PileUp_upgrade_Tracker4p0.tcl to CMS_PhaseII_140PU_conf4.tcl
  • Property mode set to 100644
File size: 108.5 KB
RevLine 
[9a16c72]1set MaxEvents 100
[4f558c7]2
3#######################################
4# Order of execution of various modules
5#######################################
6
7set ExecutionPath {
8
9 PileUpMerger
10
11 ParticlePropagator
12
13 ChargedHadronTrackingEfficiency
14 ElectronTrackingEfficiency
15 MuonTrackingEfficiency
16
17 ChargedHadronMomentumSmearing
18 ElectronEnergySmearing
19 MuonMomentumSmearing
20
21 TrackMerger
22
23 Calorimeter
24
25 TrackPileUpSubtractor
26
27 NeutralTowerMerger
28 EFlowMergerAllTracks
29 EFlowMerger
30 RunPUPPI
31
32 NeutrinoFilter
33
34 RhoGridFastJet
35 RhoPuppiGridFastJet
36
37 GenJetFinder
38 RecoJetFinder
39 PuppiJetFinder
40
41 RecoJetPileUpSubtractor
42 PuppiJetPileUpSubtractor
43
44 RecoJetFlavorAssociation
45 PuppiJetFlavorAssociation
46
47 RecoJetPileUpID
48 PuppiJetPileUpID
49
50 RecoJetBTaggingLoose
51 RecoJetBTaggingMedium
52 RecoJetBTaggingTight
53
54 PuppiJetBTaggingLoose
55 PuppiJetBTaggingMedium
56 PuppiJetBTaggingTight
57
58 GenMissingET
59 RecoMissingET
60 JetPileUpIDMissingET
61 PuppiMissingET
62
63 GenScalarHT
64 RecoScalarHT
65 JetPileUpIDScalarHT
66 PuppiScalarHT
67
68 PhotonEfficiency
69 PhotonIsolation
70
71 ElectronEfficiency
72 ElectronIsolation
73
74 MuonEfficiency
75 MuonIsolation
76
77 TauTagging
78
79 StatusPidFilter
80
81 TreeWriter
82
83}
84
85
86
87###############
88# PileUp Merger
89###############
90
91module PileUpMerger PileUpMerger {
92 set InputArray Delphes/stableParticles
93
94 set ParticleOutputArray stableParticles
95 set VertexOutputArray vertices
96
97 # pre-generated minbias input file
98 set PileUpFile MinBias.pileup
99
100 # Get rid of beam spot from http://red-gridftp11.unl.edu/Snowmass/MinBias100K_14TeV.pileup ...
101 set InputBSX 2.44
102 set InputBSY 3.39
103 # replace it with beam spot from CMSSW files
104 set OutputBSX 0.24
105 set OutputBSY 0.39
106
107 # average expected pile up
108 set MeanPileUp 140
109
110 # maximum spread in the beam direction in m
111 set ZVertexSpread 0.15
112
113 # maximum spread in time in s
114 set TVertexSpread 1.5E-09
115
116 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
117
118 #set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
119
120
121 set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) +
122 (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) +
123 (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) +
124 (abs(t) > 1.0e-09) * (abs(z) > 0.15) * (0.00)}
125
126
127}
128
129#####################################################################################
130# Propagate particles in cylinder and divide charged particles in different classes #
131#####################################################################################
132
133module ParticlePropagator ParticlePropagator {
134 ## take particles after beam spot smearing
135 set InputArray PileUpMerger/stableParticles
136 ## produce independent output collection: all particles, only charged hadrons, electrons and muons
137 set OutputArray stableParticles
138 set ChargedHadronOutputArray chargedHadrons
139 set ElectronOutputArray electrons
140 set MuonOutputArray muons
141 ## radius of the magnetic field coverage, in m
142 set Radius 1.29
143 ## half-length of the magnetic field coverage, in m
144 set HalfLength 3.00
145 ## magnetic field
146 set Bz 3.8
147}
148
149####################################
150# Charged hadron tracking efficiency
151####################################
152
153module Efficiency ChargedHadronTrackingEfficiency {
154 ## particles after propagation
155 set InputArray ParticlePropagator/chargedHadrons
156 set OutputArray chargedHadrons
157 # tracking efficiency formula for charged hadrons
158 set EfficiencyFormula {
159 (pt <= 0.2) * (0.00) + \
160 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
161 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
162 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
163 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.87) + \
164 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
165 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.82) + \
166 (abs(eta) > 4.0) * (0.00)
167 }
168}
169
170
171#####################################
172# Electron tracking efficiency - ID
173####################################
174
175module Efficiency ElectronTrackingEfficiency {
176 set InputArray ParticlePropagator/electrons
177 set OutputArray electrons
178 # tracking efficiency formula for electrons
179 set EfficiencyFormula {
180 (pt <= 0.2) * (0.00) + \
181 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
182 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
183 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
184 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 10.0) * (0.82+pt*0.01) + \
185 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 10.0) * (0.90) + \
186 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
187 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.8+pt*0.01) + \
188 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.85) + \
189 (abs(eta) > 4.0) * (0.00)
190
191 }
192}
193
194##########################
195# Muon tracking efficiency
196##########################
197
198module Efficiency MuonTrackingEfficiency {
199 set InputArray ParticlePropagator/muons
200 set OutputArray muons
201 # tracking efficiency formula for muons
202 set EfficiencyFormula {
203 (pt <= 0.2) * (0.00) + \
204 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.998) + \
205 (abs(eta) <= 1.2) * (pt > 1.0) * (0.998) + \
206 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.99) + \
207 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.99) + \
208 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
209 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
210 (abs(eta) > 4.0) * (0.00)
211
212 }
213}
214
215########################################
216# Momentum resolution for charged tracks
217########################################
218
219module MomentumSmearing ChargedHadronMomentumSmearing {
220 ## hadrons after having applied the tracking efficiency
221 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
222 set OutputArray chargedHadrons
223 # resolution formula for charged hadrons
224 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
225 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.013) + \
226 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.02) + \
227 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
228 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
229 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.015) + \
230 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
231 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05) + \
232 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 0.1 && pt <= 1.0) * (0.017) + \
233 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 1.0 && pt <= 10.0) * (0.03) + \
234 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 100.0) * (0.05) + \
235 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 100.0) * (0.30) + \
236 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
237 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 1.0 && pt <= 10.0) * (0.04) + \
238 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 10.0 && pt <= 100.0) * (0.07) + \
239 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 100.0) * (0.30) + \
240 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.025) + \
241 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.05) + \
242 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 100.0) * (0.20) + \
243 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 100.0) * (0.80)
244 }
245}
246
247#################################
248# Energy resolution for electrons
249#################################
250
251module EnergySmearing ElectronEnergySmearing {
252 set InputArray ElectronTrackingEfficiency/electrons
253 set OutputArray electrons
254 # set ResolutionFormula {resolution formula as a function of eta and energy}
255 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.02) + \
256 (abs(eta) <= 2.5) * (energy > 2.5e1) * (energy*0.028) + \
257 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (energy*0.032) + \
258 (abs(eta) > 3.0 && abs(eta) <= 5.0) * (energy*0.07)
259 }
260}
261
262###############################
263# Momentum resolution for muons
264###############################
265
266module MomentumSmearing MuonMomentumSmearing {
267 set InputArray MuonTrackingEfficiency/muons
268 set OutputArray muons
269 # resolution formula for muons
270 set ResolutionFormula {
271 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.012) + \
272 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
273 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.012) + \
274 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.030) + \
275 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.011) + \
276 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.011) + \
277 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.022) + \
278 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.030) + \
279 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 0.1 && pt <= 1.0) * (0.017) + \
280 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 1.0 && pt <= 10.0) * (0.03) + \
281 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 100.0) * (0.05) + \
282 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 100.0) * (0.30) + \
283 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
284 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 1.0 && pt <= 10.0) * (0.04) + \
285 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 10.0 && pt <= 100.0) * (0.07) + \
286 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 100.0) * (0.30) + \
287 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.025) + \
288 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.05) + \
289 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 100.0) * (0.20) + \
290 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 100.0) * (0.80)
291
292 }
293}
294
295
296##############
297# Track merger
298##############
299
300module Merger TrackMerger {
301# add InputArray InputArray
302 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
303 add InputArray ElectronEnergySmearing/electrons
304 add InputArray MuonMomentumSmearing/muons
305 set OutputArray tracks
306}
307
308
309##############################################################
310# Calorimeter : emulate calorimiter answer making caloTowers #
311##############################################################
312
313module Calorimeter Calorimeter {
314 ## particle from the propagation without any efficiency or smearing (for neutrals)
315 set ParticleInputArray ParticlePropagator/stableParticles
316 ## track after smearing and efficiency: used for charged particles
317 set TrackInputArray TrackMerger/tracks
318 ## output collections
319
320 set TowerOutputArray towers
321 set PhotonOutputArray photons
322
323 set EFlowTrackOutputArray eflowTracks
324 set EFlowPhotonOutputArray eflowPhotons
325 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
326
327 set ECalEnergyMin 0.5
328 set HCalEnergyMin 0.5
329
330 set ECalEnergySignificanceMin 2.0
331 set HCalEnergySignificanceMin 2.0
332
333 set SmearTowerCenter false
334
335 set pi [expr {acos(-1)}]
336
337 # lists of the edges of each tower in eta and phi
338 # each list starts with the lower edge of the first tower
339 # the list ends with the higher edged of the last tower
340
341 # 5 degrees towers
342 set PhiBins {}
343 for {set i -70} {$i <= 70} {incr i} {
344 add PhiBins [expr {$i * $pi/70.0}]
345 }
346
347 foreach eta {-4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 -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.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4} {
348 add EtaPhiBins $eta $PhiBins
349 }
350
351 # 10 degrees towers
352 set PhiBins {}
353 for {set i -11} {$i <= 11} {incr i} {
354 add PhiBins [expr {$i * $pi/11.0}]
355 }
356 foreach eta {-4.35 -4.175 -4 4.175 4.35 4.525} {
357 add EtaPhiBins $eta $PhiBins
358 }
359
360 # 20 degrees towers
361 set PhiBins {}
362 for {set i -9} {$i <= 9} {incr i} {
363 add PhiBins [expr {$i * $pi/9.0}]
364 }
365 foreach eta {-5 -4.7 -4.525 4.7 5} {
366 add EtaPhiBins $eta $PhiBins
367 }
368
369 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
370 add EnergyFraction {0} {0.0 1.0}
371 # energy fractions for e, gamma and pi0
372 add EnergyFraction {11} {1.0 0.0}
373 add EnergyFraction {22} {1.0 0.0}
374 add EnergyFraction {111} {1.0 0.0}
375 # energy fractions for muon, neutrinos and neutralinos
376 add EnergyFraction {12} {0.0 0.0}
377 add EnergyFraction {13} {0.0 0.0}
378 add EnergyFraction {14} {0.0 0.0}
379 add EnergyFraction {16} {0.0 0.0}
380 add EnergyFraction {1000022} {0.0 0.0}
381 add EnergyFraction {1000023} {0.0 0.0}
382 add EnergyFraction {1000025} {0.0 0.0}
383 add EnergyFraction {1000035} {0.0 0.0}
384 add EnergyFraction {1000045} {0.0 0.0}
385 # energy fractions for K0short and Lambda
386 add EnergyFraction {310} {0.3 0.7}
387 add EnergyFraction {3122} {0.3 0.7}
388
389
390 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
391 set ECalResolutionFormula {(abs(eta) <= 1.5) * sqrt(energy^2*0.005^2 + energy*0.027^2 + 0.15^2) + \
392 (abs(eta) > 1.5 && abs(eta) <= 4.0) * sqrt(energy^2*0.008^2 + energy*0.092^2 + 0.088^2) + \
393 (abs(eta) > 4.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.97^2)}
394
395
396
397 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
398 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
399 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.050^2 + energy*0.706^2) + \
400 (abs(eta) > 3.0 && abs(eta) <= 4.9) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
401
402
403
404}
405
406##########################
407# Track pile-up subtractor
408##########################
409
410module TrackPileUpSubtractor TrackPileUpSubtractor {
411# add InputArray InputArray OutputArray
412 add InputArray Calorimeter/eflowTracks eflowTracks
413 add InputArray ElectronEnergySmearing/electrons electrons
414 add InputArray MuonMomentumSmearing/muons muons
415
416 set VertexInputArray PileUpMerger/vertices
417 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
418 # Z vertex resolution in m
419 set ZVertexResolution 0.0001
420}
421
422
423####################
424# Neutral tower merger
425####################
426
427module Merger NeutralTowerMerger {
428# add InputArray InputArray
429 add InputArray Calorimeter/eflowPhotons
430 add InputArray Calorimeter/eflowNeutralHadrons
431 set OutputArray eflowTowers
432}
433
434
435##################################
436# Energy flow merger (all tracks)
437##################################
438
439module Merger EFlowMergerAllTracks {
440# add InputArray InputArray
441 add InputArray TrackMerger/tracks
442 add InputArray Calorimeter/eflowPhotons
443 add InputArray Calorimeter/eflowNeutralHadrons
444 set OutputArray eflow
445}
446
447
448####################
449# Energy flow merger
450####################
451
452module Merger EFlowMerger {
453# add InputArray InputArray
454 add InputArray TrackPileUpSubtractor/eflowTracks
455 add InputArray Calorimeter/eflowPhotons
456 add InputArray Calorimeter/eflowNeutralHadrons
457 set OutputArray eflow
458}
459
460###########################
461### Run the puppi code ###
462###########################
463
464module RunPUPPI RunPUPPI {
465 ## input information
466 set TrackInputArray TrackMerger/tracks
467 set NeutralInputArray NeutralTowerMerger/eflowTowers
468 set PVInputArray PileUpMerger/vertices
469 set MinPuppiWeight 0.05
470 set UseExp false
471
472 ## define puppi algorithm parameters (more than one for the same eta region is possible)
473 add EtaMinBin 0. 2.5 2.5 3.0 3.0
474 add EtaMaxBin 2.5 3.0 3.0 10.0 10.0
475 add PtMinBin 0. 0.5 0.5 0.5 0.5
476 add ConeSizeBin 0.25 0.25 0.25 0.25 0.25
477 add RMSPtMinBin 0.1 0.5 0.5 0.5 0.5
478 add RMSScaleFactorBin 1.0 1.0 1.0 1.0 1.0
479 add NeutralMinEBin 0.2 1.0 1.0 1.5 1.5
480 add NeutralPtSlope 0.02 0.02 0.02 0.02 0.02
481 add ApplyCHS true true true true true
482 add UseCharged true false false false false
483 add ApplyLowPUCorr true true true true true
484 add MetricId 5 5 1 5 1
485
486 ## output name
487 set OutputArray PuppiParticles
488 set OutputArrayTracks puppiTracks
489 set OutputArrayNeutrals puppiNeutrals
490}
491
492
493
494
495#####################
496# Neutrino Filter
497#####################
498
499module PdgCodeFilter NeutrinoFilter {
500
501 set InputArray Delphes/stableParticles
502 set OutputArray filteredParticles
503
504 set PTMin 0.0
505
506 add PdgCode {12}
507 add PdgCode {14}
508 add PdgCode {16}
509 add PdgCode {-12}
510 add PdgCode {-14}
511 add PdgCode {-16}
512
513}
514
515
516#################################
517# Rho pile-up with std particles
518#################################
519
520module FastJetGridMedianEstimator RhoGridFastJet {
521
522 set InputArray EFlowMerger/eflow
523 set RhoOutputArray rho
524
525 # add GridRange rapmin rapmax drap dphi
526 # rapmin - the minimum rapidity extent of the grid
527 # rapmax - the maximum rapidity extent of the grid
528 # drap - the grid spacing in rapidity
529 # dphi - the grid spacing in azimuth
530
531 add GridRange -5.0 -4.0 1.0 1.0
532 add GridRange -4.0 -2.5 1.0 1.0
533 add GridRange -2.5 2.5 1.0 1.0
534 add GridRange 2.5 4.0 1.0 1.0
535 add GridRange 4.0 5.0 1.0 1.0
536
537}
538
539###################################
540# Rho pile-up with PUPPI particles
541###################################
542
543module FastJetGridMedianEstimator RhoPuppiGridFastJet {
544
545 set InputArray RunPUPPI/PuppiParticles
546 set RhoOutputArray rho
547
548 # add GridRange rapmin rapmax drap dphi
549 # rapmin - the minimum rapidity extent of the grid
550 # rapmax - the maximum rapidity extent of the grid
551 # drap - the grid spacing in rapidity
552 # dphi - the grid spacing in azimuth
553
554 add GridRange -5.0 -4.0 1.0 1.0
555 add GridRange -4.0 -2.5 1.0 1.0
556 add GridRange -2.5 2.5 1.0 1.0
557 add GridRange 2.5 4.0 1.0 1.0
558 add GridRange 4.0 5.0 1.0 1.0
559
560}
561
562
563#####################
564# MC truth jet finder
565#####################
566
567module FastJetFinder GenJetFinder {
568 set InputArray NeutrinoFilter/filteredParticles
569
570 set OutputArray jets
571
572 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
573 set JetAlgorithm 6
574 set ParameterR 0.4
575
576 set JetPTMin 15.0
577}
578
579
580################
581## Jet finder ##
582################
583
584module FastJetFinder RecoJetFinder {
585 set InputArray EFlowMerger/eflow
586 set OutputArray jets
587 # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area
588 set AreaAlgorithm 1
589 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
590 set JetAlgorithm 6
591 set ParameterR 0.4
592 set JetPTMin 10.0
593}
594
595
596
597#####################
598## Make puppi jets ##
599#####################
600
601module FastJetFinder PuppiJetFinder {
602 set InputArray RunPUPPI/PuppiParticles
603 set OutputArray jets
604 # area algorithm: 0 Do not compute area, 1 Active area explicit ghosts, 2 One ghost passive area, 3 Passive area, 4 Voronoi, 5 Active area
605 set AreaAlgorithm 1
606 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
607 set JetAlgorithm 6
608 set ParameterR 0.4
609 set JetPTMin 10.0
610}
611
612
613
614###########################
615# Jet Pile-Up Subtraction
616###########################
617
618module JetPileUpSubtractor RecoJetPileUpSubtractor {
619 set JetInputArray RecoJetFinder/jets
620 set RhoInputArray RhoGridFastJet/rho
621
622 set OutputArray jets
623
624 set JetPTMin 20.0
625}
626
627
628###########################
629# Jet Pile-Up Subtraction
630###########################
631
632module JetPileUpSubtractor PuppiJetPileUpSubtractor {
633 set JetInputArray PuppiJetFinder/jets
634 set RhoInputArray RhoPuppiGridFastJet/rho
635
636 set OutputArray jets
637
638 set JetPTMin 20.0
639}
640
641
642
643########################
644# Jet Flavor Association
645########################
646
647module JetFlavorAssociation RecoJetFlavorAssociation {
648
649 set PartonInputArray Delphes/partons
650 set ParticleInputArray Delphes/allParticles
651 set ParticleLHEFInputArray Delphes/allParticlesLHEF
652 set JetInputArray RecoJetPileUpSubtractor/jets
653
654 set DeltaR 0.5
655 set PartonPTMin 1.0
656 set PartonEtaMax 2.5
657
658}
659
660
661
662########################
663# Jet Flavor Association
664########################
665
666module JetFlavorAssociation PuppiJetFlavorAssociation {
667
668 set PartonInputArray Delphes/partons
669 set ParticleInputArray Delphes/allParticles
670 set ParticleLHEFInputArray Delphes/allParticlesLHEF
671 set JetInputArray PuppiJetPileUpSubtractor/jets
672
673 set DeltaR 0.5
674 set PartonPTMin 1.0
675 set PartonEtaMax 2.5
676
677}
678
679
680###########################
681# Jet Pile-Up ID
682###########################
683
684module PileUpJetID RecoJetPileUpID {
685
686 set JetInputArray RecoJetPileUpSubtractor/jets
687
688 # Using constituents does not make sense with Charged hadron subtraction
689 # In 0 mode, dR cut used instead
690 set UseConstituents 0
691
692 set TrackInputArray TrackMerger/tracks
693 set NeutralInputArray NeutralTowerMerger/eflowTowers
694
695 set OutputArray jets
696 set NeutralsInPassingJets eflowTowers
697
698 set ParameterR 0.4
699
700 set MeanSqDeltaRMaxBarrel 0.07
701 set BetaMinBarrel 0.13
702 set MeanSqDeltaRMaxEndcap 0.07
703 set BetaMinEndcap 0.15
704 set MeanSqDeltaRMaxForward 0.01
705
706 set JetPTMin 10.0
707
708}
709
710###########################
711# Jet Pile-Up ID
712###########################
713
714module PileUpJetID PuppiJetPileUpID {
715
716 set JetInputArray PuppiJetPileUpSubtractor/jets
717
718 # Using constituents does not make sense with Charged hadron subtraction
719 # In 0 mode, dR cut used instead
720 set UseConstituents 0
721
722 set TrackInputArray RunPUPPI/puppiTracks
723 set NeutralInputArray RunPUPPI/puppiNeutrals
724
725 set OutputArray jets
726 set NeutralsInPassingJets eflowTowers
727
728 set ParameterR 0.4
729
730 set MeanSqDeltaRMaxBarrel 0.07
731 set BetaMinBarrel 0.13
732 set MeanSqDeltaRMaxEndcap 0.07
733 set BetaMinEndcap 0.15
734 set MeanSqDeltaRMaxForward 0.01
735
736 set JetPTMin 10.0
737
738}
739#############
740# b-tagging #
741#############
742module BTagging RecoJetBTaggingLoose {
743
744 set JetInputArray RecoJetPileUpSubtractor/jets
745
746 set BitNumber 0
747
748 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
749 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.096) + \
750 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.112) + \
751 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.0812) + \
752 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.0908) + \
753 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.0792) + \
754 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.0858) + \
755 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.0917) + \
756 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.084) + \
757 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.0906) + \
758 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.0989) + \
759 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.1022) + \
760 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.1035) + \
761 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.1095) + \
762 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.1201) + \
763 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.1348) + \
764 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.1482) + \
765 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.1629) + \
766 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.1775) + \
767 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.2002) + \
768 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.1977) + \
769 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.2084) + \
770 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.2195) + \
771 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.2424) + \
772 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.2909) + \
773 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.3457) + \
774 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
775 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.074765) + \
776 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.100053) + \
777 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.071492) + \
778 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.084796) + \
779 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.076927) + \
780 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.08424) + \
781 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.093118) + \
782 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.084629) + \
783 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.092977) + \
784 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.10206) + \
785 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.102344) + \
786 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.098435) + \
787 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.105507) + \
788 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.112841) + \
789 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.126329) + \
790 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.140759) + \
791 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.153193) + \
792 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.107869) + \
793 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.119527) + \
794 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.08688) + \
795 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.089324) + \
796 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.097172) + \
797 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.000) + \
798 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.000) + \
799 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.000) + \
800 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (0.02) +\
801 (abs(eta) > 4.0) * (0.000) \
802 }
803
804 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
805 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.387) + \
806 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.448) + \
807 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.408) + \
808 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.427) + \
809 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.408) + \
810 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.425) + \
811 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.426) + \
812 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.4) + \
813 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.415) + \
814 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.416) + \
815 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.405) + \
816 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.387) + \
817 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.39) + \
818 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.389) + \
819 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.389) + \
820 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.381) + \
821 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.381) + \
822 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.367) + \
823 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.369) + \
824 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.326) + \
825 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.335) + \
826 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.326) + \
827 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.341) + \
828 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.403) + \
829 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.47) + \
830 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
831 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.2497) + \
832 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.31891) + \
833 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.273) + \
834 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.28445) + \
835 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.28036) + \
836 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.28453) + \
837 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.29495) + \
838 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.27024) + \
839 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.28912) + \
840 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.29048) + \
841 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.27507) + \
842 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.23327) + \
843 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.2493) + \
844 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.2416) + \
845 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.26652) + \
846 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.24852) + \
847 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.26927) + \
848 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.19302) + \
849 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.19433) + \
850 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.17523) + \
851 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.14981) + \
852 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.16666) + \
853 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
854 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
855 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.) + \
856 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.29*tanh(pt*0.00997 - 0.143)) + \
857 (abs(eta) > 4.0) * (0.000)
858 }
859
860 # efficiency formula for b-jets
861 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
862 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.75) + \
863 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.827) + \
864 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.837) + \
865 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.853) + \
866 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.855) + \
867 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.862) + \
868 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.868) + \
869 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.865) + \
870 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.863) + \
871 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.857) + \
872 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.851) + \
873 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.838) + \
874 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.831) + \
875 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.817) + \
876 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.796) + \
877 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.772) + \
878 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.76) + \
879 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.743) + \
880 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.703) + \
881 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.638) + \
882 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.605) + \
883 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.572) + \
884 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.541) + \
885 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.567) + \
886 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.603) + \
887 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
888 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.6063) + \
889 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.7188) + \
890 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.72) + \
891 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.7365) + \
892 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.7462) + \
893 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.7454) + \
894 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.7415) + \
895 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.727) + \
896 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.7112) + \
897 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.7112) + \
898 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.6754) + \
899 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.6359) + \
900 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.6348) + \
901 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.6115) + \
902 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.5585) + \
903 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.5608) + \
904 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.5208) + \
905 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.456) + \
906 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.4524) + \
907 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.388) + \
908 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.3928) + \
909 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.3823) + \
910 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
911 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
912 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \
913 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \
914 (abs(eta) > 4.0) * (0.000)
915 }
916}
917
918
919#############
920# b-tagging #
921#############
922module BTagging RecoJetBTaggingMedium {
923
924 set JetInputArray RecoJetPileUpSubtractor/jets
925
926 set BitNumber 1
927
928 add EfficiencyFormula {0} { (pt <= 20.0) * (0.000) + \
929 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.00654) + \
930 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.00921) + \
931 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.00573) + \
932 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.00694) + \
933 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.0062) + \
934 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.00708) + \
935 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.00779) + \
936 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.00693) + \
937 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.00777) + \
938 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.00862) + \
939 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.01038) + \
940 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.01189) + \
941 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.01279) + \
942 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.01452) + \
943 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.01696) + \
944 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.01958) + \
945 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.02253) + \
946 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.01787) + \
947 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.02154) + \
948 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.01839) + \
949 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.01987) + \
950 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.02351) + \
951 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.02937) + \
952 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.04001) + \
953 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.0542) + \
954 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
955 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.004236) + \
956 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.006653) + \
957 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.005512) + \
958 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.00754) + \
959 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.005813) + \
960 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.006439) + \
961 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.008063) + \
962 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.00647) + \
963 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.007583) + \
964 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.008543) + \
965 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.01034) + \
966 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.011253) + \
967 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.012945) + \
968 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.014474) + \
969 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.017361) + \
970 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.020912) + \
971 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.023139) + \
972 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.010756) + \
973 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.012569) + \
974 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.006046) + \
975 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.006428) + \
976 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.00887) + \
977 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
978 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
979 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.000) + \
980 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (0.001) + \
981 (abs(eta) > 4.0) * (0.000) \
982}
983
984 # efficiency formula for c-jets (misidentification rate)
985 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
986 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.1102) + \
987 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.1344) + \
988 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.1025) + \
989 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.1025) + \
990 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.1063) + \
991 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.1087) + \
992 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.1124) + \
993 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.109) + \
994 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.111) + \
995 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.1091) + \
996 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.1087) + \
997 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.1091) + \
998 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.1107) + \
999 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.1061) + \
1000 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.1017) + \
1001 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.1) + \
1002 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.0966) + \
1003 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.0697) + \
1004 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.0679) + \
1005 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.0503) + \
1006 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.0514) + \
1007 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.0481) + \
1008 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.0667) + \
1009 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.0861) + \
1010 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.092) + \
1011 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1012 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.03331) + \
1013 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.04361) + \
1014 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.03863) + \
1015 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.04287) + \
1016 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.04431) + \
1017 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.04452) + \
1018 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.04339) + \
1019 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.0436) + \
1020 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.0456) + \
1021 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.05138) + \
1022 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.04794) + \
1023 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.04004) + \
1024 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.04713) + \
1025 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.04515) + \
1026 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.05314) + \
1027 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.05143) + \
1028 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.05936) + \
1029 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.02357) + \
1030 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.03222) + \
1031 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.01523) + \
1032 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.02621) + \
1033 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.01709) + \
1034 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1035 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1036 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.) + \
1037 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.1898*tanh(pt*0.00997 - 0.143)) + \
1038 (abs(eta) > 4.0) * (0.000)
1039 }
1040
1041 # efficiency formula for b-jets
1042 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1043 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.536) + \
1044 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.6439) + \
1045 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.6504) + \
1046 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.6716) + \
1047 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.6841) + \
1048 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.6896) + \
1049 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.6916) + \
1050 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.6882) + \
1051 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.6838) + \
1052 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.6715) + \
1053 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.6554) + \
1054 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.6366) + \
1055 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.6192) + \
1056 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.595) + \
1057 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.5551) + \
1058 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.5138) + \
1059 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.4884) + \
1060 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.4009) + \
1061 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.3459) + \
1062 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.2523) + \
1063 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.2404) + \
1064 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.2198) + \
1065 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.2263) + \
1066 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.2614) + \
1067 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.3194) + \
1068 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1069 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.3254) + \
1070 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.4339) + \
1071 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.4499) + \
1072 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.4716) + \
1073 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.4766) + \
1074 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.4788) + \
1075 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.4863) + \
1076 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.4891) + \
1077 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.462) + \
1078 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.4583) + \
1079 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.4247) + \
1080 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.3775) + \
1081 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.3734) + \
1082 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.3348) + \
1083 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.2939) + \
1084 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.285) + \
1085 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.2421) + \
1086 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.1565) + \
1087 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.1522) + \
1088 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.1231) + \
1089 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.1607) + \
1090 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.1323) + \
1091 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1092 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1093 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \
1094 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \
1095 (abs(eta) > 4.0) * (0.000)
1096 }
1097}
1098
1099#############
1100# b-tagging #
1101#############
1102module BTagging RecoJetBTaggingTight {
1103
1104 set JetInputArray RecoJetPileUpSubtractor/jets
1105
1106 set BitNumber 2
1107
1108 add EfficiencyFormula {0} { (pt <= 20.0) * (0.000) + \
1109 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.000164) + \
1110 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.000235) + \
1111 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.000266) + \
1112 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.000329) + \
1113 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.000309) + \
1114 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.000309) + \
1115 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.000546) + \
1116 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.000499) + \
1117 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.000642) + \
1118 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.000742) + \
1119 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.000928) + \
1120 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.001323) + \
1121 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.001392) + \
1122 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.00154) + \
1123 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.002094) + \
1124 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.002427) + \
1125 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.002927) + \
1126 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.001854) + \
1127 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.002355) + \
1128 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.002297) + \
1129 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.002433) + \
1130 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.002706) + \
1131 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.003602) + \
1132 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.004987) + \
1133 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.007414) + \
1134 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1135 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.000573) + \
1136 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.000574) + \
1137 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.000521) + \
1138 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.000786) + \
1139 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.000539) + \
1140 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.000673) + \
1141 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.000934) + \
1142 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.000781) + \
1143 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.000949) + \
1144 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.000977) + \
1145 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.001168) + \
1146 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.000879) + \
1147 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.000812) + \
1148 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.001215) + \
1149 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.001679) + \
1150 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.001893) + \
1151 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.002723) + \
1152 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.003555) + \
1153 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.003881) + \
1154 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.006046) + \
1155 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.005563) + \
1156 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.007611) + \
1157 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1158 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1159 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.000) + \
1160 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (0.001) + \
1161 (abs(eta) > 4.0) * (0.000) \
1162 }
1163
1164 # efficiency formula for c-jets (misidentification rate)
1165 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1166 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.00531) + \
1167 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.00567) + \
1168 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.0064) + \
1169 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.00673) + \
1170 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.00766) + \
1171 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.00729) + \
1172 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.00674) + \
1173 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.00824) + \
1174 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.00888) + \
1175 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.00919) + \
1176 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.01021) + \
1177 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.01041) + \
1178 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.01027) + \
1179 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.00999) + \
1180 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.01047) + \
1181 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.01014) + \
1182 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.01021) + \
1183 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.00601) + \
1184 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.0054) + \
1185 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.00487) + \
1186 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.00519) + \
1187 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.00469) + \
1188 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.00651) + \
1189 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.01299) + \
1190 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.00897) + \
1191 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1192 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.01014) + \
1193 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.01288) + \
1194 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.01392) + \
1195 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.01533) + \
1196 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.01508) + \
1197 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.01579) + \
1198 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.01106) + \
1199 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.01346) + \
1200 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.01315) + \
1201 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.01156) + \
1202 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.0082) + \
1203 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.00439) + \
1204 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.00744) + \
1205 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.00685) + \
1206 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.00755) + \
1207 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.00706) + \
1208 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.00428) + \
1209 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.01031) + \
1210 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.00976) + \
1211 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.01523) + \
1212 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.00374) + \
1213 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.00854) + \
1214 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1215 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1216 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \
1217 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \
1218 (abs(eta) > 4.0) * (0.000)
1219}
1220
1221 # efficiency formula for b-jets
1222 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1223 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.2426) + \
1224 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.327) + \
1225 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.3559) + \
1226 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.3704) + \
1227 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.3824) + \
1228 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.3844) + \
1229 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.3848) + \
1230 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.3862) + \
1231 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.3778) + \
1232 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.3622) + \
1233 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.3299) + \
1234 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.2889) + \
1235 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.2815) + \
1236 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.253) + \
1237 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.221) + \
1238 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.1963) + \
1239 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.1739) + \
1240 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.0992) + \
1241 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.0788) + \
1242 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.0581) + \
1243 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.0534) + \
1244 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.0521) + \
1245 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.0626) + \
1246 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.0826) + \
1247 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.1022) + \
1248 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1249 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.1562) + \
1250 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.2499) + \
1251 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.2956) + \
1252 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.3128) + \
1253 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.3147) + \
1254 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.3222) + \
1255 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.304) + \
1256 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.3051) + \
1257 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.2657) + \
1258 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.2578) + \
1259 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.2087) + \
1260 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.1634) + \
1261 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.1651) + \
1262 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.1353) + \
1263 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.109) + \
1264 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.0799) + \
1265 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.0699) + \
1266 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.054) + \
1267 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.0718) + \
1268 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.0746) + \
1269 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.0803) + \
1270 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.0882) + \
1271 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1272 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1273 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \
1274 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \
1275 (abs(eta) > 4.0) * (0.000)
1276 }
1277
1278}
1279
1280
1281
1282
1283
1284
1285#############
1286# b-tagging #
1287#############
1288module BTagging PuppiJetBTaggingLoose {
1289
1290 set JetInputArray PuppiJetPileUpSubtractor/jets
1291
1292 set BitNumber 0
1293
1294 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
1295 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.096) + \
1296 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.112) + \
1297 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.0812) + \
1298 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.0908) + \
1299 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.0792) + \
1300 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.0858) + \
1301 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.0917) + \
1302 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.084) + \
1303 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.0906) + \
1304 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.0989) + \
1305 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.1022) + \
1306 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.1035) + \
1307 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.1095) + \
1308 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.1201) + \
1309 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.1348) + \
1310 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.1482) + \
1311 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.1629) + \
1312 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.1775) + \
1313 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.2002) + \
1314 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.1977) + \
1315 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.2084) + \
1316 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.2195) + \
1317 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.2424) + \
1318 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.2909) + \
1319 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.3457) + \
1320 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1321 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.074765) + \
1322 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.100053) + \
1323 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.071492) + \
1324 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.084796) + \
1325 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.076927) + \
1326 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.08424) + \
1327 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.093118) + \
1328 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.084629) + \
1329 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.092977) + \
1330 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.10206) + \
1331 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.102344) + \
1332 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.098435) + \
1333 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.105507) + \
1334 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.112841) + \
1335 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.126329) + \
1336 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.140759) + \
1337 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.153193) + \
1338 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.107869) + \
1339 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.119527) + \
1340 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.08688) + \
1341 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.089324) + \
1342 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.097172) + \
1343 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.000) + \
1344 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.000) + \
1345 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.000) + \
1346 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (0.02) +\
1347 (abs(eta) > 4.0) * (0.000) \
1348 }
1349
1350 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1351 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.387) + \
1352 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.448) + \
1353 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.408) + \
1354 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.427) + \
1355 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.408) + \
1356 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.425) + \
1357 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.426) + \
1358 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.4) + \
1359 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.415) + \
1360 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.416) + \
1361 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.405) + \
1362 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.387) + \
1363 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.39) + \
1364 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.389) + \
1365 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.389) + \
1366 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.381) + \
1367 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.381) + \
1368 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.367) + \
1369 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.369) + \
1370 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.326) + \
1371 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.335) + \
1372 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.326) + \
1373 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.341) + \
1374 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.403) + \
1375 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.47) + \
1376 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1377 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.2497) + \
1378 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.31891) + \
1379 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.273) + \
1380 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.28445) + \
1381 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.28036) + \
1382 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.28453) + \
1383 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.29495) + \
1384 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.27024) + \
1385 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.28912) + \
1386 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.29048) + \
1387 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.27507) + \
1388 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.23327) + \
1389 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.2493) + \
1390 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.2416) + \
1391 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.26652) + \
1392 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.24852) + \
1393 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.26927) + \
1394 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.19302) + \
1395 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.19433) + \
1396 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.17523) + \
1397 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.14981) + \
1398 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.16666) + \
1399 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1400 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1401 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.) + \
1402 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.29*tanh(pt*0.00997 - 0.143)) + \
1403 (abs(eta) > 4.0) * (0.000)
1404 }
1405
1406 # efficiency formula for b-jets
1407 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1408 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.75) + \
1409 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.827) + \
1410 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.837) + \
1411 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.853) + \
1412 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.855) + \
1413 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.862) + \
1414 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.868) + \
1415 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.865) + \
1416 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.863) + \
1417 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.857) + \
1418 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.851) + \
1419 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.838) + \
1420 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.831) + \
1421 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.817) + \
1422 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.796) + \
1423 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.772) + \
1424 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.76) + \
1425 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.743) + \
1426 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.703) + \
1427 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.638) + \
1428 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.605) + \
1429 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.572) + \
1430 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.541) + \
1431 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.567) + \
1432 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.603) + \
1433 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1434 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.6063) + \
1435 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.7188) + \
1436 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.72) + \
1437 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.7365) + \
1438 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.7462) + \
1439 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.7454) + \
1440 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.7415) + \
1441 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.727) + \
1442 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.7112) + \
1443 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.7112) + \
1444 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.6754) + \
1445 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.6359) + \
1446 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.6348) + \
1447 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.6115) + \
1448 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.5585) + \
1449 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.5608) + \
1450 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.5208) + \
1451 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.456) + \
1452 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.4524) + \
1453 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.388) + \
1454 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.3928) + \
1455 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.3823) + \
1456 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1457 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1458 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \
1459 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \
1460 (abs(eta) > 4.0) * (0.000)
1461 }
1462}
1463
1464
1465#############
1466# b-tagging #
1467#############
1468module BTagging PuppiJetBTaggingMedium {
1469
1470 set JetInputArray PuppiJetPileUpSubtractor/jets
1471
1472 set BitNumber 1
1473
1474 add EfficiencyFormula {0} { (pt <= 20.0) * (0.000) + \
1475 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.00654) + \
1476 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.00921) + \
1477 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.00573) + \
1478 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.00694) + \
1479 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.0062) + \
1480 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.00708) + \
1481 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.00779) + \
1482 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.00693) + \
1483 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.00777) + \
1484 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.00862) + \
1485 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.01038) + \
1486 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.01189) + \
1487 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.01279) + \
1488 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.01452) + \
1489 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.01696) + \
1490 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.01958) + \
1491 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.02253) + \
1492 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.01787) + \
1493 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.02154) + \
1494 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.01839) + \
1495 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.01987) + \
1496 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.02351) + \
1497 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.02937) + \
1498 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.04001) + \
1499 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.0542) + \
1500 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1501 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.004236) + \
1502 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.006653) + \
1503 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.005512) + \
1504 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.00754) + \
1505 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.005813) + \
1506 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.006439) + \
1507 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.008063) + \
1508 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.00647) + \
1509 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.007583) + \
1510 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.008543) + \
1511 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.01034) + \
1512 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.011253) + \
1513 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.012945) + \
1514 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.014474) + \
1515 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.017361) + \
1516 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.020912) + \
1517 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.023139) + \
1518 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.010756) + \
1519 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.012569) + \
1520 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.006046) + \
1521 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.006428) + \
1522 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.00887) + \
1523 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1524 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1525 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.000) + \
1526 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (0.001) + \
1527 (abs(eta) > 4.0) * (0.000) \
1528}
1529
1530 # efficiency formula for c-jets (misidentification rate)
1531 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1532 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.1102) + \
1533 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.1344) + \
1534 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.1025) + \
1535 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.1025) + \
1536 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.1063) + \
1537 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.1087) + \
1538 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.1124) + \
1539 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.109) + \
1540 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.111) + \
1541 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.1091) + \
1542 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.1087) + \
1543 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.1091) + \
1544 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.1107) + \
1545 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.1061) + \
1546 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.1017) + \
1547 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.1) + \
1548 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.0966) + \
1549 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.0697) + \
1550 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.0679) + \
1551 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.0503) + \
1552 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.0514) + \
1553 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.0481) + \
1554 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.0667) + \
1555 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.0861) + \
1556 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.092) + \
1557 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1558 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.03331) + \
1559 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.04361) + \
1560 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.03863) + \
1561 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.04287) + \
1562 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.04431) + \
1563 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.04452) + \
1564 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.04339) + \
1565 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.0436) + \
1566 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.0456) + \
1567 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.05138) + \
1568 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.04794) + \
1569 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.04004) + \
1570 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.04713) + \
1571 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.04515) + \
1572 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.05314) + \
1573 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.05143) + \
1574 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.05936) + \
1575 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.02357) + \
1576 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.03222) + \
1577 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.01523) + \
1578 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.02621) + \
1579 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.01709) + \
1580 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1581 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1582 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.) + \
1583 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.1898*tanh(pt*0.00997 - 0.143)) + \
1584 (abs(eta) > 4.0) * (0.000)
1585 }
1586
1587 # efficiency formula for b-jets
1588 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1589 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.536) + \
1590 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.6439) + \
1591 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.6504) + \
1592 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.6716) + \
1593 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.6841) + \
1594 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.6896) + \
1595 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.6916) + \
1596 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.6882) + \
1597 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.6838) + \
1598 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.6715) + \
1599 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.6554) + \
1600 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.6366) + \
1601 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.6192) + \
1602 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.595) + \
1603 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.5551) + \
1604 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.5138) + \
1605 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.4884) + \
1606 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.4009) + \
1607 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.3459) + \
1608 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.2523) + \
1609 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.2404) + \
1610 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.2198) + \
1611 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.2263) + \
1612 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.2614) + \
1613 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.3194) + \
1614 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1615 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.3254) + \
1616 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.4339) + \
1617 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.4499) + \
1618 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.4716) + \
1619 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.4766) + \
1620 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.4788) + \
1621 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.4863) + \
1622 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.4891) + \
1623 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.462) + \
1624 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.4583) + \
1625 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.4247) + \
1626 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.3775) + \
1627 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.3734) + \
1628 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.3348) + \
1629 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.2939) + \
1630 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.285) + \
1631 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.2421) + \
1632 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.1565) + \
1633 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.1522) + \
1634 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.1231) + \
1635 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.1607) + \
1636 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.1323) + \
1637 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1638 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1639 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \
1640 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \
1641 (abs(eta) > 4.0) * (0.000)
1642 }
1643}
1644
1645#############
1646# b-tagging #
1647#############
1648module BTagging PuppiJetBTaggingTight {
1649
1650 set JetInputArray PuppiJetPileUpSubtractor/jets
1651
1652 set BitNumber 2
1653
1654 add EfficiencyFormula {0} { (pt <= 20.0) * (0.000) + \
1655 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.000164) + \
1656 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.000235) + \
1657 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.000266) + \
1658 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.000329) + \
1659 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.000309) + \
1660 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.000309) + \
1661 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.000546) + \
1662 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.000499) + \
1663 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.000642) + \
1664 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.000742) + \
1665 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.000928) + \
1666 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.001323) + \
1667 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.001392) + \
1668 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.00154) + \
1669 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.002094) + \
1670 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.002427) + \
1671 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.002927) + \
1672 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.001854) + \
1673 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.002355) + \
1674 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.002297) + \
1675 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.002433) + \
1676 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.002706) + \
1677 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.003602) + \
1678 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.004987) + \
1679 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.007414) + \
1680 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1681 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.000573) + \
1682 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.000574) + \
1683 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.000521) + \
1684 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.000786) + \
1685 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.000539) + \
1686 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.000673) + \
1687 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.000934) + \
1688 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.000781) + \
1689 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.000949) + \
1690 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.000977) + \
1691 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.001168) + \
1692 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.000879) + \
1693 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.000812) + \
1694 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.001215) + \
1695 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.001679) + \
1696 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.001893) + \
1697 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.002723) + \
1698 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.003555) + \
1699 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.003881) + \
1700 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.006046) + \
1701 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.005563) + \
1702 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.007611) + \
1703 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1704 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1705 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.000) + \
1706 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (0.001) + \
1707 (abs(eta) > 4.0) * (0.000) \
1708 }
1709
1710 # efficiency formula for c-jets (misidentification rate)
1711 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1712 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.00531) + \
1713 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.00567) + \
1714 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.0064) + \
1715 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.00673) + \
1716 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.00766) + \
1717 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.00729) + \
1718 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.00674) + \
1719 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.00824) + \
1720 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.00888) + \
1721 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.00919) + \
1722 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.01021) + \
1723 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.01041) + \
1724 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.01027) + \
1725 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.00999) + \
1726 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.01047) + \
1727 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.01014) + \
1728 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.01021) + \
1729 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.00601) + \
1730 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.0054) + \
1731 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.00487) + \
1732 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.00519) + \
1733 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.00469) + \
1734 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.00651) + \
1735 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.01299) + \
1736 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.00897) + \
1737 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1738 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.01014) + \
1739 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.01288) + \
1740 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.01392) + \
1741 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.01533) + \
1742 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.01508) + \
1743 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.01579) + \
1744 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.01106) + \
1745 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.01346) + \
1746 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.01315) + \
1747 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.01156) + \
1748 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.0082) + \
1749 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.00439) + \
1750 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.00744) + \
1751 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.00685) + \
1752 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.00755) + \
1753 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.00706) + \
1754 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.00428) + \
1755 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.01031) + \
1756 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.00976) + \
1757 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.01523) + \
1758 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.00374) + \
1759 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.00854) + \
1760 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1761 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1762 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \
1763 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \
1764 (abs(eta) > 4.0) * (0.000)
1765}
1766
1767 # efficiency formula for b-jets
1768 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1769 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.2426) + \
1770 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.327) + \
1771 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.3559) + \
1772 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.3704) + \
1773 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.3824) + \
1774 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.3844) + \
1775 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.3848) + \
1776 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.3862) + \
1777 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.3778) + \
1778 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.3622) + \
1779 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.3299) + \
1780 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.2889) + \
1781 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.2815) + \
1782 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.253) + \
1783 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.221) + \
1784 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.1963) + \
1785 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.1739) + \
1786 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.0992) + \
1787 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.0788) + \
1788 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.0581) + \
1789 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.0534) + \
1790 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.0521) + \
1791 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.0626) + \
1792 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.0826) + \
1793 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.1022) + \
1794 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1795 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.1562) + \
1796 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.2499) + \
1797 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.2956) + \
1798 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.3128) + \
1799 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.3147) + \
1800 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.3222) + \
1801 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.304) + \
1802 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.3051) + \
1803 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.2657) + \
1804 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.2578) + \
1805 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.2087) + \
1806 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.1634) + \
1807 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.1651) + \
1808 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.1353) + \
1809 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.109) + \
1810 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.0799) + \
1811 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.0699) + \
1812 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.054) + \
1813 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.0718) + \
1814 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.0746) + \
1815 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.0803) + \
1816 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.0882) + \
1817 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1818 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1819 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \
1820 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \
1821 (abs(eta) > 4.0) * (0.000)
1822 }
1823
1824}
1825
1826
1827
1828#####################
1829## tau tagging ##
1830#####################
1831
1832
1833module TauTagging TauTagging {
1834
1835 set ParticleInputArray Delphes/allParticles
1836 set PartonInputArray Delphes/partons
1837 set JetInputArray RecoJetPileUpSubtractor/jets
1838
1839 set DeltaR 0.4
1840 set TauPTMin 1.0
1841 set TauEtaMax 4.0
1842
1843 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1844 # default efficiency formula (misidentification rate)
1845
1846 add EfficiencyFormula {0} {(abs(eta)<1.8)*0.006+(abs(eta)>1.8)*0.015}
1847
1848 # efficiency formula for tau-jets
1849
1850 add EfficiencyFormula {15} {0.65}
1851}
1852
1853
1854#####################
1855## MissingET ##
1856#####################
1857
1858### -sum of all particles after filtering neutrinos
1859
1860module Merger GenMissingET {
1861 add InputArray NeutrinoFilter/filteredParticles
1862 set MomentumOutputArray momentum
1863}
1864
1865module Merger RecoMissingET {
1866 add InputArray EFlowMerger/eflow
1867 set MomentumOutputArray momentum
1868}
1869
1870module Merger JetPileUpIDMissingET {
1871 add InputArray TrackPileUpSubtractor/eflowTracks
1872 add InputArray RecoJetPileUpID/eflowTowers
1873 set MomentumOutputArray momentum
1874}
1875
1876module Merger PuppiMissingET {
1877 add InputArray RunPUPPI/PuppiParticles
1878 set MomentumOutputArray momentum
1879}
1880
1881
1882##############
1883## Scalar HT #
1884##############
1885
1886module Merger GenScalarHT {
1887 # add InputArray InputArray
1888 add InputArray NeutrinoFilter/filteredParticles
1889 set EnergyOutputArray energy
1890}
1891
1892module Merger RecoScalarHT {
1893 # add InputArray InputArray
1894 add InputArray EFlowMerger/eflow
1895 set EnergyOutputArray energy
1896}
1897
1898module Merger JetPileUpIDScalarHT {
1899 add InputArray TrackPileUpSubtractor/eflowTracks
1900 add InputArray RecoJetPileUpID/eflowTowers
1901 set MomentumOutputArray energy
1902}
1903
1904module Merger PuppiScalarHT {
1905 # add InputArray InputArray
1906 add InputArray RunPUPPI/PuppiParticles
1907 set EnergyOutputArray energy
1908}
1909
1910
1911#####################
1912# Photon efficiency #
1913#####################
1914
1915module Efficiency PhotonEfficiency {
1916
1917 ## input particles
1918 set InputArray Calorimeter/eflowPhotons
1919 ## output particles
1920 set OutputArray photons
1921 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1922 # efficiency formula for photons
1923 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
1924 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
1925 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
1926 (abs(eta) > 4.0) * (0.00)}
1927
1928}
1929
1930####################
1931# Photon isolation #
1932####################
1933
1934module Isolation PhotonIsolation {
1935
1936 # particle for which calculate the isolation
1937 set CandidateInputArray PhotonEfficiency/photons
1938
1939 # neutral and charged particles for the whole event (no CHS applied)
1940 set NeutralIsolationInputArray NeutralTowerMerger/eflowTowers
1941 set ChargedIsolationInputArray TrackMerger/tracks
1942
1943 # select a rho for the isolation
1944 set RhoInputArray RhoGridFastJet/rho
1945
1946 # output array
1947 set OutputArray photons
1948
1949 # isolation cone
1950 set DeltaRMax 0.3
1951
1952 # minimum pT
1953 set PTMin 0.5
1954
1955 # iso ratio to cut
1956 set PTRatioMax 9999.
1957}
1958
1959
1960#######################
1961# Electron efficiency #
1962#######################
1963
1964module Efficiency ElectronEfficiency {
1965
1966 set InputArray TrackPileUpSubtractor/electrons
1967 set OutputArray electrons
1968
1969 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1970 # efficiency formula for electrons
1971 set EfficiencyFormula {
1972 (pt <= 4.0) * (0.00) + \
1973 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
1974 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
1975 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
1976 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
1977 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
1978 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
1979 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
1980 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
1981 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
1982 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
1983 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
1984 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
1985 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
1986 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
1987 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
1988 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
1989 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1990 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1991 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1992 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1993 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
1994 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1995 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1996 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1997 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1998 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
1999 (abs(eta) > 4.0) * (0.00)
2000
2001 }
2002}
2003
2004######################
2005# Electron isolation #
2006######################
2007
2008module Isolation ElectronIsolation {
2009
2010 set CandidateInputArray ElectronEfficiency/electrons
2011
2012 set NeutralIsolationInputArray NeutralTowerMerger/eflowTowers
2013 set ChargedIsolationInputArray TrackMerger/tracks
2014
2015 set RhoInputArray RhoGridFastJet/rho
2016
2017 set OutputArray electrons
2018
2019 set DeltaRMax 0.3
2020 set PTMin 0.5
2021 set PTRatioMax 9999.
2022
2023}
2024
2025###################
2026# Muon efficiency #
2027###################
2028
2029module Efficiency MuonEfficiency {
2030
2031 set InputArray TrackPileUpSubtractor/muons
2032
2033 set OutputArray muons
2034 # set EfficiencyFormula {efficiency as a function of eta and pt}
2035 # efficiency formula for muons
2036 set EfficiencyFormula { (pt <= 2.0) * (0.00) + \
2037 (abs(eta) <= 4.00) * (pt > 2.0 && pt <= 3.0) * (0.51) + \
2038 (abs(eta) <= 4.00) * (pt > 3.0 && pt <= 4.0) * (0.85) + \
2039 (abs(eta) <= 4.00) * (pt > 4.0 && pt <= 11.0) * (0.93) + \
2040 (abs(eta) <= 4.00) * (pt > 11. && pt <= 50.) * (0.96) + \
2041 (abs(eta) <= 4.00) * (pt > 50. && pt <= 70.) * (0.98) + \
2042 (abs(eta) <= 4.00) * (pt > 70.0 ) * (1.00) + \
2043 (abs(eta) > 4.00) * (0.00)
2044 }
2045
2046}
2047
2048##################
2049# Muon isolation #
2050##################
2051
2052module Isolation MuonIsolation {
2053 set CandidateInputArray MuonEfficiency/muons
2054
2055 set NeutralIsolationInputArray NeutralTowerMerger/eflowTowers
2056 set ChargedIsolationInputArray TrackMerger/tracks
2057
2058 set RhoInputArray RhoGridFastJet/rho
2059
2060 set OutputArray muons
2061
2062 set DeltaRMax 0.3
2063 set PTMin 0.5
2064 set PTRatioMax 9999.
2065
2066}
2067
2068
2069###############################################################################################################
2070# StatusPidFilter: this module removes all generated particles except electrons, muons, taus, and status == 3 #
2071###############################################################################################################
2072
2073module StatusPidFilter StatusPidFilter {
2074 ## take the particles from Pythia8 not adding pile-up
2075 set InputArray Delphes/allParticles
2076 set OutputArray filteredParticles
2077 set PTMin 5.0
2078
2079}
2080
2081##################
2082# ROOT tree writer
2083##################
2084
2085module TreeWriter TreeWriter {
2086
2087 ## branch notation : <particle collection> <branch name> <type of object in classes/DelphesClass.h<
2088
2089 # add Branch Delphes/stableParticles GenParticles GenParticle
2090
2091 ## gen particle stuff
2092 add Branch StatusPidFilter/filteredParticles GenParticles GenParticle
2093 add Branch PileUpMerger/vertices Vertex Vertex
2094
2095 ## Detector deposits
2096
2097 #add Branch TrackMerger/tracks Track Track
2098 #add Branch Calorimeter/towers Tower Tower
2099
2100 #add Branch Calorimeter/eflowTracks EFlowTrack Track
2101 #add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
2102 #add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
2103
2104 add Branch RhoGridFastJet/rho RhoGridFastJet Rho
2105 add Branch RhoPuppiGridFastJet/rho RhoPuppiGridFastJet Rho
2106
2107 ## Jets
2108 #add Branch GenJetFinder/jets GenJet Jet
2109 #add Branch RecoJetFinder/jets RecoJet Jet
2110 #add Branch PuppiJetFinder/jets PuppiJet Jet
2111
2112
2113 add Branch RecoJetPileUpID/jets RecoJetPileUpID Jet
2114 add Branch PuppiJetPileUpID/jets PuppiJetPileUpID Jet
2115
2116 ## MET
2117 add Branch GenMissingET/momentum GenMissingET MissingET
2118 add Branch RecoMissingET/momentum RecoMissingET MissingET
2119 add Branch JetPileUpIDMissingET/momentum JetPileUpIDMissingET MissingET
2120 add Branch PuppiMissingET/momentum PuppiMissingET MissingET
2121
2122 ## HT
2123 add Branch GenScalarHT/energy GenHT ScalarHT
2124 add Branch RecoScalarHT/energy RecoHT ScalarHT
2125 add Branch JetPileUpIDScalarHT/energy JetPileUpIDHT ScalarHT
2126 add Branch PuppiScalarHT/energy PuppiHT ScalarHT
2127
2128 ## photons and leptons
2129 add Branch ElectronIsolation/electrons Electron Electron
2130 add Branch PhotonIsolation/photons Photon Photon
2131 add Branch MuonIsolation/muons Muon Muon
2132
2133}
Note: See TracBrowser for help on using the repository browser.