Fork me on GitHub

source: git/cards/delphes_card_CMS.tcl@ 94f8f5f

ImprovedOutputFile Timing
Last change on this file since 94f8f5f was dda357d, checked in by GitHub <noreply@…>, 7 years ago

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