Fork me on GitHub

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

Last change on this file since 1397 was 1397, checked in by mselvaggi, 6 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.