Fork me on GitHub

source: git/cards/FCC/FCChh.tcl@ 46936ad

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 46936ad was 46936ad, checked in by Michele Selvaggi <michele.selvaggi@…>, 6 years ago

updated FCChh card with new DTF

  • Property mode set to 100644
File size: 28.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 ExecutionPath {
15
16 ParticlePropagator
17 TrackMergerProp
18
19 DenseProp
20 DenseMergeTracks
21 DenseTrackFilter
22
23 ChargedHadronTrackingEfficiency
24 ElectronTrackingEfficiency
25 MuonTrackingEfficiency
26
27 ChargedHadronMomentumSmearing
28 ElectronMomentumSmearing
29 MuonMomentumSmearing
30
31 TrackMerger
32
33 ECal
34 HCal
35
36 Calorimeter
37 EFlowMerger
38 EFlowFilter
39
40 PhotonEfficiency
41 PhotonIsolation
42
43 ElectronFilter
44 ElectronIsolation
45
46 ChargedHadronFilter
47
48 MuonIsolation
49
50 NeutrinoFilter
51
52 MissingET
53 GenMissingET
54
55 GenJetFinder02
56 GenJetFinder04
57 GenJetFinder08
58
59 FastJetFinder02
60 FastJetFinder04
61 FastJetFinder08
62
63 CaloJetFinder02
64 CaloJetFinder04
65 CaloJetFinder08
66
67 TrackJetFinder02
68 TrackJetFinder04
69 TrackJetFinder08
70
71 JetEnergyScale
72
73 JetFlavorAssociation
74
75 BTagging
76 CTagging
77 TauTagging
78
79 ScalarHT
80
81 UniqueObjectFinder
82
83 TreeWriter
84}
85
86#####################################
87# Track propagation to calorimeters
88#####################################
89
90module ParticlePropagator ParticlePropagator {
91 set InputArray Delphes/stableParticles
92 set OutputArray stableParticles
93 set ChargedHadronOutputArray chargedHadrons
94 set ElectronOutputArray electrons
95 set MuonOutputArray muons
96
97 # radius of the magnetic field coverage, in m
98 set Radius 1.5
99 # half-length of the magnetic field coverage, in m
100 set HalfLength 5
101
102 # magnetic field
103 set Bz 4.0
104}
105
106
107##############
108# Track merger
109##############
110
111module Merger TrackMergerProp {
112# add InputArray InputArray
113 add InputArray ParticlePropagator/chargedHadrons
114 add InputArray ParticlePropagator/electrons
115 add InputArray ParticlePropagator/muons
116 set OutputArray tracks
117}
118
119
120####################################
121# Track propagation to pseudo-pixel
122####################################
123
124module ParticlePropagator DenseProp {
125
126 set InputArray TrackMergerProp/tracks
127
128 # radius of the magnetic field coverage, in m
129 set Radius 0.45
130 set RadiusMax 1.5
131 # half-length of the magnetic field coverage, in m
132 set HalfLength 0.8
133 set HalfLengthMax 5
134
135 # magnetic field
136 set Bz 4.0
137}
138
139#####################
140# Dense Track merger
141#####################
142
143module Merger DenseMergeTracks {
144# add InputArray InputArray
145 add InputArray DenseProp/chargedHadrons
146 add InputArray DenseProp/electrons
147 add InputArray DenseProp/muons
148 set OutputArray tracks
149}
150
151
152######################
153# Dense Track Filter
154######################
155
156module DenseTrackFilter DenseTrackFilter {
157
158 set TrackInputArray DenseMergeTracks/tracks
159
160 set TrackOutputArray tracks
161 set ChargedHadronOutputArray chargedHadrons
162 set ElectronOutputArray electrons
163 set MuonOutputArray muons
164
165 set EtaPhiRes 0.001
166 set EtaMax 6.0
167
168 set pi [expr {acos(-1)}]
169
170 set nbins_phi [expr {$pi/$EtaPhiRes} ]
171 set nbins_phi [expr {int($nbins_phi)} ]
172
173 set PhiBins {}
174 for {set i -$nbins_phi} {$i <= $nbins_phi} {incr i} {
175 add PhiBins [expr {$i * $pi/$nbins_phi}]
176 }
177}
178
179####################################
180# Charged hadron tracking efficiency
181####################################
182
183module Efficiency ChargedHadronTrackingEfficiency {
184 set InputArray DenseTrackFilter/chargedHadrons
185 set OutputArray chargedHadrons
186
187 # TBC (which eta_max ? which pT min?)
188
189 # tracking efficiency formula for charged hadrons
190
191 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
192(abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
193(abs(eta) <= 2.5) * (pt > 1) * (0.95) + \
194(abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
195(abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
196(abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
197(abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
198(abs(eta) > 6.0) * (0.00)}
199
200}
201
202##############################
203# Electron tracking efficiency
204##############################
205
206module Efficiency ElectronTrackingEfficiency {
207 set InputArray DenseTrackFilter/electrons
208 set OutputArray electrons
209
210# TBC (which eta_max ?)
211# putting same as charged hadrons for now...
212
213 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
214 (abs(eta) <= 2.5) * (pt > 0.5 && pt <= 1) * (0.90) + \
215 (abs(eta) <= 2.5) * (pt > 1) * (0.95) + \
216 (abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 0.5 && pt <= 1) * (0.85) + \
217 (abs(eta) > 2.5 && abs(eta) <= 4) * (pt > 1) * (0.90) + \
218 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 0.5 && pt <= 1) * (0.80) + \
219 (abs(eta) > 4 && abs(eta) <= 6) * (pt > 1.0) * (0.85) + \
220 (abs(eta) > 6.0) * (0.00)}
221
222}
223
224##########################
225# Muon tracking efficiency
226##########################
227
228module Efficiency MuonTrackingEfficiency {
229 set InputArray DenseTrackFilter/muons
230 set OutputArray muons
231
232# TBC (which eta_max ? why eff = 0 for 4 < eta < 6 ? for now put the same as central)
233# what about high pT ?
234 # tracking efficiency formula for muons
235 set EfficiencyFormula { (pt <= 0.5) * (0.00) + \
236 (abs(eta) <= 6.0) * (pt > 0.5 && pt <= 1) * (0.90) + \
237 (abs(eta) <= 6.0) * (pt > 1) * (0.99) + \
238 (abs(eta) > 6.0) * (0.00)}
239
240}
241
242########################################
243# Momentum resolution for charged tracks
244########################################
245
246module MomentumSmearing ChargedHadronMomentumSmearing {
247 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
248 set OutputArray chargedHadrons
249
250 source momentumResolutionVsP.tcl
251}
252
253
254###################################
255# Momentum resolution for electrons
256###################################
257
258module MomentumSmearing ElectronMomentumSmearing {
259 set InputArray ElectronTrackingEfficiency/electrons
260 set OutputArray electrons
261
262 source momentumResolutionVsP.tcl
263}
264
265
266###############################
267# Momentum resolution for muons
268###############################
269
270module MomentumSmearing MuonMomentumSmearing {
271 set InputArray MuonTrackingEfficiency/muons
272 set OutputArray muons
273
274 # TBC for just putting tracker resolution/ need to add improvement at high pT
275
276 source muonMomentumResolutionVsP.tcl
277}
278
279##############
280# Track merger
281##############
282
283module Merger TrackMerger {
284# add InputArray InputArray
285 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
286 add InputArray ElectronMomentumSmearing/electrons
287 add InputArray MuonMomentumSmearing/muons
288 set OutputArray tracks
289}
290
291
292#############
293# ECAL
294#############
295
296# TBC : calos seems ok, check eta max value though.
297
298module SimpleCalorimeter ECal {
299 set ParticleInputArray ParticlePropagator/stableParticles
300 set TrackInputArray TrackMerger/tracks
301
302 set TowerOutputArray ecalTowers
303 set EFlowTrackOutputArray eflowTracks
304 set EFlowTowerOutputArray eflowPhotons
305
306 set IsEcal true
307
308 set EnergyMin 0.5
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.012 rad towers up to eta = 2.5 (barrel)
320 set PhiBins {}
321 for {set i -256} {$i <= 256} {incr i} {
322 add PhiBins [expr {$i * $pi/256.0}]
323 }
324
325 # 0.01 unit in eta up to eta = 2.5 (barrel)
326 for {set i -249} {$i <= 250} {incr i} {
327 set eta [expr {$i * 0.01}]
328 add EtaPhiBins $eta $PhiBins
329 }
330
331 # 0.025 rad between 2.5 and 6.0
332 set PhiBins {}
333 for {set i -128} {$i <= 128} {incr i} {
334 add PhiBins [expr {$i * $pi/128.0}]
335 }
336
337 # 0.025 unit in eta between 2.5 and 6.0
338 for {set i 0} {$i <= 140} {incr i} {
339 set eta [expr { -6.00 + $i * 0.025}]
340 add EtaPhiBins $eta $PhiBins
341 }
342
343 for {set i 0} {$i <= 139} {incr i} {
344 set eta [expr { 2.525 + $i * 0.025}]
345 add EtaPhiBins $eta $PhiBins
346 }
347
348 # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
349
350 add EnergyFraction {0} {0.0}
351 # energy fractions for e, gamma and pi0
352 add EnergyFraction {11} {1.0}
353 add EnergyFraction {22} {1.0}
354 add EnergyFraction {111} {1.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.3}
367 # add EnergyFraction {3122} {0.3}
368
369 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
370 set ResolutionFormula { (abs(eta) <= 4.0) * sqrt(energy^2*0.007^2 + energy*0.10^2) + \
371 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.035^2 + energy*0.30^2)}
372
373
374}
375
376#############
377# HCAL
378#############
379
380module SimpleCalorimeter HCal {
381 set ParticleInputArray ParticlePropagator/stableParticles
382 set TrackInputArray ECal/eflowTracks
383
384 set TowerOutputArray hcalTowers
385 set EFlowTrackOutputArray eflowTracks
386 set EFlowTowerOutputArray eflowNeutralHadrons
387
388 set IsEcal false
389
390 set EnergyMin 1.0
391 set EnergySignificanceMin 2.0
392
393 set SmearTowerCenter true
394
395 set pi [expr {acos(-1)}]
396
397 # lists of the edges of each tower in eta and phi
398 # each list starts with the lower edge of the first tower
399 # the list ends with the higher edged of the last tower
400
401 # 0.025 rad towers up to eta = 2.5 (barrel)
402 set PhiBins {}
403 for {set i -128} {$i <= 128} {incr i} {
404 add PhiBins [expr {$i * $pi/128.0}]
405 }
406
407 # 0.025 unit in eta up to eta = 2.5 (barrel)
408 for {set i -99} {$i <= 100} {incr i} {
409 set eta [expr {$i * 0.025}]
410 add EtaPhiBins $eta $PhiBins
411 }
412
413 # 0.05 rad between 2.5 and 6.0
414 set PhiBins {}
415 for {set i -64} {$i <= 64} {incr i} {
416 add PhiBins [expr {$i * $pi/64.0}]
417 }
418
419 # 0.05 unit in eta between 2.5 and 6.0
420 for {set i 0} {$i <= 70} {incr i} {
421 set eta [expr { -6.00 + $i * 0.05}]
422 add EtaPhiBins $eta $PhiBins
423 }
424
425 for {set i 0} {$i <= 69} {incr i} {
426 set eta [expr { 2.55 + $i * 0.05}]
427 add EtaPhiBins $eta $PhiBins
428 }
429
430
431 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
432 add EnergyFraction {0} {1.0}
433 # energy fractions for e, gamma and pi0
434 add EnergyFraction {11} {0.0}
435 add EnergyFraction {22} {0.0}
436 add EnergyFraction {111} {0.0}
437 # energy fractions for muon, neutrinos and neutralinos
438 add EnergyFraction {12} {0.0}
439 add EnergyFraction {13} {0.0}
440 add EnergyFraction {14} {0.0}
441 add EnergyFraction {16} {0.0}
442 add EnergyFraction {1000022} {0.0}
443 add EnergyFraction {1000023} {0.0}
444 add EnergyFraction {1000025} {0.0}
445 add EnergyFraction {1000035} {0.0}
446 add EnergyFraction {1000045} {0.0}
447 # energy fractions for K0short and Lambda
448 # add EnergyFraction {310} {0.7}
449 # add EnergyFraction {3122} {0.7}
450
451 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
452 set ResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.03^2 + energy*0.50^2) + \
453 (abs(eta) > 1.7 && abs(eta) <= 4.0) * sqrt(energy^2*0.03^2 + energy*0.60^2) + \
454 (abs(eta) > 4.0 && abs(eta) <= 6.0) * sqrt(energy^2*0.10^2 + energy*1.00^2)}
455}
456
457#################
458# Electron filter
459#################
460
461module PdgCodeFilter ElectronFilter {
462 set InputArray HCal/eflowTracks
463 set OutputArray electrons
464 set Invert true
465 add PdgCode {11}
466 add PdgCode {-11}
467}
468
469
470######################
471# ChargedHadronFilter
472######################
473
474module PdgCodeFilter ChargedHadronFilter {
475 set InputArray HCal/eflowTracks
476 set OutputArray chargedHadrons
477
478 add PdgCode {11}
479 add PdgCode {-11}
480 add PdgCode {13}
481 add PdgCode {-13}
482}
483
484
485###################################################
486# Tower Merger (in case not using e-flow algorithm)
487###################################################
488
489module Merger Calorimeter {
490# add InputArray InputArray
491 add InputArray ECal/ecalTowers
492 add InputArray HCal/hcalTowers
493 add InputArray MuonMomentumSmearing/muons
494 set OutputArray towers
495}
496
497####################
498# Energy flow merger
499####################
500
501module Merger EFlowMerger {
502# add InputArray InputArray
503 add InputArray HCal/eflowTracks
504 add InputArray ECal/eflowPhotons
505 add InputArray HCal/eflowNeutralHadrons
506 set OutputArray eflow
507}
508
509######################
510# EFlowFilter
511######################
512
513module PdgCodeFilter EFlowFilter {
514 set InputArray EFlowMerger/eflow
515 set OutputArray eflow
516
517 add PdgCode {11}
518 add PdgCode {-11}
519 add PdgCode {13}
520 add PdgCode {-13}
521}
522
523
524###################
525# Missing ET merger
526###################
527
528module Merger MissingET {
529# add InputArray InputArray
530 add InputArray EFlowMerger/eflow
531 set MomentumOutputArray momentum
532}
533
534
535
536##################
537# Scalar HT merger
538##################
539
540module Merger ScalarHT {
541# add InputArray InputArray
542 add InputArray EFlowMerger/eflow
543 set EnergyOutputArray energy
544}
545
546#################
547# Neutrino Filter
548#################
549
550module PdgCodeFilter NeutrinoFilter {
551
552 set InputArray Delphes/stableParticles
553 set OutputArray filteredParticles
554
555 set PTMin 0.0
556
557 add PdgCode {12}
558 add PdgCode {14}
559 add PdgCode {16}
560 add PdgCode {-12}
561 add PdgCode {-14}
562 add PdgCode {-16}
563
564}
565
566
567#########################
568# Gen Missing ET merger
569########################
570
571module Merger GenMissingET {
572
573# add InputArray InputArray
574 add InputArray NeutrinoFilter/filteredParticles
575 set MomentumOutputArray momentum
576}
577
578
579#####################
580# MC truth jet finder
581#####################
582
583# TBC: is jet radius fine?
584
585module FastJetFinder GenJetFinder02 {
586 set InputArray NeutrinoFilter/filteredParticles
587
588 set OutputArray jets
589
590 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
591 set JetAlgorithm 6
592 set ParameterR 0.2
593
594 set ComputeNsubjettiness 1
595 set Beta 1.0
596 set AxisMode 4
597
598 set ComputeSoftDrop 1
599 set BetaSoftDrop 0.0
600 set SymmetryCutSoftDrop 0.1
601 set R0SoftDrop 0.2
602
603 set JetPTMin 25.0
604}
605
606
607#####################
608# MC truth jet finder
609#####################
610
611# TBC: is jet radius fine?
612
613module FastJetFinder GenJetFinder04 {
614 set InputArray NeutrinoFilter/filteredParticles
615
616 set OutputArray jets
617
618 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
619 set JetAlgorithm 6
620 set ParameterR 0.4
621
622 set ComputeNsubjettiness 1
623 set Beta 1.0
624 set AxisMode 4
625
626 set ComputeSoftDrop 1
627 set BetaSoftDrop 0.0
628 set SymmetryCutSoftDrop 0.1
629 set R0SoftDrop 0.4
630
631 set JetPTMin 25.0
632}
633#####################
634# MC truth jet finder
635#####################
636
637# TBC: is jet radius fine?
638
639module FastJetFinder GenJetFinder08 {
640 set InputArray NeutrinoFilter/filteredParticles
641
642 set OutputArray jets
643
644 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
645 set JetAlgorithm 6
646 set ParameterR 0.8
647
648 set ComputeNsubjettiness 1
649 set Beta 1.0
650 set AxisMode 4
651
652 set ComputeSoftDrop 1
653 set BetaSoftDrop 0.0
654 set SymmetryCutSoftDrop 0.1
655 set R0SoftDrop 0.8
656
657 set JetPTMin 25.0
658}
659
660
661##################
662# Fast Jet finder
663##################
664
665module FastJetFinder FastJetFinder02 {
666 set InputArray EFlowMerger/eflow
667
668 set OutputArray jets
669
670 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
671 set JetAlgorithm 6
672 set ParameterR 0.2
673
674 set ComputeNsubjettiness 1
675 set Beta 1.0
676 set AxisMode 4
677
678 set ComputeSoftDrop 1
679 set BetaSoftDrop 0.0
680 set SymmetryCutSoftDrop 0.1
681 set R0SoftDrop 0.2
682
683 set JetPTMin 25.0
684}
685
686##################
687# Fast Jet finder
688##################
689
690module FastJetFinder FastJetFinder04 {
691 set InputArray EFlowMerger/eflow
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##################
713# Fast Jet finder
714##################
715
716module FastJetFinder FastJetFinder08 {
717 set InputArray EFlowMerger/eflow
718
719 set OutputArray jets
720
721 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
722 set JetAlgorithm 6
723 set ParameterR 0.8
724
725 set ComputeNsubjettiness 1
726 set Beta 1.0
727 set AxisMode 4
728
729 set ComputeSoftDrop 1
730 set BetaSoftDrop 0.0
731 set SymmetryCutSoftDrop 0.1
732 set R0SoftDrop 0.8
733
734 set JetPTMin 25.0
735}
736
737
738
739##################
740# Fast Jet finder
741##################
742
743module FastJetFinder CaloJetFinder02 {
744 set InputArray Calorimeter/towers
745
746 set OutputArray jets
747
748 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
749 set JetAlgorithm 6
750 set ParameterR 0.2
751
752 set ComputeNsubjettiness 1
753 set Beta 1.0
754 set AxisMode 4
755
756 set ComputeSoftDrop 1
757 set BetaSoftDrop 0.0
758 set SymmetryCutSoftDrop 0.1
759 set R0SoftDrop 0.2
760
761 set JetPTMin 25.0
762}
763
764##################
765# Fast Jet finder
766##################
767
768module FastJetFinder CaloJetFinder04 {
769 set InputArray Calorimeter/towers
770
771 set OutputArray jets
772
773 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
774 set JetAlgorithm 6
775 set ParameterR 0.4
776
777 set ComputeNsubjettiness 1
778 set Beta 1.0
779 set AxisMode 4
780
781 set ComputeSoftDrop 1
782 set BetaSoftDrop 0.0
783 set SymmetryCutSoftDrop 0.1
784 set R0SoftDrop 0.4
785
786 set JetPTMin 25.0
787}
788
789
790##################
791# Fast Jet finder
792##################
793
794module FastJetFinder CaloJetFinder08 {
795 set InputArray Calorimeter/towers
796
797 set OutputArray jets
798
799 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
800 set JetAlgorithm 6
801 set ParameterR 0.8
802
803 set ComputeNsubjettiness 1
804 set Beta 1.0
805 set AxisMode 4
806
807 set ComputeSoftDrop 1
808 set BetaSoftDrop 0.0
809 set SymmetryCutSoftDrop 0.1
810 set R0SoftDrop 0.8
811
812 set JetPTMin 25.0
813}
814
815
816##################
817# Fast Jet finder
818##################
819
820module FastJetFinder TrackJetFinder02 {
821 set InputArray TrackMerger/tracks
822
823 set OutputArray jets
824
825 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
826 set JetAlgorithm 6
827 set ParameterR 0.2
828
829 set ComputeNsubjettiness 1
830 set Beta 1.0
831 set AxisMode 4
832
833 set ComputeSoftDrop 1
834 set BetaSoftDrop 0.0
835 set SymmetryCutSoftDrop 0.1
836 set R0SoftDrop 0.2
837
838 set JetPTMin 25.0
839}
840
841##################
842# Fast Jet finder
843##################
844
845module FastJetFinder TrackJetFinder04 {
846 set InputArray TrackMerger/tracks
847
848 set OutputArray jets
849
850 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
851 set JetAlgorithm 6
852 set ParameterR 0.4
853
854 set ComputeNsubjettiness 1
855 set Beta 1.0
856 set AxisMode 4
857
858 set ComputeSoftDrop 1
859 set BetaSoftDrop 0.0
860 set SymmetryCutSoftDrop 0.1
861 set R0SoftDrop 0.4
862
863 set JetPTMin 25.0
864}
865
866
867##################
868# Fast Jet finder
869##################
870
871module FastJetFinder TrackJetFinder08 {
872 set InputArray TrackMerger/tracks
873
874 set OutputArray jets
875
876 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
877 set JetAlgorithm 6
878 set ParameterR 0.8
879
880 set ComputeNsubjettiness 1
881 set Beta 1.0
882 set AxisMode 4
883
884 set ComputeSoftDrop 1
885 set BetaSoftDrop 0.0
886 set SymmetryCutSoftDrop 0.1
887 set R0SoftDrop 0.8
888
889 set JetPTMin 25.0
890}
891
892
893
894##################
895# Jet Energy Scale
896##################
897
898module EnergyScale JetEnergyScale {
899 set InputArray FastJetFinder04/jets
900 set OutputArray jets
901
902 # scale formula for jets
903 set ScaleFormula {1.00}
904}
905
906
907########################
908# Jet Flavor Association
909########################
910
911module JetFlavorAssociation JetFlavorAssociation {
912
913 set PartonInputArray Delphes/partons
914 set ParticleInputArray Delphes/allParticles
915 set ParticleLHEFInputArray Delphes/allParticlesLHEF
916 set JetInputArray JetEnergyScale/jets
917
918 set DeltaR 0.5
919 set PartonPTMin 5.0
920 set PartonEtaMax 6.0
921
922}
923
924###################
925# Photon efficiency
926###################
927
928module Efficiency PhotonEfficiency {
929 set InputArray ECal/eflowPhotons
930 set OutputArray photons
931
932 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
933
934 set EfficiencyFormula {
935 (pt <= 1.0) * (0.00) + \
936 (abs(eta) <= 2.5) * (pt > 1.0 && pt < 5.0) * (0.70) +
937 (abs(eta) <= 2.5) * (pt > 5.0 && pt < 10.0) * (0.85) +
938 (abs(eta) <= 2.5) * (pt > 10.0) * (0.95) +
939
940 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 1.0 && pt < 5.0) * (0.60) +
941 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 5.0 && pt < 10.0) * (0.80) +
942 (abs(eta) > 2.5 && abs(eta) <= 4.0) * (pt > 10.0) * (0.90) +
943
944 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 1.0 && pt < 5.0) * (0.50) + \
945 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 5.0 && pt < 10.0) * (0.70) + \
946 (abs(eta) > 4.0 && abs(eta) <= 6.0) * (pt > 10.0) * (0.80) + \
947 (abs(eta) > 6.0) * (0.00)}
948
949}
950
951##################
952# Photon isolation
953##################
954
955# TBC: check values for iso cuts
956
957module Isolation PhotonIsolation {
958 set CandidateInputArray PhotonEfficiency/photons
959 set IsolationInputArray EFlowFilter/eflow
960
961 set OutputArray photons
962
963 set DeltaRMax 0.3
964
965 set PTMin 0.5
966
967 set PTRatioMax 0.1
968}
969
970
971####################
972# Electron isolation
973####################
974
975# TBC: check values for iso cuts
976
977module Isolation ElectronIsolation {
978 set CandidateInputArray ElectronFilter/electrons
979 set IsolationInputArray EFlowFilter/eflow
980
981 set OutputArray electrons
982
983 set DeltaRMax 0.3
984
985 set PTMin 0.5
986
987 set PTRatioMax 0.1
988}
989
990
991################
992# Muon isolation
993################
994
995# TBC: check values for iso cuts
996
997module Isolation MuonIsolation {
998 set CandidateInputArray MuonMomentumSmearing/muons
999 set IsolationInputArray EFlowFilter/eflow
1000
1001 set OutputArray muons
1002
1003 set DeltaRMax 0.3
1004
1005 set PTMin 0.5
1006
1007 set PTRatioMax 0.2
1008}
1009
1010
1011
1012###########
1013# b-tagging
1014###########
1015
1016module BTagging BTagging {
1017 set JetInputArray JetEnergyScale/jets
1018
1019 set BitNumber 0
1020
1021 add EfficiencyFormula {0} {
1022
1023 (pt <= 10.0) * (0.00) +
1024 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.01) + \
1025 (abs(eta) < 2.5) * (pt > 500.0 && pt < 15000.0) * (0.01)*(1.0 - pt/15000.) + \
1026 (abs(eta) < 2.5) * (pt > 15000.0) * (0.00) + \
1027 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.0075) + \
1028 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0.0075)*(1.0 - pt/15000.) + \
1029 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 15000.0) * (0.000) + \
1030 (abs(eta) > 4.0) * (0.00)}
1031
1032 add EfficiencyFormula {4} {
1033
1034 (pt <= 10.0) * (0.00) +
1035 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.05) + \
1036 (abs(eta) < 2.5) * (pt > 500.0 && pt < 15000.0) * (0.05)*(1.0 - pt/15000.) + \
1037 (abs(eta) < 2.5) * (pt > 15000.0) * (0.000) + \
1038 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.03) + \
1039 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0.03)*(1.0 - pt/15000.) + \
1040 (abs(eta) < 2.5 && abs(eta) < 4.0) * (pt > 15000.0) * (0.000) + \
1041 (abs(eta) > 4.0) * (0.00)}
1042
1043 add EfficiencyFormula {5} {
1044
1045 (pt <= 10.0) * (0.00) +
1046 (abs(eta) < 2.5) * (pt > 10.0 && pt < 500) * (0.85) +
1047 (abs(eta) < 2.5) * (pt > 500.0 && pt < 15000.0) * (0.85)*(1.0 - pt/15000.) +
1048 (abs(eta) < 2.5) * (pt > 15000.0) * (0.000) +
1049 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 500) * (0.64) +
1050 (abs(eta) >= 2.5 && abs(eta) < 4.0) * (pt > 500.0 && pt < 15000.0) * (0.64)*(1.0 - pt/15000.) +
1051 (abs(eta) <= 2.5 && abs(eta) < 4.0) * (pt > 15000.0) * (0.000) +
1052 (abs(eta) >= 4.0) * (0.00)}
1053
1054}
1055
1056###########
1057# c-tagging
1058###########
1059
1060module BTagging CTagging {
1061 set JetInputArray JetEnergyScale/jets
1062
1063 set BitNumber 1
1064
1065 add EfficiencyFormula {0} {
1066
1067 (pt <= 10.0) * (0.00) +
1068 (abs(eta) < 4.0) * (pt > 10.0) * (0.01) + \
1069 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
1070
1071 add EfficiencyFormula {4} {
1072
1073 (pt <= 10.0) * (0.00) +
1074 (abs(eta) < 4.0) * (pt > 10.0) * (0.25) + \
1075 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
1076
1077 add EfficiencyFormula {5} {
1078
1079 (pt <= 10.0) * (0.00) +
1080 (abs(eta) < 4.0) * (pt > 10.0) * (0.03) + \
1081 (abs(eta) > 4.0) * (pt > 10.0) * (0.00)}
1082
1083}
1084
1085
1086#############
1087# tau-tagging
1088#############
1089
1090
1091module TauTagging TauTagging {
1092 set ParticleInputArray Delphes/allParticles
1093 set PartonInputArray Delphes/partons
1094 set JetInputArray JetEnergyScale/jets
1095
1096 set DeltaR 0.5
1097
1098 set TauPTMin 1.0
1099
1100 set TauEtaMax 4.0
1101
1102 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
1103 add EfficiencyFormula {0} {
1104
1105 (pt <= 10.0) * (0.00) +
1106 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.01) + \
1107 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.01) *(8./9. - pt/30000.) + \
1108 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
1109 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.0075) + \
1110 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.0075)*(8./9. - pt/30000.) + \
1111 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
1112 (abs(eta) > 4.0) * (0.00)}
1113
1114 add EfficiencyFormula {11} {
1115
1116 (pt <= 10.0) * (0.00) +
1117 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.005) + \
1118 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.005) *(8./9. - pt/30000.) + \
1119 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
1120 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.00375) + \
1121 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.00375)*(8./9. - pt/30000.) + \
1122 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
1123 (abs(eta) > 4.0) * (0.00)}
1124
1125 add EfficiencyFormula {15} {
1126
1127 (pt <= 10.0) * (0.00) +
1128 (abs(eta) < 2.5) * (pt > 10.0 && pt < 5000.0) * (0.6) + \
1129 (abs(eta) < 2.5) * (pt > 5000.0 && pt < 34000.0) * (0.6) *(8./9. - pt/30000.) + \
1130 (abs(eta) < 2.5) * (pt > 34000.0) * (0.000) + \
1131 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 10.0 && pt < 5000.0) * (0.45) + \
1132 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 5000.0 && pt < 34000.0) * (0.45)*(8./9. - pt/30000.) + \
1133 (abs(eta) > 2.5 && abs(eta) < 4.0) * (pt > 34000.0) * (0.00) + \
1134 (abs(eta) > 4.0) * (0.00)}
1135
1136}
1137
1138#####################################################
1139# Find uniquely identified photons/electrons/tau/jets
1140#####################################################
1141
1142module UniqueObjectFinder UniqueObjectFinder {
1143# earlier arrays take precedence over later ones
1144# add InputArray InputArray OutputArray
1145 add InputArray PhotonIsolation/photons photons
1146 add InputArray ElectronIsolation/electrons electrons
1147 add InputArray MuonIsolation/muons muons
1148 add InputArray JetEnergyScale/jets jets
1149
1150}
1151
1152
1153##################
1154# ROOT tree writer
1155##################
1156
1157module TreeWriter TreeWriter {
1158# add Branch InputArray BranchName BranchClass
1159 add Branch Delphes/allParticles Particle GenParticle
1160
1161 add Branch GenMissingET/momentum GenMissingET MissingET
1162 add Branch GenJetFinder02/jets GenJet Jet
1163
1164 add Branch TrackMerger/tracks Track Track
1165 add Branch Calorimeter/towers Tower Tower
1166
1167 add Branch HCal/eflowTracks EFlowTrack Track
1168 add Branch ECal/eflowPhotons EFlowPhoton Tower
1169 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
1170
1171 add Branch UniqueObjectFinder/photons Photon Photon
1172 add Branch UniqueObjectFinder/electrons Electron Electron
1173 add Branch UniqueObjectFinder/muons Muon Muon
1174 add Branch UniqueObjectFinder/jets Jet Jet
1175
1176 add Branch GenJetFinder02/jets GenJet02 Jet
1177 add Branch GenJetFinder04/jets GenJet04 Jet
1178 add Branch GenJetFinder08/jets GenJet08 Jet
1179
1180 add Branch FastJetFinder02/jets ParticleFlowJet02 Jet
1181 add Branch FastJetFinder04/jets ParticleFlowJet04 Jet
1182 add Branch FastJetFinder08/jets ParticleFlowJet08 Jet
1183
1184 add Branch CaloJetFinder02/jets CaloJet02 Jet
1185 add Branch CaloJetFinder04/jets CaloJet04 Jet
1186 add Branch CaloJetFinder08/jets CaloJet08 Jet
1187
1188 add Branch TrackJetFinder02/jets TrackJet02 Jet
1189 add Branch TrackJetFinder04/jets TrackJet04 Jet
1190 add Branch TrackJetFinder08/jets TrackJet08 Jet
1191
1192 add Branch MissingET/momentum MissingET MissingET
1193 add Branch ScalarHT/energy ScalarHT ScalarHT
1194}
1195
Note: See TracBrowser for help on using the repository browser.