Fork me on GitHub

source: git/cards/FCC/FCChh_PileUpVtx.tcl@ 0976f6a

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

test2

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