Fork me on GitHub

Ticket #207: delphes_card_ATLAS_edit.tcl

File delphes_card_ATLAS_edit.tcl, 18.5 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 GenBTagging
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 <= 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.2
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.2
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 GenBTagging {
482 set PartonInputArray Delphes/partons
483 set JetInputArray GenJetFinder/jets
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 # efficiency formula for c-jets (misidentification rate)
498 add EfficiencyFormula {4} {0.1}
499 # efficiency formula for b-jets
500 add EfficiencyFormula {5} {1.0}
501}
502
503module BTagging BTagging {
504 set PartonInputArray Delphes/partons
505 set JetInputArray FastJetFinder/jets
506
507 set DeltaR 0.5
508
509 set PartonPTMin 1.0
510
511 set PartonEtaMax 2.5
512
513 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
514 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
515 # gluon's PDG code has the lowest priority
516
517 # default efficiency formula (misidentification rate)
518 add EfficiencyFormula {0} {0.001}
519 # efficiency formula for c-jets (misidentification rate)
520 add EfficiencyFormula {4} {0.1}
521 # efficiency formula for b-jets
522 add EfficiencyFormula {5} {1.0}
523}
524
525module TauTagging TauTagging {
526 set ParticleInputArray Delphes/allParticles
527 set PartonInputArray Delphes/partons
528 set JetInputArray FastJetFinder/jets
529
530 set DeltaR 0.5
531
532 set TauPTMin 1.0
533
534 set TauEtaMax 2.5
535
536 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
537
538 # default efficiency formula (misidentification rate)
539 add EfficiencyFormula {0} {0.001}
540 # efficiency formula for tau-jets
541 add EfficiencyFormula {15} {0.4}
542}
543
544#####################################################
545# Find uniquely identified photons/electrons/tau/jets
546#####################################################
547
548module UniqueObjectFinder UniqueObjectFinder {
549# earlier arrays take precedence over later ones
550# add InputArray InputArray OutputArray
551 add InputArray PhotonIsolation/photons photons
552 add InputArray ElectronIsolation/electrons electrons
553 add InputArray FastJetFinder/jets jets
554}
555
556##################
557# ROOT tree writer
558##################
559
560module TreeWriter TreeWriter {
561# add Branch InputArray BranchName BranchClass
562 add Branch Delphes/partons Parton GenParticle
563 add Branch Delphes/allParticles Particle GenParticle
564 add Branch TrackMerger/tracks Track Track
565 add Branch Calorimeter/towers Tower Tower
566 add Branch ConstituentFilter/eflowTracks EFlowTrack Track
567 add Branch ConstituentFilter/eflowTowers EFlowTower Tower
568 add Branch ConstituentFilter/muons EFlowMuon Muon
569 add Branch GenJetFinder/jets GenJet Jet
570 add Branch UniqueObjectFinder/jets Jet Jet
571 add Branch UniqueObjectFinder/electrons Electron Electron
572 add Branch UniqueObjectFinder/photons Photon Photon
573 add Branch MuonIsolation/muons Muon Muon
574 add Branch MissingET/momentum MissingET MissingET
575 add Branch ScalarHT/energy ScalarHT ScalarHT
576}
577