PublicAnalysisDatabase: delphes_card_cms_exo_16_022.tcl

File delphes_card_cms_exo_16_022.tcl, 22.1 KB (added by BenjF, 7 months 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 MA5EfficiencyD0 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 MA5EfficiencyD0 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 && <=50)  * (0.78) +
511                                            (pt >50 && <=70)  * (0.84) +
512                                            (pt >70 && <=90)  * (0.85) +
513                                            (pt >90 && <=120)  * (0.86) +
514                                            (pt >120 && <=160)  * (0.87) +
515                                            (pt >160 && <=200)  * (0.88) +
516                                            (pt >200 && <=300)  * (0.90) +
517                                            (pt >300 && <=500)  * (0.92) +
518                                            (pt >500 && <=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 && <=50)  * (0.79) +
535                                            (pt >50 && <=70)  * (0.81) +
536                                            (pt >70 && <=90)  * (0.84) +
537                                            (pt >90 && <=120)  * (0.86) +
538                                            (pt >120 && <=160)  * (0.87) +
539                                            (pt >160 && <=200)  * (0.85) +
540                                            (pt >200 && <=300)  * (0.86) +
541                                            (pt >300 && <=500)  * (0.87) +
542                                            (pt >500 && <=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}