Fork me on GitHub

source: git/examples/delphes_card_ATLAS.tcl@ fcc959c

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

corrected constant term in fwd ATLAS calorimeter (was 10x too high)
  • Property mode set to 100644
File size: 18.6 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 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
31 NeutrinoFilter
32 GenJetFinder
33 FastJetFinder
34
35 JetEnergyScale
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
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) + \
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
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) + \
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
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) + \
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
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) + \
145 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2) * (0.05)}
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
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) + \
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
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) + \
183 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2) * (0.10)}
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
194 add InputArray MuonMomentumSmearing/muons
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
210 set EFlowPhotonOutputArray eflowPhotons
211 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
212
213 set pi [expr {acos(-1)}]
214
215 # lists of the edges of each tower in eta and phi
216 # each list starts with the lower edge of the first tower
217 # the list ends with the higher edged of the last tower
218
219 # 10 degrees towers
220 set PhiBins {}
221 for {set i -18} {$i <= 18} {incr i} {
222 add PhiBins [expr {$i * $pi/18.0}]
223 }
224 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} {
225 add EtaPhiBins $eta $PhiBins
226 }
227
228 # 20 degrees towers
229 set PhiBins {}
230 for {set i -9} {$i <= 9} {incr i} {
231 add PhiBins [expr {$i * $pi/9.0}]
232 }
233 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} {
234 add EtaPhiBins $eta $PhiBins
235 }
236
237 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
238 add EnergyFraction {0} {0.0 1.0}
239 # energy fractions for e, gamma and pi0
240 add EnergyFraction {11} {1.0 0.0}
241 add EnergyFraction {22} {1.0 0.0}
242 add EnergyFraction {111} {1.0 0.0}
243 # energy fractions for muon, neutrinos and neutralinos
244 add EnergyFraction {12} {0.0 0.0}
245 add EnergyFraction {13} {0.0 0.0}
246 add EnergyFraction {14} {0.0 0.0}
247 add EnergyFraction {16} {0.0 0.0}
248 add EnergyFraction {1000022} {0.0 0.0}
249 add EnergyFraction {1000023} {0.0 0.0}
250 add EnergyFraction {1000025} {0.0 0.0}
251 add EnergyFraction {1000035} {0.0 0.0}
252 add EnergyFraction {1000045} {0.0 0.0}
253 # energy fractions for K0short and Lambda
254 add EnergyFraction {310} {0.3 0.7}
255 add EnergyFraction {3122} {0.3 0.7}
256
257 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
258 # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
259 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
260 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
261 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + \
262 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
263
264 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
265 # http://arxiv.org/pdf/hep-ex/0004009v1
266 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
267 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
268 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \
269 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*0.075^2)}
270}
271
272####################
273# Energy flow merger
274####################
275
276module Merger EFlowMerger {
277# add InputArray InputArray
278 add InputArray Calorimeter/eflowTracks
279 add InputArray Calorimeter/eflowPhotons
280 add InputArray Calorimeter/eflowNeutralHadrons
281 set OutputArray eflow
282}
283
284###################
285# Photon efficiency
286###################
287
288module Efficiency PhotonEfficiency {
289 set InputArray Calorimeter/photons
290 set OutputArray photons
291
292 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
293
294 # efficiency formula for photons
295 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
296 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
297 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
298 (abs(eta) > 2.5) * (0.00)}
299}
300
301##################
302# Photon isolation
303##################
304
305module Isolation PhotonIsolation {
306 set CandidateInputArray PhotonEfficiency/photons
307 set IsolationInputArray EFlowMerger/eflow
308
309 set OutputArray photons
310
311 set DeltaRMax 0.5
312
313 set PTMin 0.5
314
315 set PTRatioMax 0.1
316}
317
318#####################
319# Electron efficiency
320#####################
321
322module Efficiency ElectronEfficiency {
323 set InputArray ElectronEnergySmearing/electrons
324 set OutputArray electrons
325
326 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
327
328 # efficiency formula for electrons
329 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
330 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
331 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
332 (abs(eta) > 2.5) * (0.00)}
333}
334
335####################
336# Electron isolation
337####################
338
339module Isolation ElectronIsolation {
340 set CandidateInputArray ElectronEfficiency/electrons
341 set IsolationInputArray EFlowMerger/eflow
342
343 set OutputArray electrons
344
345 set DeltaRMax 0.5
346
347 set PTMin 0.5
348
349 set PTRatioMax 0.1
350}
351
352#################
353# Muon efficiency
354#################
355
356module Efficiency MuonEfficiency {
357 set InputArray MuonMomentumSmearing/muons
358 set OutputArray muons
359
360 # set EfficiencyFormula {efficiency as a function of eta and pt}
361
362 # efficiency formula for muons
363 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
364 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
365 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) + \
366 (abs(eta) > 2.7) * (0.00)}
367}
368
369################
370# Muon isolation
371################
372
373module Isolation MuonIsolation {
374 set CandidateInputArray MuonEfficiency/muons
375 set IsolationInputArray EFlowMerger/eflow
376
377 set OutputArray muons
378
379 set DeltaRMax 0.5
380
381 set PTMin 0.5
382
383 set PTRatioMax 0.1
384}
385
386###################
387# Missing ET merger
388###################
389
390module Merger MissingET {
391# add InputArray InputArray
392 add InputArray EFlowMerger/eflow
393 set MomentumOutputArray momentum
394}
395
396##################
397# Scalar HT merger
398##################
399
400module Merger ScalarHT {
401# add InputArray InputArray
402 add InputArray UniqueObjectFinder/jets
403 add InputArray UniqueObjectFinder/electrons
404 add InputArray UniqueObjectFinder/photons
405 add InputArray UniqueObjectFinder/muons
406 set EnergyOutputArray energy
407}
408
409
410#####################
411# Neutrino Filter
412#####################
413
414module PdgCodeFilter NeutrinoFilter {
415
416 set InputArray Delphes/stableParticles
417 set OutputArray filteredParticles
418
419 set PTMin 0.0
420
421 add PdgCode {12}
422 add PdgCode {14}
423 add PdgCode {16}
424 add PdgCode {-12}
425 add PdgCode {-14}
426 add PdgCode {-16}
427
428}
429
430#####################
431# MC truth jet finder
432#####################
433
434module FastJetFinder GenJetFinder {
435 set InputArray NeutrinoFilter/filteredParticles
436
437 set OutputArray jets
438
439 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
440 set JetAlgorithm 6
441 set ParameterR 0.6
442
443 set JetPTMin 20.0
444}
445
446
447############
448# Jet finder
449############
450
451module FastJetFinder FastJetFinder {
452 set InputArray Calorimeter/towers
453
454 set OutputArray jets
455
456 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
457 set JetAlgorithm 6
458 set ParameterR 0.6
459
460 set JetPTMin 20.0
461}
462
463##################
464# Jet Energy Scale
465##################
466
467module EnergyScale JetEnergyScale {
468 set InputArray FastJetFinder/jets
469 set OutputArray jets
470
471 # scale formula for jets
472 set ScaleFormula {1.00}
473}
474
475###########
476# b-tagging
477###########
478
479module BTagging BTagging {
480 set PartonInputArray Delphes/partons
481 set JetInputArray JetEnergyScale/jets
482
483 set BitNumber 0
484
485 set DeltaR 0.5
486
487 set PartonPTMin 1.0
488
489 set PartonEtaMax 2.5
490
491 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
492 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
493 # gluon's PDG code has the lowest priority
494
495 # default efficiency formula (misidentification rate)
496 add EfficiencyFormula {0} {0.001}
497
498 # efficiency formula for c-jets (misidentification rate)
499 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
500 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
501 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
502 (abs(eta) > 2.5) * (0.000)}
503
504 # efficiency formula for b-jets
505 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \
506 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
507 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
508 (abs(eta) > 2.5) * (0.000)}
509}
510
511module TauTagging TauTagging {
512 set ParticleInputArray Delphes/allParticles
513 set PartonInputArray Delphes/partons
514 set JetInputArray JetEnergyScale/jets
515
516 set DeltaR 0.5
517
518 set TauPTMin 1.0
519
520 set TauEtaMax 2.5
521
522 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
523
524 # default efficiency formula (misidentification rate)
525 add EfficiencyFormula {0} {0.001}
526 # efficiency formula for tau-jets
527 add EfficiencyFormula {15} {0.4}
528}
529
530#####################################################
531# Find uniquely identified photons/electrons/tau/jets
532#####################################################
533
534module UniqueObjectFinder UniqueObjectFinder {
535# earlier arrays take precedence over later ones
536# add InputArray InputArray OutputArray
537 add InputArray PhotonIsolation/photons photons
538 add InputArray ElectronIsolation/electrons electrons
539 add InputArray MuonIsolation/muons muons
540 add InputArray JetEnergyScale/jets jets
541}
542
543##################
544# ROOT tree writer
545##################
546
547# tracks, towers and eflow objects are not stored by default in the output.
548# if needed (for jet constituent or other studies), uncomment the relevant
549# "add Branch ..." lines.
550
551module TreeWriter TreeWriter {
552# add Branch InputArray BranchName BranchClass
553 add Branch Delphes/allParticles Particle GenParticle
554
555 add Branch TrackMerger/tracks Track Track
556 add Branch Calorimeter/towers Tower Tower
557
558 add Branch Calorimeter/eflowTracks EFlowTrack Track
559 add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
560 add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
561
562 add Branch GenJetFinder/jets GenJet Jet
563 add Branch UniqueObjectFinder/jets Jet Jet
564 add Branch UniqueObjectFinder/electrons Electron Electron
565 add Branch UniqueObjectFinder/photons Photon Photon
566 add Branch UniqueObjectFinder/muons Muon Muon
567 add Branch MissingET/momentum MissingET MissingET
568 add Branch ScalarHT/energy ScalarHT ScalarHT
569}
570
Note: See TracBrowser for help on using the repository browser.