← Back to branch summary

~maddevelopers/mg5amcnlo/2.5.4_run_py8_at_evtgen

« back to all changes in this revision

Viewing changes to Template/Common/Cards/delphes_card_CMS.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_CMS.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
35
37
 
 
38
  JetFlavorAssociation
 
39
 
36
40
  BTagging
37
41
  TauTagging
38
42
 
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
  # based on arXiv:1405.6569
 
142
  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
 
143
                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
 
144
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
145
145
}
146
146
 
147
 
#################################
148
 
# Energy resolution for electrons
149
 
#################################
 
147
###################################
 
148
# Momentum resolution for electrons
 
149
###################################
150
150
 
151
 
module EnergySmearing ElectronEnergySmearing {
 
151
module MomentumSmearing ElectronMomentumSmearing {
152
152
  set InputArray ElectronTrackingEfficiency/electrons
153
153
  set OutputArray electrons
154
154
 
155
155
  # set ResolutionFormula {resolution formula as a function of eta and energy}
156
156
 
157
157
  # resolution formula for electrons
158
 
  set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.0225) + \
159
 
                                           (abs(eta) <= 2.5) * (energy > 2.0e1)                    * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
160
 
                         (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
161
 
                         (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
162
 
 
 
158
  # based on arXiv:1405.6569
 
159
  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
 
160
                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
 
161
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
163
162
}
164
163
 
165
164
###############################
173
172
  # set ResolutionFormula {resolution formula as a function of eta and pt}
174
173
 
175
174
  # resolution formula for muons
176
 
  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1   && pt <= 5.0)   * (0.02) + \
177
 
                                           (abs(eta) <= 0.5) * (pt > 5.0   && pt <= 1.0e2) * (0.015) + \
178
 
                                           (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + \
179
 
                                           (abs(eta) <= 0.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) + \
180
 
                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1   && pt <= 5.0)   * (0.03) + \
181
 
                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0   && pt <= 1.0e2) * (0.02) + \
182
 
                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + \
183
 
                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) + \
184
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 5.0)   * (0.04) + \
185
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0   && pt <= 1.0e2) * (0.035) + \
186
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + \
187
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4)}
 
175
  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*2.0e-4^2) +
 
176
                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.02^2 + pt^2*3.0e-4^2) +
 
177
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*6.0e-4^2)}
188
178
}
189
179
 
190
180
##############
194
184
module Merger TrackMerger {
195
185
# add InputArray InputArray
196
186
  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
197
 
  add InputArray ElectronEnergySmearing/electrons
 
187
  add InputArray ElectronMomentumSmearing/electrons
198
188
  add InputArray MuonMomentumSmearing/muons
199
189
  set OutputArray tracks
200
190
}
214
204
  set EFlowPhotonOutputArray eflowPhotons
215
205
  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
216
206
 
 
207
  set ECalEnergyMin 0.5
 
208
  set HCalEnergyMin 1.0
 
209
 
 
210
  set ECalEnergySignificanceMin 1.0
 
211
  set HCalEnergySignificanceMin 1.0
 
212
 
 
213
  set SmearTowerCenter true
 
214
 
217
215
  set pi [expr {acos(-1)}]
218
216
 
219
217
  # lists of the edges of each tower in eta and phi
268
266
  add EnergyFraction {3122} {0.3 0.7}
269
267
 
270
268
  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
271
 
  set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2)  + \
 
269
  set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2)  +
272
270
                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
273
271
 
274
272
  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
275
 
  set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \
 
273
  set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
276
274
                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
277
275
}
278
276
 
299
297
  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
300
298
 
301
299
  # efficiency formula for photons
302
 
  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
303
 
                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
304
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
 
300
  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
 
301
                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
 
302
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
305
303
                         (abs(eta) > 2.5)                                   * (0.00)}
306
304
}
307
305
 
322
320
  set PTRatioMax 0.1
323
321
}
324
322
 
 
323
#################
 
324
# Electron filter
 
325
#################
 
326
 
 
327
module PdgCodeFilter ElectronFilter {
 
328
  set InputArray Calorimeter/eflowTracks
 
329
  set OutputArray electrons
 
330
  set Invert true
 
331
  add PdgCode {11}
 
332
  add PdgCode {-11}
 
333
}
 
334
 
325
335
#####################
326
336
# Electron efficiency
327
337
#####################
328
338
 
329
339
module Efficiency ElectronEfficiency {
330
 
  set InputArray ElectronEnergySmearing/electrons
 
340
  set InputArray ElectronFilter/electrons
331
341
  set OutputArray electrons
332
342
 
333
343
  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
334
344
 
335
345
  # efficiency formula for electrons
336
 
  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
337
 
                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
338
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
 
346
  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
 
347
                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
 
348
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
339
349
                         (abs(eta) > 2.5)                                   * (0.00)}
340
350
}
341
351
 
367
377
  # set EfficiencyFormula {efficiency as a function of eta and pt}
368
378
 
369
379
  # efficiency formula for muons
370
 
  set EfficiencyFormula {                                      (pt <= 10.0)               * (0.00) + \
371
 
                                           (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
372
 
                                           (abs(eta) <= 1.5) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) + \
373
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
374
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) + \
 
380
  set EfficiencyFormula {                                      (pt <= 10.0)               * (0.00) +
 
381
                                           (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
 
382
                                           (abs(eta) <= 1.5) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) +
 
383
                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
 
384
                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) +
375
385
                         (abs(eta) > 2.4)                                                 * (0.00)}
376
386
}
377
387
 
415
425
  set EnergyOutputArray energy
416
426
}
417
427
 
 
428
 
 
429
#####################
 
430
# Neutrino Filter
 
431
#####################
 
432
 
 
433
module PdgCodeFilter NeutrinoFilter {
 
434
 
 
435
  set InputArray Delphes/stableParticles
 
436
  set OutputArray filteredParticles
 
437
 
 
438
  set PTMin 0.0
 
439
 
 
440
  add PdgCode {12}
 
441
  add PdgCode {14}
 
442
  add PdgCode {16}
 
443
  add PdgCode {-12}
 
444
  add PdgCode {-14}
 
445
  add PdgCode {-16}
 
446
 
 
447
}
 
448
 
 
449
 
418
450
#####################
419
451
# MC truth jet finder
420
452
#####################
421
453
 
422
454
module FastJetFinder GenJetFinder {
423
 
  set InputArray Delphes/stableParticles
 
455
  set InputArray NeutrinoFilter/filteredParticles
424
456
 
425
457
  set OutputArray jets
426
458
 
431
463
  set JetPTMin 20.0
432
464
}
433
465
 
 
466
 
434
467
############
435
468
# Jet finder
436
469
############
456
489
  set InputArray FastJetFinder/jets
457
490
  set OutputArray jets
458
491
 
459
 
 # scale formula for jets
460
 
  set ScaleFormula {1.00}
 
492
  # scale formula for jets
 
493
  set ScaleFormula {sqrt( (2.5 - 0.15*(abs(eta)))^2 / pt + 1.0 )}
 
494
}
 
495
 
 
496
########################
 
497
# Jet Flavor Association
 
498
########################
 
499
 
 
500
module JetFlavorAssociation JetFlavorAssociation {
 
501
 
 
502
  set PartonInputArray Delphes/partons
 
503
  set ParticleInputArray Delphes/allParticles
 
504
  set ParticleLHEFInputArray Delphes/allParticlesLHEF
 
505
  set JetInputArray JetEnergyScale/jets
 
506
 
 
507
  set DeltaR 0.5
 
508
  set PartonPTMin 1.0
 
509
  set PartonEtaMax 2.5
 
510
 
461
511
}
462
512
 
463
513
###########
465
515
###########
466
516
 
467
517
module BTagging BTagging {
468
 
  set PartonInputArray Delphes/partons
469
518
  set JetInputArray JetEnergyScale/jets
470
519
 
471
520
  set BitNumber 0
472
521
 
473
 
  set DeltaR 0.5
474
 
 
475
 
  set PartonPTMin 1.0
476
 
 
477
 
  set PartonEtaMax 2.5
478
 
 
479
522
  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
480
523
  # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
481
524
  # gluon's PDG code has the lowest priority
482
525
 
483
 
  # https://twiki.cern.ch/twiki/bin/view/CMSPublic/PhysicsResultsBTV
 
526
  # based on arXiv:1211.4462
 
527
  
484
528
  # default efficiency formula (misidentification rate)
485
 
  add EfficiencyFormula {0} {0.001}
 
529
  add EfficiencyFormula {0} {0.01+0.00038*pt}
486
530
 
487
531
  # efficiency formula for c-jets (misidentification rate)
488
 
  add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) + \
489
 
                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
490
 
                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
491
 
                              (abs(eta) > 2.5)                                  * (0.000)}
 
532
  add EfficiencyFormula {4} {0.25*tanh(0.018*pt)*(1/(1+ 0.0013*pt))}
492
533
 
493
534
  # efficiency formula for b-jets
494
 
  add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) + \
495
 
                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
496
 
                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
497
 
                              (abs(eta) > 2.5)                                  * (0.000)}
 
535
  add EfficiencyFormula {5} {0.85*tanh(0.0025*pt)*(25.0/(1+0.063*pt))}
498
536
}
499
537
 
 
538
#############
 
539
# tau-tagging
 
540
#############
 
541
 
500
542
module TauTagging TauTagging {
501
543
  set ParticleInputArray Delphes/allParticles
502
544
  set PartonInputArray Delphes/partons
511
553
  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
512
554
 
513
555
  # default efficiency formula (misidentification rate)
514
 
  add EfficiencyFormula {0} {0.001}
 
556
  add EfficiencyFormula {0} {0.01}
515
557
  # efficiency formula for tau-jets
516
 
  add EfficiencyFormula {15} {0.4}
 
558
  add EfficiencyFormula {15} {0.6}
517
559
}
518
560
 
519
561
#####################################################
556
598
  add Branch MissingET/momentum MissingET MissingET
557
599
  add Branch ScalarHT/energy ScalarHT ScalarHT
558
600
}
559