← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

Viewing changes to Template/Common/Cards/delphes_card_ATLAS.dat

  • Committer: olivier Mattelaer
  • Date: 2016-05-12 11:00:18 UTC
  • mfrom: (262.1.150 2.3.4)
  • Revision ID: olivier.mattelaer@uclouvain.be-20160512110018-sevb79f0wm4g8mpp
pass to 2.4.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/Common/Cards/delphes_card_ATLAS.dat
10
10
  MuonTrackingEfficiency
11
11
 
12
12
  ChargedHadronMomentumSmearing
13
 
  ElectronEnergySmearing
 
13
  ElectronMomentumSmearing
14
14
  MuonMomentumSmearing
15
15
 
16
16
  TrackMerger
20
20
  PhotonEfficiency
21
21
  PhotonIsolation
22
22
 
 
23
  ElectronFilter
23
24
  ElectronEfficiency
24
25
  ElectronIsolation
25
26
 
28
29
 
29
30
  MissingET
30
31
 
 
32
  NeutrinoFilter
31
33
  GenJetFinder
32
34
  FastJetFinder
33
35
 
34
36
  JetEnergyScale
 
37
  
 
38
  JetFlavorAssociation
35
39
 
36
40
  BTagging
37
41
  TauTagging
75
79
  # add EfficiencyFormula {efficiency formula as a function of eta and pt}
76
80
 
77
81
  # tracking efficiency formula for charged hadrons
78
 
  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
79
 
                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
80
 
                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) + \
81
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) + \
82
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) + \
 
82
  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
 
83
                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
 
84
                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) +
 
85
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) +
 
86
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) +
83
87
                         (abs(eta) > 2.5)                                                  * (0.00)}
84
88
}
85
89
 
94
98
  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
95
99
 
96
100
  # tracking efficiency formula for electrons
97
 
  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
98
 
                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) + \
99
 
                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) + \
100
 
                                           (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) + \
101
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) + \
102
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) + \
103
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) + \
 
101
  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
 
102
                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) +
 
103
                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) +
 
104
                                           (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) +
 
105
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) +
 
106
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) +
 
107
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) +
104
108
                         (abs(eta) > 2.5)                                                  * (0.00)}
105
109
}
106
110
 
115
119
  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
116
120
 
117
121
  # tracking efficiency formula for muons
118
 
  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
119
 
                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) + \
120
 
                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) + \
121
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
122
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) + \
 
122
  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
 
123
                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
 
124
                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
 
125
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
 
126
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) +
123
127
                         (abs(eta) > 2.5)                                                  * (0.00)}
124
128
}
125
129
 
134
138
  # set ResolutionFormula {resolution formula as a function of eta and pt}
135
139
 
136
140
  # resolution formula for charged hadrons
137
 
  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) + \
138
 
                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) + \
139
 
                                           (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
140
 
                                           (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) + \
141
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
142
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) + \
143
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
144
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
 
141
  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
 
142
                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
 
143
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
145
144
}
146
145
 
147
 
#################################
148
 
# Energy resolution for electrons
149
 
#################################
 
146
###################################
 
147
# Momentum resolution for electrons
 
148
###################################
150
149
 
151
 
module EnergySmearing ElectronEnergySmearing {
 
150
module MomentumSmearing ElectronMomentumSmearing {
152
151
  set InputArray ElectronTrackingEfficiency/electrons
153
152
  set OutputArray electrons
154
153
 
155
154
  # set ResolutionFormula {resolution formula as a function of eta and energy}
156
155
 
157
 
  set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.5e1) * (energy*0.015) + \
158
 
                                           (abs(eta) <= 2.5) * (energy > 2.5e1)                    * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
159
 
                         (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
160
 
                         (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
161
 
 
 
156
  # resolution formula for electrons
 
157
  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
 
158
                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
 
159
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
162
160
}
163
161
 
164
162
###############################
172
170
  # set ResolutionFormula {resolution formula as a function of eta and pt}
173
171
 
174
172
  # resolution formula for muons
175
 
  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
176
 
                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 5.0e1) * (0.03) + \
177
 
                                           (abs(eta) <= 1.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.04) + \
178
 
                                           (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.07) + \
179
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.04) + \
180
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 5.0e1) * (0.04) + \
181
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0e1 && pt <= 1.0e2) * (0.05) + \
182
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.10)}
 
173
  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.02^2 + pt^2*2.0e-4^2) +
 
174
                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*3.0e-4^2) +
 
175
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*6.0e-4^2)}
183
176
}
184
177
 
185
178
##############
189
182
module Merger TrackMerger {
190
183
# add InputArray InputArray
191
184
  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
192
 
  add InputArray ElectronEnergySmearing/electrons
 
185
  add InputArray ElectronMomentumSmearing/electrons
193
186
  add InputArray MuonMomentumSmearing/muons
194
187
  set OutputArray tracks
195
188
}
209
202
  set EFlowPhotonOutputArray eflowPhotons
210
203
  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
211
204
 
 
205
  set ECalEnergyMin 0.5
 
206
  set HCalEnergyMin 1.0
 
207
 
 
208
  set ECalEnergySignificanceMin 1.0
 
209
  set HCalEnergySignificanceMin 1.0
 
210
 
 
211
  set SmearTowerCenter true
 
212
 
212
213
  set pi [expr {acos(-1)}]
213
214
 
214
215
  # lists of the edges of each tower in eta and phi
257
258
  # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
258
259
  # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
259
260
  # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
260
 
  set ECalResolutionFormula {                  (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) + \
 
261
  set ECalResolutionFormula {                  (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) +
261
262
                             (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
262
263
 
263
264
  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
264
265
  # http://arxiv.org/pdf/hep-ex/0004009v1
265
266
  # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
266
 
  set HCalResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) + \
267
 
                             (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) + \
268
 
                             (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.9420^2 + energy*0.075^2)}
 
267
  set HCalResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
 
268
                             (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
 
269
                             (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
269
270
}
270
271
 
271
272
####################
285
286
###################
286
287
 
287
288
module Efficiency PhotonEfficiency {
288
 
  set InputArray Calorimeter/photons
 
289
  set InputArray Calorimeter/eflowPhotons
289
290
  set OutputArray photons
290
291
 
291
292
  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
292
293
 
293
294
  # efficiency formula for photons
294
 
  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
295
 
                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
296
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
 
295
  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
 
296
                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
 
297
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
297
298
                         (abs(eta) > 2.5)                                   * (0.00)}
298
299
}
299
300
 
314
315
  set PTRatioMax 0.1
315
316
}
316
317
 
 
318
#################
 
319
# Electron filter
 
320
#################
 
321
 
 
322
module PdgCodeFilter ElectronFilter {
 
323
  set InputArray Calorimeter/eflowTracks
 
324
  set OutputArray electrons
 
325
  set Invert true
 
326
  add PdgCode {11}
 
327
  add PdgCode {-11}
 
328
}
 
329
 
317
330
#####################
318
331
# Electron efficiency
319
332
#####################
320
333
 
321
334
module Efficiency ElectronEfficiency {
322
 
  set InputArray ElectronEnergySmearing/electrons
 
335
  set InputArray ElectronFilter/electrons
323
336
  set OutputArray electrons
324
337
 
325
338
  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
326
339
 
327
340
  # efficiency formula for electrons
328
 
  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
329
 
                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
330
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
 
341
  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
 
342
                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
 
343
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
331
344
                         (abs(eta) > 2.5)                                   * (0.00)}
332
345
}
333
346
 
359
372
  # set EfficiencyFormula {efficiency as a function of eta and pt}
360
373
 
361
374
  # efficiency formula for muons
362
 
  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
363
 
                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
364
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0)  * (0.85) + \
 
375
  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
 
376
                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
 
377
                         (abs(eta) > 1.5 && abs(eta) <= 2.7) * (pt > 10.0)  * (0.85) +
365
378
                         (abs(eta) > 2.7)                                   * (0.00)}
366
379
}
367
380
 
405
418
  set EnergyOutputArray energy
406
419
}
407
420
 
 
421
 
 
422
#####################
 
423
# Neutrino Filter
 
424
#####################
 
425
 
 
426
module PdgCodeFilter NeutrinoFilter {
 
427
 
 
428
  set InputArray Delphes/stableParticles
 
429
  set OutputArray filteredParticles
 
430
 
 
431
  set PTMin 0.0
 
432
 
 
433
  add PdgCode {12}
 
434
  add PdgCode {14}
 
435
  add PdgCode {16}
 
436
  add PdgCode {-12}
 
437
  add PdgCode {-14}
 
438
  add PdgCode {-16}
 
439
 
 
440
}
 
441
 
408
442
#####################
409
443
# MC truth jet finder
410
444
#####################
411
445
 
412
446
module FastJetFinder GenJetFinder {
413
 
  set InputArray Delphes/stableParticles
 
447
  set InputArray NeutrinoFilter/filteredParticles
414
448
 
415
449
  set OutputArray jets
416
450
 
421
455
  set JetPTMin 20.0
422
456
}
423
457
 
 
458
 
424
459
############
425
460
# Jet finder
426
461
############
445
480
  set InputArray FastJetFinder/jets
446
481
  set OutputArray jets
447
482
 
448
 
 # scale formula for jets
449
 
  set ScaleFormula {1.00}
 
483
  # scale formula for jets
 
484
  set ScaleFormula {  sqrt( (3.0 - 0.2*(abs(eta)))^2 / pt + 1.0 )  }
 
485
}
 
486
 
 
487
########################
 
488
# Jet Flavor Association
 
489
########################
 
490
 
 
491
module JetFlavorAssociation JetFlavorAssociation {
 
492
  
 
493
  set PartonInputArray Delphes/partons
 
494
  set ParticleInputArray Delphes/allParticles
 
495
  set ParticleLHEFInputArray Delphes/allParticlesLHEF
 
496
  set JetInputArray JetEnergyScale/jets
 
497
  
 
498
  set DeltaR 0.5
 
499
  set PartonPTMin 1.0
 
500
  set PartonEtaMax 2.5
 
501
 
450
502
}
451
503
 
452
504
###########
454
506
###########
455
507
 
456
508
module BTagging BTagging {
457
 
  set PartonInputArray Delphes/partons
458
509
  set JetInputArray JetEnergyScale/jets
459
510
 
460
511
  set BitNumber 0
461
512
 
462
 
  set DeltaR 0.5
463
 
 
464
 
  set PartonPTMin 1.0
465
 
 
466
 
  set PartonEtaMax 2.5
467
 
 
468
513
  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
469
514
  # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
470
515
  # gluon's PDG code has the lowest priority
471
516
 
 
517
  # based on ATL-PHYS-PUB-2015-022
 
518
 
472
519
  # default efficiency formula (misidentification rate)
473
 
  add EfficiencyFormula {0} {0.001}
 
520
  add EfficiencyFormula {0} {0.002+7.3e-06*pt}
474
521
 
475
522
  # efficiency formula for c-jets (misidentification rate)
476
 
  add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) + \
477
 
                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
478
 
                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
479
 
                              (abs(eta) > 2.5)                                  * (0.000)}
 
523
  add EfficiencyFormula {4} {0.20*tanh(0.02*pt)*(1/(1+0.0034*pt))}
480
524
 
481
525
  # efficiency formula for b-jets
482
 
  add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) + \
483
 
                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
484
 
                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
485
 
                              (abs(eta) > 2.5)                                  * (0.000)}
 
526
  add EfficiencyFormula {5} {0.80*tanh(0.003*pt)*(30/(1+0.086*pt))}
486
527
}
487
528
 
 
529
#############
 
530
# tau-tagging
 
531
#############
 
532
 
488
533
module TauTagging TauTagging {
489
534
  set ParticleInputArray Delphes/allParticles
490
535
  set PartonInputArray Delphes/partons
499
544
  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
500
545
 
501
546
  # default efficiency formula (misidentification rate)
502
 
  add EfficiencyFormula {0} {0.001}
 
547
  add EfficiencyFormula {0} {0.01}
503
548
  # efficiency formula for tau-jets
504
 
  add EfficiencyFormula {15} {0.4}
 
549
  add EfficiencyFormula {15} {0.6}
505
550
}
506
551
 
507
552
#####################################################