Fork me on GitHub

source: git/cards/delphes_card_CMS.tcl@ 93b4a8c

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

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