Fork me on GitHub

source: git/cards/FCC/FCChh_PileUpVtx.tcl@ 2c81caa

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

adapted pythia reader to rhadrons

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