Fork me on GitHub

source: git/cards/FCC/FCChh.tcl@ bdf9e2e

Last change on this file since bdf9e2e was eb332c8, checked in by Michele Selvaggi <michele.selvaggi@…>, 8 years ago

latest calo resolutions from M. Aleksa
  • Property mode set to 100644
File size: 21.7 KB
Line 
1#
2# Official Delphes card prepared by FCC-hh collaboration
3#
4# Main authors: Michele Selvaggi (CERN)
5#
6# Released on: Apr. 6th, 2017
7#
8# Configuration: FCC-hh baseline detector
9#
10#######################################
11# Order of execution of various modules
12#######################################
13
14set ExecutionPath {
15
16 ParticlePropagator
17
18 ChargedHadronTrackingEfficiency
19 ElectronTrackingEfficiency
20 MuonTrackingEfficiency
21
22 ChargedHadronMomentumSmearing
23 ElectronMomentumSmearing
24 MuonMomentumSmearing
25
26 TrackMerger
27
28 ECal
29 HCal
30
31 Calorimeter
32 EFlowMerger
33 EFlowFilter
34
35 PhotonEfficiency
36 PhotonIsolation
37
38 ElectronFilter
39 ElectronIsolation
40
41 ChargedHadronFilter
42
43 MuonIsolation
44
45 NeutrinoFilter
46
47 MissingET
48 GenMissingET
49
50 GenJetFinder
51 FastJetFinder
52 FatJetFinder
53
54 JetEnergyScale
55
56 JetFlavorAssociation
57
58 BTagging
59 CTagging
60 TauTagging
61
62 ScalarHT
63
64 UniqueObjectFinder
65
66 TreeWriter
67}
68
69#################################
70# Propagate particles in cylinder
71#################################
72
73module ParticlePropagator ParticlePropagator {
74 set InputArray Delphes/stableParticles
75
76 set OutputArray stableParticles
77 set ChargedHadronOutputArray chargedHadrons
78 set ElectronOutputArray electrons
79 set MuonOutputArray muons
80
81 # radius of the magnetic field coverage, in m
82 set Radius 1.5
83 # half-length of the magnetic field coverage, in m
84 set HalfLength 5
85
86 # magnetic field
87 set Bz 4.0
88}
89
90####################################
91# Charged hadron tracking efficiency
92####################################
93
94module Efficiency ChargedHadronTrackingEfficiency {
95 set InputArray ParticlePropagator/chargedHadrons
96 set OutputArray chargedHadrons
97
98 # TBC (which eta_max ? which pT min?)
99
100 # tracking efficiency formula for charged hadrons
101
102 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
103(abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
104(abs(eta) <= 2.5) * (pt > 1) * (0.95) + \
105(abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
106(abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
107(abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
108(abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
109(abs(eta) > 6.0) * (0.00)}
110
111}
112
113##############################
114# Electron tracking efficiency
115##############################
116
117module Efficiency ElectronTrackingEfficiency {
118 set InputArray ParticlePropagator/electrons
119 set OutputArray electrons
120
121# TBC (which eta_max ?)
122# putting same as charged hadrons for now...
123
124 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
125 (abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
126 (abs(eta) <= 2.5) * (pt > 1) * (0.95) + \
127 (abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
128 (abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
129 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
130 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
131 (abs(eta) > 6.0) * (0.00)}
132
133}
134##########################
135# Muon tracking efficiency
136##########################
137
138module Efficiency MuonTrackingEfficiency {
139 set InputArray ParticlePropagator/muons
140 set OutputArray muons
141
142# TBC (which eta_max ? why eff = 0 for 4 < eta < 6 ? for now put the same as central)
143# what about high pT ?
144 # tracking efficiency formula for muons
145 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
146 (abs(eta) <= 6.0) * (pt > 0.5 && pt <= 1) * (0.90) + \
147 (abs(eta) <= 6.0) * (pt > 1) * (0.99) + \
148 (abs(eta) > 6.0) * (0.00)}
149
150}
151
152########################################
153# Momentum resolution for charged tracks
154########################################
155
156module MomentumSmearing ChargedHadronMomentumSmearing {
157 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
158 set OutputArray chargedHadrons
159
160 source momentumResolutionVsP.tcl
161}
162
163
164###################################
165# Momentum resolution for electrons
166###################################
167
168module MomentumSmearing ElectronMomentumSmearing {
169 set InputArray ElectronTrackingEfficiency/electrons
170 set OutputArray electrons
171
172 source momentumResolutionVsP.tcl
173}
174
175
176###############################
177# Momentum resolution for muons
178###############################
179
180module MomentumSmearing MuonMomentumSmearing {
181 set InputArray MuonTrackingEfficiency/muons
182 set OutputArray muons
183
184 # TBC for just putting tracker resolution/ need to add improvement at high pT
185
186 source muonMomentumResolutionVsP.tcl
187}
188
189##############
190# Track merger
191##############
192
193module Merger TrackMerger {
194# add InputArray InputArray
195 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
196 add InputArray ElectronMomentumSmearing/electrons
197 add InputArray MuonMomentumSmearing/muons
198 set OutputArray tracks
199}
200
201#############
202# ECAL
203#############
204
205# TBC : calos seems ok, check eta max value though.
206
207module SimpleCalorimeter ECal {
208 set ParticleInputArray ParticlePropagator/stableParticles
209 set TrackInputArray TrackMerger/tracks
210
211 set TowerOutputArray ecalTowers
212 set EFlowTrackOutputArray eflowTracks
213 set EFlowTowerOutputArray eflowPhotons
214
215 set IsEcal true
216
217 set EnergyMin 0.5
218 set EnergySignificanceMin 2.0
219
220 set SmearTowerCenter true
221
222 set pi [expr {acos(-1)}]
223
224 # lists of the edges of each tower in eta and phi
225 # each list starts with the lower edge of the first tower
226 # the list ends with the higher edged of the last tower
227
228 # 0.012 rad towers up to eta = 2.5 (barrel)
229 set PhiBins {}
230 for {set i -256} {$i <= 256} {incr i} {
231 add PhiBins [expr {$i * $pi/256.0}]
232 }
233
234 # 0.01 unit in eta up to eta = 2.5 (barrel)
235 for {set i -249} {$i <= 250} {incr i} {
236 set eta [expr {$i * 0.01}]
237 add EtaPhiBins $eta $PhiBins
238 }
239
240 # 0.025 rad between 2.5 and 6.0
241 set PhiBins {}
242 for {set i -128} {$i <= 128} {incr i} {
243 add PhiBins [expr {$i * $pi/128.0}]
244 }
245
246 # 0.025 unit in eta between 2.5 and 6.0
247 for {set i 0} {$i <= 140} {incr i} {
248 set eta [expr { -6.00 + $i * 0.025}]
249 add EtaPhiBins $eta $PhiBins
250 }
251
252 for {set i 0} {$i <= 139} {incr i} {
253 set eta [expr { 2.525 + $i * 0.025}]
254 add EtaPhiBins $eta $PhiBins
255 }
256
257 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
258
259 add EnergyFraction {0} {0.0}
260 # energy fractions for e, gamma and pi0
261 add EnergyFraction {11} {1.0}
262 add EnergyFraction {22} {1.0}
263 add EnergyFraction {111} {1.0}
264 # energy fractions for muon, neutrinos and neutralinos
265 add EnergyFraction {12} {0.0}
266 add EnergyFraction {13} {0.0}
267 add EnergyFraction {14} {0.0}
268 add EnergyFraction {16} {0.0}
269 add EnergyFraction {1000022} {0.0}
270 add EnergyFraction {1000023} {0.0}
271 add EnergyFraction {1000025} {0.0}
272 add EnergyFraction {1000035} {0.0}
273 add EnergyFraction {1000045} {0.0}
274 # energy fractions for K0short and Lambda
275 # add EnergyFraction {310} {0.3}
276 # add EnergyFraction {3122} {0.3}
277
278 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
279 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.007^2 + energy*0.10^2) + \
280 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.035^2 + energy*0.30^2)}
281
282
283}
284
285#############
286# HCAL
287#############
288
289module SimpleCalorimeter HCal {
290 set ParticleInputArray ParticlePropagator/stableParticles
291 set TrackInputArray ECal/eflowTracks
292
293 set TowerOutputArray hcalTowers
294 set EFlowTrackOutputArray eflowTracks
295 set EFlowTowerOutputArray eflowNeutralHadrons
296
297 set IsEcal false
298
299 set EnergyMin 1.0
300 set EnergySignificanceMin 2.0
301
302 set SmearTowerCenter true
303
304 set pi [expr {acos(-1)}]
305
306 # lists of the edges of each tower in eta and phi
307 # each list starts with the lower edge of the first tower
308 # the list ends with the higher edged of the last tower
309
310 # 0.025 rad towers up to eta = 2.5 (barrel)
311 set PhiBins {}
312 for {set i -128} {$i <= 128} {incr i} {
313 add PhiBins [expr {$i * $pi/128.0}]
314 }
315
316 # 0.025 unit in eta up to eta = 2.5 (barrel)
317 for {set i -99} {$i <= 100} {incr i} {
318 set eta [expr {$i * 0.025}]
319 add EtaPhiBins $eta $PhiBins
320 }
321
322 # 0.05 rad between 2.5 and 6.0
323 set PhiBins {}
324 for {set i -64} {$i <= 64} {incr i} {
325 add PhiBins [expr {$i * $pi/64.0}]
326 }
327
328 # 0.05 unit in eta between 2.5 and 6.0
329 for {set i 0} {$i <= 70} {incr i} {
330 set eta [expr { -6.00 + $i * 0.05}]
331 add EtaPhiBins $eta $PhiBins
332 }
333
334 for {set i 0} {$i <= 69} {incr i} {
335 set eta [expr { 2.55 + $i * 0.05}]
336 add EtaPhiBins $eta $PhiBins
337 }
338
339
340 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
341 add EnergyFraction {0} {1.0}
342 # energy fractions for e, gamma and pi0
343 add EnergyFraction {11} {0.0}
344 add EnergyFraction {22} {0.0}
345 add EnergyFraction {111} {0.0}
346 # energy fractions for muon, neutrinos and neutralinos
347 add EnergyFraction {12} {0.0}
348 add EnergyFraction {13} {0.0}
349 add EnergyFraction {14} {0.0}
350 add EnergyFraction {16} {0.0}
351 add EnergyFraction {1000022} {0.0}
352 add EnergyFraction {1000023} {0.0}
353 add EnergyFraction {1000025} {0.0}
354 add EnergyFraction {1000035} {0.0}
355 add EnergyFraction {1000045} {0.0}
356 # energy fractions for K0short and Lambda
357 # add EnergyFraction {310} {0.7}
358 # add EnergyFraction {3122} {0.7}
359
360 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
361 set ResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.03^2 + energy*0.50^2) + \
362 (abs(eta) > 1.7 && abs(eta) <= 4.0) * sqrt(energy^2*0.03^2 + energy*0.60^2) + \
363 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.10^2 + energy*1.00^2)}
364}
365
366#################
367# Electron filter
368#################
369
370module PdgCodeFilter ElectronFilter {
371 set InputArray HCal/eflowTracks
372 set OutputArray electrons
373 set Invert true
374 add PdgCode {11}
375 add PdgCode {-11}
376}
377
378
379######################
380# ChargedHadronFilter
381######################
382
383module PdgCodeFilter ChargedHadronFilter {
384 set InputArray HCal/eflowTracks
385 set OutputArray chargedHadrons
386
387 add PdgCode {11}
388 add PdgCode {-11}
389 add PdgCode {13}
390 add PdgCode {-13}
391}
392
393
394###################################################
395# Tower Merger (in case not using e-flow algorithm)
396###################################################
397
398module Merger Calorimeter {
399# add InputArray InputArray
400 add InputArray ECal/ecalTowers
401 add InputArray HCal/hcalTowers
402 set OutputArray towers
403}
404
405####################
406# Energy flow merger
407####################
408
409module Merger EFlowMerger {
410# add InputArray InputArray
411 add InputArray HCal/eflowTracks
412 add InputArray ECal/eflowPhotons
413 add InputArray HCal/eflowNeutralHadrons
414 set OutputArray eflow
415}
416
417######################
418# EFlowFilter
419######################
420
421module PdgCodeFilter EFlowFilter {
422 set InputArray EFlowMerger/eflow
423 set OutputArray eflow
424
425 add PdgCode {11}
426 add PdgCode {-11}
427 add PdgCode {13}
428 add PdgCode {-13}
429}
430
431
432###################
433# Missing ET merger
434###################
435
436module Merger MissingET {
437# add InputArray InputArray
438 add InputArray EFlowMerger/eflow
439 set MomentumOutputArray momentum
440}
441
442
443
444##################
445# Scalar HT merger
446##################
447
448module Merger ScalarHT {
449# add InputArray InputArray
450 add InputArray EFlowMerger/eflow
451 set EnergyOutputArray energy
452}
453
454#################
455# Neutrino Filter
456#################
457
458module PdgCodeFilter NeutrinoFilter {
459
460 set InputArray Delphes/stableParticles
461 set OutputArray filteredParticles
462
463 set PTMin 0.0
464
465 add PdgCode {12}
466 add PdgCode {14}
467 add PdgCode {16}
468 add PdgCode {-12}
469 add PdgCode {-14}
470 add PdgCode {-16}
471
472}
473
474
475#####################
476# MC truth jet finder
477#####################
478
479# TBC: is jet radius fine?
480
481module FastJetFinder GenJetFinder {
482# set InputArray NeutrinoFilter/filteredParticles
483 set InputArray Delphes/stableParticles
484
485 set OutputArray jets
486
487 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
488 set JetAlgorithm 6
489 set ParameterR 0.4
490
491 set JetPTMin 5.0
492}
493
494#########################
495# Gen Missing ET merger
496########################
497
498module Merger GenMissingET {
499
500# add InputArray InputArray
501 add InputArray NeutrinoFilter/filteredParticles
502 set MomentumOutputArray momentum
503}
504
505
506
507############
508# Jet finder
509############
510
511# TBC need to include jet substructure variables
512# TBC is jet radius fine?
513
514module FastJetFinder FastJetFinder {
515# set InputArray Calorimeter/towers
516 set InputArray EFlowMerger/eflow
517
518 set OutputArray jets
519
520 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
521 # 7: anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness
522
523 set JetAlgorithm 6
524 set ParameterR 0.4
525
526 set JetPTMin 30.0
527}
528
529##################
530# Fat Jet finder
531##################
532
533module FastJetFinder FatJetFinder {
534 set InputArray EFlowMerger/eflow
535
536 set OutputArray jets
537
538 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
539 set JetAlgorithm 6
540 set ParameterR 0.8
541
542 set ComputeNsubjettiness 1
543 set Beta 1.0
544 set AxisMode 4
545
546 set ComputeTrimming 1
547 set RTrim 0.2
548 set PtFracTrim 0.05
549
550 set ComputePruning 1
551 set ZcutPrun 0.1
552 set RcutPrun 0.5
553 set RPrun 0.8
554
555 set ComputeSoftDrop 1
556 set BetaSoftDrop 0.0
557 set SymmetryCutSoftDrop 0.1
558 set R0SoftDrop 0.8
559
560 set JetPTMin 200.0
561}
562
563
564
565##################
566# Jet Energy Scale
567##################
568
569module EnergyScale JetEnergyScale {
570 set InputArray FastJetFinder/jets
571 set OutputArray jets
572
573 # scale formula for jets
574 set ScaleFormula {1.00}
575}
576
577
578########################
579# Jet Flavor Association
580########################
581
582module JetFlavorAssociation JetFlavorAssociation {
583
584 set PartonInputArray Delphes/partons
585 set ParticleInputArray Delphes/allParticles
586 set ParticleLHEFInputArray Delphes/allParticlesLHEF
587 set JetInputArray JetEnergyScale/jets
588
589 set DeltaR 0.5
590 set PartonPTMin 5.0
591 set PartonEtaMax 6.0
592
593}
594
595###################
596# Photon efficiency
597###################
598
599module Efficiency PhotonEfficiency {
600 set InputArray ECal/eflowPhotons
601 set OutputArray photons
602
603 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
604
605 set EfficiencyFormula {
606 (pt <= 1.0) * (0.00) + \
607 (abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0) * (0.70) +
608 (abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) +
609 (abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
610
611 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0) * (0.60) +
612 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 5.0 && pt < 10.0) * (0.80) +
613 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) +
614
615 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 1.0 && pt < 5.0) * (0.50) + \
616 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 5.0 && pt < 10.0) * (0.70) + \
617 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 10.0) * (0.80) + \
618 (abs(eta) > 6.0) * (0.00)}
619
620}
621
622##################
623# Photon isolation
624##################
625
626# TBC: check values for iso cuts
627
628module Isolation PhotonIsolation {
629 set CandidateInputArray PhotonEfficiency/photons
630 set IsolationInputArray EFlowFilter/eflow
631
632 set OutputArray photons
633
634 set DeltaRMax 0.3
635
636 set PTMin 0.5
637
638 set PTRatioMax 0.1
639}
640
641
642####################
643# Electron isolation
644####################
645
646# TBC: check values for iso cuts
647
648module Isolation ElectronIsolation {
649 set CandidateInputArray ElectronFilter/electrons
650 set IsolationInputArray EFlowFilter/eflow
651
652 set OutputArray electrons
653
654 set DeltaRMax 0.3
655
656 set PTMin 0.5
657
658 set PTRatioMax 0.1
659}
660
661
662################
663# Muon isolation
664################
665
666# TBC: check values for iso cuts
667
668module Isolation MuonIsolation {
669 set CandidateInputArray MuonMomentumSmearing/muons
670 set IsolationInputArray EFlowFilter/eflow
671
672 set OutputArray muons
673
674 set DeltaRMax 0.3
675
676 set PTMin 0.5
677
678 set PTRatioMax 0.2
679}
680
681
682###########
683# b-tagging
684###########
685
686module BTagging BTagging {
687 set JetInputArray JetEnergyScale/jets
688
689 set BitNumber 0
690
691 add EfficiencyFormula {0} {
692
693 (pt <= 10.0) * (0.00) +
694 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.001) + \
695 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.001)*(1.0 - pt/20000.) + \
696 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \
697 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.00075) + \
698 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.00075)*(1.0 - pt/20000.) + \
699 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \
700 (abs(eta) > 4.0) * (0.00)}
701
702 add EfficiencyFormula {4} {
703
704 (pt <= 10.0) * (0.00) +
705 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.04) + \
706 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.04)*(1.0 - pt/20000.) + \
707 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) + \
708 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.03) + \
709 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.03)*(1.0 - pt/20000.) + \
710 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) + \
711 (abs(eta) > 4.0) * (0.00)}
712
713 add EfficiencyFormula {5} {
714
715 (pt <= 10.0) * (0.00) +
716 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.85) +
717 (abs(eta) < 2.5) * (pt > 500.0 && pt < 20000.0) * (0.85)*(1.0 - pt/20000.) +
718 (abs(eta) < 2.5) * (pt > 20000.0) * (0.000) +
719 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.64) +
720 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 20000.0) * (0.64)*(1.0 - pt/20000.) +
721 (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 20000.0) * (0.000) +
722 (abs(eta) >= 4.0) * (0.00)}
723
724}
725
726###########
727# c-tagging
728###########
729
730module BTagging CTagging {
731 set JetInputArray JetEnergyScale/jets
732
733 set BitNumber 1
734
735 add EfficiencyFormula {0} {
736
737 (pt <= 10.0) * (0.00) +
738 (abs(eta) < 4.0) * (pt > 10.0) * (0.01) + \
739 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
740
741 add EfficiencyFormula {4} {
742
743 (pt <= 10.0) * (0.00) +
744 (abs(eta) < 4.0) * (pt > 10.0) * (0.10) + \
745 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
746
747 add EfficiencyFormula {5} {
748
749 (pt <= 10.0) * (0.00) +
750 (abs(eta) < 4.0) * (pt > 10.0) * (0.25) + \
751 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
752
753}
754
755
756#############
757# tau-tagging
758#############
759
760
761module TauTagging TauTagging {
762 set ParticleInputArray Delphes/allParticles
763 set PartonInputArray Delphes/partons
764 set JetInputArray JetEnergyScale/jets
765
766 set DeltaR 0.5
767
768 set TauPTMin 1.0
769
770 set TauEtaMax 4.0
771
772 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
773 add EfficiencyFormula {0} {
774
775 (pt <= 10.0) * (0.00) +
776 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.01) + \
777 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.01) *(8./9. - pt/30000.) + \
778 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
779 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.0075) + \
780 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.0075)*(8./9. - pt/30000.) + \
781 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
782 (abs(eta) > 4.0) * (0.00)}
783
784 add EfficiencyFormula {11} {
785
786 (pt <= 10.0) * (0.00) +
787 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.005) + \
788 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.005) *(8./9. - pt/30000.) + \
789 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
790 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.00375) + \
791 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.00375)*(8./9. - pt/30000.) + \
792 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
793 (abs(eta) > 4.0) * (0.00)}
794
795 add EfficiencyFormula {15} {
796
797 (pt <= 10.0) * (0.00) +
798 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.6) + \
799 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.6) *(8./9. - pt/30000.) + \
800 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
801 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.45) + \
802 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.45)*(8./9. - pt/30000.) + \
803 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
804 (abs(eta) > 4.0) * (0.00)}
805
806}
807
808#####################################################
809# Find uniquely identified photons/electrons/tau/jets
810#####################################################
811
812module UniqueObjectFinder UniqueObjectFinder {
813# earlier arrays take precedence over later ones
814# add InputArray InputArray OutputArray
815 add InputArray PhotonIsolation/photons photons
816 add InputArray ElectronIsolation/electrons electrons
817 add InputArray MuonIsolation/muons muons
818 add InputArray JetEnergyScale/jets jets
819
820}
821
822
823##################
824# ROOT tree writer
825##################
826
827module TreeWriter TreeWriter {
828# add Branch InputArray BranchName BranchClass
829 add Branch Delphes/allParticles Particle GenParticle
830
831 add Branch GenJetFinder/jets GenJet Jet
832 add Branch GenMissingET/momentum GenMissingET MissingET
833
834 add Branch TrackMerger/tracks Track Track
835 add Branch Calorimeter/towers Tower Tower
836
837 add Branch HCal/eflowTracks EFlowTrack Track
838 add Branch ECal/eflowPhotons EFlowPhoton Tower
839 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
840
841 add Branch UniqueObjectFinder/photons Photon Photon
842 add Branch UniqueObjectFinder/electrons Electron Electron
843 add Branch UniqueObjectFinder/muons Muon Muon
844 add Branch UniqueObjectFinder/jets Jet Jet
845
846 add Branch FatJetFinder/jets FatJet Jet
847
848 add Branch MissingET/momentum MissingET MissingET
849 add Branch ScalarHT/energy ScalarHT ScalarHT
850}
851
Note: See TracBrowser for help on using the repository browser.