Fork me on GitHub

source: git/cards/delphes_card_ATLAS.tcl@ 20dc57b

Last change on this file since 20dc57b was a2983ec, checked in by Pavel Demin <pavel.demin@…>, 10 years ago

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