← Back to branch summary

~maddevelopers/mg5amcnlo/2.6.3_rwgt

« back to all changes in this revision

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

  • Committer: olivier-mattelaer
  • Date: 2017-05-26 07:48:55 UTC
  • mfrom: (271.1.33 2.5.5)
  • Revision ID: olivier-mattelaer-20170526074855-r463wfxlom110fiu
pass the 2.5.5

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/Common/Cards/delphes_card_default.dat
14
14
  MuonMomentumSmearing
15
15
 
16
16
  TrackMerger
 
17
 
 
18
  ECal
 
19
  HCal
 
20
 
17
21
  Calorimeter
18
22
  EFlowMerger
19
 
 
 
23
  EFlowFilter
 
24
  
20
25
  PhotonEfficiency
21
26
  PhotonIsolation
22
27
 
24
29
  ElectronEfficiency
25
30
  ElectronIsolation
26
31
 
 
32
  ChargedHadronFilter
 
33
 
27
34
  MuonEfficiency
28
35
  MuonIsolation
29
36
 
34
41
  GenMissingET
35
42
  
36
43
  FastJetFinder
 
44
  FatJetFinder
37
45
 
38
46
  JetEnergyScale
39
47
 
123
131
  # tracking efficiency formula for muons
124
132
  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
125
133
                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
126
 
                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
 
134
                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e3) * (0.99) +
 
135
                                           (abs(eta) <= 1.5) * (pt > 1.0e3 )               * (0.99 * exp(0.5 - pt*5.0e-4)) +
 
136
 
127
137
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
128
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) +
 
138
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e3) * (0.98) +
 
139
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3)                * (0.98 * exp(0.5 - pt*5.0e-4)) +
129
140
                         (abs(eta) > 2.5)                                                  * (0.00)}
130
141
}
131
142
 
141
152
 
142
153
  # resolution formula for charged hadrons
143
154
  # based on arXiv:1405.6569
144
 
  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.5e-4^2) +
145
 
                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*2.5e-4^2) +
146
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*5.5e-4^2)}
 
155
  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
 
156
                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
 
157
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
147
158
}
148
159
 
149
160
###################################
191
202
  set OutputArray tracks
192
203
}
193
204
 
194
 
#############
195
 
# Calorimeter
196
 
#############
197
 
 
198
 
module Calorimeter Calorimeter {
 
205
 
 
206
 
 
207
#############
 
208
#   ECAL
 
209
#############
 
210
 
 
211
module SimpleCalorimeter ECal {
199
212
  set ParticleInputArray ParticlePropagator/stableParticles
200
213
  set TrackInputArray TrackMerger/tracks
201
214
 
202
 
  set TowerOutputArray towers
203
 
  set PhotonOutputArray photons
204
 
 
205
 
  set EFlowTrackOutputArray eflowTracks
206
 
  set EFlowPhotonOutputArray eflowPhotons
207
 
  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
208
 
 
209
 
  set ECalEnergyMin 0.5
210
 
  set HCalEnergyMin 1.0
211
 
 
212
 
  set ECalEnergySignificanceMin 1.0
213
 
  set HCalEnergySignificanceMin 1.0
 
215
  set TowerOutputArray ecalTowers
 
216
  set EFlowTrackOutputArray eflowTracks
 
217
  set EFlowTowerOutputArray eflowPhotons
 
218
 
 
219
  set IsEcal true
 
220
 
 
221
  set EnergyMin 0.5
 
222
  set EnergySignificanceMin 2.0
 
223
 
 
224
  set SmearTowerCenter true
 
225
 
 
226
  set pi [expr {acos(-1)}]
 
227
 
 
228
  # lists of the edges of each tower in eta and phi
 
229
  # each list starts with the lower edge of the first tower
 
230
  # the list ends with the higher edged of the last tower
 
231
 
 
232
  # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
 
233
 
 
234
  set PhiBins {}
 
235
  for {set i -180} {$i <= 180} {incr i} {
 
236
    add PhiBins [expr {$i * $pi/180.0}]
 
237
  }
 
238
 
 
239
  # 0.02 unit in eta up to eta = 1.5 (barrel)
 
240
  for {set i -85} {$i <= 86} {incr i} {
 
241
    set eta [expr {$i * 0.0174}]
 
242
    add EtaPhiBins $eta $PhiBins
 
243
  }
 
244
 
 
245
  # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
 
246
 
 
247
  set PhiBins {}
 
248
  for {set i -180} {$i <= 180} {incr i} {
 
249
    add PhiBins [expr {$i * $pi/180.0}]
 
250
  }
 
251
 
 
252
  # 0.02 unit in eta up to eta = 3
 
253
  for {set i 1} {$i <= 84} {incr i} {
 
254
    set eta [expr { -2.958 + $i * 0.0174}]
 
255
    add EtaPhiBins $eta $PhiBins
 
256
  }
 
257
 
 
258
  for {set i 1} {$i <= 84} {incr i} {
 
259
    set eta [expr { 1.4964 + $i * 0.0174}]
 
260
    add EtaPhiBins $eta $PhiBins
 
261
  }
 
262
 
 
263
  # take present CMS granularity for HF
 
264
 
 
265
  # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
 
266
  set PhiBins {}
 
267
  for {set i -18} {$i <= 18} {incr i} {
 
268
    add PhiBins [expr {$i * $pi/18.0}]
 
269
  }
 
270
 
 
271
  foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.958 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
 
272
    add EtaPhiBins $eta $PhiBins
 
273
  }
 
274
 
 
275
 
 
276
  add EnergyFraction {0} {0.0}
 
277
  # energy fractions for e, gamma and pi0
 
278
  add EnergyFraction {11} {1.0}
 
279
  add EnergyFraction {22} {1.0}
 
280
  add EnergyFraction {111} {1.0}
 
281
  # energy fractions for muon, neutrinos and neutralinos
 
282
  add EnergyFraction {12} {0.0}
 
283
  add EnergyFraction {13} {0.0}
 
284
  add EnergyFraction {14} {0.0}
 
285
  add EnergyFraction {16} {0.0}
 
286
  add EnergyFraction {1000022} {0.0}
 
287
  add EnergyFraction {1000023} {0.0}
 
288
  add EnergyFraction {1000025} {0.0}
 
289
  add EnergyFraction {1000035} {0.0}
 
290
  add EnergyFraction {1000045} {0.0}
 
291
  # energy fractions for K0short and Lambda
 
292
  add EnergyFraction {310} {0.3}
 
293
  add EnergyFraction {3122} {0.3}
 
294
 
 
295
  # set ResolutionFormula {resolution formula as a function of eta and energy}
 
296
 
 
297
  # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
 
298
 
 
299
  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
 
300
  # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
 
301
  set ResolutionFormula {                      (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
 
302
                             (abs(eta) > 1.5 && abs(eta) <= 2.5) * (2.16 + 5.6*(abs(eta)-2)^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
 
303
                             (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
 
304
 
 
305
}
 
306
 
 
307
 
 
308
#############
 
309
#   HCAL
 
310
#############
 
311
 
 
312
module SimpleCalorimeter HCal {
 
313
  set ParticleInputArray ParticlePropagator/stableParticles
 
314
  set TrackInputArray ECal/eflowTracks
 
315
 
 
316
  set TowerOutputArray hcalTowers
 
317
  set EFlowTrackOutputArray eflowTracks
 
318
  set EFlowTowerOutputArray eflowNeutralHadrons
 
319
 
 
320
  set IsEcal false
 
321
 
 
322
  set EnergyMin 1.0
 
323
  set EnergySignificanceMin 1.0
214
324
 
215
325
  set SmearTowerCenter true
216
326
 
248
358
  }
249
359
 
250
360
  # default energy fractions {abs(PDG code)} {Fecal Fhcal}
251
 
  add EnergyFraction {0} {0.0 1.0}
 
361
  add EnergyFraction {0} {1.0}
252
362
  # energy fractions for e, gamma and pi0
253
 
  add EnergyFraction {11} {1.0 0.0}
254
 
  add EnergyFraction {22} {1.0 0.0}
255
 
  add EnergyFraction {111} {1.0 0.0}
 
363
  add EnergyFraction {11} {0.0}
 
364
  add EnergyFraction {22} {0.0}
 
365
  add EnergyFraction {111} {0.0}
256
366
  # energy fractions for muon, neutrinos and neutralinos
257
 
  add EnergyFraction {12} {0.0 0.0}
258
 
  add EnergyFraction {13} {0.0 0.0}
259
 
  add EnergyFraction {14} {0.0 0.0}
260
 
  add EnergyFraction {16} {0.0 0.0}
261
 
  add EnergyFraction {1000022} {0.0 0.0}
262
 
  add EnergyFraction {1000023} {0.0 0.0}
263
 
  add EnergyFraction {1000025} {0.0 0.0}
264
 
  add EnergyFraction {1000035} {0.0 0.0}
265
 
  add EnergyFraction {1000045} {0.0 0.0}
 
367
  add EnergyFraction {12} {0.0}
 
368
  add EnergyFraction {13} {0.0}
 
369
  add EnergyFraction {14} {0.0}
 
370
  add EnergyFraction {16} {0.0}
 
371
  add EnergyFraction {1000022} {0.0}
 
372
  add EnergyFraction {1000023} {0.0}
 
373
  add EnergyFraction {1000025} {0.0}
 
374
  add EnergyFraction {1000035} {0.0}
 
375
  add EnergyFraction {1000045} {0.0}
266
376
  # energy fractions for K0short and Lambda
267
 
  add EnergyFraction {310} {0.3 0.7}
268
 
  add EnergyFraction {3122} {0.3 0.7}
269
 
 
270
 
  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
271
 
  # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
272
 
  set ECalResolutionFormula {                  (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
273
 
                             (abs(eta) > 1.5 && abs(eta) <= 2.5) * (2.16 + 5.6*(abs(eta)-2)^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
274
 
                             (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
 
377
  add EnergyFraction {310} {0.7}
 
378
  add EnergyFraction {3122} {0.7}
275
379
 
276
380
  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
277
 
  set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
 
381
  set ResolutionFormula {                      (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
278
382
                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
279
 
}
 
383
 
 
384
}
 
385
 
 
386
 
 
387
#################
 
388
# Electron filter
 
389
#################
 
390
 
 
391
module PdgCodeFilter ElectronFilter {
 
392
  set InputArray HCal/eflowTracks
 
393
  set OutputArray electrons
 
394
  set Invert true
 
395
  add PdgCode {11}
 
396
  add PdgCode {-11}
 
397
}
 
398
 
 
399
######################
 
400
# ChargedHadronFilter
 
401
######################
 
402
 
 
403
module PdgCodeFilter ChargedHadronFilter {
 
404
  set InputArray HCal/eflowTracks
 
405
  set OutputArray chargedHadrons
 
406
  
 
407
  add PdgCode {11}
 
408
  add PdgCode {-11}
 
409
  add PdgCode {13}
 
410
  add PdgCode {-13}
 
411
}
 
412
 
 
413
 
 
414
###################################################
 
415
# Tower Merger (in case not using e-flow algorithm)
 
416
###################################################
 
417
 
 
418
module Merger Calorimeter {
 
419
# add InputArray InputArray
 
420
  add InputArray ECal/ecalTowers
 
421
  add InputArray HCal/hcalTowers
 
422
  set OutputArray towers
 
423
}
 
424
 
 
425
 
280
426
 
281
427
####################
282
428
# Energy flow merger
284
430
 
285
431
module Merger EFlowMerger {
286
432
# add InputArray InputArray
287
 
  add InputArray Calorimeter/eflowTracks
288
 
  add InputArray Calorimeter/eflowPhotons
289
 
  add InputArray Calorimeter/eflowNeutralHadrons
290
 
  set OutputArray eflow
291
 
}
 
433
  add InputArray HCal/eflowTracks
 
434
  add InputArray ECal/eflowPhotons
 
435
  add InputArray HCal/eflowNeutralHadrons
 
436
  set OutputArray eflow
 
437
}
 
438
 
 
439
######################
 
440
# EFlowFilter
 
441
######################
 
442
 
 
443
module PdgCodeFilter EFlowFilter {
 
444
  set InputArray EFlowMerger/eflow
 
445
  set OutputArray eflow
 
446
  
 
447
  add PdgCode {11}
 
448
  add PdgCode {-11}
 
449
  add PdgCode {13}
 
450
  add PdgCode {-13}
 
451
}
 
452
 
292
453
 
293
454
###################
294
455
# Photon efficiency
295
456
###################
296
457
 
297
458
module Efficiency PhotonEfficiency {
298
 
  set InputArray Calorimeter/eflowPhotons
 
459
  set InputArray ECal/eflowPhotons
299
460
  set OutputArray photons
300
461
 
301
462
  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
313
474
 
314
475
module Isolation PhotonIsolation {
315
476
  set CandidateInputArray PhotonEfficiency/photons
316
 
  set IsolationInputArray EFlowMerger/eflow
 
477
  set IsolationInputArray EFlowFilter/eflow
317
478
 
318
479
  set OutputArray photons
319
480
 
324
485
  set PTRatioMax 0.12
325
486
}
326
487
 
327
 
#################
328
 
# Electron filter
329
 
#################
330
 
 
331
 
module PdgCodeFilter ElectronFilter {
332
 
  set InputArray Calorimeter/eflowTracks
333
 
  set OutputArray electrons
334
 
  set Invert true
335
 
  add PdgCode {11}
336
 
  add PdgCode {-11}
337
 
}
338
488
 
339
489
#####################
340
490
# Electron efficiency
359
509
 
360
510
module Isolation ElectronIsolation {
361
511
  set CandidateInputArray ElectronEfficiency/electrons
362
 
  set IsolationInputArray EFlowMerger/eflow
 
512
  set IsolationInputArray EFlowFilter/eflow
363
513
 
364
514
  set OutputArray electrons
365
515
 
381
531
  # set EfficiencyFormula {efficiency as a function of eta and pt}
382
532
 
383
533
  # efficiency formula for muons
384
 
  set EfficiencyFormula {                                      (pt <= 10.0)               * (0.00) +
385
 
                                           (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
386
 
                                           (abs(eta) <= 1.5) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) +
387
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
388
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) +
 
534
  set EfficiencyFormula {                                     (pt <= 10.0)                * (0.00) +
 
535
                                           (abs(eta) <= 1.5) * (pt > 10.0)                * (0.95) +
 
536
                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0)                * (0.95) +
389
537
                         (abs(eta) > 2.4)                                                 * (0.00)}
390
538
}
391
539
 
395
543
 
396
544
module Isolation MuonIsolation {
397
545
  set CandidateInputArray MuonEfficiency/muons
398
 
  set IsolationInputArray EFlowMerger/eflow
 
546
  set IsolationInputArray EFlowFilter/eflow
399
547
 
400
548
  set OutputArray muons
401
549
 
497
645
}
498
646
 
499
647
##################
 
648
# Fat Jet finder
 
649
##################
 
650
 
 
651
module FastJetFinder FatJetFinder {
 
652
  set InputArray EFlowMerger/eflow
 
653
 
 
654
  set OutputArray jets
 
655
 
 
656
  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
 
657
  set JetAlgorithm 6
 
658
  set ParameterR 0.8
 
659
 
 
660
  set ComputeNsubjettiness 1
 
661
  set Beta 1.0
 
662
  set AxisMode 4
 
663
 
 
664
  set ComputeTrimming 1
 
665
  set RTrim 0.2
 
666
  set PtFracTrim 0.05
 
667
 
 
668
  set ComputePruning 1
 
669
  set ZcutPrun 0.1
 
670
  set RcutPrun 0.5
 
671
  set RPrun 0.8
 
672
 
 
673
  set ComputeSoftDrop 1
 
674
  set BetaSoftDrop 0.0
 
675
  set SymmetryCutSoftDrop 0.1
 
676
  set R0SoftDrop 0.8
 
677
 
 
678
  set JetPTMin 200.0
 
679
}
 
680
 
 
681
 
 
682
 
 
683
 
 
684
##################
500
685
# Jet Energy Scale
501
686
##################
502
687
 
601
786
  add Branch TrackMerger/tracks Track Track
602
787
  add Branch Calorimeter/towers Tower Tower
603
788
 
604
 
  add Branch Calorimeter/eflowTracks EFlowTrack Track
605
 
  add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
606
 
  add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
 
789
  add Branch HCal/eflowTracks EFlowTrack Track
 
790
  add Branch ECal/eflowPhotons EFlowPhoton Tower
 
791
  add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
607
792
 
608
793
  add Branch GenJetFinder/jets GenJet Jet
609
794
  add Branch GenMissingET/momentum GenMissingET MissingET
610
 
 
 
795
 
611
796
  add Branch UniqueObjectFinder/jets Jet Jet
612
797
  add Branch UniqueObjectFinder/electrons Electron Electron
613
798
  add Branch UniqueObjectFinder/photons Photon Photon
614
799
  add Branch UniqueObjectFinder/muons Muon Muon
 
800
 
 
801
  add Branch FatJetFinder/jets FatJet Jet
 
802
 
615
803
  add Branch MissingET/momentum MissingET MissingET
616
804
  add Branch ScalarHT/energy ScalarHT ScalarHT
617
805
}