Fork me on GitHub

Changeset a2983ec in git for cards/delphes_card_CMS.tcl


Ignore:
Timestamp:
Dec 22, 2014, 11:38:33 AM (10 years ago)
Author:
Pavel Demin <pavel.demin@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
8c72f09
Parents:
4e09c3a
Message:

add parameter SmearTowerCenter to all cards

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/delphes_card_CMS.tcl

    r4e09c3a ra2983ec  
    7777
    7878  # tracking efficiency formula for charged hadrons
    79   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    80                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    81                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) + \
    82                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) + \
    83                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) + \
     79  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     80                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     81                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) +
     82                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) +
     83                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) +
    8484                         (abs(eta) > 2.5)                                                  * (0.00)}
    8585}
    … …  
    9696
    9797  # tracking efficiency formula for electrons
    98   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    99                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) + \
    100                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) + \
    101                                            (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) + \
    102                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) + \
    103                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) + \
    104                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) + \
     98  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     99                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) +
     100                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) +
     101                                           (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) +
     102                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) +
     103                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) +
     104                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) +
    105105                         (abs(eta) > 2.5)                                                  * (0.00)}
    106106}
    … …  
    117117
    118118  # tracking efficiency formula for muons
    119   set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) + \
    120                                            (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) + \
    121                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) + \
    122                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) + \
    123                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) + \
     119  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     120                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
     121                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
     122                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     123                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) +
    124124                         (abs(eta) > 2.5)                                                  * (0.00)}
    125125}
    … …  
    136136
    137137  # resolution formula for charged hadrons
    138   set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) + \
    139                                            (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) + \
    140                                            (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) + \
    141                                            (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) + \
    142                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) + \
    143                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) + \
    144                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) + \
     138  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) +
     139                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) +
     140                                           (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) +
     141                                           (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) +
     142                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
     143                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) +
     144                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) +
    145145                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
    146146}
    … …  
    157157
    158158  # resolution formula for electrons
    159   set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.0225) + \
    160                                            (abs(eta) <= 2.5) * (energy > 2.0e1)                    * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
    161                          (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) + \
     159  set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.0225) +
     160                                           (abs(eta) <= 2.5) * (energy > 2.0e1)                    * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
     161                         (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
    162162                         (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    163163
    … …  
    175175
    176176  # resolution formula for muons
    177   set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1   && pt <= 5.0)   * (0.02) + \
    178                                            (abs(eta) <= 0.5) * (pt > 5.0   && pt <= 1.0e2) * (0.015) + \
    179                                            (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) + \
    180                                            (abs(eta) <= 0.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) + \
    181                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1   && pt <= 5.0)   * (0.03) + \
    182                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0   && pt <= 1.0e2) * (0.02) + \
    183                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) + \
    184                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) + \
    185                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 5.0)   * (0.04) + \
    186                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0   && pt <= 1.0e2) * (0.035) + \
    187                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) + \
     177  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1   && pt <= 5.0)   * (0.02) +
     178                                           (abs(eta) <= 0.5) * (pt > 5.0   && pt <= 1.0e2) * (0.015) +
     179                                           (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) +
     180                                           (abs(eta) <= 0.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) +
     181                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1   && pt <= 5.0)   * (0.03) +
     182                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0   && pt <= 1.0e2) * (0.02) +
     183                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) +
     184                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) +
     185                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 5.0)   * (0.04) +
     186                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0   && pt <= 1.0e2) * (0.035) +
     187                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) +
    188188                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4)}
    189189}
    … …  
    212212  set PhotonOutputArray photons
    213213
    214   set ECalEnergyMin 0.5
    215   set HCalEnergyMin 1.0
    216 
    217   set ECalEnergySignificanceMin 1.0
    218   set HCalEnergySignificanceMin 1.0
    219    
    220214  set EFlowTrackOutputArray eflowTracks
    221215  set EFlowPhotonOutputArray eflowPhotons
    222216  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
     217
     218  set ECalEnergyMin 0.5
     219  set HCalEnergyMin 1.0
     220
     221  set ECalEnergySignificanceMin 1.0
     222  set HCalEnergySignificanceMin 1.0
     223
     224  set SmearTowerCenter true
    223225
    224226  set pi [expr {acos(-1)}]
    … …  
    276278
    277279  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    278   set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2)  + \
     280  set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2)  +
    279281                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    280282
    281283  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    282   set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) + \
     284  set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
    283285                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
    284286}
    … …  
    307309
    308310  # efficiency formula for photons
    309   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
    310                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
    311                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
     311  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     312                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     313                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
    312314                         (abs(eta) > 2.5)                                   * (0.00)}
    313315}
    … …  
    341343
    342344  # efficiency formula for electrons
    343   set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) + \
    344                                            (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) + \
    345                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) + \
     345  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     346                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     347                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
    346348                         (abs(eta) > 2.5)                                   * (0.00)}
    347349}
    … …  
    375377
    376378  # efficiency formula for muons
    377   set EfficiencyFormula {                                      (pt <= 10.0)               * (0.00) + \
    378                                            (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
    379                                            (abs(eta) <= 1.5) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) + \
    380                          (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) + \
    381                          (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) + \
     379  set EfficiencyFormula {                                      (pt <= 10.0)               * (0.00) +
     380                                           (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
     381                                           (abs(eta) <= 1.5) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) +
     382                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
     383                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) +
    382384                         (abs(eta) > 2.4)                                                 * (0.00)}
    383385}
    … …  
    429431
    430432module PdgCodeFilter NeutrinoFilter {
    431  
     433
    432434  set InputArray Delphes/stableParticles
    433435  set OutputArray filteredParticles
    434436
    435437  set PTMin 0.0
    436  
     438
    437439  add PdgCode {12}
    438440  add PdgCode {14}
    … …  
    518520
    519521  # efficiency formula for c-jets (misidentification rate)
    520   add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) + \
    521                                                 (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) + \
    522                               (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) + \
     522  add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) +
     523                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) +
     524                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) +
    523525                              (abs(eta) > 2.5)                                  * (0.000)}
    524526
    525527  # efficiency formula for b-jets
    526   add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) + \
    527                                                 (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) + \
    528                               (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) + \
     528  add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) +
     529                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) +
     530                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) +
    529531                              (abs(eta) > 2.5)                                  * (0.000)}
    530532}
    … …  
    589591  add Branch ScalarHT/energy ScalarHT ScalarHT
    590592}
    591 
Note: See TracChangeset for help on using the changeset viewer.