Fork me on GitHub

source: git/cards/delphes_card_CircularEE.tcl@ 0d65492

Last change on this file since 0d65492 was 201679f, checked in by Michele Selvaggi <michele.selvaggi@…>, 4 years ago

added 2nd tau tagging WP

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