Fork me on GitHub

source: git/cards/delphes_card_ATLAS.tcl@ db316cd

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

added muons in tower collection #1118
  • Property mode set to 100644
File size: 21.0 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.15
73 # half-length of the magnetic field coverage, in m
74 set HalfLength 3.51
75
76 # magnetic field
77 set Bz 2.0
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) * (0.99) +
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) * (0.98) +
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 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
151 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
152 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
153}
154
155###################################
156# Momentum resolution for electrons
157###################################
158
159module MomentumSmearing ElectronMomentumSmearing {
160 set InputArray ElectronTrackingEfficiency/electrons
161 set OutputArray electrons
162
163 # set ResolutionFormula {resolution formula as a function of eta and energy}
164
165 # resolution formula for electrons
166 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
167 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
168 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
169}
170
171###############################
172# Momentum resolution for muons
173###############################
174
175module MomentumSmearing MuonMomentumSmearing {
176 set InputArray MuonTrackingEfficiency/muons
177 set OutputArray muons
178
179 # set ResolutionFormula {resolution formula as a function of eta and pt}
180 # resolution formula for muons
181 set ResolutionFormula { (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
182 (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
183 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
184}
185
186##############
187# Track merger
188##############
189
190module Merger TrackMerger {
191# add InputArray InputArray
192 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
193 add InputArray ElectronMomentumSmearing/electrons
194 add InputArray MuonMomentumSmearing/muons
195 set OutputArray tracks
196}
197
198
199#############
200# ECAL
201#############
202
203module SimpleCalorimeter ECal {
204 set ParticleInputArray ParticlePropagator/stableParticles
205 set TrackInputArray TrackMerger/tracks
206
207 set TowerOutputArray ecalTowers
208 set EFlowTrackOutputArray eflowTracks
209 set EFlowTowerOutputArray eflowPhotons
210
211 set IsEcal true
212
213 set EnergyMin 0.5
214 set EnergySignificanceMin 2.0
215
216 set SmearTowerCenter true
217
218 set pi [expr {acos(-1)}]
219
220 # lists of the edges of each tower in eta and phi
221 # each list starts with the lower edge of the first tower
222 # the list ends with the higher edged of the last tower
223
224 # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
225
226 set PhiBins {}
227 for {set i -180} {$i <= 180} {incr i} {
228 add PhiBins [expr {$i * $pi/180.0}]
229 }
230
231 # 0.02 unit in eta up to eta = 1.5 (barrel)
232 for {set i -85} {$i <= 86} {incr i} {
233 set eta [expr {$i * 0.0174}]
234 add EtaPhiBins $eta $PhiBins
235 }
236
237 # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0
238 set PhiBins {}
239 for {set i -180} {$i <= 180} {incr i} {
240 add PhiBins [expr {$i * $pi/180.0}]
241 }
242
243 # 0.02 unit in eta up to eta = 3
244 for {set i 1} {$i <= 84} {incr i} {
245 set eta [expr { -2.958 + $i * 0.0174}]
246 add EtaPhiBins $eta $PhiBins
247 }
248
249 for {set i 1} {$i <= 84} {incr i} {
250 set eta [expr { 1.4964 + $i * 0.0174}]
251 add EtaPhiBins $eta $PhiBins
252 }
253
254 # take present CMS granularity for HF
255
256 # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
257 set PhiBins {}
258 for {set i -18} {$i <= 18} {incr i} {
259 add PhiBins [expr {$i * $pi/18.0}]
260 }
261
262 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} {
263 add EtaPhiBins $eta $PhiBins
264 }
265
266
267 add EnergyFraction {0} {0.0}
268 # energy fractions for e, gamma and pi0
269 add EnergyFraction {11} {1.0}
270 add EnergyFraction {22} {1.0}
271 add EnergyFraction {111} {1.0}
272 # energy fractions for muon, neutrinos and neutralinos
273 add EnergyFraction {12} {0.0}
274 add EnergyFraction {13} {0.0}
275 add EnergyFraction {14} {0.0}
276 add EnergyFraction {16} {0.0}
277 add EnergyFraction {1000022} {0.0}
278 add EnergyFraction {1000023} {0.0}
279 add EnergyFraction {1000025} {0.0}
280 add EnergyFraction {1000035} {0.0}
281 add EnergyFraction {1000045} {0.0}
282 # energy fractions for K0short and Lambda
283 add EnergyFraction {310} {0.3}
284 add EnergyFraction {3122} {0.3}
285
286 # set ResolutionFormula {resolution formula as a function of eta and energy}
287
288 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
289 # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
290 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
291 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
292 set ResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) +
293 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
294
295
296}
297
298
299
300#############
301# HCAL
302#############
303
304module SimpleCalorimeter HCal {
305 set ParticleInputArray ParticlePropagator/stableParticles
306 set TrackInputArray ECal/eflowTracks
307
308 set TowerOutputArray hcalTowers
309 set EFlowTrackOutputArray eflowTracks
310 set EFlowTowerOutputArray eflowNeutralHadrons
311
312 set IsEcal false
313
314 set EnergyMin 1.0
315 set EnergySignificanceMin 2.0
316
317 set SmearTowerCenter true
318
319 set pi [expr {acos(-1)}]
320
321 # lists of the edges of each tower in eta and phi
322 # each list starts with the lower edge of the first tower
323 # the list ends with the higher edged of the last tower
324
325 # 10 degrees towers
326 set PhiBins {}
327 for {set i -18} {$i <= 18} {incr i} {
328 add PhiBins [expr {$i * $pi/18.0}]
329 }
330 foreach eta {-3.2 -2.5 -2.4 -2.3 -2.2 -2.1 -2 -1.9 -1.8 -1.7 -1.6 -1.5 -1.4 -1.3 -1.2 -1.1 -1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 2.6 3.3} {
331 add EtaPhiBins $eta $PhiBins
332 }
333
334 # 20 degrees towers
335 set PhiBins {}
336 for {set i -9} {$i <= 9} {incr i} {
337 add PhiBins [expr {$i * $pi/9.0}]
338 }
339 foreach eta {-4.9 -4.7 -4.5 -4.3 -4.1 -3.9 -3.7 -3.5 -3.3 -3 -2.8 -2.6 2.8 3 3.2 3.5 3.7 3.9 4.1 4.3 4.5 4.7 4.9} {
340 add EtaPhiBins $eta $PhiBins
341 }
342
343 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
344 add EnergyFraction {0} {1.0}
345 # energy fractions for e, gamma and pi0
346 add EnergyFraction {11} {0.0}
347 add EnergyFraction {22} {0.0}
348 add EnergyFraction {111} {0.0}
349 # energy fractions for muon, neutrinos and neutralinos
350 add EnergyFraction {12} {0.0}
351 add EnergyFraction {13} {0.0}
352 add EnergyFraction {14} {0.0}
353 add EnergyFraction {16} {0.0}
354 add EnergyFraction {1000022} {0.0}
355 add EnergyFraction {1000023} {0.0}
356 add EnergyFraction {1000025} {0.0}
357 add EnergyFraction {1000035} {0.0}
358 add EnergyFraction {1000045} {0.0}
359 # energy fractions for K0short and Lambda
360 add EnergyFraction {310} {0.7}
361 add EnergyFraction {3122} {0.7}
362
363 # http://arxiv.org/pdf/hep-ex/0004009v1
364 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
365 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
366 set ResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
367 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
368 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
369}
370
371
372#################
373# Electron filter
374#################
375
376module PdgCodeFilter ElectronFilter {
377 set InputArray HCal/eflowTracks
378 set OutputArray electrons
379 set Invert true
380 add PdgCode {11}
381 add PdgCode {-11}
382}
383
384######################
385# ChargedHadronFilter
386######################
387
388module PdgCodeFilter ChargedHadronFilter {
389 set InputArray HCal/eflowTracks
390 set OutputArray chargedHadrons
391
392 add PdgCode {11}
393 add PdgCode {-11}
394 add PdgCode {13}
395 add PdgCode {-13}
396}
397
398
399
400###################################################
401# Tower Merger (in case not using e-flow algorithm)
402###################################################
403
404module Merger Calorimeter {
405# add InputArray InputArray
406 add InputArray ECal/ecalTowers
407 add InputArray HCal/hcalTowers
408 add InputArray MuonMomentumSmearing/muons
409 set OutputArray towers
410}
411
412####################
413# Energy flow merger
414####################
415
416module Merger EFlowMerger {
417# add InputArray InputArray
418 add InputArray HCal/eflowTracks
419 add InputArray ECal/eflowPhotons
420 add InputArray HCal/eflowNeutralHadrons
421 set OutputArray eflow
422}
423
424######################
425# EFlowFilter
426######################
427
428module PdgCodeFilter EFlowFilter {
429 set InputArray EFlowMerger/eflow
430 set OutputArray eflow
431
432 add PdgCode {11}
433 add PdgCode {-11}
434 add PdgCode {13}
435 add PdgCode {-13}
436}
437
438###################
439# Photon efficiency
440###################
441
442module Efficiency PhotonEfficiency {
443 set InputArray ECal/eflowPhotons
444 set OutputArray photons
445
446 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
447
448 # efficiency formula for photons
449 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
450 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
451 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
452 (abs(eta) > 2.5) * (0.00)}
453}
454
455##################
456# Photon isolation
457##################
458
459module Isolation PhotonIsolation {
460 set CandidateInputArray PhotonEfficiency/photons
461 set IsolationInputArray EFlowFilter/eflow
462
463 set OutputArray photons
464
465 set DeltaRMax 0.5
466
467 set PTMin 0.5
468
469 set PTRatioMax 0.12
470}
471
472
473#####################
474# Electron efficiency
475#####################
476
477module Efficiency ElectronEfficiency {
478 set InputArray ElectronFilter/electrons
479 set OutputArray electrons
480
481 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
482
483 # efficiency formula for electrons
484 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
485 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
486 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) +
487 (abs(eta) > 2.5) * (0.00)}
488}
489
490####################
491# Electron isolation
492####################
493
494module Isolation ElectronIsolation {
495 set CandidateInputArray ElectronEfficiency/electrons
496 set IsolationInputArray EFlowFilter/eflow
497
498 set OutputArray electrons
499
500 set DeltaRMax 0.5
501
502 set PTMin 0.5
503
504 set PTRatioMax 0.12
505}
506
507#################
508# Muon efficiency
509#################
510
511module Efficiency MuonEfficiency {
512 set InputArray MuonMomentumSmearing/muons
513 set OutputArray muons
514
515 # set EfficiencyFormula {efficiency as a function of eta and pt}
516
517 # efficiency formula for muons
518 set EfficiencyFormula { (pt <= 10.0) * (0.00) +
519 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) +
520 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) +
521 (abs(eta) > 2.7) * (0.00)}
522}
523
524################
525# Muon isolation
526################
527
528module Isolation MuonIsolation {
529 set CandidateInputArray MuonEfficiency/muons
530 set IsolationInputArray EFlowFilter/eflow
531
532 set OutputArray muons
533
534 set DeltaRMax 0.5
535
536 set PTMin 0.5
537
538 set PTRatioMax 0.25
539}
540
541###################
542# Missing ET merger
543###################
544
545module Merger MissingET {
546# add InputArray InputArray
547 add InputArray Calorimeter/towers
548 set MomentumOutputArray momentum
549}
550
551##################
552# Scalar HT merger
553##################
554
555module Merger ScalarHT {
556# add InputArray InputArray
557 add InputArray UniqueObjectFinder/jets
558 add InputArray UniqueObjectFinder/electrons
559 add InputArray UniqueObjectFinder/photons
560 add InputArray UniqueObjectFinder/muons
561 set EnergyOutputArray energy
562}
563
564
565#####################
566# Neutrino Filter
567#####################
568
569module PdgCodeFilter NeutrinoFilter {
570
571 set InputArray Delphes/stableParticles
572 set OutputArray filteredParticles
573
574 set PTMin 0.0
575
576 add PdgCode {12}
577 add PdgCode {14}
578 add PdgCode {16}
579 add PdgCode {-12}
580 add PdgCode {-14}
581 add PdgCode {-16}
582
583}
584
585#####################
586# MC truth jet finder
587#####################
588
589module FastJetFinder GenJetFinder {
590 set InputArray NeutrinoFilter/filteredParticles
591
592 set OutputArray jets
593
594 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
595 set JetAlgorithm 6
596 set ParameterR 0.6
597
598 set JetPTMin 20.0
599}
600
601
602#########################
603# Gen Missing ET merger
604########################
605
606module Merger GenMissingET {
607# add InputArray InputArray
608 add InputArray NeutrinoFilter/filteredParticles
609 set MomentumOutputArray momentum
610}
611
612
613
614############
615# Jet finder
616############
617
618module FastJetFinder FastJetFinder {
619 set InputArray Calorimeter/towers
620
621 set OutputArray jets
622
623 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
624 set JetAlgorithm 6
625 set ParameterR 0.6
626
627 set JetPTMin 20.0
628}
629
630##################
631# Jet Energy Scale
632##################
633
634module EnergyScale JetEnergyScale {
635 set InputArray FastJetFinder/jets
636 set OutputArray jets
637
638 # scale formula for jets
639 set ScaleFormula { sqrt( (3.0 - 0.2*(abs(eta)))^2 / pt + 1.0 ) }
640}
641
642########################
643# Jet Flavor Association
644########################
645
646module JetFlavorAssociation JetFlavorAssociation {
647
648 set PartonInputArray Delphes/partons
649 set ParticleInputArray Delphes/allParticles
650 set ParticleLHEFInputArray Delphes/allParticlesLHEF
651 set JetInputArray JetEnergyScale/jets
652
653 set DeltaR 0.5
654 set PartonPTMin 1.0
655 set PartonEtaMax 2.5
656
657}
658
659###########
660# b-tagging
661###########
662
663module BTagging BTagging {
664 set JetInputArray JetEnergyScale/jets
665
666 set BitNumber 0
667
668 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
669 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
670 # gluon's PDG code has the lowest priority
671
672 # based on ATL-PHYS-PUB-2015-022
673
674 # default efficiency formula (misidentification rate)
675 add EfficiencyFormula {0} {0.002+7.3e-06*pt}
676
677 # efficiency formula for c-jets (misidentification rate)
678 add EfficiencyFormula {4} {0.20*tanh(0.02*pt)*(1/(1+0.0034*pt))}
679
680 # efficiency formula for b-jets
681 add EfficiencyFormula {5} {0.80*tanh(0.003*pt)*(30/(1+0.086*pt))}
682}
683
684#############
685# tau-tagging
686#############
687
688module TrackCountingTauTagging TauTagging {
689
690 set ParticleInputArray Delphes/allParticles
691 set PartonInputArray Delphes/partons
692 set TrackInputArray TrackMerger/tracks
693 set JetInputArray JetEnergyScale/jets
694
695 set DeltaR 0.2
696 set DeltaRTrack 0.2
697
698 set TrackPTMin 1.0
699
700 set TauPTMin 1.0
701 set TauEtaMax 2.5
702
703 # instructions: {n-prongs} {eff}
704
705 # 1 - one prong efficiency
706 # 2 - two or more efficiency
707 # -1 - one prong mistag rate
708 # -2 - two or more mistag rate
709
710 set BitNumber 0
711
712 # taken from ATL-PHYS-PUB-2015-045 (medium working point)
713 add EfficiencyFormula {1} {0.70}
714 add EfficiencyFormula {2} {0.60}
715 add EfficiencyFormula {-1} {0.02}
716 add EfficiencyFormula {-2} {0.01}
717
718}
719
720#####################################################
721# Find uniquely identified photons/electrons/tau/jets
722#####################################################
723
724module UniqueObjectFinder UniqueObjectFinder {
725# earlier arrays take precedence over later ones
726# add InputArray InputArray OutputArray
727 add InputArray PhotonIsolation/photons photons
728 add InputArray ElectronIsolation/electrons electrons
729 add InputArray MuonIsolation/muons muons
730 add InputArray JetEnergyScale/jets jets
731}
732
733##################
734# ROOT tree writer
735##################
736
737# tracks, towers and eflow objects are not stored by default in the output.
738# if needed (for jet constituent or other studies), uncomment the relevant
739# "add Branch ..." lines.
740
741module TreeWriter TreeWriter {
742# add Branch InputArray BranchName BranchClass
743 add Branch Delphes/allParticles Particle GenParticle
744
745 add Branch TrackMerger/tracks Track Track
746 add Branch Calorimeter/towers Tower Tower
747
748 add Branch HCal/eflowTracks EFlowTrack Track
749 add Branch ECal/eflowPhotons EFlowPhoton Tower
750 add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
751
752 add Branch GenJetFinder/jets GenJet Jet
753 add Branch GenMissingET/momentum GenMissingET MissingET
754
755 add Branch UniqueObjectFinder/jets Jet Jet
756 add Branch UniqueObjectFinder/electrons Electron Electron
757 add Branch UniqueObjectFinder/photons Photon Photon
758 add Branch UniqueObjectFinder/muons Muon Muon
759 add Branch MissingET/momentum MissingET MissingET
760 add Branch ScalarHT/energy ScalarHT ScalarHT
761}
762
Note: See TracBrowser for help on using the repository browser.