Fork me on GitHub

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

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

revert to pre 3.0.10 charged hadron resolution for ATLAS
File size: 17.9 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
193 set OutputArray tracks
194}
195
[881]196#############
197# Calorimeter
198#############
199
200module Calorimeter Calorimeter {
[941]201 set ParticleInputArray ParticlePropagator/stableParticles
202 set TrackInputArray TrackMerger/tracks
[881]203
204 set TowerOutputArray towers
205 set PhotonOutputArray photons
206
[941]207 set EFlowTrackOutputArray eflowTracks
208 set EFlowTowerOutputArray eflowTowers
209
[881]210 set pi [expr {acos(-1)}]
[1032]211
[883]212 # lists of the edges of each tower in eta and phi
213 # each list starts with the lower edge of the first tower
214 # the list ends with the higher edged of the last tower
215
[881]216 # 10 degrees towers
[883]217 set PhiBins {}
[1032]218 for {set i -18} {$i <= 18} {incr i} {
[883]219 add PhiBins [expr {$i * $pi/18.0}]
[881]220 }
[1032]221 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]222 add EtaPhiBins $eta $PhiBins
[881]223 }
[1032]224
[881]225 # 20 degrees towers
[883]226 set PhiBins {}
[1032]227 for {set i -9} {$i <= 9} {incr i} {
[883]228 add PhiBins [expr {$i * $pi/9.0}]
[1032]229 }
[1083]230 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]231 add EtaPhiBins $eta $PhiBins
[881]232 }
233
234 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
235 add EnergyFraction {0} {0.0 1.0}
236 # energy fractions for e, gamma and pi0
237 add EnergyFraction {11} {1.0 0.0}
238 add EnergyFraction {22} {1.0 0.0}
239 add EnergyFraction {111} {1.0 0.0}
[1076]240 # energy fractions for muon, neutrinos and neutralinos
[881]241 add EnergyFraction {12} {0.0 0.0}
242 add EnergyFraction {13} {0.0 0.0}
243 add EnergyFraction {14} {0.0 0.0}
244 add EnergyFraction {16} {0.0 0.0}
[1076]245 add EnergyFraction {1000022} {0.0 0.0}
246 add EnergyFraction {1000023} {0.0 0.0}
247 add EnergyFraction {1000025} {0.0 0.0}
248 add EnergyFraction {1000035} {0.0 0.0}
249 add EnergyFraction {1000045} {0.0 0.0}
[881]250 # energy fractions for K0short and Lambda
251 add EnergyFraction {310} {0.3 0.7}
252 add EnergyFraction {3122} {0.3 0.7}
[1032]253
[881]254 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
255 # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
[1064]256 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
[1063]257 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
[881]258 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + \
[1063]259 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
[881]260
261 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
262 # http://arxiv.org/pdf/hep-ex/0004009v1
263 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
264 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
265 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \
266 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.9420^2 + energy*0.075^2)}
267}
268
[941]269####################
270# Energy flow merger
271####################
272
273module Merger EFlowMerger {
274# add InputArray InputArray
275 add InputArray Calorimeter/eflowTracks
276 add InputArray Calorimeter/eflowTowers
277 add InputArray MuonMomentumSmearing/muons
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
389 add InputArray EFlowMerger/eflow
[881]390 set MomentumOutputArray momentum
391}
392
[894]393##################
394# Scalar HT merger
395##################
396
397module Merger ScalarHT {
398# add InputArray InputArray
399 add InputArray UniqueObjectFinder/jets
400 add InputArray UniqueObjectFinder/electrons
401 add InputArray UniqueObjectFinder/photons
[1258]402 add InputArray UniqueObjectFinder/muons
[894]403 set EnergyOutputArray energy
404}
405
[1076]406#####################
407# MC truth jet finder
408#####################
409
410module FastJetFinder GenJetFinder {
411 set InputArray Delphes/stableParticles
412
413 set OutputArray jets
414
415 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
416 set JetAlgorithm 6
[1258]417 set ParameterR 0.6
[1076]418
419 set JetPTMin 20.0
420}
421
[881]422############
423# Jet finder
424############
425
426module FastJetFinder FastJetFinder {
[1258]427 set InputArray Calorimeter/towers
[1032]428
[881]429 set OutputArray jets
430
431 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
432 set JetAlgorithm 6
[1258]433 set ParameterR 0.6
[881]434
435 set JetPTMin 20.0
436}
437
[1258]438##################
439# Jet Energy Scale
440##################
[1076]441
[1258]442module EnergyScale JetEnergyScale {
443 set InputArray FastJetFinder/jets
444 set OutputArray jets
[1076]445
[1258]446 # scale formula for jets
447 set ScaleFormula {1.08}
[1076]448}
449
[881]450###########
451# b-tagging
452###########
453
454module BTagging BTagging {
455 set PartonInputArray Delphes/partons
[1259]456 set JetInputArray JetEnergyScale/jets
[881]457
[1099]458 set BitNumber 0
459
[881]460 set DeltaR 0.5
461
462 set PartonPTMin 1.0
463
464 set PartonEtaMax 2.5
465
466 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
467 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
468 # gluon's PDG code has the lowest priority
469
470 # default efficiency formula (misidentification rate)
471 add EfficiencyFormula {0} {0.001}
[1258]472
[881]473 # efficiency formula for c-jets (misidentification rate)
[1258]474 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
475 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
476 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
477 (abs(eta) > 2.5) * (0.000)}
478
[881]479 # efficiency formula for b-jets
[1258]480 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \
481 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
482 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
483 (abs(eta) > 2.5) * (0.000)}
[881]484}
485
486module TauTagging TauTagging {
[941]487 set ParticleInputArray Delphes/allParticles
[881]488 set PartonInputArray Delphes/partons
[1259]489 set JetInputArray JetEnergyScale/jets
[881]490
491 set DeltaR 0.5
492
493 set TauPTMin 1.0
494
495 set TauEtaMax 2.5
496
497 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
498
499 # default efficiency formula (misidentification rate)
500 add EfficiencyFormula {0} {0.001}
501 # efficiency formula for tau-jets
502 add EfficiencyFormula {15} {0.4}
503}
504
505#####################################################
506# Find uniquely identified photons/electrons/tau/jets
507#####################################################
508
509module UniqueObjectFinder UniqueObjectFinder {
510# earlier arrays take precedence over later ones
511# add InputArray InputArray OutputArray
512 add InputArray PhotonIsolation/photons photons
513 add InputArray ElectronIsolation/electrons electrons
[1258]514 add InputArray MuonIsolation/muons muons
[1259]515 add InputArray JetEnergyScale/jets jets
[881]516}
517
518##################
519# ROOT tree writer
520##################
521
522module TreeWriter TreeWriter {
523# add Branch InputArray BranchName BranchClass
[941]524 add Branch Delphes/allParticles Particle GenParticle
[951]525 add Branch TrackMerger/tracks Track Track
526 add Branch Calorimeter/towers Tower Tower
[1258]527# add Branch Calorimeter/eflowTracks EFlowTrack Track
528# add Branch Calorimeter/eflowTowers EFlowTower Tower
529# add Branch MuonMomentumSmearing/muons EFlowMuon Muon
[1076]530 add Branch GenJetFinder/jets GenJet Jet
[954]531 add Branch UniqueObjectFinder/jets Jet Jet
532 add Branch UniqueObjectFinder/electrons Electron Electron
[881]533 add Branch UniqueObjectFinder/photons Photon Photon
[1258]534 add Branch UniqueObjectFinder/muons Muon Muon
[881]535 add Branch MissingET/momentum MissingET MissingET
[894]536 add Branch ScalarHT/energy ScalarHT ScalarHT
[881]537}
[941]538
Note: See TracBrowser for help on using the repository browser.