Fork me on GitHub

source: git/cards/FCC/FCChh_PileUpVtx_new.tcl@ 7939c6c

Timing
Last change on this file since 7939c6c was 7939c6c, checked in by Kaan Yüksel Oyulmaz <kaanyukseloyulmaz@…>, 5 years ago

TimeSmearing Module is updated for neutral particles

  • Property mode set to 100644
File size: 25.3 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 MaxEvents 100
15set RandomSeed 123
16
17set ExecutionPath {
18
19 PileUpMerger
20 ParticlePropagator
21
22 ChargedHadronTrackingEfficiency
23 ElectronTrackingEfficiency
24 MuonTrackingEfficiency
25
26 ChargedHadronMomentumSmearing
27 ElectronMomentumSmearing
28 MuonMomentumSmearing
29
30 TrackMerger
31
32 TrackSmearing
33 TimeSmearing
34
35 VertexFinderDA4D
36
37 TrackTimingPileUpSubtractor
38
39 ECal
40 HCal
41
42 Calorimeter
43 EFlowMerger
44 EFlowFilter
45
46 PhotonEfficiency
47 PhotonIsolation
48
49 ElectronFilter
50 ElectronIsolation
51
52 ChargedHadronFilter
53
54 MuonIsolation
55
56 NeutrinoFilter
57
58 MissingET
59 GenMissingET
60
61 GenJetFinder04
62 FastJetFinder04
63 CaloJetFinder04
64 TrackJetFinder04
65
66 JetEnergyScale
67
68 JetFlavorAssociation
69
70 BTagging
71 CTagging
72 TauTagging
73
74 ScalarHT
75
76 UniqueObjectFinder
77
78 TreeWriter
79}
80
81###############
82# PileUp Merger
83###############
84
85module PileUpMerger PileUpMerger {
86 set InputArray Delphes/stableParticles
87
88 set ParticleOutputArray stableParticles
89 set VertexOutputArray vertices
90
91 # pre-generated minbias input file
92 set PileUpFile MinBias.pileup
93
94 # average expected pile up
95 set MeanPileUp 10
96
97 # 0-poisson, 1-uniform, 2-delta
98 set PileUpDistribution 2
99
100 # maximum spread in the beam direction in m
101 set ZVertexSpread 0.25
102
103 # maximum spread in time in s
104 set TVertexSpread 800E-12
105
106 # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
107
108 #set VertexDistributionFormula {exp(-(t^2/(2*(0.063/2.99792458E8*exp(-(z^2/(2*(0.063)^2))))^2)))}
109 set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
110
111 # taking 5.3 cm x 160 ps
112
113 #set VertexDistributionFormula { (abs(t) <= 160e-12) * (abs(z) <= 0.053) * (1.00) +
114 # (abs(t) > 160e-12) * (abs(z) <= 0.053) * (0.00) +
115 # (abs(t) <= 160e-12) * (abs(z) > 0.053) * (0.00) +
116 # (abs(t) > 160e-12) * (abs(z) > 0.053) * (0.00)}
117
118}
119
120#################################
121# Propagate particles in cylinder
122#################################
123
124module ParticlePropagator ParticlePropagator {
125 set InputArray PileUpMerger/stableParticles
126 set OutputArray stableParticles
127 set ChargedHadronOutputArray chargedHadrons
128 set ElectronOutputArray electrons
129 set MuonOutputArray muons
130
131 # radius of the magnetic field coverage, in m
132 set Radius 1.5
133 # half-length of the magnetic field coverage, in m
134 set HalfLength 5
135
136 # magnetic field
137 set Bz 4.0
138}
139
140
141
142####################################
143# Charged hadron tracking efficiency
144####################################
145
146module Efficiency ChargedHadronTrackingEfficiency {
147 set InputArray ParticlePropagator/chargedHadrons
148 set OutputArray chargedHadrons
149
150 # TBC (which eta_max ? which pT min?)
151
152 # tracking efficiency formula for charged hadrons
153
154 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
155(abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
156(abs(eta) <= 2.5) * (pt > 1) * (0.95) + \
157(abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
158(abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
159(abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
160(abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
161(abs(eta) > 6.0) * (0.00)}
162
163}
164
165##############################
166# Electron tracking efficiency
167##############################
168
169module Efficiency ElectronTrackingEfficiency {
170 set InputArray ParticlePropagator/electrons
171 set OutputArray electrons
172
173# TBC (which eta_max ?)
174# putting same as charged hadrons for now...
175
176 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
177 (abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
178 (abs(eta) <= 2.5) * (pt > 1) * (0.95) + \
179 (abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
180 (abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
181 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
182 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
183 (abs(eta) > 6.0) * (0.00)}
184
185}
186
187##########################
188# Muon tracking efficiency
189##########################
190
191module Efficiency MuonTrackingEfficiency {
192 set InputArray ParticlePropagator/muons
193 set OutputArray muons
194
195# TBC (which eta_max ? why eff = 0 for 4 < eta < 6 ? for now put the same as central)
196# what about high pT ?
197 # tracking efficiency formula for muons
198 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
199 (abs(eta) <= 6.0) * (pt > 0.5 && pt <= 1) * (0.90) + \
200 (abs(eta) <= 6.0) * (pt > 1) * (0.99) + \
201 (abs(eta) > 6.0) * (0.00)}
202
203}
204
205########################################
206# Momentum resolution for charged tracks
207########################################
208
209module MomentumSmearing ChargedHadronMomentumSmearing {
210 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
211 set OutputArray chargedHadrons
212
213 source momentumResolutionVsP.tcl
214}
215
216
217###################################
218# Momentum resolution for electrons
219###################################
220
221module MomentumSmearing ElectronMomentumSmearing {
222 set InputArray ElectronTrackingEfficiency/electrons
223 set OutputArray electrons
224
225 source momentumResolutionVsP.tcl
226}
227
228
229###############################
230# Momentum resolution for muons
231###############################
232
233module MomentumSmearing MuonMomentumSmearing {
234 set InputArray MuonTrackingEfficiency/muons
235 set OutputArray muons
236
237 # TBC for just putting tracker resolution/ need to add improvement at high pT
238
239 source muonMomentumResolutionVsP.tcl
240}
241
242##############
243# Track merger
244##############
245
246module Merger TrackMerger {
247# add InputArray InputArray
248 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
249 add InputArray ElectronMomentumSmearing/electrons
250 add InputArray MuonMomentumSmearing/muons
251 set OutputArray tracks
252}
253
254
255########################################
256# Smear tracks
257########################################
258
259module TrackSmearing TrackSmearing {
260 set InputArray TrackMerger/tracks
261 set OutputArray tracks
262 set ApplyToPileUp true
263
264 # from http://mersi.web.cern.ch/mersi/layouts/.private/Baseline_tilted_200_Pixel_1_1_1/index.html
265 source trackResolutionCMS.tcl
266 # FIXME !!!! we need to add track resolution of FCC-hh baseline detector !!!!!
267}
268
269########################################
270# Time Smearing
271########################################
272
273module TimeSmearing TimeSmearing {
274 set InputArray TrackSmearing/tracks
275 set OutputArray tracks
276
277 # assume 20 ps resolution for now
278 set TimeResolution {20E-12}
279}
280
281##################################
282# Primary vertex reconstruction
283##################################
284
285
286module VertexFinderDA4D VertexFinderDA4D {
287 set InputArray TimeSmearing/tracks
288
289 set OutputArray tracks
290 set VertexOutputArray vertices
291
292 set Verbose 0
293 set MinPT 1.0
294
295 # in mm
296 set VertexSpaceSize 0.5
297
298 # in s
299 set VertexTimeSize 10E-12
300
301 set UseTc 1
302 set BetaMax 0.1
303 set BetaStop 1.0
304 set CoolingFactor 0.8
305 set MaxIterations 100
306
307 # in mm
308 set DzCutOff 40
309 set D0CutOff 30
310
311}
312
313##########################
314# Track pile-up subtractor
315##########################
316
317module TrackTimingPileUpSubtractor TrackTimingPileUpSubtractor {
318# add InputArray InputArray OutputArray
319
320 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
321 add InputArray ElectronMomentumSmearing/electrons
322 add InputArray MuonMomentumSmearing/muons
323
324 set VertexInputArray VertexFinderDA4D/vertices
325 # assume perfect pile-up subtraction for tracks with |z| > fZVertexResolution
326 # Z vertex resolution in m
327 set ZVertexResolution {0.0001}
328}
329
330
331#############
332# ECAL
333#############
334
335# TBC : calos seems ok, check eta max value though.
336
337module SimpleCalorimeter ECal {
338 set ParticleInputArray ParticlePropagator/stableParticles
339 set TrackInputArray TimeSmearing/tracks
340
341 set TowerOutputArray ecalTowers
342 set EFlowTrackOutputArray eflowTracks
343 set EFlowTowerOutputArray eflowPhotons
344
345 set IsEcal true
346
347 set EnergyMin 0.5
348 set EnergySignificanceMin 2.0
349
350 set SmearTowerCenter true
351
352 set pi [expr {acos(-1)}]
353
354 # lists of the edges of each tower in eta and phi
355 # each list starts with the lower edge of the first tower
356 # the list ends with the higher edged of the last tower
357
358 # 0.012 rad towers up to eta = 2.5 (barrel)
359 set PhiBins {}
360 for {set i -256} {$i <= 256} {incr i} {
361 add PhiBins [expr {$i * $pi/256.0}]
362 }
363
364 # 0.01 unit in eta up to eta = 2.5 (barrel)
365 for {set i -249} {$i <= 250} {incr i} {
366 set eta [expr {$i * 0.01}]
367 add EtaPhiBins $eta $PhiBins
368 }
369
370 # 0.025 rad between 2.5 and 6.0
371 set PhiBins {}
372 for {set i -128} {$i <= 128} {incr i} {
373 add PhiBins [expr {$i * $pi/128.0}]
374 }
375
376 # 0.025 unit in eta between 2.5 and 6.0
377 for {set i 0} {$i <= 140} {incr i} {
378 set eta [expr { -6.00 + $i * 0.025}]
379 add EtaPhiBins $eta $PhiBins
380 }
381
382 for {set i 0} {$i <= 139} {incr i} {
383 set eta [expr { 2.525 + $i * 0.025}]
384 add EtaPhiBins $eta $PhiBins
385 }
386
387 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
388
389 add EnergyFraction {0} {0.0}
390 # energy fractions for e, gamma and pi0
391 add EnergyFraction {11} {1.0}
392 add EnergyFraction {22} {1.0}
393 add EnergyFraction {111} {1.0}
394 # energy fractions for muon, neutrinos and neutralinos
395 add EnergyFraction {12} {0.0}
396 add EnergyFraction {13} {0.0}
397 add EnergyFraction {14} {0.0}
398 add EnergyFraction {16} {0.0}
399 add EnergyFraction {1000022} {0.0}
400 add EnergyFraction {1000023} {0.0}
401 add EnergyFraction {1000025} {0.0}
402 add EnergyFraction {1000035} {0.0}
403 add EnergyFraction {1000045} {0.0}
404 # energy fractions for K0short and Lambda
405 # add EnergyFraction {310} {0.3}
406 # add EnergyFraction {3122} {0.3}
407
408 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
409 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.007^2 + energy*0.10^2) + \
410 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.035^2 + energy*0.30^2)}
411
412
413}
414
415#############
416# HCAL
417#############
418
419module SimpleCalorimeter HCal {
420 set ParticleInputArray ParticlePropagator/stableParticles
421 set TrackInputArray ECal/eflowTracks
422
423 set TowerOutputArray hcalTowers
424 set EFlowTrackOutputArray eflowTracks
425 set EFlowTowerOutputArray eflowNeutralHadrons
426
427 set IsEcal false
428
429 set EnergyMin 1.0
430 set EnergySignificanceMin 2.0
431
432 set SmearTowerCenter true
433
434 set pi [expr {acos(-1)}]
435
436 # lists of the edges of each tower in eta and phi
437 # each list starts with the lower edge of the first tower
438 # the list ends with the higher edged of the last tower
439
440 # 0.025 rad towers up to eta = 2.5 (barrel)
441 set PhiBins {}
442 for {set i -128} {$i <= 128} {incr i} {
443 add PhiBins [expr {$i * $pi/128.0}]
444 }
445
446 # 0.025 unit in eta up to eta = 2.5 (barrel)
447 for {set i -99} {$i <= 100} {incr i} {
448 set eta [expr {$i * 0.025}]
449 add EtaPhiBins $eta $PhiBins
450 }
451
452 # 0.05 rad between 2.5 and 6.0
453 set PhiBins {}
454 for {set i -64} {$i <= 64} {incr i} {
455 add PhiBins [expr {$i * $pi/64.0}]
456 }
457
458 # 0.05 unit in eta between 2.5 and 6.0
459 for {set i 0} {$i <= 70} {incr i} {
460 set eta [expr { -6.00 + $i * 0.05}]
461 add EtaPhiBins $eta $PhiBins
462 }
463
464 for {set i 0} {$i <= 69} {incr i} {
465 set eta [expr { 2.55 + $i * 0.05}]
466 add EtaPhiBins $eta $PhiBins
467 }
468
469
470 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
471 add EnergyFraction {0} {1.0}
472 # energy fractions for e, gamma and pi0
473 add EnergyFraction {11} {0.0}
474 add EnergyFraction {22} {0.0}
475 add EnergyFraction {111} {0.0}
476 # energy fractions for muon, neutrinos and neutralinos
477 add EnergyFraction {12} {0.0}
478 add EnergyFraction {13} {0.0}
479 add EnergyFraction {14} {0.0}
480 add EnergyFraction {16} {0.0}
481 add EnergyFraction {1000022} {0.0}
482 add EnergyFraction {1000023} {0.0}
483 add EnergyFraction {1000025} {0.0}
484 add EnergyFraction {1000035} {0.0}
485 add EnergyFraction {1000045} {0.0}
486 # energy fractions for K0short and Lambda
487 # add EnergyFraction {310} {0.7}
488 # add EnergyFraction {3122} {0.7}
489
490 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
491 set ResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.03^2 + energy*0.50^2) + \
492 (abs(eta) > 1.7 && abs(eta) <= 4.0) * sqrt(energy^2*0.03^2 + energy*0.60^2) + \
493 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.10^2 + energy*1.00^2)}
494}
495
496
497#################
498# Electron filter
499#################
500
501module PdgCodeFilter ElectronFilter {
502 set InputArray HCal/eflowTracks
503 set OutputArray electrons
504 set Invert true
505 add PdgCode {11}
506 add PdgCode {-11}
507}
508
509
510######################
511# ChargedHadronFilter
512######################
513
514module PdgCodeFilter ChargedHadronFilter {
515 set InputArray HCal/eflowTracks
516 set OutputArray chargedHadrons
517
518 add PdgCode {11}
519 add PdgCode {-11}
520 add PdgCode {13}
521 add PdgCode {-13}
522}
523
524
525###################################################
526# Tower Merger (in case not using e-flow algorithm)
527###################################################
528
529module Merger Calorimeter {
530# add InputArray InputArray
531 add InputArray ECal/ecalTowers
532 add InputArray HCal/hcalTowers
533 add InputArray MuonMomentumSmearing/muons
534 set OutputArray towers
535}
536
537####################
538# Energy flow merger
539####################
540
541module Merger EFlowMerger {
542# add InputArray InputArray
543 add InputArray HCal/eflowTracks
544 add InputArray ECal/eflowPhotons
545 add InputArray HCal/eflowNeutralHadrons
546 set OutputArray eflow
547}
548
549######################
550# EFlowFilter
551######################
552
553module PdgCodeFilter EFlowFilter {
554 set InputArray EFlowMerger/eflow
555 set OutputArray eflow
556
557 add PdgCode {11}
558 add PdgCode {-11}
559 add PdgCode {13}
560 add PdgCode {-13}
561}
562
563
564###################
565# Missing ET merger
566###################
567
568module Merger MissingET {
569# add InputArray InputArray
570 add InputArray EFlowMerger/eflow
571 set MomentumOutputArray momentum
572}
573
574
575
576##################
577# Scalar HT merger
578##################
579
580module Merger ScalarHT {
581# add InputArray InputArray
582 add InputArray EFlowMerger/eflow
583 set EnergyOutputArray energy
584}
585
586#################
587# Neutrino Filter
588#################
589
590module PdgCodeFilter NeutrinoFilter {
591
592 set InputArray Delphes/stableParticles
593 set OutputArray filteredParticles
594
595 set PTMin 0.0
596
597 add PdgCode {12}
598 add PdgCode {14}
599 add PdgCode {16}
600 add PdgCode {-12}
601 add PdgCode {-14}
602 add PdgCode {-16}
603
604}
605
606
607#########################
608# Gen Missing ET merger
609########################
610
611module Merger GenMissingET {
612
613# add InputArray InputArray
614 add InputArray NeutrinoFilter/filteredParticles
615 set MomentumOutputArray momentum
616}
617
618
619#####################
620# MC truth jet finder
621#####################
622
623# TBC: is jet radius fine?
624
625module FastJetFinder GenJetFinder04 {
626 set InputArray NeutrinoFilter/filteredParticles
627
628 set OutputArray jets
629
630 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
631 set JetAlgorithm 6
632 set ParameterR 0.4
633
634 set ComputeNsubjettiness 1
635 set Beta 1.0
636 set AxisMode 4
637
638 set ComputeSoftDrop 1
639 set BetaSoftDrop 0.0
640 set SymmetryCutSoftDrop 0.1
641 set R0SoftDrop 0.4
642
643 set JetPTMin 25.0
644}
645
646
647
648##################
649# Fast Jet finder
650##################
651
652module FastJetFinder FastJetFinder04 {
653 set InputArray EFlowMerger/eflow
654
655 set OutputArray jets
656
657 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
658 set JetAlgorithm 6
659 set ParameterR 0.4
660
661 set ComputeNsubjettiness 1
662 set Beta 1.0
663 set AxisMode 4
664
665 set ComputeSoftDrop 1
666 set BetaSoftDrop 0.0
667 set SymmetryCutSoftDrop 0.1
668 set R0SoftDrop 0.4
669
670 set JetPTMin 25.0
671}
672
673
674##################
675# Fast Jet finder
676##################
677
678module FastJetFinder CaloJetFinder04 {
679 set InputArray Calorimeter/towers
680
681 set OutputArray jets
682
683 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
684 set JetAlgorithm 6
685 set ParameterR 0.4
686
687 set ComputeNsubjettiness 1
688 set Beta 1.0
689 set AxisMode 4
690
691 set ComputeSoftDrop 1
692 set BetaSoftDrop 0.0
693 set SymmetryCutSoftDrop 0.1
694 set R0SoftDrop 0.4
695
696 set JetPTMin 25.0
697}
698
699##################
700# Fast Jet finder
701##################
702
703module FastJetFinder TrackJetFinder04 {
704 set InputArray TrackMerger/tracks
705
706 set OutputArray jets
707
708 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
709 set JetAlgorithm 6
710 set ParameterR 0.4
711
712 set ComputeNsubjettiness 1
713 set Beta 1.0
714 set AxisMode 4
715
716 set ComputeSoftDrop 1
717 set BetaSoftDrop 0.0
718 set SymmetryCutSoftDrop 0.1
719 set R0SoftDrop 0.4
720
721 set JetPTMin 25.0
722}
723
724
725##################
726# Jet Energy Scale
727##################
728
729module EnergyScale JetEnergyScale {
730 set InputArray FastJetFinder04/jets
731 set OutputArray jets
732
733 # scale formula for jets
734 set ScaleFormula {1.00}
735}
736
737
738########################
739# Jet Flavor Association
740########################
741
742module JetFlavorAssociation JetFlavorAssociation {
743
744 set PartonInputArray Delphes/partons
745 set ParticleInputArray Delphes/allParticles
746 set ParticleLHEFInputArray Delphes/allParticlesLHEF
747 set JetInputArray JetEnergyScale/jets
748
749 set DeltaR 0.5
750 set PartonPTMin 5.0
751 set PartonEtaMax 6.0
752
753}
754
755###################
756# Photon efficiency
757###################
758
759module Efficiency PhotonEfficiency {
760 set InputArray ECal/eflowPhotons
761 set OutputArray photons
762
763 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
764
765 set EfficiencyFormula {
766 (pt <= 1.0) * (0.00) + \
767 (abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0) * (0.70) +
768 (abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) +
769 (abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
770
771 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0) * (0.60) +
772 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 5.0 && pt < 10.0) * (0.80) +
773 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) +
774
775 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 1.0 && pt < 5.0) * (0.50) + \
776 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 5.0 && pt < 10.0) * (0.70) + \
777 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 10.0) * (0.80) + \
778 (abs(eta) > 6.0) * (0.00)}
779
780}
781
782##################
783# Photon isolation
784##################
785
786# TBC: check values for iso cuts
787
788module Isolation PhotonIsolation {
789 set CandidateInputArray PhotonEfficiency/photons
790 set IsolationInputArray EFlowFilter/eflow
791
792 set OutputArray photons
793
794 set DeltaRMax 0.3
795
796 set PTMin 0.5
797
798 set PTRatioMax 0.1
799}
800
801
802####################
803# Electron isolation
804####################
805
806# TBC: check values for iso cuts
807
808module Isolation ElectronIsolation {
809 set CandidateInputArray ElectronFilter/electrons
810 set IsolationInputArray EFlowFilter/eflow
811
812 set OutputArray electrons
813
814 set DeltaRMax 0.3
815
816 set PTMin 0.5
817
818 set PTRatioMax 0.1
819}
820
821
822################
823# Muon isolation
824################
825
826# TBC: check values for iso cuts
827
828module Isolation MuonIsolation {
829 set CandidateInputArray MuonMomentumSmearing/muons
830 set IsolationInputArray EFlowFilter/eflow
831
832 set OutputArray muons
833
834 set DeltaRMax 0.3
835
836 set PTMin 0.5
837
838 set PTRatioMax 0.2
839}
840
841
842
843###########
844# b-tagging
845###########
846
847module BTagging BTagging {
848 set JetInputArray JetEnergyScale/jets
849
850 set BitNumber 0
851
852 add EfficiencyFormula {0} {
853
854 (pt <= 10.0) * (0.00) +
855 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.01) + \
856 (abs(eta) < 2.5) * (pt > 500.0 && pt < 15000.0) * (0.01)*(1.0 - pt/15000.) + \
857 (abs(eta) < 2.5) * (pt > 15000.0) * (0.00) + \
858 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.0075) + \
859 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0.0075)*(1.0 - pt/15000.) + \
860 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 15000.0) * (0.000) + \
861 (abs(eta) > 4.0) * (0.00)}
862
863 add EfficiencyFormula {4} {
864
865 (pt <= 10.0) * (0.00) +
866 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.05) + \
867 (abs(eta) < 2.5) * (pt > 500.0 && pt < 15000.0) * (0.05)*(1.0 - pt/15000.) + \
868 (abs(eta) < 2.5) * (pt > 15000.0) * (0.000) + \
869 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.03) + \
870 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0.03)*(1.0 - pt/15000.) + \
871 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 15000.0) * (0.000) + \
872 (abs(eta) > 4.0) * (0.00)}
873
874 add EfficiencyFormula {5} {
875
876 (pt <= 10.0) * (0.00) +
877 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.85) +
878 (abs(eta) < 2.5) * (pt > 500.0 && pt < 15000.0) * (0.85)*(1.0 - pt/15000.) +
879 (abs(eta) < 2.5) * (pt > 15000.0) * (0.000) +
880 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.64) +
881 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0.64)*(1.0 - pt/15000.) +
882 (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 15000.0) * (0.000) +
883 (abs(eta) >= 4.0) * (0.00)}
884
885}
886
887###########
888# c-tagging
889###########
890
891module BTagging CTagging {
892 set JetInputArray JetEnergyScale/jets
893
894 set BitNumber 1
895
896 add EfficiencyFormula {0} {
897
898 (pt <= 10.0) * (0.00) +
899 (abs(eta) < 4.0) * (pt > 10.0) * (0.01) + \
900 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
901
902 add EfficiencyFormula {4} {
903
904 (pt <= 10.0) * (0.00) +
905 (abs(eta) < 4.0) * (pt > 10.0) * (0.25) + \
906 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
907
908 add EfficiencyFormula {5} {
909
910 (pt <= 10.0) * (0.00) +
911 (abs(eta) < 4.0) * (pt > 10.0) * (0.03) + \
912 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
913
914}
915
916
917#############
918# tau-tagging
919#############
920
921
922module TauTagging TauTagging {
923 set ParticleInputArray Delphes/allParticles
924 set PartonInputArray Delphes/partons
925 set JetInputArray JetEnergyScale/jets
926
927 set DeltaR 0.5
928
929 set TauPTMin 1.0
930
931 set TauEtaMax 4.0
932
933 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
934 add EfficiencyFormula {0} {
935
936 (pt <= 10.0) * (0.00) +
937 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.01) + \
938 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.01) *(8./9. - pt/30000.) + \
939 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
940 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.0075) + \
941 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.0075)*(8./9. - pt/30000.) + \
942 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
943 (abs(eta) > 4.0) * (0.00)}
944
945 add EfficiencyFormula {11} {
946
947 (pt <= 10.0) * (0.00) +
948 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.005) + \
949 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.005) *(8./9. - pt/30000.) + \
950 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
951 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.00375) + \
952 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.00375)*(8./9. - pt/30000.) + \
953 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
954 (abs(eta) > 4.0) * (0.00)}
955
956 add EfficiencyFormula {15} {
957
958 (pt <= 10.0) * (0.00) +
959 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.8) + \
960 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.8) *(8./9. - pt/30000.) + \
961 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
962 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.65) + \
963 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.65)*(8./9. - pt/30000.) + \
964 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
965 (abs(eta) > 4.0) * (0.00)}
966
967}
968
969#####################################################
970# Find uniquely identified photons/electrons/tau/jets
971#####################################################
972
973module UniqueObjectFinder UniqueObjectFinder {
974# earlier arrays take precedence over later ones
975# add InputArray InputArray OutputArray
976 add InputArray PhotonIsolation/photons photons
977 add InputArray ElectronIsolation/electrons electrons
978 add InputArray MuonIsolation/muons muons
979 add InputArray JetEnergyScale/jets jets
980
981}
982
983
984##################
985# ROOT tree writer
986##################
987
988module TreeWriter TreeWriter {
989# add Branch InputArray BranchName BranchClass
990 add Branch Delphes/allParticles Particle GenParticle
991
992 add Branch GenMissingET/momentum GenMissingET MissingET
993
994 add Branch TimeSmearing/tracks Track Track
995 add Branch Calorimeter/towers Tower Tower
996
997 add Branch HCal/eflowTracks EFlowTrack Track
998 add Branch ECal/eflowPhotons EFlowPhoton Tower
999 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1000
1001 add Branch UniqueObjectFinder/photons Photon Photon
1002 add Branch UniqueObjectFinder/electrons Electron Electron
1003 add Branch UniqueObjectFinder/muons Muon Muon
1004 add Branch UniqueObjectFinder/jets Jet Jet
1005
1006 add Branch GenJetFinder04/jets GenJet04 Jet
1007 add Branch FastJetFinder04/jets ParticleFlowJet04 Jet
1008 add Branch CaloJetFinder04/jets CaloJet04 Jet
1009 add Branch TrackJetFinder04/jets TrackJet04 Jet
1010
1011 add Branch MissingET/momentum MissingET MissingET
1012 add Branch ScalarHT/energy ScalarHT ScalarHT
1013 add Branch VertexFinderDA4D/vertices Vertex4D Vertex
1014}
1015
Note: See TracBrowser for help on using the repository browser.