Fork me on GitHub

source: svn/trunk/examples/delphes_card_ATLAS.tcl@ 1323

Last change on this file since 1323 was 1322, checked in by Pavel Demin, 11 years ago

fix energy flow for muons in ATLAS cards
File size: 18.0 KB
RevLine 
[881]1#######################################
2# Order of execution of various modules
3#######################################
4
5set ExecutionPath {
6 ParticlePropagator
[941]7
8 ChargedHadronTrackingEfficiency
9 ElectronTrackingEfficiency
10 MuonTrackingEfficiency
11
12 ChargedHadronMomentumSmearing
[881]13 ElectronEnergySmearing
[941]14 MuonMomentumSmearing
15
16 TrackMerger
[881]17 Calorimeter
[941]18 EFlowMerger
19
[881]20 PhotonEfficiency
21 PhotonIsolation
[1032]22
[881]23 ElectronEfficiency
24 ElectronIsolation
[1032]25
[881]26 MuonEfficiency
27 MuonIsolation
[1032]28
[881]29 MissingET
[1032]30
[1076]31 GenJetFinder
[881]32 FastJetFinder
[1076]33
[1258]34 JetEnergyScale
[1076]35
[881]36 BTagging
37 TauTagging
[1032]38
[881]39 UniqueObjectFinder
[1032]40
[894]41 ScalarHT
[1032]42
[881]43 TreeWriter
44}
45
46#################################
47# Propagate particles in cylinder
48#################################
49
50module ParticlePropagator ParticlePropagator {
[941]51 set InputArray Delphes/stableParticles
[881]52
[941]53 set OutputArray stableParticles
54 set ChargedHadronOutputArray chargedHadrons
55 set ElectronOutputArray electrons
56 set MuonOutputArray muons
[881]57
58 # radius of the magnetic field coverage, in m
59 set Radius 1.15
60 # half-length of the magnetic field coverage, in m
61 set HalfLength 3.51
62
63 # magnetic field
64 set Bz 2.0
65}
66
[941]67####################################
68# Charged hadron tracking efficiency
69####################################
[881]70
[941]71module Efficiency ChargedHadronTrackingEfficiency {
72 set InputArray ParticlePropagator/chargedHadrons
73 set OutputArray chargedHadrons
[881]74
[941]75 # add EfficiencyFormula {efficiency formula as a function of eta and pt}
[881]76
[941]77 # tracking efficiency formula for charged hadrons
78 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
79 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
80 (abs(eta) <= 1.5) * (pt > 1.0) * (0.95) + \
81 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.60) + \
82 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.85) + \
83 (abs(eta) > 2.5) * (0.00)}
84}
[881]85
[941]86##############################
87# Electron tracking efficiency
88##############################
[881]89
[941]90module Efficiency ElectronTrackingEfficiency {
91 set InputArray ParticlePropagator/electrons
92 set OutputArray electrons
93
94 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
[1032]95
[941]96 # tracking efficiency formula for electrons
97 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
98 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.73) + \
99 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e2) * (0.95) + \
100 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.99) + \
101 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.50) + \
102 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e2) * (0.83) + \
103 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.90) + \
104 (abs(eta) > 2.5) * (0.00)}
[881]105}
106
[1032]107##########################
[941]108# Muon tracking efficiency
[1032]109##########################
[881]110
[941]111module Efficiency MuonTrackingEfficiency {
112 set InputArray ParticlePropagator/muons
113 set OutputArray muons
[881]114
[941]115 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
[1032]116
[941]117 # tracking efficiency formula for muons
118 set EfficiencyFormula { (pt <= 0.1) * (0.00) + \
119 (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.75) + \
120 (abs(eta) <= 1.5) * (pt > 1.0) * (0.99) + \
121 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.70) + \
122 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0) * (0.98) + \
123 (abs(eta) > 2.5) * (0.00)}
124}
[881]125
[941]126########################################
127# Momentum resolution for charged tracks
128########################################
[881]129
[941]130module MomentumSmearing ChargedHadronMomentumSmearing {
131 set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
132 set OutputArray chargedHadrons
133
134 # set ResolutionFormula {resolution formula as a function of eta and pt}
135
136 # resolution formula for charged hadrons
[1279]137 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.02) + \
138 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 1.0e1) * (0.01) + \
139 (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
140 (abs(eta) <= 1.5) * (pt > 2.0e2) * (0.05) + \
141 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
142 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 1.0e1) * (0.02) + \
143 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
144 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
[881]145}
146
[941]147#################################
148# Energy resolution for electrons
149#################################
[881]150
151module EnergySmearing ElectronEnergySmearing {
[941]152 set InputArray ElectronTrackingEfficiency/electrons
[881]153 set OutputArray electrons
154
[941]155 # set ResolutionFormula {resolution formula as a function of eta and energy}
[881]156
[941]157 set ResolutionFormula { (abs(eta) <= 2.5) * (energy > 0.1 && energy <= 2.5e1) * (energy*0.015) + \
158 (abs(eta) <= 2.5) * (energy > 2.5e1) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
159 (abs(eta) > 2.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
160 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
[881]161
162}
163
[941]164###############################
165# Momentum resolution for muons
166###############################
167
168module MomentumSmearing MuonMomentumSmearing {
169 set InputArray MuonTrackingEfficiency/muons
170 set OutputArray muons
171
172 # set ResolutionFormula {resolution formula as a function of eta and pt}
173
174 # resolution formula for muons
175 set ResolutionFormula { (abs(eta) <= 1.5) * (pt > 0.1 && pt <= 1.0) * (0.03) + \
[1258]176 (abs(eta) <= 1.5) * (pt > 1.0 && pt <= 5.0e1) * (0.03) + \
177 (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) + \
178 (abs(eta) <= 1.5) * (pt > 1.0e2) * (0.07) + \
[941]179 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1 && pt <= 1.0) * (0.04) + \
[1258]180 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0 && pt <= 5.0e1) * (0.04) + \
181 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) + \
182 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.10)}
[941]183}
184
185##############
186# Track merger
187##############
188
189module Merger TrackMerger {
190# add InputArray InputArray
191 add InputArray ChargedHadronMomentumSmearing/chargedHadrons
192 add InputArray ElectronEnergySmearing/electrons
[1322]193 add InputArray MuonMomentumSmearing/muons
[941]194 set OutputArray tracks
195}
196
[881]197#############
198# Calorimeter
199#############
200
201module Calorimeter Calorimeter {
[941]202 set ParticleInputArray ParticlePropagator/stableParticles
203 set TrackInputArray TrackMerger/tracks
[881]204
205 set TowerOutputArray towers
206 set PhotonOutputArray photons
207
[941]208 set EFlowTrackOutputArray eflowTracks
209 set EFlowTowerOutputArray eflowTowers
210
[881]211 set pi [expr {acos(-1)}]
[1032]212
[883]213 # lists of the edges of each tower in eta and phi
214 # each list starts with the lower edge of the first tower
215 # the list ends with the higher edged of the last tower
216
[881]217 # 10 degrees towers
[883]218 set PhiBins {}
[1032]219 for {set i -18} {$i <= 18} {incr i} {
[883]220 add PhiBins [expr {$i * $pi/18.0}]
[881]221 }
[1032]222 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} {
[883]223 add EtaPhiBins $eta $PhiBins
[881]224 }
[1032]225
[881]226 # 20 degrees towers
[883]227 set PhiBins {}
[1032]228 for {set i -9} {$i <= 9} {incr i} {
[883]229 add PhiBins [expr {$i * $pi/9.0}]
[1032]230 }
[1083]231 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} {
[883]232 add EtaPhiBins $eta $PhiBins
[881]233 }
234
235 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
236 add EnergyFraction {0} {0.0 1.0}
237 # energy fractions for e, gamma and pi0
238 add EnergyFraction {11} {1.0 0.0}
239 add EnergyFraction {22} {1.0 0.0}
240 add EnergyFraction {111} {1.0 0.0}
[1076]241 # energy fractions for muon, neutrinos and neutralinos
[881]242 add EnergyFraction {12} {0.0 0.0}
243 add EnergyFraction {13} {0.0 0.0}
244 add EnergyFraction {14} {0.0 0.0}
245 add EnergyFraction {16} {0.0 0.0}
[1076]246 add EnergyFraction {1000022} {0.0 0.0}
247 add EnergyFraction {1000023} {0.0 0.0}
248 add EnergyFraction {1000025} {0.0 0.0}
249 add EnergyFraction {1000035} {0.0 0.0}
250 add EnergyFraction {1000045} {0.0 0.0}
[881]251 # energy fractions for K0short and Lambda
252 add EnergyFraction {310} {0.3 0.7}
253 add EnergyFraction {3122} {0.3 0.7}
[1032]254
[881]255 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
256 # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
[1064]257 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
[1063]258 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
[881]259 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + \
[1063]260 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
[881]261
262 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
263 # http://arxiv.org/pdf/hep-ex/0004009v1
264 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
265 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
266 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \
267 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.9420^2 + energy*0.075^2)}
268}
269
[941]270####################
271# Energy flow merger
272####################
273
274module Merger EFlowMerger {
275# add InputArray InputArray
276 add InputArray Calorimeter/eflowTracks
277 add InputArray Calorimeter/eflowTowers
278 set OutputArray eflow
279}
280
[881]281###################
282# Photon efficiency
283###################
284
285module Efficiency PhotonEfficiency {
286 set InputArray Calorimeter/photons
287 set OutputArray photons
288
[941]289 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
[881]290
[941]291 # efficiency formula for photons
292 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
293 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
294 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
295 (abs(eta) > 2.5) * (0.00)}
[881]296}
297
298##################
299# Photon isolation
300##################
301
302module Isolation PhotonIsolation {
303 set CandidateInputArray PhotonEfficiency/photons
[1259]304 set IsolationInputArray EFlowMerger/eflow
[1032]305
[881]306 set OutputArray photons
307
308 set DeltaRMax 0.5
309
310 set PTMin 0.5
311
312 set PTRatioMax 0.1
313}
314
315#####################
316# Electron efficiency
317#####################
318
319module Efficiency ElectronEfficiency {
320 set InputArray ElectronEnergySmearing/electrons
321 set OutputArray electrons
322
[941]323 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
[881]324
325 # efficiency formula for electrons
[941]326 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
327 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
328 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
329 (abs(eta) > 2.5) * (0.00)}
[881]330}
331
332####################
333# Electron isolation
334####################
335
336module Isolation ElectronIsolation {
337 set CandidateInputArray ElectronEfficiency/electrons
[1259]338 set IsolationInputArray EFlowMerger/eflow
[1032]339
[881]340 set OutputArray electrons
341
342 set DeltaRMax 0.5
343
344 set PTMin 0.5
345
[894]346 set PTRatioMax 0.1
[881]347}
348
349#################
350# Muon efficiency
351#################
352
353module Efficiency MuonEfficiency {
[941]354 set InputArray MuonMomentumSmearing/muons
[881]355 set OutputArray muons
356
[941]357 # set EfficiencyFormula {efficiency as a function of eta and pt}
[881]358
359 # efficiency formula for muons
[941]360 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
361 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
362 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) + \
363 (abs(eta) > 2.7) * (0.00)}
[881]364}
365
366################
367# Muon isolation
368################
369
370module Isolation MuonIsolation {
371 set CandidateInputArray MuonEfficiency/muons
[1259]372 set IsolationInputArray EFlowMerger/eflow
[1032]373
[881]374 set OutputArray muons
375
376 set DeltaRMax 0.5
377
378 set PTMin 0.5
379
[894]380 set PTRatioMax 0.1
[881]381}
382
383###################
384# Missing ET merger
385###################
386
387module Merger MissingET {
[941]388# add InputArray InputArray
[1322]389 add InputArray Calorimeter/eflowTracks
390 add InputArray Calorimeter/eflowTowers
[881]391 set MomentumOutputArray momentum
392}
393
[894]394##################
395# Scalar HT merger
396##################
397
398module Merger ScalarHT {
399# add InputArray InputArray
400 add InputArray UniqueObjectFinder/jets
401 add InputArray UniqueObjectFinder/electrons
402 add InputArray UniqueObjectFinder/photons
[1258]403 add InputArray UniqueObjectFinder/muons
[894]404 set EnergyOutputArray energy
405}
406
[1076]407#####################
408# MC truth jet finder
409#####################
410
411module FastJetFinder GenJetFinder {
412 set InputArray Delphes/stableParticles
413
414 set OutputArray jets
415
416 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
417 set JetAlgorithm 6
[1258]418 set ParameterR 0.6
[1076]419
420 set JetPTMin 20.0
421}
422
[881]423############
424# Jet finder
425############
426
427module FastJetFinder FastJetFinder {
[1258]428 set InputArray Calorimeter/towers
[1032]429
[881]430 set OutputArray jets
431
432 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
433 set JetAlgorithm 6
[1258]434 set ParameterR 0.6
[881]435
436 set JetPTMin 20.0
437}
438
[1258]439##################
440# Jet Energy Scale
441##################
[1076]442
[1258]443module EnergyScale JetEnergyScale {
444 set InputArray FastJetFinder/jets
445 set OutputArray jets
[1076]446
[1258]447 # scale formula for jets
448 set ScaleFormula {1.08}
[1076]449}
450
[881]451###########
452# b-tagging
453###########
454
455module BTagging BTagging {
456 set PartonInputArray Delphes/partons
[1259]457 set JetInputArray JetEnergyScale/jets
[881]458
[1099]459 set BitNumber 0
460
[881]461 set DeltaR 0.5
462
463 set PartonPTMin 1.0
464
465 set PartonEtaMax 2.5
466
467 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
468 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
469 # gluon's PDG code has the lowest priority
470
471 # default efficiency formula (misidentification rate)
472 add EfficiencyFormula {0} {0.001}
[1258]473
[881]474 # efficiency formula for c-jets (misidentification rate)
[1258]475 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
476 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
477 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
478 (abs(eta) > 2.5) * (0.000)}
479
[881]480 # efficiency formula for b-jets
[1258]481 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \
482 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
483 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
484 (abs(eta) > 2.5) * (0.000)}
[881]485}
486
487module TauTagging TauTagging {
[941]488 set ParticleInputArray Delphes/allParticles
[881]489 set PartonInputArray Delphes/partons
[1259]490 set JetInputArray JetEnergyScale/jets
[881]491
492 set DeltaR 0.5
493
494 set TauPTMin 1.0
495
496 set TauEtaMax 2.5
497
498 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
499
500 # default efficiency formula (misidentification rate)
501 add EfficiencyFormula {0} {0.001}
502 # efficiency formula for tau-jets
503 add EfficiencyFormula {15} {0.4}
504}
505
506#####################################################
507# Find uniquely identified photons/electrons/tau/jets
508#####################################################
509
510module UniqueObjectFinder UniqueObjectFinder {
511# earlier arrays take precedence over later ones
512# add InputArray InputArray OutputArray
513 add InputArray PhotonIsolation/photons photons
514 add InputArray ElectronIsolation/electrons electrons
[1258]515 add InputArray MuonIsolation/muons muons
[1259]516 add InputArray JetEnergyScale/jets jets
[881]517}
518
519##################
520# ROOT tree writer
521##################
522
523module TreeWriter TreeWriter {
524# add Branch InputArray BranchName BranchClass
[941]525 add Branch Delphes/allParticles Particle GenParticle
[951]526 add Branch TrackMerger/tracks Track Track
527 add Branch Calorimeter/towers Tower Tower
[1258]528# add Branch Calorimeter/eflowTracks EFlowTrack Track
529# add Branch Calorimeter/eflowTowers EFlowTower Tower
530# add Branch MuonMomentumSmearing/muons EFlowMuon Muon
[1076]531 add Branch GenJetFinder/jets GenJet Jet
[954]532 add Branch UniqueObjectFinder/jets Jet Jet
533 add Branch UniqueObjectFinder/electrons Electron Electron
[881]534 add Branch UniqueObjectFinder/photons Photon Photon
[1258]535 add Branch UniqueObjectFinder/muons Muon Muon
[881]536 add Branch MissingET/momentum MissingET MissingET
[894]537 add Branch ScalarHT/energy ScalarHT ScalarHT
[881]538}
[941]539
Note: See TracBrowser for help on using the repository browser.