← Back to branch summary

~madteam/mg5amcnlo/series2.0

« back to all changes in this revision

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

  • Committer: Marco Zaro
  • Date: 2014-06-10 17:31:23 UTC
  • mfrom: (253.1.38 2.1.2)
  • mto: (253.17.22 2.2.0)
  • mto: This revision was merged to the branch mainline in revision 255.
  • Revision ID: mzaro@ingrid-ui1.cism.ucl.ac.be-20140610173123-773ewpd60dvj7ehv
merged with 2.1.2 r 291

Show diffs side-by-side

added added

removed removed

Lines of Context:
Template/Common/Cards/delphes_card_CMS.dat
19
19
 
20
20
  PhotonEfficiency
21
21
  PhotonIsolation
22
 
  
 
22
 
23
23
  ElectronEfficiency
24
24
  ElectronIsolation
25
 
  
 
25
 
26
26
  MuonEfficiency
27
27
  MuonIsolation
28
 
  
 
28
 
29
29
  MissingET
30
 
  
 
30
 
 
31
  GenJetFinder
31
32
  FastJetFinder
 
33
 
 
34
  JetEnergyScale
 
35
 
32
36
  BTagging
33
37
  TauTagging
34
 
  
 
38
 
35
39
  UniqueObjectFinder
36
 
  
 
40
 
37
41
  ScalarHT
38
 
  
 
42
 
39
43
  TreeWriter
40
44
}
41
45
 
88
92
  set OutputArray electrons
89
93
 
90
94
  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
91
 
              
 
95
 
92
96
  # tracking efficiency formula for electrons
93
97
  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
94
98
                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) + \
100
104
                         (abs(eta) > 2.5)                                                  * (0.00)}
101
105
}
102
106
 
103
 
##############################
 
107
##########################
104
108
# Muon tracking efficiency
105
 
##############################
 
109
##########################
106
110
 
107
111
module Efficiency MuonTrackingEfficiency {
108
112
  set InputArray ParticlePropagator/muons
109
113
  set OutputArray muons
110
114
 
111
115
  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
112
 
              
 
116
 
113
117
  # tracking efficiency formula for muons
114
118
  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
115
119
                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) + \
150
154
 
151
155
  # set ResolutionFormula {resolution formula as a function of eta and energy}
152
156
 
153
 
  set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.5e1) * (energy*0.015) + \
154
 
                                           (abs(eta) <= 2.5) * (energy > 2.5e1)                    * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
155
 
                         (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
 
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) + \
156
161
                         (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
157
162
 
158
163
}
168
173
  # set ResolutionFormula {resolution formula as a function of eta and pt}
169
174
 
170
175
  # resolution formula for muons
171
 
  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
172
 
                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) + \
173
 
                                           (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
174
 
                                           (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) + \
175
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.04) + \
176
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.03) + \
177
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
178
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
 
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)}
179
188
}
180
189
 
181
190
##############
186
195
# add InputArray InputArray
187
196
  add InputArray ChargedHadronMomentumSmearing/chargedHadrons
188
197
  add InputArray ElectronEnergySmearing/electrons
 
198
  add InputArray MuonMomentumSmearing/muons
189
199
  set OutputArray tracks
190
200
}
191
201
 
201
211
  set PhotonOutputArray photons
202
212
 
203
213
  set EFlowTrackOutputArray eflowTracks
204
 
  set EFlowTowerOutputArray eflowTowers
 
214
  set EFlowPhotonOutputArray eflowPhotons
 
215
  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
205
216
 
206
217
  set pi [expr {acos(-1)}]
207
 
  
 
218
 
208
219
  # lists of the edges of each tower in eta and phi
209
220
  # each list starts with the lower edge of the first tower
210
221
  # the list ends with the higher edged of the last tower
211
222
 
212
223
  # 5 degrees towers
213
224
  set PhiBins {}
214
 
  for {set i -36} {$i <= 36} {incr i} { 
 
225
  for {set i -36} {$i <= 36} {incr i} {
215
226
    add PhiBins [expr {$i * $pi/36.0}]
216
227
  }
217
 
  foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} { 
 
228
  foreach eta {-1.566 -1.479 -1.392 -1.305 -1.218 -1.131 -1.044 -0.957 -0.87 -0.783 -0.696 -0.609 -0.522 -0.435 -0.348 -0.261 -0.174 -0.087 0 0.087 0.174 0.261 0.348 0.435 0.522 0.609 0.696 0.783 0.87 0.957 1.044 1.131 1.218 1.305 1.392 1.479 1.566 1.653} {
218
229
    add EtaPhiBins $eta $PhiBins
219
230
  }
220
231
 
221
232
  # 10 degrees towers
222
233
  set PhiBins {}
223
 
  for {set i -18} {$i <= 18} {incr i} { 
 
234
  for {set i -18} {$i <= 18} {incr i} {
224
235
    add PhiBins [expr {$i * $pi/18.0}]
225
236
  }
226
 
  foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} { 
 
237
  foreach eta {-4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -2.95 -2.868 -2.65 -2.5 -2.322 -2.172 -2.043 -1.93 -1.83 -1.74 -1.653 1.74 1.83 1.93 2.043 2.172 2.322 2.5 2.65 2.868 2.95 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525} {
227
238
    add EtaPhiBins $eta $PhiBins
228
239
  }
229
 
 
 
240
 
230
241
  # 20 degrees towers
231
242
  set PhiBins {}
232
 
  for {set i -9} {$i <= 9} {incr i} { 
 
243
  for {set i -9} {$i <= 9} {incr i} {
233
244
    add PhiBins [expr {$i * $pi/9.0}]
234
 
  }  
235
 
  foreach eta {-5 -4.7 -4.525 4.7 5} { 
 
245
  }
 
246
  foreach eta {-5 -4.7 -4.525 4.7 5} {
236
247
    add EtaPhiBins $eta $PhiBins
237
248
  }
238
249
 
242
253
  add EnergyFraction {11} {1.0 0.0}
243
254
  add EnergyFraction {22} {1.0 0.0}
244
255
  add EnergyFraction {111} {1.0 0.0}
245
 
  # energy fractions for muon and neutrinos
 
256
  # energy fractions for muon, neutrinos and neutralinos
246
257
  add EnergyFraction {12} {0.0 0.0}
247
258
  add EnergyFraction {13} {0.0 0.0}
248
259
  add EnergyFraction {14} {0.0 0.0}
249
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}
250
266
  # energy fractions for K0short and Lambda
251
267
  add EnergyFraction {310} {0.3 0.7}
252
268
  add EnergyFraction {3122} {0.3 0.7}
253
 
  
 
269
 
254
270
  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
255
 
  set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.005^2 + energy*0.05^2 + 0.25^2) + \
 
271
  set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2)  + \
256
272
                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
257
273
 
258
274
  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
267
283
module Merger EFlowMerger {
268
284
# add InputArray InputArray
269
285
  add InputArray Calorimeter/eflowTracks
270
 
  add InputArray Calorimeter/eflowTowers
271
 
  add InputArray MuonMomentumSmearing/muons
 
286
  add InputArray Calorimeter/eflowPhotons
 
287
  add InputArray Calorimeter/eflowNeutralHadrons
272
288
  set OutputArray eflow
273
289
}
274
290
 
277
293
###################
278
294
 
279
295
module Efficiency PhotonEfficiency {
280
 
  set InputArray Calorimeter/photons
 
296
  set InputArray Calorimeter/eflowPhotons
281
297
  set OutputArray photons
282
298
 
283
299
  # set EfficiencyFormula {efficiency formula as a function of eta and pt}
296
312
module Isolation PhotonIsolation {
297
313
  set CandidateInputArray PhotonEfficiency/photons
298
314
  set IsolationInputArray EFlowMerger/eflow
299
 
  
 
315
 
300
316
  set OutputArray photons
301
317
 
302
318
  set DeltaRMax 0.5
330
346
module Isolation ElectronIsolation {
331
347
  set CandidateInputArray ElectronEfficiency/electrons
332
348
  set IsolationInputArray EFlowMerger/eflow
333
 
  
 
349
 
334
350
  set OutputArray electrons
335
351
 
336
352
  set DeltaRMax 0.5
351
367
  # set EfficiencyFormula {efficiency as a function of eta and pt}
352
368
 
353
369
  # efficiency formula for muons
354
 
  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
355
 
                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
356
 
                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0)  * (0.85) + \
357
 
                         (abs(eta) > 2.4)                                   * (0.00)}
 
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)) + \
 
375
                         (abs(eta) > 2.4)                                                 * (0.00)}
358
376
}
359
377
 
360
378
################
364
382
module Isolation MuonIsolation {
365
383
  set CandidateInputArray MuonEfficiency/muons
366
384
  set IsolationInputArray EFlowMerger/eflow
367
 
  
 
385
 
368
386
  set OutputArray muons
369
387
 
370
388
  set DeltaRMax 0.5
393
411
  add InputArray UniqueObjectFinder/jets
394
412
  add InputArray UniqueObjectFinder/electrons
395
413
  add InputArray UniqueObjectFinder/photons
396
 
  add InputArray MuonIsolation/muons
 
414
  add InputArray UniqueObjectFinder/muons
397
415
  set EnergyOutputArray energy
398
416
}
399
417
 
 
418
#####################
 
419
# MC truth jet finder
 
420
#####################
 
421
 
 
422
module FastJetFinder GenJetFinder {
 
423
  set InputArray Delphes/stableParticles
 
424
 
 
425
  set OutputArray jets
 
426
 
 
427
  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
 
428
  set JetAlgorithm 6
 
429
  set ParameterR 0.5
 
430
 
 
431
  set JetPTMin 20.0
 
432
}
 
433
 
400
434
############
401
435
# Jet finder
402
436
############
404
438
module FastJetFinder FastJetFinder {
405
439
#  set InputArray Calorimeter/towers
406
440
  set InputArray EFlowMerger/eflow
407
 
  
 
441
 
408
442
  set OutputArray jets
409
443
 
410
444
  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
411
445
  set JetAlgorithm 6
412
 
  set ParameterR 0.7
413
 
 
414
 
  set ConeRadius 0.5
415
 
  set SeedThreshold 1.0
416
 
  set ConeAreaFraction 1.0
417
 
  set AdjacencyCut 2.0
418
 
  set OverlapThreshold 0.75
419
 
 
420
 
  set MaxIterations 100
421
 
  set MaxPairSize 2
422
 
  set Iratch 1
423
 
  
 
446
  set ParameterR 0.5
 
447
 
424
448
  set JetPTMin 20.0
425
449
}
426
450
 
 
451
##################
 
452
# Jet Energy Scale
 
453
##################
 
454
 
 
455
module EnergyScale JetEnergyScale {
 
456
  set InputArray FastJetFinder/jets
 
457
  set OutputArray jets
 
458
 
 
459
 # scale formula for jets
 
460
  set ScaleFormula {1.00}
 
461
}
 
462
 
427
463
###########
428
464
# b-tagging
429
465
###########
430
466
 
431
467
module BTagging BTagging {
432
468
  set PartonInputArray Delphes/partons
433
 
  set JetInputArray FastJetFinder/jets
 
469
  set JetInputArray JetEnergyScale/jets
 
470
 
 
471
  set BitNumber 0
434
472
 
435
473
  set DeltaR 0.5
436
474
 
450
488
  add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) + \
451
489
                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
452
490
                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
453
 
                              (abs(eta) > 2.5)                                  * (0.000)}  
454
 
  
 
491
                              (abs(eta) > 2.5)                                  * (0.000)}
 
492
 
455
493
  # efficiency formula for b-jets
456
494
  add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) + \
457
495
                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
462
500
module TauTagging TauTagging {
463
501
  set ParticleInputArray Delphes/allParticles
464
502
  set PartonInputArray Delphes/partons
465
 
  set JetInputArray FastJetFinder/jets
 
503
  set JetInputArray JetEnergyScale/jets
466
504
 
467
505
  set DeltaR 0.5
468
506
 
487
525
# add InputArray InputArray OutputArray
488
526
  add InputArray PhotonIsolation/photons photons
489
527
  add InputArray ElectronIsolation/electrons electrons
490
 
  add InputArray FastJetFinder/jets jets
 
528
  add InputArray MuonIsolation/muons muons
 
529
  add InputArray JetEnergyScale/jets jets
491
530
}
492
531
 
493
532
##################
494
533
# ROOT tree writer
495
534
##################
496
535
 
 
536
# tracks, towers and eflow objects are not stored by default in the output.
 
537
# if needed (for jet constituent or other studies), uncomment the relevant
 
538
# "add Branch ..." lines.
 
539
 
497
540
module TreeWriter TreeWriter {
498
541
# add Branch InputArray BranchName BranchClass
499
542
  add Branch Delphes/allParticles Particle GenParticle
 
543
 
500
544
  add Branch TrackMerger/tracks Track Track
501
545
  add Branch Calorimeter/towers Tower Tower
 
546
 
502
547
  add Branch Calorimeter/eflowTracks EFlowTrack Track
503
 
  add Branch Calorimeter/eflowTowers EFlowTower Tower
504
 
  add Branch MuonMomentumSmearing/muons EFlowMuon Muon
 
548
  add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
 
549
  add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
 
550
 
 
551
  add Branch GenJetFinder/jets GenJet Jet
505
552
  add Branch UniqueObjectFinder/jets Jet Jet
506
553
  add Branch UniqueObjectFinder/electrons Electron Electron
507
554
  add Branch UniqueObjectFinder/photons Photon Photon
508
 
  add Branch MuonIsolation/muons Muon Muon
 
555
  add Branch UniqueObjectFinder/muons Muon Muon
509
556
  add Branch MissingET/momentum MissingET MissingET
510
557
  add Branch ScalarHT/energy ScalarHT ScalarHT
511
558
}