Fork me on GitHub

source: git/cards/delphes_card_CircularEE.tcl@ 8dd735c

Last change on this file since 8dd735c was e7cccc2, checked in by Michele Selvaggi <michele.selvaggi@…>, 4 years ago

added alpha version of common card for circular ee studies for Snowmass

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