Fork me on GitHub

source: git/cards/delphes_card_CMS.tcl@ feada94

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

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