Fork me on GitHub

Ticket #215: delphes_card_ATLAS2.tcl

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