Fork me on GitHub

source: git/cards/CMS_PhaseII_140PU_conf4.tcl@ 7442f9b

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

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