Fork me on GitHub

source: git/cards/FCC/FCChh_PileUpVtx.tcl@ 1641926

Timing
Last change on this file since 1641926 was aa0b10f, checked in by Michele Selvaggi <michele.selvaggi@…>, 5 years ago

first dummy card for test pu vertexing in FCChh

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