MA5SandBox: delphes_card_cms_sus_19_006.tcl

File delphes_card_cms_sus_19_006.tcl, 23.3 KB (added by Benjamin Fuks, 6 years ago)
Line 
1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6 ParticlePropagator
7
8 ChargedHadronTrackingEfficiency
9 ElectronTrackingEfficiency
10 MuonTrackingEfficiency
11
12 ChargedHadronMomentumSmearing
13 ElectronMomentumSmearing
14 MuonMomentumSmearing
15
16 ElectronTrackingEfficiencyD0
17 MuonTrackingEfficiencyD0
18
19 TrackMerger
20
21 ECal
22 HCal
23
24 Calorimeter
25 EFlowMerger
26 EFlowFilter
27
28 PhotonEfficiency
29
30 ElectronFilter
31 ElectronEfficiency
32
33 ChargedHadronFilter
34
35 MuonEfficiency
36
37 MissingET
38
39 NeutrinoFilter
40 GenJetFinder
41 GenMissingET
42
43 FastJetFinder
44 FatJetFinder
45
46 JetEnergyScale
47
48 JetFlavorAssociation
49
50 BTagging
51 TauTagging
52
53 TreeWriter
54}
55
56#################################
57# Propagate particles in cylinder
58#################################
59
60module ParticlePropagator ParticlePropagator {
61 set InputArray Delphes/stableParticles
62
63 set OutputArray stableParticles
64 set ChargedHadronOutputArray chargedHadrons
65 set ElectronOutputArray electrons
66 set MuonOutputArray muons
67
68 # radius of the magnetic field coverage, in m
69 set Radius 1.29
70 # half-length of the magnetic field coverage, in m
71 set HalfLength 3.00
72
73 # magnetic field
74 set Bz 3.8
75}
76
77####################################
78# Charged hadron tracking efficiency
79####################################
80
81module Efficiency ChargedHadronTrackingEfficiency {
82 set InputArray ParticlePropagator/chargedHadrons
83 set OutputArray chargedHadrons
84
85 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
86
87 # tracking efficiency formula for charged hadrons
88 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
89 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
90 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
91 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
92 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
93 (abs(eta) > 2.5) * (0.00)}
94}
95
96##############################
97# Electron tracking efficiency
98##############################
99
100module Efficiency ElectronTrackingEfficiency {
101 set InputArray ParticlePropagator/electrons
102 set OutputArray electrons
103
104 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
105
106 # tracking efficiency formula for electrons
107 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
108 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
109 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
110 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
111 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
112 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
113 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
114 (abs(eta) > 2.5) * (0.00)}
115}
116
117##########################
118# Muon tracking efficiency
119##########################
120
121module Efficiency MuonTrackingEfficiency {
122 set InputArray ParticlePropagator/muons
123 set OutputArray muons
124
125 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
126
127 # tracking efficiency formula for muons
128 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
129 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
130 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e3) * (0.99) +
131 (abs(eta) <= 1.5) * (pt > 1.0e3 ) * (0.99 * exp(0.5 - pt*5.0e-4)) +
132
133 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
134 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e3) * (0.98) +
135 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3) * (0.98 * exp(0.5 - pt*5.0e-4)) +
136 (abs(eta) > 2.5) * (0.00)}
137}
138
139########################################
140# Momentum resolution for charged tracks
141########################################
142
143module MomentumSmearing ChargedHadronMomentumSmearing {
144 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
145 set OutputArray chargedHadrons
146
147 # set ResolutionFormula {resolution formula as a function of eta and pt}
148
149 # resolution formula for charged hadrons
150 # based on arXiv:1405.6569
151 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
152 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
153 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
154}
155
156###################################
157# Momentum resolution for electrons
158###################################
159
160module MomentumSmearing ElectronMomentumSmearing {
161 set InputArray ElectronTrackingEfficiency/electrons
162 set OutputArray electrons
163
164 # set ResolutionFormula {resolution formula as a function of eta and energy}
165
166 # resolution formula for electrons
167 # based on arXiv:1405.6569
168 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
169 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
170 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
171}
172
173###############################
174# Momentum resolution for muons
175###############################
176
177module MomentumSmearing MuonMomentumSmearing {
178 set InputArray MuonTrackingEfficiency/muons
179 set OutputArray muons
180
181 # set ResolutionFormula {resolution formula as a function of eta and pt}
182
183 # resolution formula for muons
184 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
185 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
186 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
187}
188
189##############################
190# Electron tracking efficiency
191##############################
192
193module Efficiency ElectronTrackingEfficiencyD0 {
194 set InputArray ElectronMomentumSmearing/electrons
195 set OutputArray electrons
196
197 # standard-algo (d0 in mm)
198 # set EfficiencyFormula { (d0 <=600) * (7.48324e-16 * d0**6 - 1.51117e-12 * d0**5 + 1.19454e-09 * d0**4 - 4.69901e-07 * d0**3 + 9.79836E-05 * d0**2 - 0.0115465 * d0 + 1) +
199 # (d0 >=600) * (0.00) }
200
201 # displaced-algo for 13 TeV (d0 in mm)
202 set EfficiencyFormula { (d0 <=600) * (1.13477e-16 * d0**6 - 2.54327e-13 * d0**5 + 2.28586e-10 * d0**4 - 1.05512E-07 * d0**3 + 2.62108E-05 * d0**2 - 0.00375463 * d0 + 1) +
203 (d0 >=600) * (0.00) }
204}
205
206
207
208##############################
209# Muon tracking efficiency
210##############################
211
212module Efficiency MuonTrackingEfficiencyD0 {
213 set InputArray MuonMomentumSmearing/muons
214 set OutputArray muons
215
216 # standard-algo (d0 in mm)
217 # set EfficiencyFormula { (d0 <=600) * (7.48324e-16 * d0**6 - 1.51117e-12 * d0**5 + 1.19454e-09 * d0**4 - 4.69901e-07 * d0**3 + 9.79836E-05 * d0**2 - 0.0115465 * d0 + 1) +
218 # (d0 >=600) * (0.00) }
219
220 # displaced-algo for 13 TeV (d0 in mm)
221 set EfficiencyFormula { (d0 <=600) * (1.13477e-16 * d0**6 - 2.54327e-13 * d0**5 + 2.28586e-10 * d0**4 - 1.05512E-07 * d0**3 + 2.62108E-05 * d0**2 - 0.00375463 * d0 + 1) +
222 (d0 >=600) * (0.00) }
223
224}
225
226
227
228##############
229# Track merger
230##############
231
232module Merger TrackMerger {
233# add InputArray InputArray
234 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
235 add InputArray ElectronTrackingEfficiencyD0/electrons
236 add InputArray MuonTrackingEfficiencyD0/muons
237 set OutputArray tracks
238}
239
240
241
242#############
243# ECAL
244#############
245
246module SimpleCalorimeter ECal {
247 set ParticleInputArray ParticlePropagator/stableParticles
248 set TrackInputArray TrackMerger/tracks
249
250 set TowerOutputArray ecalTowers
251 set EFlowTrackOutputArray eflowTracks
252 set EFlowTowerOutputArray eflowPhotons
253
254 set IsEcal true
255
256 set EnergyMin 0.5
257 set EnergySignificanceMin 2.0
258
259 set SmearTowerCenter true
260
261 set pi [expr {acos(-1)}]
262
263 # lists of the edges of each tower in eta and phi
264 # each list starts with the lower edge of the first tower
265 # the list ends with the higher edged of the last tower
266
267 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
268
269 set PhiBins {}
270 for {set i -180} {$i <= 180} {incr i} {
271 add PhiBins [expr {$i * $pi/180.0}]
272 }
273
274 # 0.02 unit in eta up to eta = 1.5 (barrel)
275 for {set i -85} {$i <= 86} {incr i} {
276 set eta [expr {$i * 0.0174}]
277 add EtaPhiBins $eta $PhiBins
278 }
279
280 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
281
282 set PhiBins {}
283 for {set i -180} {$i <= 180} {incr i} {
284 add PhiBins [expr {$i * $pi/180.0}]
285 }
286
287 # 0.02 unit in eta up to eta = 3
288 for {set i 1} {$i <= 84} {incr i} {
289 set eta [expr { -2.958 + $i * 0.0174}]
290 add EtaPhiBins $eta $PhiBins
291 }
292
293 for {set i 1} {$i <= 84} {incr i} {
294 set eta [expr { 1.4964 + $i * 0.0174}]
295 add EtaPhiBins $eta $PhiBins
296 }
297
298 # take present CMS granularity for HF
299
300 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
301 set PhiBins {}
302 for {set i -18} {$i <= 18} {incr i} {
303 add PhiBins [expr {$i * $pi/18.0}]
304 }
305
306 foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
307 add EtaPhiBins $eta $PhiBins
308 }
309
310
311 add EnergyFraction {0} {0.0}
312 # energy fractions for e, gamma and pi0
313 add EnergyFraction {11} {1.0}
314 add EnergyFraction {22} {1.0}
315 add EnergyFraction {111} {1.0}
316 # energy fractions for muon, neutrinos and neutralinos
317 add EnergyFraction {12} {0.0}
318 add EnergyFraction {13} {0.0}
319 add EnergyFraction {14} {0.0}
320 add EnergyFraction {16} {0.0}
321 add EnergyFraction {1000022} {0.0}
322 add EnergyFraction {1000023} {0.0}
323 add EnergyFraction {1000025} {0.0}
324 add EnergyFraction {1000035} {0.0}
325 add EnergyFraction {1000045} {0.0}
326 # energy fractions for K0short and Lambda
327 add EnergyFraction {310} {0.3}
328 add EnergyFraction {3122} {0.3}
329
330 # set ResolutionFormula {resolution formula as a function of eta and energy}
331
332 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
333
334 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
335 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
336 set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
337 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (2.16 + 5.6*(abs(eta)-2)^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
338 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
339
340}
341
342
343#############
344# HCAL
345#############
346
347module SimpleCalorimeter HCal {
348 set ParticleInputArray ParticlePropagator/stableParticles
349 set TrackInputArray ECal/eflowTracks
350
351 set TowerOutputArray hcalTowers
352 set EFlowTrackOutputArray eflowTracks
353 set EFlowTowerOutputArray eflowNeutralHadrons
354
355 set IsEcal false
356
357 set EnergyMin 1.0
358 set EnergySignificanceMin 1.0
359
360 set SmearTowerCenter true
361
362 set pi [expr {acos(-1)}]
363
364 # lists of the edges of each tower in eta and phi
365 # each list starts with the lower edge of the first tower
366 # the list ends with the higher edged of the last tower
367
368 # 5 degrees towers
369 set PhiBins {}
370 for {set i -36} {$i <= 36} {incr i} {
371 add PhiBins [expr {$i * $pi/36.0}]
372 }
373 foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
374 add EtaPhiBins $eta $PhiBins
375 }
376
377 # 10 degrees towers
378 set PhiBins {}
379 for {set i -18} {$i <= 18} {incr i} {
380 add PhiBins [expr {$i * $pi/18.0}]
381 }
382 foreach eta {-4.35 -4.175 -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.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 4.175 4.35 4.525} {
383 add EtaPhiBins $eta $PhiBins
384 }
385
386 # 20 degrees towers
387 set PhiBins {}
388 for {set i -9} {$i <= 9} {incr i} {
389 add PhiBins [expr {$i * $pi/9.0}]
390 }
391 foreach eta {-5 -4.7 -4.525 4.7 5} {
392 add EtaPhiBins $eta $PhiBins
393 }
394
395 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
396 add EnergyFraction {0} {1.0}
397 # energy fractions for e, gamma and pi0
398 add EnergyFraction {11} {0.0}
399 add EnergyFraction {22} {0.0}
400 add EnergyFraction {111} {0.0}
401 # energy fractions for muon, neutrinos and neutralinos
402 add EnergyFraction {12} {0.0}
403 add EnergyFraction {13} {0.0}
404 add EnergyFraction {14} {0.0}
405 add EnergyFraction {16} {0.0}
406 add EnergyFraction {1000022} {0.0}
407 add EnergyFraction {1000023} {0.0}
408 add EnergyFraction {1000025} {0.0}
409 add EnergyFraction {1000035} {0.0}
410 add EnergyFraction {1000045} {0.0}
411 # energy fractions for K0short and Lambda
412 add EnergyFraction {310} {0.7}
413 add EnergyFraction {3122} {0.7}
414
415 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
416 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
417 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
418
419}
420
421
422#################
423# Electron filter
424#################
425
426module PdgCodeFilter ElectronFilter {
427 set InputArray HCal/eflowTracks
428 set OutputArray electrons
429 set Invert true
430 add PdgCode {11}
431 add PdgCode {-11}
432}
433
434######################
435# ChargedHadronFilter
436######################
437
438module PdgCodeFilter ChargedHadronFilter {
439 set InputArray HCal/eflowTracks
440 set OutputArray chargedHadrons
441
442 add PdgCode {11}
443 add PdgCode {-11}
444 add PdgCode {13}
445 add PdgCode {-13}
446}
447
448
449###################################################
450# Tower Merger (in case not using e-flow algorithm)
451###################################################
452
453module Merger Calorimeter {
454# add InputArray InputArray
455 add InputArray ECal/ecalTowers
456 add InputArray HCal/hcalTowers
457 set OutputArray towers
458}
459
460
461
462####################
463# Energy flow merger
464####################
465
466module Merger EFlowMerger {
467# add InputArray InputArray
468 add InputArray HCal/eflowTracks
469 add InputArray ECal/eflowPhotons
470 add InputArray HCal/eflowNeutralHadrons
471 set OutputArray eflow
472}
473
474######################
475# EFlowFilter
476######################
477
478module PdgCodeFilter EFlowFilter {
479 set InputArray EFlowMerger/eflow
480 set OutputArray eflow
481
482 add PdgCode {11}
483 add PdgCode {-11}
484 add PdgCode {13}
485 add PdgCode {-13}
486}
487
488
489###################
490# Photon efficiency
491###################
492
493module Efficiency PhotonEfficiency {
494 set InputArray ECal/eflowPhotons
495 set OutputArray photons
496
497 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
498
499 # efficiency formula for photons
500 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
501 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
502 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
503 (abs(eta) > 2.5) * (0.00)}
504}
505
506
507#####################
508# Electron efficiency
509#####################
510
511module Efficiency ElectronEfficiency {
512 set InputArray ElectronFilter/electrons
513 set OutputArray electrons
514
515 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
516
517 # efficiency formula for electrons
518 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
519 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
520 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
521 (abs(eta) > 2.5) * (0.00)}
522}
523
524#################
525# Muon efficiency
526#################
527
528module Efficiency MuonEfficiency {
529 set InputArray MuonTrackingEfficiencyD0/muons
530 set OutputArray muons
531
532 # set EfficiencyFormula {efficiency as a function of eta and pt}
533
534 # efficiency formula for muons
535 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
536 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
537 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0) * (0.95) +
538 (abs(eta) > 2.4) * (0.00)}
539}
540
541
542###################
543# Missing ET merger
544###################
545
546module Merger MissingET {
547# add InputArray InputArray
548 add InputArray EFlowMerger/eflow
549 set MomentumOutputArray momentum
550}
551
552
553
554#####################
555# Neutrino Filter
556#####################
557
558module PdgCodeFilter NeutrinoFilter {
559
560 set InputArray Delphes/stableParticles
561 set OutputArray filteredParticles
562
563 set PTMin 0.0
564
565 add PdgCode {12}
566 add PdgCode {14}
567 add PdgCode {16}
568 add PdgCode {-12}
569 add PdgCode {-14}
570 add PdgCode {-16}
571
572}
573
574
575#####################
576# MC truth jet finder
577#####################
578
579module FastJetFinder GenJetFinder {
580 set InputArray NeutrinoFilter/filteredParticles
581
582 set OutputArray jets
583
584 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
585 set JetAlgorithm 6
586 set ParameterR 0.4
587
588 set JetPTMin 20.0
589}
590
591#########################
592# Gen Missing ET merger
593########################
594
595module Merger GenMissingET {
596# add InputArray InputArray
597 add InputArray NeutrinoFilter/filteredParticles
598 set MomentumOutputArray momentum
599}
600
601
602
603############
604# Jet finder
605############
606
607module FastJetFinder FastJetFinder {
608# set InputArray Calorimeter/towers
609 set InputArray EFlowMerger/eflow
610
611 set OutputArray jets
612
613 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
614 set JetAlgorithm 6
615 set ParameterR 0.4
616
617 set JetPTMin 20.0
618}
619
620##################
621# Fat Jet finder
622##################
623
624module FastJetFinder FatJetFinder {
625 set InputArray EFlowMerger/eflow
626
627 set OutputArray jets
628
629 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
630 set JetAlgorithm 6
631 set ParameterR 0.8
632
633 set ComputeNsubjettiness 1
634 set Beta 1.0
635 set AxisMode 4
636
637 set ComputeTrimming 1
638 set RTrim 0.2
639 set PtFracTrim 0.05
640
641 set ComputePruning 1
642 set ZcutPrun 0.1
643 set RcutPrun 0.5
644 set RPrun 0.8
645
646 set ComputeSoftDrop 1
647 set BetaSoftDrop 0.0
648 set SymmetryCutSoftDrop 0.1
649 set R0SoftDrop 0.8
650
651 set JetPTMin 200.0
652}
653
654
655
656
657##################
658# Jet Energy Scale
659##################
660
661module EnergyScale JetEnergyScale {
662 set InputArray FastJetFinder/jets
663 set OutputArray jets
664
665 # scale formula for jets
666# set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.05 )}
667 set ScaleFormula {1.00}
668}
669
670########################
671# Jet Flavor Association
672########################
673
674module JetFlavorAssociation JetFlavorAssociation {
675
676 set PartonInputArray Delphes/partons
677 set ParticleInputArray Delphes/allParticles
678 set ParticleLHEFInputArray Delphes/allParticlesLHEF
679 set JetInputArray JetEnergyScale/jets
680
681 set DeltaR 0.4
682 set PartonPTMin 1.0
683 set PartonEtaMax 2.5
684
685}
686
687###########
688# b-tagging
689###########
690
691module BTagging BTagging {
692 set JetInputArray JetEnergyScale/jets
693
694 set BitNumber 0
695
696 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
697 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
698 # gluon's PDG code has the lowest priority
699
700 # based on arXiv:1211.4462
701
702 # default efficiency formula (misidentification rate)
703# add EfficiencyFormula {0} {0.01+0.000038*pt}
704 add EfficiencyFormula{0} {(((pt>20) && (pt<=225)) * (0.019 \
705 -0.00031*pt \
706 +3.39*pow(10,-6)*pow(pt,2)\
707 -1.47*pow(10,-8)*pow(pt,3)\
708 +2.92*pow(10,-11)*pow(pt,4)\
709 -2.12*pow(10,-14)*pow(pt,5))\
710 +\
711 ((pt>225) && (pt<=1000))*(0.0218\
712 +2.46*pow(10,-5)*pt\
713 +4.7*pow(10,-9)*pow(pt,2))\
714 +\
715 ((pt>1000)*(0.0511)))}
716 # efficiency formula for c-jets (misidentification rate)
717# add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
718 add EfficiencyFormula {4} {(((pt>20) && (pt<=200))*(0.136 \
719 -0.000639*pt\
720 +6.188*pow(10,-6)*pow(pt,2)\
721 -2.26*pow(10,-8)*pow(pt,3)\
722 +3.61*pow(10,-11)*pow(pt,4)\
723 +2.09*pow(10,-14)*pow(pt,5))\
724 +\
725 ((pt>200) && (pt<=1000))*(0.103\
726 +0.00014*pt\
727 -1.15*pow(10,-7)*pow(pt,2))\
728 +\
729 ((pt>1000)*(0.128)))}
730
731 # efficiency formula for b-jets
732# add EfficiencyFormula {5} {0.9*(0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt)))}
733 add EfficiencyFormula {5} {((((pt > 20) && (pt<=50))*(0.194\
734 +0.0211*pt\
735 -0.000348*pow(pt,2)\
736 +pow(10,-6)*2.761*pow(pt,3)\
737 -pow(10,-8)*1.044*pow(pt,4)\
738 +pow(10,-11)*1.499*pow(pt,5)))\
739 +\
740 (((pt > 50) && (pt<=250))*(0.577+\
741 0.003417*pt -\
742 pow(10,-5)*3.26*pow(pt,2)+\
743 pow(10,-7)*1.506*pow(pt,3)-\
744 pow(10,-10)*3.63*pow(pt,4)+\
745 pow(10,-13)*3.522*pow(pt,5)))\
746 +\
747 (((pt > 250) && (pt<=1300))*(0.768-\
748 0.00045*pt))\
749
750 +\
751 ((pt>1300)*(0.183))) }
752
753
754}
755
756#############
757# tau-tagging
758#############
759
760module TauTagging TauTagging {
761 set ParticleInputArray Delphes/allParticles
762 set PartonInputArray Delphes/partons
763 set JetInputArray JetEnergyScale/jets
764
765 set DeltaR 0.5
766
767 set TauPTMin 1.0
768
769 set TauEtaMax 2.5
770
771 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
772
773 # default efficiency formula (misidentification rate)
774 add EfficiencyFormula {0} {0.01}
775 # efficiency formula for tau-jets
776 add EfficiencyFormula {15} {0.6}
777}
778
779
780##################
781# ROOT tree writer
782##################
783
784# tracks, towers and eflow objects are not stored by default in the output.
785# if needed (for jet constituent or other studies), uncomment the relevant
786# "add Branch ..." lines.
787
788module TreeWriter TreeWriter {
789# add Branch InputArray BranchName BranchClass
790 add Branch Delphes/allParticles Particle GenParticle
791
792 add Branch TrackMerger/tracks Track Track
793 add Branch Calorimeter/towers Tower Tower
794
795 add Branch HCal/eflowTracks EFlowTrack Track
796 add Branch ECal/eflowPhotons EFlowPhoton Tower
797 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
798
799 add Branch GenJetFinder/jets GenJet Jet
800 add Branch GenMissingET/momentum GenMissingET MissingET
801
802 add Branch JetEnergyScale/jets JetMA5 Jet
803 add Branch ElectronEfficiency/electrons ElectronMA5 Electron
804 add Branch PhotonEfficiency/photons PhotonMA5 Photon
805 add Branch MuonEfficiency/muons MuonMA5 Muon
806
807 add Branch FatJetFinder/jets FatJet Jet
808
809 add Branch MissingET/momentum MissingET MissingET
810}