Fork me on GitHub

source: git/examples/delphes_card_ATLAS.tcl@ 9e991f8

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 9e991f8 was 35b9204, checked in by Michele <michele.selvaggi@…>, 10 years ago

Updated Jet Energy scale values for delphes_card_ATLAS.tcl and delphes_card_CMS.tcl cards.
Added Example4.C in the examples directory as an example on how to compute JES.
  • Property mode set to 100644
File size: 18.7 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 EcalTowerMinEnergy 0.50
210 set HcalTowerMinEnergy 1.00
211
212 set EcalTowerMinSignificance 1.0
213 set HcalTowerMinSignificance 1.0
214
215 set EFlowTrackOutputArray eflowTracks
216 set EFlowPhotonOutputArray eflowPhotons
217 set EFlowNeutralHadronOutputArray eflowNeutralHadrons
218
219 set pi [expr {acos(-1)}]
220
221 # lists of the edges of each tower in eta and phi
222 # each list starts with the lower edge of the first tower
223 # the list ends with the higher edged of the last tower
224
225 # 10 degrees towers
226 set PhiBins {}
227 for {set i -18} {$i <= 18} {incr i} {
228 add PhiBins [expr {$i * $pi/18.0}]
229 }
230 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} {
231 add EtaPhiBins $eta $PhiBins
232 }
233
234 # 20 degrees towers
235 set PhiBins {}
236 for {set i -9} {$i <= 9} {incr i} {
237 add PhiBins [expr {$i * $pi/9.0}]
238 }
239 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} {
240 add EtaPhiBins $eta $PhiBins
241 }
242
243 # default energy fractions {abs(PDG code)} {Fecal Fhcal}
244 add EnergyFraction {0} {0.0 1.0}
245 # energy fractions for e, gamma and pi0
246 add EnergyFraction {11} {1.0 0.0}
247 add EnergyFraction {22} {1.0 0.0}
248 add EnergyFraction {111} {1.0 0.0}
249 # energy fractions for muon, neutrinos and neutralinos
250 add EnergyFraction {12} {0.0 0.0}
251 add EnergyFraction {13} {0.0 0.0}
252 add EnergyFraction {14} {0.0 0.0}
253 add EnergyFraction {16} {0.0 0.0}
254 add EnergyFraction {1000022} {0.0 0.0}
255 add EnergyFraction {1000023} {0.0 0.0}
256 add EnergyFraction {1000025} {0.0 0.0}
257 add EnergyFraction {1000035} {0.0 0.0}
258 add EnergyFraction {1000045} {0.0 0.0}
259 # energy fractions for K0short and Lambda
260 add EnergyFraction {310} {0.3 0.7}
261 add EnergyFraction {3122} {0.3 0.7}
262
263 # set ECalResolutionFormula {resolution formula as a function of eta and energy}
264 # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
265 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
266 # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
267 set ECalResolutionFormula { (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + \
268 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
269
270 # set HCalResolutionFormula {resolution formula as a function of eta and energy}
271 # http://arxiv.org/pdf/hep-ex/0004009v1
272 # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
273 set HCalResolutionFormula { (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
274 (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \
275 (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
276}
277
278####################
279# Energy flow merger
280####################
281
282module Merger EFlowMerger {
283# add InputArray InputArray
284 add InputArray Calorimeter/eflowTracks
285 add InputArray Calorimeter/eflowPhotons
286 add InputArray Calorimeter/eflowNeutralHadrons
287 set OutputArray eflow
288}
289
290###################
291# Photon efficiency
292###################
293
294module Efficiency PhotonEfficiency {
295 set InputArray Calorimeter/photons
296 set OutputArray photons
297
298 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
299
300 # efficiency formula for photons
301 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
302 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
303 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
304 (abs(eta) > 2.5) * (0.00)}
305}
306
307##################
308# Photon isolation
309##################
310
311module Isolation PhotonIsolation {
312 set CandidateInputArray PhotonEfficiency/photons
313 set IsolationInputArray EFlowMerger/eflow
314
315 set OutputArray photons
316
317 set DeltaRMax 0.5
318
319 set PTMin 0.5
320
321 set PTRatioMax 0.1
322}
323
324#####################
325# Electron efficiency
326#####################
327
328module Efficiency ElectronEfficiency {
329 set InputArray ElectronEnergySmearing/electrons
330 set OutputArray electrons
331
332 # set EfficiencyFormula {efficiency formula as a function of eta and pt}
333
334 # efficiency formula for electrons
335 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
336 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
337 (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0) * (0.85) + \
338 (abs(eta) > 2.5) * (0.00)}
339}
340
341####################
342# Electron isolation
343####################
344
345module Isolation ElectronIsolation {
346 set CandidateInputArray ElectronEfficiency/electrons
347 set IsolationInputArray EFlowMerger/eflow
348
349 set OutputArray electrons
350
351 set DeltaRMax 0.5
352
353 set PTMin 0.5
354
355 set PTRatioMax 0.1
356}
357
358#################
359# Muon efficiency
360#################
361
362module Efficiency MuonEfficiency {
363 set InputArray MuonMomentumSmearing/muons
364 set OutputArray muons
365
366 # set EfficiencyFormula {efficiency as a function of eta and pt}
367
368 # efficiency formula for muons
369 set EfficiencyFormula { (pt <= 10.0) * (0.00) + \
370 (abs(eta) <= 1.5) * (pt > 10.0) * (0.95) + \
371 (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0) * (0.85) + \
372 (abs(eta) > 2.7) * (0.00)}
373}
374
375################
376# Muon isolation
377################
378
379module Isolation MuonIsolation {
380 set CandidateInputArray MuonEfficiency/muons
381 set IsolationInputArray EFlowMerger/eflow
382
383 set OutputArray muons
384
385 set DeltaRMax 0.5
386
387 set PTMin 0.5
388
389 set PTRatioMax 0.1
390}
391
392###################
393# Missing ET merger
394###################
395
396module Merger MissingET {
397# add InputArray InputArray
398 add InputArray EFlowMerger/eflow
399 set MomentumOutputArray momentum
400}
401
402##################
403# Scalar HT merger
404##################
405
406module Merger ScalarHT {
407# add InputArray InputArray
408 add InputArray UniqueObjectFinder/jets
409 add InputArray UniqueObjectFinder/electrons
410 add InputArray UniqueObjectFinder/photons
411 add InputArray UniqueObjectFinder/muons
412 set EnergyOutputArray energy
413}
414
415
416#####################
417# Neutrino Filter
418#####################
419
420module PdgCodeFilter NeutrinoFilter {
421
422 set InputArray Delphes/stableParticles
423 set OutputArray filteredParticles
424
425 set PTMin 0.0
426
427 add PdgCode {12}
428 add PdgCode {14}
429 add PdgCode {16}
430 add PdgCode {-12}
431 add PdgCode {-14}
432 add PdgCode {-16}
433
434}
435
436#####################
437# MC truth jet finder
438#####################
439
440module FastJetFinder GenJetFinder {
441 set InputArray NeutrinoFilter/filteredParticles
442
443 set OutputArray jets
444
445 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
446 set JetAlgorithm 6
447 set ParameterR 0.6
448
449 set JetPTMin 20.0
450}
451
452
453############
454# Jet finder
455############
456
457module FastJetFinder FastJetFinder {
458 set InputArray Calorimeter/towers
459
460 set OutputArray jets
461
462 # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
463 set JetAlgorithm 6
464 set ParameterR 0.6
465
466 set JetPTMin 20.0
467}
468
469##################
470# Jet Energy Scale
471##################
472
473module EnergyScale JetEnergyScale {
474 set InputArray FastJetFinder/jets
475 set OutputArray jets
476
477 # scale formula for jets
478 set ScaleFormula { sqrt( (3.0 - 0.2*(abs(eta)))^2 / pt + 1.0 ) }
479}
480
481###########
482# b-tagging
483###########
484
485module BTagging BTagging {
486 set PartonInputArray Delphes/partons
487 set JetInputArray JetEnergyScale/jets
488
489 set BitNumber 0
490
491 set DeltaR 0.5
492
493 set PartonPTMin 1.0
494
495 set PartonEtaMax 2.5
496
497 # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
498 # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
499 # gluon's PDG code has the lowest priority
500
501 # default efficiency formula (misidentification rate)
502 add EfficiencyFormula {0} {0.001}
503
504 # efficiency formula for c-jets (misidentification rate)
505 add EfficiencyFormula {4} { (pt <= 15.0) * (0.000) + \
506 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
507 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
508 (abs(eta) > 2.5) * (0.000)}
509
510 # efficiency formula for b-jets
511 add EfficiencyFormula {5} { (pt <= 15.0) * (0.000) + \
512 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
513 (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
514 (abs(eta) > 2.5) * (0.000)}
515}
516
517module TauTagging TauTagging {
518 set ParticleInputArray Delphes/allParticles
519 set PartonInputArray Delphes/partons
520 set JetInputArray JetEnergyScale/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 MuonIsolation/muons muons
546 add InputArray JetEnergyScale/jets jets
547}
548
549##################
550# ROOT tree writer
551##################
552
553# tracks, towers and eflow objects are not stored by default in the output.
554# if needed (for jet constituent or other studies), uncomment the relevant
555# "add Branch ..." lines.
556
557module TreeWriter TreeWriter {
558# add Branch InputArray BranchName BranchClass
559 add Branch Delphes/allParticles Particle GenParticle
560
561 add Branch TrackMerger/tracks Track Track
562 add Branch Calorimeter/towers Tower Tower
563
564 add Branch Calorimeter/eflowTracks EFlowTrack Track
565 add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
566 add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
567
568 add Branch GenJetFinder/jets GenJet Jet
569 add Branch UniqueObjectFinder/jets Jet Jet
570 add Branch UniqueObjectFinder/electrons Electron Electron
571 add Branch UniqueObjectFinder/photons Photon Photon
572 add Branch UniqueObjectFinder/muons Muon Muon
573 add Branch MissingET/momentum MissingET MissingET
574 add Branch ScalarHT/energy ScalarHT ScalarHT
575}
576
Note: See TracBrowser for help on using the repository browser.