Fork me on GitHub

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

Last change on this file since 1357 was 1357, checked in by Michele Selvaggi, 10 years ago

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