Fork me on GitHub

source: git/cards/delphes_card_ATLAS.tcl@ 8d200a6

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 8d200a6 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
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 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
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
269 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) +
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
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) +
277 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
278}
279
280####################
281# Energy flow merger
282####################
283
284module Merger EFlowMerger {
285# add InputArray InputArray
286 add InputArray Calorimeter/eflowTracks
287 add InputArray Calorimeter/eflowPhotons
288 add InputArray Calorimeter/eflowNeutralHadrons
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
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) +
306 (abs(eta) > 2.5) * (0.00)}
307}
308
309##################
310# Photon isolation
311##################
312
313module Isolation PhotonIsolation {
314 set CandidateInputArray PhotonEfficiency/photons
315 set IsolationInputArray EFlowMerger/eflow
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
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) +
340 (abs(eta) > 2.5) * (0.00)}
341}
342
343####################
344# Electron isolation
345####################
346
347module Isolation ElectronIsolation {
348 set CandidateInputArray ElectronEfficiency/electrons
349 set IsolationInputArray EFlowMerger/eflow
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
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) +
374 (abs(eta) > 2.7) * (0.00)}
375}
376
377################
378# Muon isolation
379################
380
381module Isolation MuonIsolation {
382 set CandidateInputArray MuonEfficiency/muons
383 set IsolationInputArray EFlowMerger/eflow
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
400 add InputArray EFlowMerger/eflow
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
413 add InputArray UniqueObjectFinder/muons
414 set EnergyOutputArray energy
415}
416
417
418#####################
419# Neutrino Filter
420#####################
421
422module PdgCodeFilter NeutrinoFilter {
423
424 set InputArray Delphes/stableParticles
425 set OutputArray filteredParticles
426
427 set PTMin 0.0
428
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
438#####################
439# MC truth jet finder
440#####################
441
442module FastJetFinder GenJetFinder {
443 set InputArray NeutrinoFilter/filteredParticles
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
449 set ParameterR 0.6
450
451 set JetPTMin 20.0
452}
453
454
455############
456# Jet finder
457############
458
459module FastJetFinder FastJetFinder {
460 set InputArray Calorimeter/towers
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
466 set ParameterR 0.6
467
468 set JetPTMin 20.0
469}
470
471##################
472# Jet Energy Scale
473##################
474
475module EnergyScale JetEnergyScale {
476 set InputArray FastJetFinder/jets
477 set OutputArray jets
478
479 # scale formula for jets
480 set ScaleFormula { sqrt( (3.0 - 0.2*(abs(eta)))^2 / pt + 1.0 ) }
481}
482
483###########
484# b-tagging
485###########
486
487module BTagging BTagging {
488 set PartonInputArray Delphes/partons
489 set JetInputArray JetEnergyScale/jets
490
491 set BitNumber 0
492
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}
505
506 # efficiency formula for c-jets (misidentification rate)
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)) +
510 (abs(eta) > 2.5) * (0.000)}
511
512 # efficiency formula for b-jets
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)) +
516 (abs(eta) > 2.5) * (0.000)}
517}
518
519module TauTagging TauTagging {
520 set ParticleInputArray Delphes/allParticles
521 set PartonInputArray Delphes/partons
522 set JetInputArray JetEnergyScale/jets
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
547 add InputArray MuonIsolation/muons muons
548 add InputArray JetEnergyScale/jets jets
549}
550
551##################
552# ROOT tree writer
553##################
554
555# tracks, towers and eflow objects are not stored by default in the output.
556# if needed (for jet constituent or other studies), uncomment the relevant
557# "add Branch ..." lines.
558
559module TreeWriter TreeWriter {
560# add Branch InputArray BranchName BranchClass
561 add Branch Delphes/allParticles Particle GenParticle
562
563 add Branch TrackMerger/tracks Track Track
564 add Branch Calorimeter/towers Tower Tower
565
566 add Branch Calorimeter/eflowTracks EFlowTrack Track
567 add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
568 add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
569
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
574 add Branch UniqueObjectFinder/muons Muon Muon
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.