Fork me on GitHub

source: git/cards/delphes_card_CMS.tcl

Last change on this file was 46b3e01, checked in by youngkwon jo <cccpy98@…>, 3 years ago

add code for keeping the BHadrons

  • Property mode set to 100644
File size: 22.1 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 InputArray Delphes/allParticles
589 set OutputArray filteredParticles
590
591 set PTMin 0.0
592 set RequireStatus true
593 set RequireKeepGhostBHadron true
594
595 add PdgCode {12}
596 add PdgCode {14}
597 add PdgCode {16}
598 add PdgCode {-12}
599 add PdgCode {-14}
600 add PdgCode {-16}
601
602}
603
604
605#####################
606# MC truth jet finder
607#####################
608
609module FastJetFinder GenJetFinder {
610 set InputArray NeutrinoFilter/filteredParticles
611
612 set OutputArray jets
613
614 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
615 set JetAlgorithm 6
616 set ParameterR 0.5
617
618 set JetPTMin 20.0
619}
620
621#########################
622# Gen Missing ET merger
623########################
624
625module Merger GenMissingET {
626# add InputArray InputArray
627 add InputArray NeutrinoFilter/filteredParticles
628 set MomentumOutputArray momentum
629}
630
631
632
633############
634# Jet finder
635############
636
637module FastJetFinder FastJetFinder {
638# set InputArray Calorimeter/towers
639 set InputArray EFlowMerger/eflow
640
641 set OutputArray jets
642
643 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
644 set JetAlgorithm 6
645 set ParameterR 0.5
646
647 set JetPTMin 20.0
648}
649
650##################
651# Fat Jet finder
652##################
653
654module FastJetFinder FatJetFinder {
655 set InputArray EFlowMerger/eflow
656
657 set OutputArray jets
658
659 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
660 set JetAlgorithm 6
661 set ParameterR 0.8
662
663 set ComputeNsubjettiness 1
664 set Beta 1.0
665 set AxisMode 4
666
667 set ComputeTrimming 1
668 set RTrim 0.2
669 set PtFracTrim 0.05
670
671 set ComputePruning 1
672 set ZcutPrun 0.1
673 set RcutPrun 0.5
674 set RPrun 0.8
675
676 set ComputeSoftDrop 1
677 set BetaSoftDrop 0.0
678 set SymmetryCutSoftDrop 0.1
679 set R0SoftDrop 0.8
680
681 set JetPTMin 200.0
682}
683
684
685
686
687##################
688# Jet Energy Scale
689##################
690
691module EnergyScale JetEnergyScale {
692 set InputArray FastJetFinder/jets
693 set OutputArray jets
694
695 # scale formula for jets
696 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
697}
698
699########################
700# Jet Flavor Association
701########################
702
703module JetFlavorAssociation JetFlavorAssociation {
704
705 set PartonInputArray Delphes/partons
706 set ParticleInputArray Delphes/allParticles
707 set ParticleLHEFInputArray Delphes/allParticlesLHEF
708 set JetInputArray JetEnergyScale/jets
709
710 set DeltaR 0.5
711 set PartonPTMin 1.0
712 set PartonEtaMax 2.5
713
714}
715
716###########
717# b-tagging
718###########
719
720module BTagging BTagging {
721 set JetInputArray JetEnergyScale/jets
722
723 set BitNumber 0
724
725 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
726 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
727 # gluon's PDG code has the lowest priority
728
729 # based on arXiv:1211.4462
730
731 # default efficiency formula (misidentification rate)
732 add EfficiencyFormula {0} {0.01+0.000038*pt}
733
734 # efficiency formula for c-jets (misidentification rate)
735 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
736
737 # efficiency formula for b-jets
738 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
739}
740
741#############
742# tau-tagging
743#############
744
745module TauTagging TauTagging {
746 set ParticleInputArray Delphes/allParticles
747 set PartonInputArray Delphes/partons
748 set JetInputArray JetEnergyScale/jets
749
750 set DeltaR 0.5
751
752 set TauPTMin 1.0
753
754 set TauEtaMax 2.5
755
756 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
757
758 # default efficiency formula (misidentification rate)
759 add EfficiencyFormula {0} {0.01}
760 # efficiency formula for tau-jets
761 add EfficiencyFormula {15} {0.6}
762}
763
764#####################################################
765# Find uniquely identified photons/electrons/tau/jets
766#####################################################
767
768module UniqueObjectFinder UniqueObjectFinder {
769# earlier arrays take precedence over later ones
770# add InputArray InputArray OutputArray
771 add InputArray PhotonIsolation/photons photons
772 add InputArray ElectronIsolation/electrons electrons
773 add InputArray MuonIsolation/muons muons
774 add InputArray JetEnergyScale/jets jets
775}
776
777##################
778# ROOT tree writer
779##################
780
781# tracks, towers and eflow objects are not stored by default in the output.
782# if needed (for jet constituent or other studies), uncomment the relevant
783# "add Branch ..." lines.
784
785module TreeWriter TreeWriter {
786# add Branch InputArray BranchName BranchClass
787 add Branch Delphes/allParticles Particle GenParticle
788
789 add Branch TrackMerger/tracks Track Track
790 add Branch Calorimeter/towers Tower Tower
791
792 add Branch HCal/eflowTracks EFlowTrack Track
793 add Branch ECal/eflowPhotons EFlowPhoton Tower
794 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
795
796 add Branch GenJetFinder/jets GenJet Jet
797 add Branch GenMissingET/momentum GenMissingET MissingET
798
799 add Branch UniqueObjectFinder/jets Jet Jet
800 add Branch UniqueObjectFinder/electrons Electron Electron
801 add Branch UniqueObjectFinder/photons Photon Photon
802 add Branch UniqueObjectFinder/muons Muon Muon
803
804 add Branch FatJetFinder/jets FatJet Jet
805
806 add Branch MissingET/momentum MissingET MissingET
807 add Branch ScalarHT/energy ScalarHT ScalarHT
808}
Note: See TracBrowser for help on using the repository browser.