Fork me on GitHub

source: git/cards/delphes_card_CMS.tcl@ 66b8143

Last change on this file since 66b8143 was 15772b7, checked in by Michele Selvaggi <michele.selvaggi@…>, 8 years ago

make cards compliant with the FCCSW configuration
  • Property mode set to 100644
File size: 21.4 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 MissingET
38
39 NeutrinoFilter
40 GenJetFinder
41 GenMissingET
42
43 FastJetFinder
44
45 JetEnergyScale
46
47 JetFlavorAssociation
48
49 BTagging
50 TauTagging
51
52 UniqueObjectFinder
53
54 ScalarHT
55
56 TreeWriter
57}
58
59#################################
60# Propagate particles in cylinder
61#################################
62
63module ParticlePropagator ParticlePropagator {
64 set InputArray Delphes/stableParticles
65
66 set OutputArray stableParticles
67 set ChargedHadronOutputArray chargedHadrons
68 set ElectronOutputArray electrons
69 set MuonOutputArray muons
70
71 # radius of the magnetic field coverage, in m
72 set Radius 1.29
73 # half-length of the magnetic field coverage, in m
74 set HalfLength 3.00
75
76 # magnetic field
77 set Bz 3.8
78}
79
80####################################
81# Charged hadron tracking efficiency
82####################################
83
84module Efficiency ChargedHadronTrackingEfficiency {
85 set InputArray ParticlePropagator/chargedHadrons
86 set OutputArray chargedHadrons
87
88 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
89
90 # tracking efficiency formula for charged hadrons
91 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
92 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
93 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) +
94 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) +
95 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) +
96 (abs(eta) > 2.5) * (0.00)}
97}
98
99##############################
100# Electron tracking efficiency
101##############################
102
103module Efficiency ElectronTrackingEfficiency {
104 set InputArray ParticlePropagator/electrons
105 set OutputArray electrons
106
107 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
108
109 # tracking efficiency formula for electrons
110 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
111 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) +
112 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) +
113 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) +
114 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) +
115 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) +
116 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) +
117 (abs(eta) > 2.5) * (0.00)}
118}
119
120##########################
121# Muon tracking efficiency
122##########################
123
124module Efficiency MuonTrackingEfficiency {
125 set InputArray ParticlePropagator/muons
126 set OutputArray muons
127
128 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
129
130 # tracking efficiency formula for muons
131 set EfficiencyFormula { (pt <= 0.1) * (0.00) +
132 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) +
133 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e3) * (0.99) +
134 (abs(eta) <= 1.5) * (pt > 1.0e3 ) * (0.99 * exp(0.5 - pt*5.0e-4)) +
135
136 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) +
137 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e3) * (0.98) +
138 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3) * (0.98 * exp(0.5 - pt*5.0e-4)) +
139 (abs(eta) > 2.5) * (0.00)}
140}
141
142########################################
143# Momentum resolution for charged tracks
144########################################
145
146module MomentumSmearing ChargedHadronMomentumSmearing {
147 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
148 set OutputArray chargedHadrons
149
150 # set ResolutionFormula {resolution formula as a function of eta and pt}
151
152 # resolution formula for charged hadrons
153 # based on arXiv:1405.6569
154 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
155 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
156 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
157}
158
159###################################
160# Momentum resolution for electrons
161###################################
162
163module MomentumSmearing ElectronMomentumSmearing {
164 set InputArray ElectronTrackingEfficiency/electrons
165 set OutputArray electrons
166
167 # set ResolutionFormula {resolution formula as a function of eta and energy}
168
169 # resolution formula for electrons
170 # based on arXiv:1405.6569
171 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
172 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
173 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^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 muons
187 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
188 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
189 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
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
206#############
207# ECAL
208#############
209
210module SimpleCalorimeter ECal {
211 set ParticleInputArray ParticlePropagator/stableParticles
212 set TrackInputArray TrackMerger/tracks
213
214 set TowerOutputArray ecalTowers
215 set EFlowTrackOutputArray eflowTracks
216 set EFlowTowerOutputArray eflowPhotons
217
218 set IsEcal true
219
220 set EnergyMin 0.5
221 set EnergySignificanceMin 2.0
222
223 set SmearTowerCenter true
224
225 set pi [expr {acos(-1)}]
226
227 # lists of the edges of each tower in eta and phi
228 # each list starts with the lower edge of the first tower
229 # the list ends with the higher edged of the last tower
230
231 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
232
233 set PhiBins {}
234 for {set i -180} {$i <= 180} {incr i} {
235 add PhiBins [expr {$i * $pi/180.0}]
236 }
237
238 # 0.02 unit in eta up to eta = 1.5 (barrel)
239 for {set i -85} {$i <= 86} {incr i} {
240 set eta [expr {$i * 0.0174}]
241 add EtaPhiBins $eta $PhiBins
242 }
243
244 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
245
246 set PhiBins {}
247 for {set i -180} {$i <= 180} {incr i} {
248 add PhiBins [expr {$i * $pi/180.0}]
249 }
250
251 # 0.02 unit in eta up to eta = 3
252 for {set i 1} {$i <= 84} {incr i} {
253 set eta [expr { -2.958 + $i * 0.0174}]
254 add EtaPhiBins $eta $PhiBins
255 }
256
257 for {set i 1} {$i <= 84} {incr i} {
258 set eta [expr { 1.4964 + $i * 0.0174}]
259 add EtaPhiBins $eta $PhiBins
260 }
261
262 # take present CMS granularity for HF
263
264 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
265 set PhiBins {}
266 for {set i -18} {$i <= 18} {incr i} {
267 add PhiBins [expr {$i * $pi/18.0}]
268 }
269
270 foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
271 add EtaPhiBins $eta $PhiBins
272 }
273
274
275 add EnergyFraction {0} {0.0}
276 # energy fractions for e, gamma and pi0
277 add EnergyFraction {11} {1.0}
278 add EnergyFraction {22} {1.0}
279 add EnergyFraction {111} {1.0}
280 # energy fractions for muon, neutrinos and neutralinos
281 add EnergyFraction {12} {0.0}
282 add EnergyFraction {13} {0.0}
283 add EnergyFraction {14} {0.0}
284 add EnergyFraction {16} {0.0}
285 add EnergyFraction {1000022} {0.0}
286 add EnergyFraction {1000023} {0.0}
287 add EnergyFraction {1000025} {0.0}
288 add EnergyFraction {1000035} {0.0}
289 add EnergyFraction {1000045} {0.0}
290 # energy fractions for K0short and Lambda
291 add EnergyFraction {310} {0.3}
292 add EnergyFraction {3122} {0.3}
293
294 # set ResolutionFormula {resolution formula as a function of eta and energy}
295
296 # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
297
298 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
299 # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
300 set ResolutionFormula { (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
301 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (2.16 + 5.6*(abs(eta)-2)^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
302 (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
303
304}
305
306
307#############
308# HCAL
309#############
310
311module SimpleCalorimeter HCal {
312 set ParticleInputArray ParticlePropagator/stableParticles
313 set TrackInputArray ECal/eflowTracks
314
315 set TowerOutputArray hcalTowers
316 set EFlowTrackOutputArray eflowTracks
317 set EFlowTowerOutputArray eflowNeutralHadrons
318
319 set IsEcal false
320
321 set EnergyMin 1.0
322 set EnergySignificanceMin 1.0
323
324 set SmearTowerCenter true
325
326 set pi [expr {acos(-1)}]
327
328 # lists of the edges of each tower in eta and phi
329 # each list starts with the lower edge of the first tower
330 # the list ends with the higher edged of the last tower
331
332 # 5 degrees towers
333 set PhiBins {}
334 for {set i -36} {$i <= 36} {incr i} {
335 add PhiBins [expr {$i * $pi/36.0}]
336 }
337 foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
338 add EtaPhiBins $eta $PhiBins
339 }
340
341 # 10 degrees towers
342 set PhiBins {}
343 for {set i -18} {$i <= 18} {incr i} {
344 add PhiBins [expr {$i * $pi/18.0}]
345 }
346 foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
347 add EtaPhiBins $eta $PhiBins
348 }
349
350 # 20 degrees towers
351 set PhiBins {}
352 for {set i -9} {$i <= 9} {incr i} {
353 add PhiBins [expr {$i * $pi/9.0}]
354 }
355 foreach eta {-5 -4.7 -4.525 4.7 5} {
356 add EtaPhiBins $eta $PhiBins
357 }
358
359 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
360 add EnergyFraction {0} {1.0}
361 # energy fractions for e, gamma and pi0
362 add EnergyFraction {11} {0.0}
363 add EnergyFraction {22} {0.0}
364 add EnergyFraction {111} {0.0}
365 # energy fractions for muon, neutrinos and neutralinos
366 add EnergyFraction {12} {0.0}
367 add EnergyFraction {13} {0.0}
368 add EnergyFraction {14} {0.0}
369 add EnergyFraction {16} {0.0}
370 add EnergyFraction {1000022} {0.0}
371 add EnergyFraction {1000023} {0.0}
372 add EnergyFraction {1000025} {0.0}
373 add EnergyFraction {1000035} {0.0}
374 add EnergyFraction {1000045} {0.0}
375 # energy fractions for K0short and Lambda
376 add EnergyFraction {310} {0.7}
377 add EnergyFraction {3122} {0.7}
378
379 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
380 set ResolutionFormula { (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
381 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
382
383}
384
385
386#################
387# Electron filter
388#################
389
390module PdgCodeFilter ElectronFilter {
391 set InputArray HCal/eflowTracks
392 set OutputArray electrons
393 set Invert true
394 add PdgCode {11}
395 add PdgCode {-11}
396}
397
398######################
399# ChargedHadronFilter
400######################
401
402module PdgCodeFilter ChargedHadronFilter {
403 set InputArray HCal/eflowTracks
404 set OutputArray chargedHadrons
405
406 add PdgCode {11}
407 add PdgCode {-11}
408 add PdgCode {13}
409 add PdgCode {-13}
410}
411
412
413###################################################
414# Tower Merger (in case not using e-flow algorithm)
415###################################################
416
417module Merger Calorimeter {
418# add InputArray InputArray
419 add InputArray ECal/ecalTowers
420 add InputArray HCal/hcalTowers
421 set OutputArray towers
422}
423
424
425
426####################
427# Energy flow merger
428####################
429
430module Merger EFlowMerger {
431# add InputArray InputArray
432 add InputArray HCal/eflowTracks
433 add InputArray ECal/eflowPhotons
434 add InputArray HCal/eflowNeutralHadrons
435 set OutputArray eflow
436}
437
438######################
439# EFlowFilter
440######################
441
442module PdgCodeFilter EFlowFilter {
443 set InputArray EFlowMerger/eflow
444 set OutputArray eflow
445
446 add PdgCode {11}
447 add PdgCode {-11}
448 add PdgCode {13}
449 add PdgCode {-13}
450}
451
452
453###################
454# Photon efficiency
455###################
456
457module Efficiency PhotonEfficiency {
458 set InputArray ECal/eflowPhotons
459 set OutputArray photons
460
461 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
462
463 # efficiency formula for photons
464 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
465 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
466 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
467 (abs(eta) > 2.5) * (0.00)}
468}
469
470##################
471# Photon isolation
472##################
473
474module Isolation PhotonIsolation {
475 set CandidateInputArray PhotonEfficiency/photons
476 set IsolationInputArray EFlowFilter/eflow
477
478 set OutputArray photons
479
480 set DeltaRMax 0.5
481
482 set PTMin 0.5
483
484 set PTRatioMax 0.12
485}
486
487
488#####################
489# Electron efficiency
490#####################
491
492module Efficiency ElectronEfficiency {
493 set InputArray ElectronFilter/electrons
494 set OutputArray electrons
495
496 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
497
498 # efficiency formula for electrons
499 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
500 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
501 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
502 (abs(eta) > 2.5) * (0.00)}
503}
504
505####################
506# Electron isolation
507####################
508
509module Isolation ElectronIsolation {
510 set CandidateInputArray ElectronEfficiency/electrons
511 set IsolationInputArray EFlowFilter/eflow
512
513 set OutputArray electrons
514
515 set DeltaRMax 0.5
516
517 set PTMin 0.5
518
519 set PTRatioMax 0.12
520}
521
522#################
523# Muon efficiency
524#################
525
526module Efficiency MuonEfficiency {
527 set InputArray MuonMomentumSmearing/muons
528 set OutputArray muons
529
530 # set EfficiencyFormula {efficiency as a function of eta and pt}
531
532 # efficiency formula for muons
533 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
534 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
535 (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0) * (0.95) +
536 (abs(eta) > 2.4) * (0.00)}
537}
538
539################
540# Muon isolation
541################
542
543module Isolation MuonIsolation {
544 set CandidateInputArray MuonEfficiency/muons
545 set IsolationInputArray EFlowFilter/eflow
546
547 set OutputArray muons
548
549 set DeltaRMax 0.5
550
551 set PTMin 0.5
552
553 set PTRatioMax 0.25
554}
555
556###################
557# Missing ET merger
558###################
559
560module Merger MissingET {
561# add InputArray InputArray
562 add InputArray EFlowMerger/eflow
563 set MomentumOutputArray momentum
564}
565
566##################
567# Scalar HT merger
568##################
569
570module Merger ScalarHT {
571# add InputArray InputArray
572 add InputArray UniqueObjectFinder/jets
573 add InputArray UniqueObjectFinder/electrons
574 add InputArray UniqueObjectFinder/photons
575 add InputArray UniqueObjectFinder/muons
576 set EnergyOutputArray energy
577}
578
579
580#####################
581# Neutrino Filter
582#####################
583
584module PdgCodeFilter NeutrinoFilter {
585
586 set InputArray Delphes/stableParticles
587 set OutputArray filteredParticles
588
589 set PTMin 0.0
590
591 add PdgCode {12}
592 add PdgCode {14}
593 add PdgCode {16}
594 add PdgCode {-12}
595 add PdgCode {-14}
596 add PdgCode {-16}
597
598}
599
600
601#####################
602# MC truth jet finder
603#####################
604
605module FastJetFinder GenJetFinder {
606 set InputArray NeutrinoFilter/filteredParticles
607
608 set OutputArray jets
609
610 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
611 set JetAlgorithm 6
612 set ParameterR 0.5
613
614 set JetPTMin 20.0
615}
616
617#########################
618# Gen Missing ET merger
619########################
620
621module Merger GenMissingET {
622# add InputArray InputArray
623 add InputArray NeutrinoFilter/filteredParticles
624 set MomentumOutputArray momentum
625}
626
627
628
629############
630# Jet finder
631############
632
633module FastJetFinder FastJetFinder {
634# set InputArray Calorimeter/towers
635 set InputArray EFlowMerger/eflow
636
637 set OutputArray jets
638
639 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
640 set JetAlgorithm 6
641 set ParameterR 0.5
642
643 set JetPTMin 20.0
644}
645
646##################
647# Jet Energy Scale
648##################
649
650module EnergyScale JetEnergyScale {
651 set InputArray FastJetFinder/jets
652 set OutputArray jets
653
654 # scale formula for jets
655 set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
656}
657
658########################
659# Jet Flavor Association
660########################
661
662module JetFlavorAssociation JetFlavorAssociation {
663
664 set PartonInputArray Delphes/partons
665 set ParticleInputArray Delphes/allParticles
666 set ParticleLHEFInputArray Delphes/allParticlesLHEF
667 set JetInputArray JetEnergyScale/jets
668
669 set DeltaR 0.5
670 set PartonPTMin 1.0
671 set PartonEtaMax 2.5
672
673}
674
675###########
676# b-tagging
677###########
678
679module BTagging BTagging {
680 set JetInputArray JetEnergyScale/jets
681
682 set BitNumber 0
683
684 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
685 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
686 # gluon's PDG code has the lowest priority
687
688 # based on arXiv:1211.4462
689
690 # default efficiency formula (misidentification rate)
691 add EfficiencyFormula {0} {0.01+0.000038*pt}
692
693 # efficiency formula for c-jets (misidentification rate)
694 add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
695
696 # efficiency formula for b-jets
697 add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
698}
699
700#############
701# tau-tagging
702#############
703
704module TauTagging TauTagging {
705 set ParticleInputArray Delphes/allParticles
706 set PartonInputArray Delphes/partons
707 set JetInputArray JetEnergyScale/jets
708
709 set DeltaR 0.5
710
711 set TauPTMin 1.0
712
713 set TauEtaMax 2.5
714
715 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
716
717 # default efficiency formula (misidentification rate)
718 add EfficiencyFormula {0} {0.01}
719 # efficiency formula for tau-jets
720 add EfficiencyFormula {15} {0.6}
721}
722
723#####################################################
724# Find uniquely identified photons/electrons/tau/jets
725#####################################################
726
727module UniqueObjectFinder UniqueObjectFinder {
728# earlier arrays take precedence over later ones
729# add InputArray InputArray OutputArray
730 add InputArray PhotonIsolation/photons photons
731 add InputArray ElectronIsolation/electrons electrons
732 add InputArray MuonIsolation/muons muons
733 add InputArray JetEnergyScale/jets jets
734}
735
736##################
737# ROOT tree writer
738##################
739
740# tracks, towers and eflow objects are not stored by default in the output.
741# if needed (for jet constituent or other studies), uncomment the relevant
742# "add Branch ..." lines.
743
744module TreeWriter TreeWriter {
745# add Branch InputArray BranchName BranchClass
746 add Branch Delphes/allParticles Particle GenParticle
747
748 add Branch TrackMerger/tracks Track Track
749 add Branch Calorimeter/towers Tower Tower
750
751 add Branch HCal/eflowTracks EFlowTrack Track
752 add Branch ECal/eflowPhotons EFlowPhoton Tower
753 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
754
755 add Branch GenJetFinder/jets GenJet Jet
756 add Branch GenMissingET/momentum GenMissingET MissingET
757
758 add Branch UniqueObjectFinder/jets Jet Jet
759 add Branch UniqueObjectFinder/electrons Electron Electron
760 add Branch UniqueObjectFinder/photons Photon Photon
761 add Branch UniqueObjectFinder/muons Muon Muon
762 add Branch MissingET/momentum MissingET MissingET
763 add Branch ScalarHT/energy ScalarHT ScalarHT
764}
Note: See TracBrowser for help on using the repository browser.