Fork me on GitHub

source: git/cards/CMS_PhaseII_140PU_conf4.tcl@ fc2ef70

Last change on this file since fc2ef70 was ff9f67fb, checked in by Michele Selvaggi <michele.selvaggi@…>, 9 years ago

cards/CMS_PhaseII_0PU.tcl
  • Property mode set to 100644
File size: 67.2 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 ElectronFilter
24
25 TrackPileUpSubtractor
26
27 NeutralTowerMerger
28 EFlowMergerAllTracks
29 EFlowMerger
30
31 NeutrinoFilter
32
33 Rho
34
35 GenJetFinder
36 RecoJetFinder
37
38 JetPileUpSubtractor
39
40 JetFlavorAssociation
41
42 RecoJetBTaggingLoose
43 RecoJetBTaggingMedium
44 RecoJetBTaggingTight
45
46 GenMissingET
47 RecoMissingET
48
49 GenScalarHT
50 RecoScalarHT
51
52 PhotonEfficiency
53 PhotonIsolation
54
55 ElectronEfficiency
56 ElectronIsolation
57
58 MuonEfficiency
59 MuonIsolation
60
61 TauTagging
62
63 StatusPidFilter
64
65 TreeWriter
66
67}
68
69
70
71###############
72# PileUp Merger
73###############
74
75module PileUpMerger PileUpMerger {
76 set InputArray Delphes/stableParticles
77
78 set ParticleOutputArray stableParticles
79 set VertexOutputArray vertices
80
81 # pre-generated minbias input file
82 set PileUpFile MinBias.pileup
83
84 # Get rid of beam spot from http://red-gridftp11.unl.edu/Snowmass/MinBias100K_14TeV.pileup ...
85 set InputBSX 2.44
86 set InputBSY 3.39
87 # replace it with beam spot from CMSSW files
88 set OutputBSX 0.24
89 set OutputBSY 0.39
90
91 # average expected pile up
92 set MeanPileUp 140
93
94 # maximum spread in the beam direction in m
95 set ZVertexSpread 0.15
96
97 # maximum spread in time in s
98 set TVertexSpread 1.5E-09
99
100 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
101
102 #set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
103
104
105 set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) +
106 (abs(t) > 1.0e-09) * (abs(z) <= 0.15) * (0.00) +
107 (abs(t) <= 1.0e-09) * (abs(z) > 0.15) * (0.00) +
108 (abs(t) > 1.0e-09) * (abs(z) > 0.15) * (0.00)}
109
110
111}
112
113#####################################################################################
114# Propagate particles in cylinder and divide charged particles in different classes #
115#####################################################################################
116
117module ParticlePropagator ParticlePropagator {
118 ## take particles after beam spot smearing
119 set InputArray PileUpMerger/stableParticles
120 ## produce independent output collection: all particles, only charged hadrons, electrons and muons
121 set OutputArray stableParticles
122 set ChargedHadronOutputArray chargedHadrons
123 set ElectronOutputArray electrons
124 set MuonOutputArray muons
125 ## radius of the magnetic field coverage, in m
126 set Radius 1.29
127 ## half-length of the magnetic field coverage, in m
128 set HalfLength 3.00
129 ## magnetic field
130 set Bz 3.8
131}
132
133####################################
134# Charged hadron tracking efficiency
135####################################
136
137module Efficiency ChargedHadronTrackingEfficiency {
138 ## particles after propagation
139 set InputArray ParticlePropagator/chargedHadrons
140 set OutputArray chargedHadrons
141 # tracking efficiency formula for charged hadrons
142 set EfficiencyFormula {
143 (pt <= 0.2) * (0.00) + \
144 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
145 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
146 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
147 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.87) + \
148 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
149 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.82) + \
150 (abs(eta) > 4.0) * (0.00)
151 }
152}
153
154
155#####################################
156# Electron tracking efficiency - ID
157####################################
158
159module Efficiency ElectronTrackingEfficiency {
160 set InputArray ParticlePropagator/electrons
161 set OutputArray electrons
162 # tracking efficiency formula for electrons
163 set EfficiencyFormula {
164 (pt <= 0.2) * (0.00) + \
165 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.96) + \
166 (abs(eta) <= 1.2) * (pt > 1.0) * (0.97) + \
167 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.85) + \
168 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 10.0) * (0.82+pt*0.01) + \
169 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 10.0) * (0.90) + \
170 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.8) + \
171 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.8+pt*0.01) + \
172 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.85) + \
173 (abs(eta) > 4.0) * (0.00)
174
175 }
176}
177
178##########################
179# Muon tracking efficiency
180##########################
181
182module Efficiency MuonTrackingEfficiency {
183 set InputArray ParticlePropagator/muons
184 set OutputArray muons
185 # tracking efficiency formula for muons
186 set EfficiencyFormula {
187 (pt <= 0.2) * (0.00) + \
188 (abs(eta) <= 1.2) * (pt > 0.2 && pt <= 1.0) * (pt * 0.998) + \
189 (abs(eta) <= 1.2) * (pt > 1.0) * (0.998) + \
190 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 0.2 && pt <= 1.0) * (pt*0.99) + \
191 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 1.0) * (0.99) + \
192 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 0.2 && pt <= 1.0) * (pt*0.95) + \
193 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0) * (0.95) + \
194 (abs(eta) > 4.0) * (0.00)
195
196 }
197}
198
199########################################
200# Momentum resolution for charged tracks
201########################################
202
203module MomentumSmearing ChargedHadronMomentumSmearing {
204 ## hadrons after having applied the tracking efficiency
205 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
206 set OutputArray chargedHadrons
207 # resolution formula for charged hadrons
208 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
209 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.013) + \
210 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.02) + \
211 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
212 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.015) + \
213 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.015) + \
214 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
215 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05) + \
216 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 0.1 && pt <= 1.0) * (0.017) + \
217 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 1.0 && pt <= 10.0) * (0.03) + \
218 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 100.0) * (0.05) + \
219 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 100.0) * (0.30) + \
220 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
221 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 1.0 && pt <= 10.0) * (0.04) + \
222 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 10.0 && pt <= 100.0) * (0.07) + \
223 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 100.0) * (0.30) + \
224 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.025) + \
225 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.05) + \
226 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 100.0) * (0.20) + \
227 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 100.0) * (0.80)
228 }
229}
230
231#################################
232# Energy resolution for electrons
233#################################
234
235module EnergySmearing ElectronEnergySmearing {
236 set InputArray ElectronTrackingEfficiency/electrons
237 set OutputArray electrons
238 # set ResolutionFormula {resolution formula as a function of eta and energy}
239 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.02) + \
240 (abs(eta) <= 2.5) * (energy > 2.5e1) * (energy*0.028) + \
241 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (energy*0.032) + \
242 (abs(eta) > 3.0 && abs(eta) <= 5.0) * (energy*0.07)
243 }
244}
245
246###############################
247# Momentum resolution for muons
248###############################
249
250module MomentumSmearing MuonMomentumSmearing {
251 set InputArray MuonTrackingEfficiency/muons
252 set OutputArray muons
253 # resolution formula for muons
254 set ResolutionFormula {
255 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.012) + \
256 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
257 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.012) + \
258 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.030) + \
259 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.011) + \
260 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.011) + \
261 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.022) + \
262 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.030) + \
263 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 0.1 && pt <= 1.0) * (0.017) + \
264 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 1.0 && pt <= 10.0) * (0.03) + \
265 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 10.0 && pt <= 100.0) * (0.05) + \
266 (abs(eta) > 2.5 && abs(eta) <= 3.0) * (pt > 100.0) * (0.30) + \
267 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
268 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 1.0 && pt <= 10.0) * (0.04) + \
269 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 10.0 && pt <= 100.0) * (0.07) + \
270 (abs(eta) > 3.0 && abs(eta) <= 3.5) * (pt > 100.0) * (0.30) + \
271 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 0.1 && pt <= 1.0) * (0.025) + \
272 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt <= 10.0) * (0.05) + \
273 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 10.0 && pt <= 100.0) * (0.20) + \
274 (abs(eta) > 3.5 && abs(eta) <= 4.0) * (pt > 100.0) * (0.80)
275
276 }
277}
278
279
280##############
281# Track merger
282##############
283
284module Merger TrackMerger {
285# add InputArray InputArray
286 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
287 add InputArray ElectronEnergySmearing/electrons
288 add InputArray MuonMomentumSmearing/muons
289 set OutputArray tracks
290}
291
292
293##############################################################
294# Calorimeter : emulate calorimiter answer making caloTowers #
295##############################################################
296
297module Calorimeter Calorimeter {
298 ## particle from the propagation without any efficiency or smearing (for neutrals)
299 set ParticleInputArray ParticlePropagator/stableParticles
300 ## track after smearing and efficiency: used for charged particles
301 set TrackInputArray TrackMerger/tracks
302 ## output collections
303
304 set TowerOutputArray towers
305 set PhotonOutputArray photons
306
307 set EFlowTrackOutputArray eflowTracks
308 set EFlowPhotonOutputArray eflowPhotons
309 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
310
311 set ECalEnergyMin 0.0
312 set HCalEnergyMin 0.0
313
314 set ECalEnergySignificanceMin 0.0
315 set HCalEnergySignificanceMin 0.0
316
317 set SmearTowerCenter false
318
319 set pi [expr {acos(-1)}]
320
321 # lists of the edges of each tower in eta and phi
322 # each list starts with the lower edge of the first tower
323 # the list ends with the higher edged of the last tower
324
325 # 5 degrees towers
326 set PhiBins {}
327 for {set i -70} {$i <= 70} {incr i} {
328 add PhiBins [expr {$i * $pi/70.0}]
329 }
330
331 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} {
332 add EtaPhiBins $eta $PhiBins
333 }
334
335 # 10 degrees towers
336 set PhiBins {}
337 for {set i -11} {$i <= 11} {incr i} {
338 add PhiBins [expr {$i * $pi/11.0}]
339 }
340 foreach eta {-4.35 -4.175 -4 4.175 4.35 4.525} {
341 add EtaPhiBins $eta $PhiBins
342 }
343
344 # 20 degrees towers
345 set PhiBins {}
346 for {set i -9} {$i <= 9} {incr i} {
347 add PhiBins [expr {$i * $pi/9.0}]
348 }
349 foreach eta {-5 -4.7 -4.525 4.7 5} {
350 add EtaPhiBins $eta $PhiBins
351 }
352
353 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
354 add EnergyFraction {0} {0.0 1.0}
355 # energy fractions for e, gamma and pi0
356 add EnergyFraction {11} {1.0 0.0}
357 add EnergyFraction {22} {1.0 0.0}
358 add EnergyFraction {111} {1.0 0.0}
359 # energy fractions for muon, neutrinos and neutralinos
360 add EnergyFraction {12} {0.0 0.0}
361 add EnergyFraction {13} {0.0 0.0}
362 add EnergyFraction {14} {0.0 0.0}
363 add EnergyFraction {16} {0.0 0.0}
364 add EnergyFraction {1000022} {0.0 0.0}
365 add EnergyFraction {1000023} {0.0 0.0}
366 add EnergyFraction {1000025} {0.0 0.0}
367 add EnergyFraction {1000035} {0.0 0.0}
368 add EnergyFraction {1000045} {0.0 0.0}
369 # energy fractions for K0short and Lambda
370 add EnergyFraction {310} {0.3 0.7}
371 add EnergyFraction {3122} {0.3 0.7}
372
373
374 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
375 set ECalResolutionFormula {(abs(eta) <= 1.5) * sqrt(energy^2*0.005^2 + energy*0.027^2 + 0.15^2) + \
376 (abs(eta) > 1.5 && abs(eta) <= 4.0) * sqrt(energy^2*0.008^2 + energy*0.092^2 + 0.088^2) + \
377 (abs(eta) > 4.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.97^2)}
378
379
380
381 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
382 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
383 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.050^2 + energy*0.706^2) + \
384 (abs(eta) > 3.0 && abs(eta) <= 4.9) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
385
386
387
388}
389
390
391#################
392# Electron filter
393#################
394
395module PdgCodeFilter ElectronFilter {
396 set InputArray Calorimeter/eflowTracks
397 set OutputArray electrons
398 set Invert true
399 add PdgCode {11}
400 add PdgCode {-11}
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 ElectronFilter/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# Rho pile-up
515#############
516
517module FastJetFinder Rho {
518# set InputArray Calorimeter/towers
519 set InputArray EFlowMerger/eflow
520
521 set ComputeRho true
522 set RhoOutputArray rho
523
524 # 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
525 set AreaAlgorithm 5
526
527 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
528 set JetAlgorithm 4
529 set ParameterR 0.4
530 set GhostEtaMax 5.0
531 set RhoEtaMax 5.0
532
533 add RhoEtaRange -5.0 -4.0
534 add RhoEtaRange -4.0 -2.5
535 add RhoEtaRange -2.5 2.5
536 add RhoEtaRange 2.5 4.0
537 add RhoEtaRange 4.0 5.0
538
539 set JetPTMin 0.0
540}
541
542
543#################################
544# Rho pile-up with std particles
545#################################
546
547module FastJetGridMedianEstimator RhoGridFastJet {
548
549 set InputArray EFlowMerger/eflow
550 set RhoOutputArray rho
551
552 # add GridRange rapmin rapmax drap dphi
553 # rapmin - the minimum rapidity extent of the grid
554 # rapmax - the maximum rapidity extent of the grid
555 # drap - the grid spacing in rapidity
556 # dphi - the grid spacing in azimuth
557
558 add GridRange -5.0 -4.0 1.0 1.0
559 add GridRange -4.0 -2.5 1.0 1.0
560 add GridRange -2.5 2.5 1.0 1.0
561 add GridRange 2.5 4.0 1.0 1.0
562 add GridRange 4.0 5.0 1.0 1.0
563
564}
565
566###################################
567# Rho pile-up with PUPPI particles
568###################################
569
570module FastJetGridMedianEstimator RhoPuppiGridFastJet {
571
572 set InputArray RunPUPPI/PuppiParticles
573 set RhoOutputArray rho
574
575 # add GridRange rapmin rapmax drap dphi
576 # rapmin - the minimum rapidity extent of the grid
577 # rapmax - the maximum rapidity extent of the grid
578 # drap - the grid spacing in rapidity
579 # dphi - the grid spacing in azimuth
580
581 add GridRange -5.0 -4.0 1.0 1.0
582 add GridRange -4.0 -2.5 1.0 1.0
583 add GridRange -2.5 2.5 1.0 1.0
584 add GridRange 2.5 4.0 1.0 1.0
585 add GridRange 4.0 5.0 1.0 1.0
586
587}
588
589
590#####################
591# MC truth jet finder
592#####################
593
594module FastJetFinder GenJetFinder {
595 set InputArray NeutrinoFilter/filteredParticles
596
597 set OutputArray jets
598
599 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
600 set JetAlgorithm 6
601 set ParameterR 0.4
602
603 set JetPTMin 15.0
604}
605
606
607
608############
609# Jet finder
610############
611
612module FastJetFinder RecoJetFinder {
613# set InputArray Calorimeter/towers
614 set InputArray EFlowMerger/eflow
615
616 set OutputArray jets
617
618 # 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
619 set AreaAlgorithm 5
620
621 # jet algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
622 set JetAlgorithm 6
623 set ParameterR 0.4
624
625 set JetPTMin 5.0
626}
627
628
629###########################
630# Jet Pile-Up Subtraction
631###########################
632
633module JetPileUpSubtractor JetPileUpSubtractor {
634 set JetInputArray RecoJetFinder/jets
635 set RhoInputArray Rho/rho
636
637 set OutputArray jets
638
639 set JetPTMin 5.0
640}
641
642
643########################
644# Jet Flavor Association
645########################
646
647module JetFlavorAssociation JetFlavorAssociation {
648
649 set PartonInputArray Delphes/partons
650 set ParticleInputArray Delphes/allParticles
651 set ParticleLHEFInputArray Delphes/allParticlesLHEF
652 set JetInputArray JetPileUpSubtractor/jets
653
654 set DeltaR 0.5
655 set PartonPTMin 1.0
656 set PartonEtaMax 4.0
657
658}
659
660
661#############
662# b-tagging #
663#############
664module BTagging RecoJetBTaggingLoose {
665
666 set JetInputArray JetPileUpSubtractor/jets
667
668 set BitNumber 0
669
670 add EfficiencyFormula {0} {(pt <= 20.0) * (0.000) + \
671 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.096) + \
672 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.112) + \
673 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.0812) + \
674 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.0908) + \
675 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.0792) + \
676 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.0858) + \
677 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.0917) + \
678 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.084) + \
679 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.0906) + \
680 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.0989) + \
681 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.1022) + \
682 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.1035) + \
683 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.1095) + \
684 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.1201) + \
685 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.1348) + \
686 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.1482) + \
687 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.1629) + \
688 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.1775) + \
689 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.2002) + \
690 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.1977) + \
691 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.2084) + \
692 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.2195) + \
693 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.2424) + \
694 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.2909) + \
695 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.3457) + \
696 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
697 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.074765) + \
698 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.100053) + \
699 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.071492) + \
700 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.084796) + \
701 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.076927) + \
702 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.08424) + \
703 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.093118) + \
704 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.084629) + \
705 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.092977) + \
706 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.10206) + \
707 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.102344) + \
708 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.098435) + \
709 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.105507) + \
710 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.112841) + \
711 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.126329) + \
712 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.140759) + \
713 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.153193) + \
714 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.107869) + \
715 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.119527) + \
716 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.08688) + \
717 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.089324) + \
718 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.097172) + \
719 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.000) + \
720 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.000) + \
721 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.000) + \
722 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (0.02) +\
723 (abs(eta) > 4.0) * (0.000) \
724 }
725
726 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
727 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.387) + \
728 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.448) + \
729 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.408) + \
730 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.427) + \
731 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.408) + \
732 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.425) + \
733 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.426) + \
734 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.4) + \
735 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.415) + \
736 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.416) + \
737 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.405) + \
738 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.387) + \
739 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.39) + \
740 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.389) + \
741 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.389) + \
742 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.381) + \
743 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.381) + \
744 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.367) + \
745 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.369) + \
746 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.326) + \
747 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.335) + \
748 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.326) + \
749 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.341) + \
750 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.403) + \
751 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.47) + \
752 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
753 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.2497) + \
754 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.31891) + \
755 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.273) + \
756 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.28445) + \
757 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.28036) + \
758 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.28453) + \
759 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.29495) + \
760 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.27024) + \
761 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.28912) + \
762 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.29048) + \
763 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.27507) + \
764 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.23327) + \
765 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.2493) + \
766 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.2416) + \
767 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.26652) + \
768 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.24852) + \
769 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.26927) + \
770 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.19302) + \
771 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.19433) + \
772 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.17523) + \
773 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.14981) + \
774 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.16666) + \
775 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
776 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
777 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.) + \
778 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.29*tanh(pt*0.00997 - 0.143)) + \
779 (abs(eta) > 4.0) * (0.000)
780 }
781
782 # efficiency formula for b-jets
783 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
784 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.75) + \
785 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.827) + \
786 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.837) + \
787 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.853) + \
788 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.855) + \
789 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.862) + \
790 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.868) + \
791 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.865) + \
792 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.863) + \
793 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.857) + \
794 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.851) + \
795 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.838) + \
796 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.831) + \
797 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.817) + \
798 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.796) + \
799 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.772) + \
800 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.76) + \
801 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.743) + \
802 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.703) + \
803 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.638) + \
804 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.605) + \
805 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.572) + \
806 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.541) + \
807 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.567) + \
808 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.603) + \
809 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
810 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.6063) + \
811 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.7188) + \
812 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.72) + \
813 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.7365) + \
814 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.7462) + \
815 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.7454) + \
816 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.7415) + \
817 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.727) + \
818 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.7112) + \
819 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.7112) + \
820 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.6754) + \
821 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.6359) + \
822 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.6348) + \
823 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.6115) + \
824 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.5585) + \
825 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.5608) + \
826 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.5208) + \
827 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.456) + \
828 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.4524) + \
829 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.388) + \
830 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.3928) + \
831 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.3823) + \
832 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
833 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
834 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \
835 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \
836 (abs(eta) > 4.0) * (0.000)
837 }
838}
839
840
841#############
842# b-tagging #
843#############
844module BTagging RecoJetBTaggingMedium {
845
846 set JetInputArray JetPileUpSubtractor/jets
847
848 set BitNumber 1
849
850 add EfficiencyFormula {0} { (pt <= 20.0) * (0.000) + \
851 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.00654) + \
852 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.00921) + \
853 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.00573) + \
854 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.00694) + \
855 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.0062) + \
856 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.00708) + \
857 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.00779) + \
858 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.00693) + \
859 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.00777) + \
860 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.00862) + \
861 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.01038) + \
862 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.01189) + \
863 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.01279) + \
864 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.01452) + \
865 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.01696) + \
866 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.01958) + \
867 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.02253) + \
868 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.01787) + \
869 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.02154) + \
870 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.01839) + \
871 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.01987) + \
872 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.02351) + \
873 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.02937) + \
874 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.04001) + \
875 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.0542) + \
876 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
877 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.004236) + \
878 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.006653) + \
879 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.005512) + \
880 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.00754) + \
881 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.005813) + \
882 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.006439) + \
883 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.008063) + \
884 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.00647) + \
885 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.007583) + \
886 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.008543) + \
887 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.01034) + \
888 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.011253) + \
889 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.012945) + \
890 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.014474) + \
891 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.017361) + \
892 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.020912) + \
893 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.023139) + \
894 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.010756) + \
895 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.012569) + \
896 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.006046) + \
897 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.006428) + \
898 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.00887) + \
899 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
900 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
901 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.000) + \
902 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (0.001) + \
903 (abs(eta) > 4.0) * (0.000) \
904}
905
906 # efficiency formula for c-jets (misidentification rate)
907 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
908 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.1102) + \
909 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.1344) + \
910 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.1025) + \
911 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.1025) + \
912 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.1063) + \
913 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.1087) + \
914 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.1124) + \
915 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.109) + \
916 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.111) + \
917 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.1091) + \
918 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.1087) + \
919 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.1091) + \
920 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.1107) + \
921 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.1061) + \
922 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.1017) + \
923 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.1) + \
924 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.0966) + \
925 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.0697) + \
926 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.0679) + \
927 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.0503) + \
928 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.0514) + \
929 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.0481) + \
930 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.0667) + \
931 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.0861) + \
932 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.092) + \
933 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
934 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.03331) + \
935 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.04361) + \
936 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.03863) + \
937 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.04287) + \
938 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.04431) + \
939 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.04452) + \
940 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.04339) + \
941 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.0436) + \
942 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.0456) + \
943 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.05138) + \
944 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.04794) + \
945 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.04004) + \
946 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.04713) + \
947 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.04515) + \
948 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.05314) + \
949 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.05143) + \
950 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.05936) + \
951 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.02357) + \
952 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.03222) + \
953 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.01523) + \
954 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.02621) + \
955 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.01709) + \
956 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
957 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
958 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.) + \
959 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.1898*tanh(pt*0.00997 - 0.143)) + \
960 (abs(eta) > 4.0) * (0.000)
961 }
962
963 # efficiency formula for b-jets
964 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
965 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.536) + \
966 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.6439) + \
967 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.6504) + \
968 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.6716) + \
969 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.6841) + \
970 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.6896) + \
971 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.6916) + \
972 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.6882) + \
973 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.6838) + \
974 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.6715) + \
975 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.6554) + \
976 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.6366) + \
977 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.6192) + \
978 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.595) + \
979 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.5551) + \
980 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.5138) + \
981 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.4884) + \
982 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.4009) + \
983 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.3459) + \
984 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.2523) + \
985 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.2404) + \
986 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.2198) + \
987 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.2263) + \
988 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.2614) + \
989 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.3194) + \
990 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
991 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.3254) + \
992 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.4339) + \
993 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.4499) + \
994 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.4716) + \
995 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.4766) + \
996 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.4788) + \
997 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.4863) + \
998 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.4891) + \
999 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.462) + \
1000 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.4583) + \
1001 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.4247) + \
1002 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.3775) + \
1003 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.3734) + \
1004 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.3348) + \
1005 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.2939) + \
1006 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.285) + \
1007 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.2421) + \
1008 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.1565) + \
1009 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.1522) + \
1010 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.1231) + \
1011 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.1607) + \
1012 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.1323) + \
1013 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1014 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1015 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \
1016 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \
1017 (abs(eta) > 4.0) * (0.000)
1018 }
1019}
1020
1021#############
1022# b-tagging #
1023#############
1024module BTagging RecoJetBTaggingTight {
1025
1026 set JetInputArray JetPileUpSubtractor/jets
1027
1028 set BitNumber 2
1029
1030 add EfficiencyFormula {0} { (pt <= 20.0) * (0.000) + \
1031 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.000164) + \
1032 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.000235) + \
1033 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.000266) + \
1034 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.000329) + \
1035 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.000309) + \
1036 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.000309) + \
1037 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.000546) + \
1038 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.000499) + \
1039 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.000642) + \
1040 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.000742) + \
1041 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.000928) + \
1042 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.001323) + \
1043 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.001392) + \
1044 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.00154) + \
1045 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.002094) + \
1046 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.002427) + \
1047 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.002927) + \
1048 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.001854) + \
1049 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.002355) + \
1050 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.002297) + \
1051 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.002433) + \
1052 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.002706) + \
1053 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.003602) + \
1054 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.004987) + \
1055 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.007414) + \
1056 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1057 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.000573) + \
1058 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.000574) + \
1059 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.000521) + \
1060 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.000786) + \
1061 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.000539) + \
1062 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.000673) + \
1063 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.000934) + \
1064 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.000781) + \
1065 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.000949) + \
1066 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.000977) + \
1067 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.001168) + \
1068 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.000879) + \
1069 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.000812) + \
1070 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.001215) + \
1071 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.001679) + \
1072 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.001893) + \
1073 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.002723) + \
1074 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.003555) + \
1075 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.003881) + \
1076 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.006046) + \
1077 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.005563) + \
1078 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.007611) + \
1079 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1080 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1081 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.000) + \
1082 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (0.001) + \
1083 (abs(eta) > 4.0) * (0.000) \
1084 }
1085
1086 # efficiency formula for c-jets (misidentification rate)
1087 add EfficiencyFormula {4} { (pt <= 20.0) * (0.000) + \
1088 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.00531) + \
1089 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.00567) + \
1090 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.0064) + \
1091 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.00673) + \
1092 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.00766) + \
1093 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.00729) + \
1094 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.00674) + \
1095 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.00824) + \
1096 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.00888) + \
1097 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.00919) + \
1098 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.01021) + \
1099 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.01041) + \
1100 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.01027) + \
1101 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.00999) + \
1102 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.01047) + \
1103 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.01014) + \
1104 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.01021) + \
1105 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.00601) + \
1106 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.0054) + \
1107 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.00487) + \
1108 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.00519) + \
1109 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.00469) + \
1110 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.00651) + \
1111 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.01299) + \
1112 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.00897) + \
1113 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1114 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.01014) + \
1115 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.01288) + \
1116 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.01392) + \
1117 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.01533) + \
1118 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.01508) + \
1119 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.01579) + \
1120 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.01106) + \
1121 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.01346) + \
1122 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.01315) + \
1123 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.01156) + \
1124 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.0082) + \
1125 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.00439) + \
1126 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.00744) + \
1127 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.00685) + \
1128 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.00755) + \
1129 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.00706) + \
1130 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.00428) + \
1131 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.01031) + \
1132 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.00976) + \
1133 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.01523) + \
1134 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.00374) + \
1135 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.00854) + \
1136 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1137 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1138 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \
1139 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \
1140 (abs(eta) > 4.0) * (0.000)
1141}
1142
1143 # efficiency formula for b-jets
1144 add EfficiencyFormula {5} { (pt <= 20.0) * (0.000) + \
1145 (abs(eta) <= 1.8) * (pt > 20.0 && pt <= 30) * (0.2426) + \
1146 (abs(eta) <= 1.8) * (pt > 30.0 && pt <= 40) * (0.327) + \
1147 (abs(eta) <= 1.8) * (pt > 40.0 && pt <= 50) * (0.3559) + \
1148 (abs(eta) <= 1.8) * (pt > 50.0 && pt <= 60) * (0.3704) + \
1149 (abs(eta) <= 1.8) * (pt > 60.0 && pt <= 70) * (0.3824) + \
1150 (abs(eta) <= 1.8) * (pt > 70.0 && pt <= 80) * (0.3844) + \
1151 (abs(eta) <= 1.8) * (pt > 80.0 && pt <= 90) * (0.3848) + \
1152 (abs(eta) <= 1.8) * (pt > 90.0 && pt <= 100) * (0.3862) + \
1153 (abs(eta) <= 1.8) * (pt > 100.0 && pt <= 120) * (0.3778) + \
1154 (abs(eta) <= 1.8) * (pt > 120.0 && pt <= 140) * (0.3622) + \
1155 (abs(eta) <= 1.8) * (pt > 140.0 && pt <= 160) * (0.3299) + \
1156 (abs(eta) <= 1.8) * (pt > 160.0 && pt <= 180) * (0.2889) + \
1157 (abs(eta) <= 1.8) * (pt > 180.0 && pt <= 200) * (0.2815) + \
1158 (abs(eta) <= 1.8) * (pt > 200.0 && pt <= 250) * (0.253) + \
1159 (abs(eta) <= 1.8) * (pt > 250.0 && pt <= 300) * (0.221) + \
1160 (abs(eta) <= 1.8) * (pt > 300.0 && pt <= 350) * (0.1963) + \
1161 (abs(eta) <= 1.8) * (pt > 350.0 && pt <= 400) * (0.1739) + \
1162 (abs(eta) <= 1.8) * (pt > 400.0 && pt <= 500) * (0.0992) + \
1163 (abs(eta) <= 1.8) * (pt > 500.0 && pt <= 600) * (0.0788) + \
1164 (abs(eta) <= 1.8) * (pt > 600.0 && pt <= 700) * (0.0581) + \
1165 (abs(eta) <= 1.8) * (pt > 700.0 && pt <= 800) * (0.0534) + \
1166 (abs(eta) <= 1.8) * (pt > 800.0 && pt <= 1000) * (0.0521) + \
1167 (abs(eta) <= 1.8) * (pt > 1000.0 && pt <= 1400) * (0.0626) + \
1168 (abs(eta) <= 1.8) * (pt > 1400.0 && pt <= 2000) * (0.0826) + \
1169 (abs(eta) <= 1.8) * (pt > 2000.0) * (0.1022) + \
1170 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt <= 20.0) * (0.000) + \
1171 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 20.0 && pt <= 30) * (0.1562) + \
1172 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 30.0 && pt <= 40) * (0.2499) + \
1173 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 40.0 && pt <= 50) * (0.2956) + \
1174 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 50.0 && pt <= 60) * (0.3128) + \
1175 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 60.0 && pt <= 70) * (0.3147) + \
1176 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 70.0 && pt <= 80) * (0.3222) + \
1177 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 80.0 && pt <= 90) * (0.304) + \
1178 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 90.0 && pt <= 100) * (0.3051) + \
1179 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 100.0 && pt <= 120) * (0.2657) + \
1180 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 120.0 && pt <= 140) * (0.2578) + \
1181 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 140.0 && pt <= 160) * (0.2087) + \
1182 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 160.0 && pt <= 180) * (0.1634) + \
1183 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 180.0 && pt <= 200) * (0.1651) + \
1184 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 200.0 && pt <= 250) * (0.1353) + \
1185 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 250.0 && pt <= 300) * (0.109) + \
1186 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 300.0 && pt <= 350) * (0.0799) + \
1187 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 350.0 && pt <= 400) * (0.0699) + \
1188 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 400.0 && pt <= 500) * (0.054) + \
1189 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 500.0 && pt <= 600) * (0.0718) + \
1190 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 600.0 && pt <= 700) * (0.0746) + \
1191 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 700.0 && pt <= 800) * (0.0803) + \
1192 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 800.0 && pt <= 1000) * (0.0882) + \
1193 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1000.0 && pt <= 1400) * (0.0) + \
1194 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 1400.0 && pt <= 2000) * (0.0) + \
1195 (abs(eta) > 1.8 && abs(eta) <= 2.4) * (pt > 2000.0) * (0.0) + \
1196 (abs(eta) > 2.4 && abs(eta) <= 4.0) * (pt > 15.0) * (0.584522*tanh(pt*0.0144387 + 0.397034)) + \
1197 (abs(eta) > 4.0) * (0.000)
1198 }
1199
1200}
1201
1202
1203
1204
1205
1206
1207
1208#####################
1209## tau tagging ##
1210#####################
1211
1212
1213module TauTagging TauTagging {
1214
1215 set ParticleInputArray Delphes/allParticles
1216 set PartonInputArray Delphes/partons
1217 set JetInputArray JetPileUpSubtractor/jets
1218
1219 set DeltaR 0.4
1220 set TauPTMin 1.0
1221 set TauEtaMax 4.0
1222
1223 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1224 # default efficiency formula (misidentification rate)
1225
1226 add EfficiencyFormula {0} {(abs(eta)<1.8)*0.006+(abs(eta)>1.8)*0.015}
1227
1228 # efficiency formula for tau-jets
1229
1230 add EfficiencyFormula {15} {0.65}
1231}
1232
1233
1234#####################
1235## MissingET ##
1236#####################
1237
1238### -sum of all particles after filtering neutrinos
1239
1240module Merger GenMissingET {
1241 add InputArray NeutrinoFilter/filteredParticles
1242 set MomentumOutputArray momentum
1243}
1244
1245module Merger RecoMissingET {
1246 add InputArray EFlowMerger/eflow
1247 set MomentumOutputArray momentum
1248}
1249
1250
1251
1252##############
1253## Scalar HT #
1254##############
1255
1256module Merger GenScalarHT {
1257 # add InputArray InputArray
1258 add InputArray NeutrinoFilter/filteredParticles
1259 set EnergyOutputArray energy
1260}
1261
1262module Merger RecoScalarHT {
1263 # add InputArray InputArray
1264 add InputArray EFlowMerger/eflow
1265 set EnergyOutputArray energy
1266}
1267
1268#####################
1269# Photon efficiency #
1270#####################
1271
1272module Efficiency PhotonEfficiency {
1273
1274 ## input particles
1275 set InputArray Calorimeter/eflowPhotons
1276 ## output particles
1277 set OutputArray photons
1278 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1279 # efficiency formula for photons
1280 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
1281 (abs(eta) <= 1.5) * (pt > 10.0) * (0.9635) + \
1282 (abs(eta) > 1.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.9624) + \
1283 (abs(eta) > 4.0) * (0.00)}
1284
1285}
1286
1287####################
1288# Photon isolation #
1289####################
1290
1291module Isolation PhotonIsolation {
1292
1293 # particle for which calculate the isolation
1294 set CandidateInputArray PhotonEfficiency/photons
1295
1296 # isolation collection
1297 set IsolationInputArray EFlowMergerAllTracks/eflow
1298
1299 # select a rho for the isolation
1300 set RhoInputArray Rho/rho
1301
1302 # output array
1303 set OutputArray photons
1304
1305 # isolation cone
1306 set DeltaRMax 0.3
1307
1308 # minimum pT
1309 set PTMin 1.0
1310
1311 # iso ratio to cut
1312 set PTRatioMax 9999.
1313}
1314
1315
1316#######################
1317# Electron efficiency #
1318#######################
1319
1320module Efficiency ElectronEfficiency {
1321
1322 set InputArray TrackPileUpSubtractor/electrons
1323 set OutputArray electrons
1324
1325 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
1326 # efficiency formula for electrons
1327 set EfficiencyFormula {
1328 (pt <= 4.0) * (0.00) + \
1329 (abs(eta) <= 1.45 ) * (pt > 4.0 && pt <= 6.0) * (0.50) + \
1330 (abs(eta) <= 1.45 ) * (pt > 6.0 && pt <= 8.0) * (0.70) + \
1331 (abs(eta) <= 1.45 ) * (pt > 8.0 && pt <= 10.0) * (0.85) + \
1332 (abs(eta) <= 1.45 ) * (pt > 10.0 && pt <= 30.0) * (0.94) + \
1333 (abs(eta) <= 1.45 ) * (pt > 30.0 && pt <= 50.0) * (0.97) + \
1334 (abs(eta) <= 1.45 ) * (pt > 50.0 && pt <= 70.0) * (0.98) + \
1335 (abs(eta) <= 1.45 ) * (pt > 70.0 ) * (1.0) + \
1336 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 4.0 && pt <= 10.0) * (0.35) + \
1337 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 10.0 && pt <= 30.0) * (0.40) + \
1338 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 30.0 && pt <= 70.0) * (0.45) + \
1339 (abs(eta) > 1.45 && abs(eta) <= 1.55) * (pt > 70.0 ) * (0.55) + \
1340 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 4.0 && pt <= 10.0) * (0.75) + \
1341 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 10.0 && pt <= 30.0) * (0.85) + \
1342 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 30.0 && pt <= 50.0) * (0.95) + \
1343 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 50.0 && pt <= 70.0) * (0.95) + \
1344 (abs(eta) >= 1.55 && abs(eta) <= 2.0 ) * (pt > 70.0 ) * (1.0) + \
1345 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1346 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1347 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1348 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1349 (abs(eta) >= 2.0 && abs(eta) <= 2.5 ) * (pt > 70.0 ) * (0.90) + \
1350 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 4.0 && pt <= 10.0) * (0.65) + \
1351 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 10.0 && pt <= 30.0) * (0.75) + \
1352 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 30.0 && pt <= 50.0) * (0.90) + \
1353 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 50.0 && pt <= 70.0) * (0.90) + \
1354 (abs(eta) > 2.5 && abs(eta) <= 4.0 ) * (pt > 70.0 ) * (0.90) + \
1355 (abs(eta) > 4.0) * (0.00)
1356
1357 }
1358}
1359
1360######################
1361# Electron isolation #
1362######################
1363
1364module Isolation ElectronIsolation {
1365
1366 set CandidateInputArray ElectronEfficiency/electrons
1367
1368 # isolation collection
1369 set IsolationInputArray EFlowMergerAllTracks/eflow
1370
1371 set RhoInputArray Rho/rho
1372
1373 set OutputArray electrons
1374
1375 set DeltaRMax 0.3
1376 set PTMin 1.0
1377 set PTRatioMax 9999.
1378
1379}
1380
1381
1382
1383###################
1384# Muon efficiency #
1385###################
1386
1387module Efficiency MuonEfficiency {
1388
1389 set InputArray TrackPileUpSubtractor/muons
1390
1391 set OutputArray muons
1392 # set EfficiencyFormula {efficiency as a function of eta and pt}
1393 # efficiency formula for muons
1394 set EfficiencyFormula { (pt <= 2.0) * (0.00) + \
1395 (abs(eta) <= 4.00) * (pt > 2.0 && pt <= 3.0) * (0.51) + \
1396 (abs(eta) <= 4.00) * (pt > 3.0 && pt <= 4.0) * (0.85) + \
1397 (abs(eta) <= 4.00) * (pt > 4.0 && pt <= 11.0) * (0.93) + \
1398 (abs(eta) <= 4.00) * (pt > 11. && pt <= 50.) * (0.96) + \
1399 (abs(eta) <= 4.00) * (pt > 50. && pt <= 70.) * (0.98) + \
1400 (abs(eta) <= 4.00) * (pt > 70.0 ) * (1.00) + \
1401 (abs(eta) > 4.00) * (0.00)
1402 }
1403
1404}
1405
1406
1407
1408##################
1409# Muon isolation #
1410##################
1411
1412module Isolation MuonIsolation {
1413 set CandidateInputArray MuonEfficiency/muons
1414
1415 # isolation collection
1416 set IsolationInputArray EFlowMergerAllTracks/eflow
1417
1418 set RhoInputArray Rho/rho
1419
1420 set OutputArray muons
1421
1422 set DeltaRMax 0.3
1423 set PTMin 1.0
1424 set PTRatioMax 9999.
1425
1426}
1427
1428
1429###############################################################################################################
1430# StatusPidFilter: this module removes all generated particles except electrons, muons, taus, and status == 3 #
1431###############################################################################################################
1432
1433module StatusPidFilter StatusPidFilter {
1434 ## take the particles from Pythia8 not adding pile-up
1435 set InputArray Delphes/allParticles
1436 set OutputArray filteredParticles
1437 set PTMin 5.0
1438
1439}
1440
1441##################
1442# ROOT tree writer
1443##################
1444
1445module TreeWriter TreeWriter {
1446
1447 ## branch notation : <particle collection> <branch name> <type of object in classes/DelphesClass.h<
1448
1449 # add Branch Delphes/stableParticles GenParticles GenParticle
1450
1451 ## gen particle stuff
1452 add Branch StatusPidFilter/filteredParticles GenParticles GenParticle
1453 add Branch PileUpMerger/vertices Vertex Vertex
1454
1455 ## Detector deposits
1456
1457 #add Branch TrackMerger/tracks Track Track
1458 #add Branch Calorimeter/towers Tower Tower
1459
1460 #add Branch Calorimeter/eflowTracks EFlowTrack Track
1461 #add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
1462 #add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
1463
1464 add Branch Rho/rho Rho Rho
1465
1466 ## Jets
1467 add Branch GenJetFinder/jets GenJet Jet
1468 add Branch JetPileUpSubtractor/jets Jet Jet
1469
1470 ## MET
1471 add Branch GenMissingET/momentum GenMissingET MissingET
1472 add Branch RecoMissingET/momentum MissingET MissingET
1473
1474 ## HT
1475 add Branch GenScalarHT/energy GenHT ScalarHT
1476 add Branch RecoScalarHT/energy RecoHT ScalarHT
1477
1478 ## photons and leptons
1479 add Branch ElectronIsolation/electrons Electron Electron
1480 add Branch PhotonIsolation/photons Photon Photon
1481 add Branch MuonIsolation/muons Muon Muon
1482
1483
1484
1485}
Note: See TracBrowser for help on using the repository browser.