Fork me on GitHub

source: svn/trunk/examples/delphes_card_ATLAS.tcl

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

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