MA5SandBox: delphes_card_cms_exo_16_022.tcl

File delphes_card_cms_exo_16_022.tcl, 22.2 KB (added by Benjamin Fuks, 5 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 set EfficiencyFormula { (d0<=20) * (-5.06107e-7 * d0**6 + 0.0000272756 * d0**5 - 0.00049321 * d0**4 + 0.00287189 * d0**3 + 0.00522007 * d0**2 -0.0917957 * d0 + 0.924921) +
198 (d0>20) * (0.00) }
199
200}
201
202
203
204##############################
205# Muon tracking efficiency
206##############################
207
208module Efficiency MuonTrackingEfficiencyD0 {
209 set InputArray MuonMomentumSmearing/muons
210 set OutputArray muons
211
212 set EfficiencyFormula { (d0<=20) * (-1.11279e-7 * d0**6 + 4.23266e-6 * d0**5 - 0.0000260292 * d0**4 - 0.000611108 * d0**3 + 0.00743217 * d0**2 - 0.0271852 * d0 + 0.99067) +
213 (d0>20) * (0.00) }
214
215}
216
217
218
219##############
220# Track merger
221##############
222
223module Merger TrackMerger {
224# add InputArray InputArray
225 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
226 add InputArray ElectronTrackingEfficiencyD0/electrons
227 add InputArray MuonTrackingEfficiencyD0/muons
228 set OutputArray tracks
229}
230
231
232
233#############
234# ECAL
235#############
236
237module SimpleCalorimeter ECal {
238 set ParticleInputArray ParticlePropagator/stableParticles
239 set TrackInputArray TrackMerger/tracks
240
241 set TowerOutputArray ecalTowers
242 set EFlowTrackOutputArray eflowTracks
243 set EFlowTowerOutputArray eflowPhotons
244
245 set IsEcal true
246
247 set EnergyMin 0.5
248 set EnergySignificanceMin 2.0
249
250 set SmearTowerCenter true
251
252 set pi [expr {acos(-1)}]
253
254 # lists of the edges of each tower in eta and phi
255 # each list starts with the lower edge of the first tower
256 # the list ends with the higher edged of the last tower
257
258 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
259
260 set PhiBins {}
261 for {set i -180} {$i <= 180} {incr i} {
262 add PhiBins [expr {$i * $pi/180.0}]
263 }
264
265 # 0.02 unit in eta up to eta = 1.5 (barrel)
266 for {set i -85} {$i <= 86} {incr i} {
267 set eta [expr {$i * 0.0174}]
268 add EtaPhiBins $eta $PhiBins
269 }
270
271 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
272
273 set PhiBins {}
274 for {set i -180} {$i <= 180} {incr i} {
275 add PhiBins [expr {$i * $pi/180.0}]
276 }
277
278 # 0.02 unit in eta up to eta = 3
279 for {set i 1} {$i <= 84} {incr i} {
280 set eta [expr { -2.958 + $i * 0.0174}]
281 add EtaPhiBins $eta $PhiBins
282 }
283
284 for {set i 1} {$i <= 84} {incr i} {
285 set eta [expr { 1.4964 + $i * 0.0174}]
286 add EtaPhiBins $eta $PhiBins
287 }
288
289 # take present CMS granularity for HF
290
291 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
292 set PhiBins {}
293 for {set i -18} {$i <= 18} {incr i} {
294 add PhiBins [expr {$i * $pi/18.0}]
295 }
296
297 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} {
298 add EtaPhiBins $eta $PhiBins
299 }
300
301
302 add EnergyFraction {0} {0.0}
303 # energy fractions for e, gamma and pi0
304 add EnergyFraction {11} {1.0}
305 add EnergyFraction {22} {1.0}
306 add EnergyFraction {111} {1.0}
307 # energy fractions for muon, neutrinos and neutralinos
308 add EnergyFraction {12} {0.0}
309 add EnergyFraction {13} {0.0}
310 add EnergyFraction {14} {0.0}
311 add EnergyFraction {16} {0.0}
312 add EnergyFraction {1000022} {0.0}
313 add EnergyFraction {1000023} {0.0}
314 add EnergyFraction {1000025} {0.0}
315 add EnergyFraction {1000035} {0.0}
316 add EnergyFraction {1000045} {0.0}
317 # energy fractions for K0short and Lambda
318 add EnergyFraction {310} {0.3}
319 add EnergyFraction {3122} {0.3}
320
321 # set ResolutionFormula {resolution formula as a function of eta and energy}
322
323 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
324
325 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
326 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
327 set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
328 (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) +
329 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
330
331}
332
333
334#############
335# HCAL
336#############
337
338module SimpleCalorimeter HCal {
339 set ParticleInputArray ParticlePropagator/stableParticles
340 set TrackInputArray ECal/eflowTracks
341
342 set TowerOutputArray hcalTowers
343 set EFlowTrackOutputArray eflowTracks
344 set EFlowTowerOutputArray eflowNeutralHadrons
345
346 set IsEcal false
347
348 set EnergyMin 1.0
349 set EnergySignificanceMin 1.0
350
351 set SmearTowerCenter true
352
353 set pi [expr {acos(-1)}]
354
355 # lists of the edges of each tower in eta and phi
356 # each list starts with the lower edge of the first tower
357 # the list ends with the higher edged of the last tower
358
359 # 5 degrees towers
360 set PhiBins {}
361 for {set i -36} {$i <= 36} {incr i} {
362 add PhiBins [expr {$i * $pi/36.0}]
363 }
364 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} {
365 add EtaPhiBins $eta $PhiBins
366 }
367
368 # 10 degrees towers
369 set PhiBins {}
370 for {set i -18} {$i <= 18} {incr i} {
371 add PhiBins [expr {$i * $pi/18.0}]
372 }
373 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} {
374 add EtaPhiBins $eta $PhiBins
375 }
376
377 # 20 degrees towers
378 set PhiBins {}
379 for {set i -9} {$i <= 9} {incr i} {
380 add PhiBins [expr {$i * $pi/9.0}]
381 }
382 foreach eta {-5 -4.7 -4.525 4.7 5} {
383 add EtaPhiBins $eta $PhiBins
384 }
385
386 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
387 add EnergyFraction {0} {1.0}
388 # energy fractions for e, gamma and pi0
389 add EnergyFraction {11} {0.0}
390 add EnergyFraction {22} {0.0}
391 add EnergyFraction {111} {0.0}
392 # energy fractions for muon, neutrinos and neutralinos
393 add EnergyFraction {12} {0.0}
394 add EnergyFraction {13} {0.0}
395 add EnergyFraction {14} {0.0}
396 add EnergyFraction {16} {0.0}
397 add EnergyFraction {1000022} {0.0}
398 add EnergyFraction {1000023} {0.0}
399 add EnergyFraction {1000025} {0.0}
400 add EnergyFraction {1000035} {0.0}
401 add EnergyFraction {1000045} {0.0}
402 # energy fractions for K0short and Lambda
403 add EnergyFraction {310} {0.7}
404 add EnergyFraction {3122} {0.7}
405
406 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
407 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
408 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
409
410}
411
412
413#################
414# Electron filter
415#################
416
417module PdgCodeFilter ElectronFilter {
418 set InputArray HCal/eflowTracks
419 set OutputArray electrons
420 set Invert true
421 add PdgCode {11}
422 add PdgCode {-11}
423}
424
425######################
426# ChargedHadronFilter
427######################
428
429module PdgCodeFilter ChargedHadronFilter {
430 set InputArray HCal/eflowTracks
431 set OutputArray chargedHadrons
432
433 add PdgCode {11}
434 add PdgCode {-11}
435 add PdgCode {13}
436 add PdgCode {-13}
437}
438
439
440###################################################
441# Tower Merger (in case not using e-flow algorithm)
442###################################################
443
444module Merger Calorimeter {
445# add InputArray InputArray
446 add InputArray ECal/ecalTowers
447 add InputArray HCal/hcalTowers
448 set OutputArray towers
449}
450
451
452
453####################
454# Energy flow merger
455####################
456
457module Merger EFlowMerger {
458# add InputArray InputArray
459 add InputArray HCal/eflowTracks
460 add InputArray ECal/eflowPhotons
461 add InputArray HCal/eflowNeutralHadrons
462 set OutputArray eflow
463}
464
465######################
466# EFlowFilter
467######################
468
469module PdgCodeFilter EFlowFilter {
470 set InputArray EFlowMerger/eflow
471 set OutputArray eflow
472
473 add PdgCode {11}
474 add PdgCode {-11}
475 add PdgCode {13}
476 add PdgCode {-13}
477}
478
479
480###################
481# Photon efficiency
482###################
483
484module Efficiency PhotonEfficiency {
485 set InputArray ECal/eflowPhotons
486 set OutputArray photons
487
488 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
489
490 # efficiency formula for photons
491 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
492 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
493 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
494 (abs(eta) > 2.5) * (0.00)}
495}
496
497
498#####################
499# Electron efficiency
500#####################
501
502module Efficiency ElectronEfficiency {
503 set InputArray ElectronFilter/electrons
504 set OutputArray electrons
505
506 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
507
508 # efficiency formula for electrons
509 set EfficiencyFormula { (pt < 38.0) * (0.00) +
510 (pt >=38 && pt<=50) * (0.78) +
511 (pt >50 && pt<=70) * (0.84) +
512 (pt >70 && pt<=90) * (0.85) +
513 (pt >90 && pt<=120) * (0.86) +
514 (pt >120 && pt<=160) * (0.87) +
515 (pt >160 && pt<=200) * (0.88) +
516 (pt >200 && pt<=300) * (0.90) +
517 (pt >300 && pt<=500) * (0.92) +
518 (pt >500 && pt<=1000) * (0.94) +
519 (abs(eta) > 2.5) * (0.00)}
520}
521
522#################
523# Muon efficiency
524#################
525
526module Efficiency MuonEfficiency {
527 set InputArray MuonTrackingEfficiencyD0/muons
528 set OutputArray muons
529
530 # set EfficiencyFormula {efficiency as a function of eta and pt}
531
532 # efficiency formula for muons
533 set EfficiencyFormula { (pt < 38.0) * (0.00) +
534 (pt >=38 && pt<=50) * (0.79) +
535 (pt >50 && pt<=70) * (0.81) +
536 (pt >70 && pt<=90) * (0.84) +
537 (pt >90 && pt<=120) * (0.86) +
538 (pt >120 && pt<=160) * (0.87) +
539 (pt >160 && pt<=200) * (0.85) +
540 (pt >200 && pt<=300) * (0.86) +
541 (pt >300 && pt<=500) * (0.87) +
542 (pt >500 && pt<=1000) * (0.9) +
543 (abs(eta) > 2.4) * (0.00)}
544}
545
546
547###################
548# Missing ET merger
549###################
550
551module Merger MissingET {
552# add InputArray InputArray
553 add InputArray EFlowMerger/eflow
554 set MomentumOutputArray momentum
555}
556
557
558
559#####################
560# Neutrino Filter
561#####################
562
563module PdgCodeFilter NeutrinoFilter {
564
565 set InputArray Delphes/stableParticles
566 set OutputArray filteredParticles
567
568 set PTMin 0.0
569
570 add PdgCode {12}
571 add PdgCode {14}
572 add PdgCode {16}
573 add PdgCode {-12}
574 add PdgCode {-14}
575 add PdgCode {-16}
576
577}
578
579
580#####################
581# MC truth jet finder
582#####################
583
584module FastJetFinder GenJetFinder {
585 set InputArray NeutrinoFilter/filteredParticles
586
587 set OutputArray jets
588
589 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
590 set JetAlgorithm 6
591 set ParameterR 0.5
592
593 set JetPTMin 20.0
594}
595
596#########################
597# Gen Missing ET merger
598########################
599
600module Merger GenMissingET {
601# add InputArray InputArray
602 add InputArray NeutrinoFilter/filteredParticles
603 set MomentumOutputArray momentum
604}
605
606
607
608############
609# Jet finder
610############
611
612module FastJetFinder FastJetFinder {
613# set InputArray Calorimeter/towers
614 set InputArray EFlowMerger/eflow
615
616 set OutputArray jets
617
618 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
619 set JetAlgorithm 6
620 set ParameterR 0.5
621
622 set JetPTMin 20.0
623}
624
625##################
626# Fat Jet finder
627##################
628
629module FastJetFinder FatJetFinder {
630 set InputArray EFlowMerger/eflow
631
632 set OutputArray jets
633
634 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
635 set JetAlgorithm 6
636 set ParameterR 0.8
637
638 set ComputeNsubjettiness 1
639 set Beta 1.0
640 set AxisMode 4
641
642 set ComputeTrimming 1
643 set RTrim 0.2
644 set PtFracTrim 0.05
645
646 set ComputePruning 1
647 set ZcutPrun 0.1
648 set RcutPrun 0.5
649 set RPrun 0.8
650
651 set ComputeSoftDrop 1
652 set BetaSoftDrop 0.0
653 set SymmetryCutSoftDrop 0.1
654 set R0SoftDrop 0.8
655
656 set JetPTMin 200.0
657}
658
659
660
661
662##################
663# Jet Energy Scale
664##################
665
666module EnergyScale JetEnergyScale {
667 set InputArray FastJetFinder/jets
668 set OutputArray jets
669
670 # scale formula for jets
671 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
672}
673
674########################
675# Jet Flavor Association
676########################
677
678module JetFlavorAssociation JetFlavorAssociation {
679
680 set PartonInputArray Delphes/partons
681 set ParticleInputArray Delphes/allParticles
682 set ParticleLHEFInputArray Delphes/allParticlesLHEF
683 set JetInputArray JetEnergyScale/jets
684
685 set DeltaR 0.5
686 set PartonPTMin 1.0
687 set PartonEtaMax 2.5
688
689}
690
691###########
692# b-tagging
693###########
694
695module BTagging BTagging {
696 set JetInputArray JetEnergyScale/jets
697
698 set BitNumber 0
699
700 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
701 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
702 # gluon's PDG code has the lowest priority
703
704 # based on arXiv:1211.4462
705
706 # default efficiency formula (misidentification rate)
707 add EfficiencyFormula {0} {0.01+0.000038*pt}
708
709 # efficiency formula for c-jets (misidentification rate)
710 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
711
712 # efficiency formula for b-jets
713 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
714}
715
716#############
717# tau-tagging
718#############
719
720module TauTagging TauTagging {
721 set ParticleInputArray Delphes/allParticles
722 set PartonInputArray Delphes/partons
723 set JetInputArray JetEnergyScale/jets
724
725 set DeltaR 0.5
726
727 set TauPTMin 1.0
728
729 set TauEtaMax 2.5
730
731 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
732
733 # default efficiency formula (misidentification rate)
734 add EfficiencyFormula {0} {0.01}
735 # efficiency formula for tau-jets
736 add EfficiencyFormula {15} {0.6}
737}
738
739
740##################
741# ROOT tree writer
742##################
743
744# tracks, towers and eflow objects are not stored by default in the output.
745# if needed (for jet constituent or other studies), uncomment the relevant
746# "add Branch ..." lines.
747
748module TreeWriter TreeWriter {
749# add Branch InputArray BranchName BranchClass
750 add Branch Delphes/allParticles Particle GenParticle
751
752 add Branch TrackMerger/tracks Track Track
753 add Branch Calorimeter/towers Tower Tower
754
755 add Branch HCal/eflowTracks EFlowTrack Track
756 add Branch ECal/eflowPhotons EFlowPhoton Tower
757 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
758
759 add Branch GenJetFinder/jets GenJet Jet
760 add Branch GenMissingET/momentum GenMissingET MissingET
761
762 add Branch JetEnergyScale/jets JetMA5 Jet
763 add Branch ElectronEfficiency/electrons ElectronMA5 Electron
764 add Branch PhotonEfficiency/photons PhotonMA5 Photon
765 add Branch MuonEfficiency/muons MuonMA5 Muon
766
767 add Branch FatJetFinder/jets FatJet Jet
768
769 add Branch MissingET/momentum MissingET MissingET
770}