Fork me on GitHub

Changeset d77b51d in git for cards/delphes_card_CMS.tcl


Ignore:
Timestamp:
Sep 29, 2015, 2:08:10 PM (9 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
a98c7ef
Parents:
d870fc5 (diff), 06ec139 (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.
Message:

Merge remote-tracking branch 'upstream/master'

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/delphes_card_CMS.tcl

    rd870fc5 rd77b51d  
    3434
    3535  JetEnergyScale
     36
     37  JetFlavorAssociation
    3638
    3739  BTagging
    … …  
    7779
    7880  # 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) + \
     81  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     82                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     83                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.95) +
     84                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.60) +
     85                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.85) +
    8486                         (abs(eta) > 2.5)                                                  * (0.00)}
    8587}
    … …  
    9698
    9799  # 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) + \
     100  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     101                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.73) +
     102                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e2) * (0.95) +
     103                                           (abs(eta) <= 1.5) * (pt > 1.0e2)                * (0.99) +
     104                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.50) +
     105                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e2) * (0.83) +
     106                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2)                * (0.90) +
    105107                         (abs(eta) > 2.5)                                                  * (0.00)}
    106108}
    … …  
    117119
    118120  # 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) + \
     121  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
     122                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
     123                                           (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
     124                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
     125                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) +
    124126                         (abs(eta) > 2.5)                                                  * (0.00)}
    125127}
    … …  
    136138
    137139  # 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) + \
     140  set ResolutionFormula {                  (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.02) +
     141                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e1) * (0.01) +
     142                                           (abs(eta) <= 1.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.03) +
     143                                           (abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05) +
     144                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.03) +
     145                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e1) * (0.02) +
     146                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e1 && pt <= 2.0e2) * (0.04) +
    145147                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05)}
    146148}
    … …  
    157159
    158160  # 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) + \
     161  set ResolutionFormula {                  (abs(eta) <= 2.5) * (energy > 0.1   && energy <= 2.0e1) * (energy*0.0225) +
     162                                           (abs(eta) <= 2.5) * (energy > 2.0e1)                    * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
     163                         (abs(eta) > 2.5 && abs(eta) <= 3.0)                                       * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2) +
    162164                         (abs(eta) > 3.0 && abs(eta) <= 5.0)                                       * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    163165
    … …  
    175177
    176178  # 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) + \
     179  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1   && pt <= 5.0)   * (0.02) +
     180                                           (abs(eta) <= 0.5) * (pt > 5.0   && pt <= 1.0e2) * (0.015) +
     181                                           (abs(eta) <= 0.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.03) +
     182                                           (abs(eta) <= 0.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) +
     183                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1   && pt <= 5.0)   * (0.03) +
     184                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 5.0   && pt <= 1.0e2) * (0.02) +
     185                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.04) +
     186                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4) +
     187                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 5.0)   * (0.04) +
     188                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 5.0   && pt <= 1.0e2) * (0.035) +
     189                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e2 && pt <= 2.0e2) * (0.05) +
    188190                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 2.0e2)                * (0.05 + pt*1.e-4)}
    189191}
    … …  
    212214  set PhotonOutputArray photons
    213215
    214   set ECalEnergyMin 0.5
    215   set HCalEnergyMin 1.0
    216 
    217   set ECalEnergySignificanceMin 1.0
    218   set HCalEnergySignificanceMin 1.0
    219    
    220216  set EFlowTrackOutputArray eflowTracks
    221217  set EFlowPhotonOutputArray eflowPhotons
    222218  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
     219
     220  set ECalEnergyMin 0.5
     221  set HCalEnergyMin 1.0
     222
     223  set ECalEnergySignificanceMin 1.0
     224  set HCalEnergySignificanceMin 1.0
     225
     226  set SmearTowerCenter true
    223227
    224228  set pi [expr {acos(-1)}]
    … …  
    276280
    277281  # 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)  + \
     282  set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2)  +
    279283                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    280284
    281285  # 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) + \
     286  set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
    283287                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
    284288}
    … …  
    307311
    308312  # 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) + \
     313  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     314                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     315                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
    312316                         (abs(eta) > 2.5)                                   * (0.00)}
    313317}
    … …  
    341345
    342346  # 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) + \
     347  set EfficiencyFormula {                                      (pt <= 10.0) * (0.00) +
     348                                           (abs(eta) <= 1.5) * (pt > 10.0)  * (0.95) +
     349                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 10.0)  * (0.85) +
    346350                         (abs(eta) > 2.5)                                   * (0.00)}
    347351}
    … …  
    375379
    376380  # 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)) + \
     381  set EfficiencyFormula {                                      (pt <= 10.0)               * (0.00) +
     382                                           (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
     383                                           (abs(eta) <= 1.5) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) +
     384                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
     385                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) +
    382386                         (abs(eta) > 2.4)                                                 * (0.00)}
    383387}
    … …  
    429433
    430434module PdgCodeFilter NeutrinoFilter {
    431  
     435
    432436  set InputArray Delphes/stableParticles
    433437  set OutputArray filteredParticles
    434438
    435439  set PTMin 0.0
    436  
     440
    437441  add PdgCode {12}
    438442  add PdgCode {14}
    … …  
    462466
    463467
    464 
    465 
    466468############
    467469# Jet finder
    … …  
    493495}
    494496
     497########################
     498# Jet Flavor Association
     499########################
     500
     501module JetFlavorAssociation JetFlavorAssociation {
     502 
     503  set PartonInputArray Delphes/partons
     504  set ParticleInputArray Delphes/allParticles
     505  set ParticleLHEFInputArray Delphes/allParticlesLHEF
     506  set JetInputArray JetEnergyScale/jets
     507 
     508  set DeltaR 0.5
     509  set PartonPTMin 1.0
     510  set PartonEtaMax 2.5
     511
     512}
     513
    495514###########
    496515# b-tagging
    … …  
    498517
    499518module BTagging BTagging {
    500   set PartonInputArray Delphes/partons
    501519  set JetInputArray JetEnergyScale/jets
    502520
    503521  set BitNumber 0
    504 
    505   set DeltaR 0.5
    506 
    507   set PartonPTMin 1.0
    508 
    509   set PartonEtaMax 2.5
    510522
    511523  # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
    … …  
    518530
    519531  # 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)) + \
     532  add EfficiencyFormula {4} {                                      (pt <= 15.0) * (0.000) +
     533                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.2*tanh(pt*0.03 - 0.4)) +
     534                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.1*tanh(pt*0.03 - 0.4)) +
    523535                              (abs(eta) > 2.5)                                  * (0.000)}
    524536
    525537  # 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)) + \
     538  add EfficiencyFormula {5} {                                      (pt <= 15.0) * (0.000) +
     539                                                (abs(eta) <= 1.2) * (pt > 15.0) * (0.5*tanh(pt*0.03 - 0.4)) +
     540                              (abs(eta) > 1.2 && abs(eta) <= 2.5) * (pt > 15.0) * (0.4*tanh(pt*0.03 - 0.4)) +
    529541                              (abs(eta) > 2.5)                                  * (0.000)}
    530542}
     543
     544#############
     545# tau-tagging
     546#############
    531547
    532548module TauTagging TauTagging {
    … …  
    589605  add Branch ScalarHT/energy ScalarHT ScalarHT
    590606}
    591 
Note: See TracChangeset for help on using the changeset viewer.