Fork me on GitHub

Changeset b9ae4c3 in git for cards


Ignore:
Timestamp:
Sep 2, 2016, 3:46:14 PM (8 years ago)
Author:
GitHub <noreply@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
fa068d3
Parents:
ec5e04b (diff), 23389ff (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.
git-author:
Pavel Demin <pavel-demin@…> (09/02/16 15:46:14)
git-committer:
GitHub <noreply@…> (09/02/16 15:46:14)
Message:

Merge pull request #32 from delphes/dev_01

Dev 01

Location:
cards
Files:
2 added
6 edited

Legend:

Unmodified
Added
Removed

  • cards/converter_card.tcl

    rec5e04b rb9ae4c3  
    1313module TreeWriter TreeWriter {
    1414# add Branch InputArray BranchName BranchClass
    15   add Branch Delphes/allParticles Particle GenParticle
     15  add Branch Delphes/stableParticles Particle GenParticle
    1616}
    1717

  • cards/delphes_card_ATLAS.tcl

    rec5e04b rb9ae4c3  
    141141
    142142  # resolution formula for charged hadrons
    143   set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.5e-4^2) +
    144                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*2.5e-4^2) +
    145                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*5.5e-4^2)}
     143  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
     144                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
     145                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
    146146}
    147147

  • cards/delphes_card_ATLAS_PileUp.tcl

    rec5e04b rb9ae4c3  
    177177  # resolution formula for charged hadrons
    178178  # based on arXiv:1405.6569
    179   set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.5e-4^2) +
    180                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*2.5e-4^2) +
    181                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*5.5e-4^2)}
     179  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
     180                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
     181                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
    182182}
    183183

  • cards/delphes_card_CMS.tcl

    rec5e04b rb9ae4c3  
    1515
    1616  TrackMerger
     17 
     18  ECal
     19  HCal
     20 
    1721  Calorimeter
    1822  EFlowMerger
    … …  
    124128  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
    125129                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
    126                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
     130                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e3) * (0.99) +
     131                                           (abs(eta) <= 1.5) * (pt > 1.0e3 )               * (0.99 * exp(0.5 - pt*5.0e-4)) +
     132
    127133                         (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) +
     134                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e3) * (0.98) +
     135                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3)                * (0.98 * exp(0.5 - pt*5.0e-4)) +
    129136                         (abs(eta) > 2.5)                                                  * (0.00)}
    130137}
    … …  
    142149  # resolution formula for charged hadrons
    143150  # 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)}
     151  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
     152                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
     153                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
    147154}
    148155
    … …  
    192199}
    193200
     201
     202
    194203#############
    195 # Calorimeter
     204#   ECAL
    196205#############
    197206
    198 module Calorimeter Calorimeter {
     207module SimpleCalorimeter ECal {
    199208  set ParticleInputArray ParticlePropagator/stableParticles
    200209  set TrackInputArray TrackMerger/tracks
    201210
    202   set TowerOutputArray towers
    203   set PhotonOutputArray photons
    204 
     211  set TowerOutputArray ecalTowers
    205212  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
     213  set EFlowTowerOutputArray eflowPhotons
     214
     215  set IsEcal true
     216
     217  set EnergyMin 0.5
     218  set EnergySignificanceMin 2.0
     219
     220  set SmearTowerCenter true
     221
     222  set pi [expr {acos(-1)}]
     223
     224  # lists of the edges of each tower in eta and phi
     225  # each list starts with the lower edge of the first tower
     226  # the list ends with the higher edged of the last tower
     227
     228  # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
     229
     230  set PhiBins {}
     231  for {set i -180} {$i <= 180} {incr i} {
     232    add PhiBins [expr {$i * $pi/180.0}]
     233  }
     234
     235  # 0.02 unit in eta up to eta = 1.5 (barrel)
     236  for {set i -85} {$i <= 86} {incr i} {
     237    set eta [expr {$i * 0.0174}]
     238    add EtaPhiBins $eta $PhiBins
     239  }
     240
     241  # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
     242
     243  set PhiBins {}
     244  for {set i -180} {$i <= 180} {incr i} {
     245    add PhiBins [expr {$i * $pi/180.0}]
     246  }
     247
     248  # 0.02 unit in eta up to eta = 3
     249  for {set i 1} {$i <= 84} {incr i} {
     250    set eta [expr { -2.958 + $i * 0.0174}]
     251    add EtaPhiBins $eta $PhiBins
     252  }
     253
     254  for {set i 1} {$i <= 84} {incr i} {
     255    set eta [expr { 1.4964 + $i * 0.0174}]
     256    add EtaPhiBins $eta $PhiBins
     257  }
     258
     259  # take present CMS granularity for HF
     260
     261  # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
     262  set PhiBins {}
     263  for {set i -18} {$i <= 18} {incr i} {
     264    add PhiBins [expr {$i * $pi/18.0}]
     265  }
     266
     267  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} {
     268    add EtaPhiBins $eta $PhiBins
     269  }
     270
     271
     272  add EnergyFraction {0} {0.0}
     273  # energy fractions for e, gamma and pi0
     274  add EnergyFraction {11} {1.0}
     275  add EnergyFraction {22} {1.0}
     276  add EnergyFraction {111} {1.0}
     277  # energy fractions for muon, neutrinos and neutralinos
     278  add EnergyFraction {12} {0.0}
     279  add EnergyFraction {13} {0.0}
     280  add EnergyFraction {14} {0.0}
     281  add EnergyFraction {16} {0.0}
     282  add EnergyFraction {1000022} {0.0}
     283  add EnergyFraction {1000023} {0.0}
     284  add EnergyFraction {1000025} {0.0}
     285  add EnergyFraction {1000035} {0.0}
     286  add EnergyFraction {1000045} {0.0}
     287  # energy fractions for K0short and Lambda
     288  add EnergyFraction {310} {0.3}
     289  add EnergyFraction {3122} {0.3}
     290
     291  # set ResolutionFormula {resolution formula as a function of eta and energy}
     292
     293  # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
     294
     295  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
     296  # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
     297  set ResolutionFormula {                      (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
     298                             (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) +
     299                             (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
     300
     301}
     302
     303
     304#############
     305#   HCAL
     306#############
     307
     308module SimpleCalorimeter HCal {
     309  set ParticleInputArray ParticlePropagator/stableParticles
     310  set TrackInputArray ECal/eflowTracks
     311
     312  set TowerOutputArray hcalTowers
     313  set EFlowTrackOutputArray eflowTracks
     314  set EFlowTowerOutputArray eflowNeutralHadrons
     315
     316  set IsEcal false
     317
     318  set EnergyMin 1.0
     319  set EnergySignificanceMin 2.0
    214320
    215321  set SmearTowerCenter true
    … …  
    249355
    250356  # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    251   add EnergyFraction {0} {0.0 1.0}
     357  add EnergyFraction {0} {1.0}
    252358  # 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}
     359  add EnergyFraction {11} {0.0}
     360  add EnergyFraction {22} {0.0}
     361  add EnergyFraction {111} {0.0}
    256362  # 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}
     363  add EnergyFraction {12} {0.0}
     364  add EnergyFraction {13} {0.0}
     365  add EnergyFraction {14} {0.0}
     366  add EnergyFraction {16} {0.0}
     367  add EnergyFraction {1000022} {0.0}
     368  add EnergyFraction {1000023} {0.0}
     369  add EnergyFraction {1000025} {0.0}
     370  add EnergyFraction {1000035} {0.0}
     371  add EnergyFraction {1000045} {0.0}
    266372  # 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)}
     373  add EnergyFraction {310} {0.7}
     374  add EnergyFraction {3122} {0.7}
    275375
    276376  # 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) +
     377  set ResolutionFormula {                      (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
    278378                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
    279 }
     379
     380}
     381
     382
     383#################
     384# Electron filter
     385#################
     386
     387module PdgCodeFilter ElectronFilter {
     388  set InputArray HCal/eflowTracks
     389  set OutputArray electrons
     390  set Invert true
     391  add PdgCode {11}
     392  add PdgCode {-11}
     393}
     394
     395###################################################
     396# Tower Merger (in case not using e-flow algorithm)
     397###################################################
     398
     399module Merger Calorimeter {
     400# add InputArray InputArray
     401  add InputArray ECal/ecalTowers
     402  add InputArray HCal/hcalTowers
     403  set OutputArray towers
     404}
     405
     406
    280407
    281408####################
    … …  
    285412module Merger EFlowMerger {
    286413# add InputArray InputArray
    287   add InputArray Calorimeter/eflowTracks
    288   add InputArray Calorimeter/eflowPhotons
    289   add InputArray Calorimeter/eflowNeutralHadrons
     414  add InputArray HCal/eflowTracks
     415  add InputArray ECal/eflowPhotons
     416  add InputArray HCal/eflowNeutralHadrons
    290417  set OutputArray eflow
    291418}
    … …  
    296423
    297424module Efficiency PhotonEfficiency {
    298   set InputArray Calorimeter/eflowPhotons
     425  set InputArray ECal/eflowPhotons
    299426  set OutputArray photons
    300427
    … …  
    325452}
    326453
    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 }
    338454
    339455#####################
    … …  
    382498
    383499  # 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)) +
     500  set EfficiencyFormula {                                     (pt <= 10.0)                * (0.00) +
     501                                           (abs(eta) <= 1.5) * (pt > 10.0)                * (0.95) +
     502                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0)                * (0.95) +
    389503                         (abs(eta) > 2.4)                                                 * (0.00)}
    390504}
    … …  
    602716  add Branch Calorimeter/towers Tower Tower
    603717
    604   add Branch Calorimeter/eflowTracks EFlowTrack Track
    605   add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
    606   add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
     718  add Branch HCal/eflowTracks EFlowTrack Track
     719  add Branch ECal/eflowPhotons EFlowPhoton Tower
     720  add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
    607721
    608722  add Branch GenJetFinder/jets GenJet Jet
    609723  add Branch GenMissingET/momentum GenMissingET MissingET
    610 
     724 
    611725  add Branch UniqueObjectFinder/jets Jet Jet
    612726  add Branch UniqueObjectFinder/electrons Electron Electron

  • cards/delphes_card_CMS_NoFastJet.tcl

    rec5e04b rb9ae4c3  
    1515
    1616  TrackMerger
     17
     18  ECal
     19  HCal
     20
    1721  Calorimeter
    1822  EFlowMerger
     23
    1924}
    2025
    … …  
    9398  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
    9499                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
    95                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
     100                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e3) * (0.99) +
     101                                           (abs(eta) <= 1.5) * (pt > 1.0e3 )               * (0.99 * exp(0.5 - pt*5.0e-4)) +
     102
    96103                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
    97                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) +
     104                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e3) * (0.98) +
     105                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3)                * (0.98 * exp(0.5 - pt*5.0e-4)) +
    98106                         (abs(eta) > 2.5)                                                  * (0.00)}
     107
    99108}
    100109
    … …  
    128137  # resolution formula for electrons
    129138  # based on arXiv:1405.6569
    130   set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
    131                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
    132                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
     139  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.03^2 + pt^2*1.3e-3^2) +
     140                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*1.7e-3^2) +
     141                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.15^2 + pt^2*3.1e-3^2)}
    133142}
    134143
    … …  
    144153
    145154  # resolution formula for muons
    146   set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*2.0e-4^2) +
    147                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.02^2 + pt^2*3.0e-4^2) +
    148                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.05^2 + pt^2*6.0e-4^2)}
    149 }
     155  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.0e-4^2) +
     156                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*1.5e-4^2) +
     157                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*3.5e-4^2)}
     158}
     159
    150160
    151161##############
    … …  
    162172
    163173#############
    164 # Calorimeter
     174#   ECAL
    165175#############
    166176
    167 module Calorimeter Calorimeter {
     177module SimpleCalorimeter ECal {
    168178  set ParticleInputArray ParticlePropagator/stableParticles
    169179  set TrackInputArray TrackMerger/tracks
    170180
    171   set TowerOutputArray towers
    172   set PhotonOutputArray photons
    173 
     181  set TowerOutputArray ecalTowers
    174182  set EFlowTrackOutputArray eflowTracks
    175   set EFlowPhotonOutputArray eflowPhotons
    176   set EFlowNeutralHadronOutputArray eflowNeutralHadrons
    177 
    178   set ECalEnergyMin 0.5
    179   set HCalEnergyMin 1.0
    180 
    181   set ECalEnergySignificanceMin 1.0
    182   set HCalEnergySignificanceMin 1.0
     183  set EFlowTowerOutputArray eflowPhotons
     184
     185  set IsEcal true
     186
     187  set EnergyMin 0.5
     188  set EnergySignificanceMin 2.0
    183189
    184190  set SmearTowerCenter true
    … …  
    190196  # the list ends with the higher edged of the last tower
    191197
     198  # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
     199
     200  set PhiBins {}
     201  for {set i -180} {$i <= 180} {incr i} {
     202    add PhiBins [expr {$i * $pi/180.0}]
     203  }
     204
     205  # 0.02 unit in eta up to eta = 1.5 (barrel)
     206  for {set i -85} {$i <= 86} {incr i} {
     207    set eta [expr {$i * 0.0174}]
     208    add EtaPhiBins $eta $PhiBins
     209  }
     210
     211  # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
     212
     213  set PhiBins {}
     214  for {set i -180} {$i <= 180} {incr i} {
     215    add PhiBins [expr {$i * $pi/180.0}]
     216  }
     217
     218  # 0.02 unit in eta up to eta = 3
     219  for {set i 1} {$i <= 84} {incr i} {
     220    set eta [expr { -2.958 + $i * 0.0174}]
     221    add EtaPhiBins $eta $PhiBins
     222  }
     223
     224  for {set i 1} {$i <= 84} {incr i} {
     225    set eta [expr { 1.4964 + $i * 0.0174}]
     226    add EtaPhiBins $eta $PhiBins
     227  }
     228
     229  # take present CMS granularity for HF
     230
     231  # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
     232  set PhiBins {}
     233  for {set i -18} {$i <= 18} {incr i} {
     234    add PhiBins [expr {$i * $pi/18.0}]
     235  }
     236
     237  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} {
     238    add EtaPhiBins $eta $PhiBins
     239  }
     240
     241
     242  add EnergyFraction {0} {0.0}
     243  # energy fractions for e, gamma and pi0
     244  add EnergyFraction {11} {1.0}
     245  add EnergyFraction {22} {1.0}
     246  add EnergyFraction {111} {1.0}
     247  # energy fractions for muon, neutrinos and neutralinos
     248  add EnergyFraction {12} {0.0}
     249  add EnergyFraction {13} {0.0}
     250  add EnergyFraction {14} {0.0}
     251  add EnergyFraction {16} {0.0}
     252  add EnergyFraction {1000022} {0.0}
     253  add EnergyFraction {1000023} {0.0}
     254  add EnergyFraction {1000025} {0.0}
     255  add EnergyFraction {1000035} {0.0}
     256  add EnergyFraction {1000045} {0.0}
     257  # energy fractions for K0short and Lambda
     258  add EnergyFraction {310} {0.3}
     259  add EnergyFraction {3122} {0.3}
     260
     261  # set ResolutionFormula {resolution formula as a function of eta and energy}
     262
     263  # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
     264
     265  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
     266  # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
     267  set ResolutionFormula {                      (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
     268                             (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) +
     269                             (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
     270
     271}
     272
     273
     274#############
     275#   HCAL
     276#############
     277
     278module SimpleCalorimeter HCal {
     279  set ParticleInputArray ParticlePropagator/stableParticles
     280  set TrackInputArray ECal/eflowTracks
     281
     282  set TowerOutputArray hcalTowers
     283  set EFlowTrackOutputArray eflowTracks
     284  set EFlowTowerOutputArray eflowNeutralHadrons
     285
     286  set IsEcal false
     287
     288  set EnergyMin 1.0
     289  set EnergySignificanceMin 2.0
     290
     291  set SmearTowerCenter true
     292
     293  set pi [expr {acos(-1)}]
     294
     295  # lists of the edges of each tower in eta and phi
     296  # each list starts with the lower edge of the first tower
     297  # the list ends with the higher edged of the last tower
     298
    192299  # 5 degrees towers
    193300  set PhiBins {}
    … …  
    218325
    219326  # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    220   add EnergyFraction {0} {0.0 1.0}
     327  add EnergyFraction {0} {1.0}
    221328  # energy fractions for e, gamma and pi0
    222   add EnergyFraction {11} {1.0 0.0}
    223   add EnergyFraction {22} {1.0 0.0}
    224   add EnergyFraction {111} {1.0 0.0}
     329  add EnergyFraction {11} {0.0}
     330  add EnergyFraction {22} {0.0}
     331  add EnergyFraction {111} {0.0}
    225332  # energy fractions for muon, neutrinos and neutralinos
    226   add EnergyFraction {12} {0.0 0.0}
    227   add EnergyFraction {13} {0.0 0.0}
    228   add EnergyFraction {14} {0.0 0.0}
    229   add EnergyFraction {16} {0.0 0.0}
    230   add EnergyFraction {1000022} {0.0 0.0}
    231   add EnergyFraction {1000023} {0.0 0.0}
    232   add EnergyFraction {1000025} {0.0 0.0}
    233   add EnergyFraction {1000035} {0.0 0.0}
    234   add EnergyFraction {1000045} {0.0 0.0}
     333  add EnergyFraction {12} {0.0}
     334  add EnergyFraction {13} {0.0}
     335  add EnergyFraction {14} {0.0}
     336  add EnergyFraction {16} {0.0}
     337  add EnergyFraction {1000022} {0.0}
     338  add EnergyFraction {1000023} {0.0}
     339  add EnergyFraction {1000025} {0.0}
     340  add EnergyFraction {1000035} {0.0}
     341  add EnergyFraction {1000045} {0.0}
    235342  # energy fractions for K0short and Lambda
    236   add EnergyFraction {310} {0.3 0.7}
    237   add EnergyFraction {3122} {0.3 0.7}
    238 
    239   # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    240   set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2)  +
    241                              (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
     343  add EnergyFraction {310} {0.7}
     344  add EnergyFraction {3122} {0.7}
    242345
    243346  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    244   set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
     347  set ResolutionFormula {                      (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
    245348                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
    246 }
     349
     350}
     351
     352###################################################
     353# Tower Merger (in case not using e-flow algorithm)
     354###################################################
     355
     356module Merger Calorimeter {
     357# add InputArray InputArray
     358  add InputArray ECal/ecalTowers
     359  add InputArray HCal/hcalTowers
     360  set OutputArray towers
     361}
     362
    247363
    248364####################
    … …  
    252368module Merger EFlowMerger {
    253369# add InputArray InputArray
    254   add InputArray Calorimeter/eflowTracks
    255   add InputArray Calorimeter/eflowPhotons
    256   add InputArray Calorimeter/eflowNeutralHadrons
     370  add InputArray HCal/eflowTracks
     371  add InputArray ECal/eflowPhotons
     372  add InputArray HCal/eflowNeutralHadrons
    257373  set OutputArray eflow
    258374}
     375

  • cards/delphes_card_CMS_PileUp.tcl

    rec5e04b rb9ae4c3  
    1717
    1818  TrackMerger
    19   Calorimeter
     19
     20  ECal
     21  HCal
     22
    2023  ElectronFilter
    2124  TrackPileUpSubtractor
    … …  
    4548
    4649  MissingET
    47  
     50
    4851  JetFlavorAssociation
    4952
    … …  
    8285  # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
    8386  set VertexDistributionFormula {exp(-(t^2/160e-12^2/2))*exp(-(z^2/0.053^2/2))}
    84  
     87
    8588
    8689}
    … …  
    160163  set EfficiencyFormula {                                                    (pt <= 0.1)   * (0.00) +
    161164                                           (abs(eta) <= 1.5) * (pt > 0.1   && pt <= 1.0)   * (0.75) +
    162                                            (abs(eta) <= 1.5) * (pt > 1.0)                  * (0.99) +
     165                                           (abs(eta) <= 1.5) * (pt > 1.0   && pt <= 1.0e3) * (0.99) +
     166                                           (abs(eta) <= 1.5) * (pt > 1.0e3 )               * (0.99 * exp(0.5 - pt*5.0e-4)) +
     167
    163168                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1   && pt <= 1.0)   * (0.70) +
    164                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0)                  * (0.98) +
     169                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0   && pt <= 1.0e3) * (0.98) +
     170                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 1.0e3)                * (0.98 * exp(0.5 - pt*5.0e-4)) +
    165171                         (abs(eta) > 2.5)                                                  * (0.00)}
    166172}
    … …  
    178184  # resolution formula for charged hadrons
    179185  # based on arXiv:1405.6569
    180   set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.01^2 + pt^2*1.5e-4^2) +
    181                          (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.015^2 + pt^2*2.5e-4^2) +
    182                          (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.025^2 + pt^2*5.5e-4^2)}
     186  set ResolutionFormula {                  (abs(eta) <= 0.5) * (pt > 0.1) * sqrt(0.06^2 + pt^2*1.3e-3^2) +
     187                         (abs(eta) > 0.5 && abs(eta) <= 1.5) * (pt > 0.1) * sqrt(0.10^2 + pt^2*1.7e-3^2) +
     188                         (abs(eta) > 1.5 && abs(eta) <= 2.5) * (pt > 0.1) * sqrt(0.25^2 + pt^2*3.1e-3^2)}
    183189}
    184190
    … …  
    229235
    230236#############
    231 # Calorimeter
     237#   ECAL
    232238#############
    233239
    234 module Calorimeter Calorimeter {
     240module SimpleCalorimeter ECal {
    235241  set ParticleInputArray ParticlePropagator/stableParticles
    236242  set TrackInputArray TrackMerger/tracks
    237243
    238   set TowerOutputArray towers
    239   set PhotonOutputArray photons
    240 
     244  set TowerOutputArray ecalTowers
    241245  set EFlowTrackOutputArray eflowTracks
    242   set EFlowPhotonOutputArray eflowPhotons
    243   set EFlowNeutralHadronOutputArray eflowNeutralHadrons
    244 
    245   set ECalEnergyMin 0.5
    246   set HCalEnergyMin 1.0
    247 
    248   set ECalEnergySignificanceMin 1.0
    249   set HCalEnergySignificanceMin 1.0
     246  set EFlowTowerOutputArray eflowPhotons
     247
     248  set IsEcal true
     249
     250  set EnergyMin 0.5
     251  set EnergySignificanceMin 2.0
     252
     253  set SmearTowerCenter true
     254
     255  set pi [expr {acos(-1)}]
     256
     257  # lists of the edges of each tower in eta and phi
     258  # each list starts with the lower edge of the first tower
     259  # the list ends with the higher edged of the last tower
     260
     261  # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
     262
     263  set PhiBins {}
     264  for {set i -180} {$i <= 180} {incr i} {
     265    add PhiBins [expr {$i * $pi/180.0}]
     266  }
     267
     268  # 0.02 unit in eta up to eta = 1.5 (barrel)
     269  for {set i -85} {$i <= 86} {incr i} {
     270    set eta [expr {$i * 0.0174}]
     271    add EtaPhiBins $eta $PhiBins
     272  }
     273
     274  # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
     275
     276  set PhiBins {}
     277  for {set i -180} {$i <= 180} {incr i} {
     278    add PhiBins [expr {$i * $pi/180.0}]
     279  }
     280
     281  # 0.02 unit in eta up to eta = 3
     282  for {set i 1} {$i <= 84} {incr i} {
     283    set eta [expr { -2.958 + $i * 0.0174}]
     284    add EtaPhiBins $eta $PhiBins
     285  }
     286
     287  for {set i 1} {$i <= 84} {incr i} {
     288    set eta [expr { 1.4964 + $i * 0.0174}]
     289    add EtaPhiBins $eta $PhiBins
     290  }
     291
     292  # take present CMS granularity for HF
     293
     294  # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
     295  set PhiBins {}
     296  for {set i -18} {$i <= 18} {incr i} {
     297    add PhiBins [expr {$i * $pi/18.0}]
     298  }
     299
     300  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} {
     301    add EtaPhiBins $eta $PhiBins
     302  }
     303
     304
     305  add EnergyFraction {0} {0.0}
     306  # energy fractions for e, gamma and pi0
     307  add EnergyFraction {11} {1.0}
     308  add EnergyFraction {22} {1.0}
     309  add EnergyFraction {111} {1.0}
     310  # energy fractions for muon, neutrinos and neutralinos
     311  add EnergyFraction {12} {0.0}
     312  add EnergyFraction {13} {0.0}
     313  add EnergyFraction {14} {0.0}
     314  add EnergyFraction {16} {0.0}
     315  add EnergyFraction {1000022} {0.0}
     316  add EnergyFraction {1000023} {0.0}
     317  add EnergyFraction {1000025} {0.0}
     318  add EnergyFraction {1000035} {0.0}
     319  add EnergyFraction {1000045} {0.0}
     320  # energy fractions for K0short and Lambda
     321  add EnergyFraction {310} {0.3}
     322  add EnergyFraction {3122} {0.3}
     323
     324  # set ResolutionFormula {resolution formula as a function of eta and energy}
     325
     326  # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
     327
     328  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
     329  # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
     330  set ResolutionFormula {                      (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
     331                             (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) +
     332                             (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
     333
     334}
     335
     336
     337#############
     338#   HCAL
     339#############
     340
     341module SimpleCalorimeter HCal {
     342  set ParticleInputArray ParticlePropagator/stableParticles
     343  set TrackInputArray ECal/eflowTracks
     344
     345  set TowerOutputArray hcalTowers
     346  set EFlowTrackOutputArray eflowTracks
     347  set EFlowTowerOutputArray eflowNeutralHadrons
     348
     349  set IsEcal false
     350
     351  set EnergyMin 1.0
     352  set EnergySignificanceMin 2.0
    250353
    251354  set SmearTowerCenter true
    … …  
    285388
    286389  # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    287   add EnergyFraction {0} {0.0 1.0}
     390  add EnergyFraction {0} {1.0}
    288391  # energy fractions for e, gamma and pi0
    289   add EnergyFraction {11} {1.0 0.0}
    290   add EnergyFraction {22} {1.0 0.0}
    291   add EnergyFraction {111} {1.0 0.0}
     392  add EnergyFraction {11} {0.0}
     393  add EnergyFraction {22} {0.0}
     394  add EnergyFraction {111} {0.0}
    292395  # energy fractions for muon, neutrinos and neutralinos
    293   add EnergyFraction {12} {0.0 0.0}
    294   add EnergyFraction {13} {0.0 0.0}
    295   add EnergyFraction {14} {0.0 0.0}
    296   add EnergyFraction {16} {0.0 0.0}
    297   add EnergyFraction {1000022} {0.0 0.0}
    298   add EnergyFraction {1000023} {0.0 0.0}
    299   add EnergyFraction {1000025} {0.0 0.0}
    300   add EnergyFraction {1000035} {0.0 0.0}
    301   add EnergyFraction {1000045} {0.0 0.0}
     396  add EnergyFraction {12} {0.0}
     397  add EnergyFraction {13} {0.0}
     398  add EnergyFraction {14} {0.0}
     399  add EnergyFraction {16} {0.0}
     400  add EnergyFraction {1000022} {0.0}
     401  add EnergyFraction {1000023} {0.0}
     402  add EnergyFraction {1000025} {0.0}
     403  add EnergyFraction {1000035} {0.0}
     404  add EnergyFraction {1000045} {0.0}
    302405  # energy fractions for K0short and Lambda
    303   add EnergyFraction {310} {0.3 0.7}
    304   add EnergyFraction {3122} {0.3 0.7}
    305 
    306   # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    307   # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
    308   set ECalResolutionFormula {                  (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
    309                              (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) +
    310                              (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
     406  add EnergyFraction {310} {0.7}
     407  add EnergyFraction {3122} {0.7}
    311408
    312409  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    313   set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
     410  set ResolutionFormula {                      (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
    314411                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
     412
    315413}
    316414
    … …  
    320418
    321419module PdgCodeFilter ElectronFilter {
    322   set InputArray Calorimeter/eflowTracks
     420  set InputArray HCal/eflowTracks
    323421  set OutputArray electrons
    324422  set Invert true
    … …  
    327425}
    328426
     427###################################################
     428# Tower Merger (in case not using e-flow algorithm)
     429###################################################
     430
     431module Merger Calorimeter {
     432# add InputArray InputArray
     433  add InputArray ECal/ecalTowers
     434  add InputArray HCal/hcalTowers
     435  set OutputArray towers
     436}
     437
     438
    329439##########################
    330440# Track pile-up subtractor
    … …  
    333443module TrackPileUpSubtractor TrackPileUpSubtractor {
    334444# add InputArray InputArray OutputArray
    335   add InputArray Calorimeter/eflowTracks eflowTracks
     445  add InputArray HCal/eflowTracks eflowTracks
    336446  add InputArray ElectronFilter/electrons electrons
    337447  add InputArray MuonMomentumSmearing/muons muons
    … …  
    343453}
    344454
     455
    345456####################
    346 # Neutral tower merger
     457# Neutral Tower merger
    347458####################
    348459
    349460module Merger NeutralTowerMerger {
    350461# add InputArray InputArray
    351   add InputArray Calorimeter/eflowPhotons
    352   add InputArray Calorimeter/eflowNeutralHadrons
    353   set OutputArray eflowTowers
    354 }
    355 
    356 ##################################
    357 # Energy flow merger (all tracks)
    358 ##################################
    359 
    360 module Merger EFlowMergerAllTracks {
    361 # add InputArray InputArray
    362   add InputArray TrackMerger/tracks
    363   add InputArray Calorimeter/eflowPhotons
    364   add InputArray Calorimeter/eflowNeutralHadrons
    365   set OutputArray eflow
     462  add InputArray ECal/eflowPhotons
     463  add InputArray HCal/eflowNeutralHadrons
     464  set OutputArray towers
    366465}
    367466
    … …  
    371470####################
    372471
     472module Merger EFlowMergerAllTracks {
     473# add InputArray InputArray
     474  add InputArray HCal/eflowTracks
     475  add InputArray ECal/eflowPhotons
     476  add InputArray HCal/eflowNeutralHadrons
     477  set OutputArray eflow
     478}
     479
     480
     481
     482
     483####################
     484# Energy flow merger
     485####################
     486
    373487module Merger EFlowMerger {
    374488# add InputArray InputArray
    375489  add InputArray TrackPileUpSubtractor/eflowTracks
    376   add InputArray Calorimeter/eflowPhotons
    377   add InputArray Calorimeter/eflowNeutralHadrons
     490  add InputArray ECal/eflowPhotons
     491  add InputArray HCal/eflowNeutralHadrons
    378492  set OutputArray eflow
    379493}
    380 
    381 
    382494
    383495#############
    … …  
    450562}
    451563
    452 
    453 
    454564############
    455565# Jet finder
    … …  
    478588module PileUpJetID PileUpJetID {
    479589  set JetInputArray FastJetFinder/jets
    480   set TrackInputArray Calorimeter/eflowTracks
    481   set NeutralInputArray NeutralTowerMerger/eflowTowers
     590  set TrackInputArray HCal/eflowTracks
     591  set NeutralInputArray NeutralTowerMerger/towers
    482592
    483593  set VertexInputArray PileUpMerger/vertices
    … …  
    524634
    525635module Efficiency PhotonEfficiency {
    526   set InputArray Calorimeter/eflowPhotons
     636  set InputArray ECal/eflowPhotons
    527637  set OutputArray photons
    528638
    … …  
    601711
    602712  # efficiency formula for muons
    603   set EfficiencyFormula {                                      (pt <= 10.0)               * (0.00) +
    604                                            (abs(eta) <= 1.5) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
    605                                            (abs(eta) <= 1.5) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) +
    606                          (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0 && pt <= 1.0e3) * (0.95) +
    607                          (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 1.0e3)               * (0.95 * exp(0.5 - pt*5.0e-4)) +
     713  set EfficiencyFormula {                                     (pt <= 10.0)                * (0.00) +
     714                                           (abs(eta) <= 1.5) * (pt > 10.0)                * (0.95) +
     715                         (abs(eta) > 1.5 && abs(eta) <= 2.4) * (pt > 10.0)                * (0.95) +
    608716                         (abs(eta) > 2.4)                                                 * (0.00)}
     717
    609718}
    610719
    … …  
    657766
    658767module JetFlavorAssociation JetFlavorAssociation {
    659  
     768
    660769  set PartonInputArray Delphes/partons
    661770  set ParticleInputArray Delphes/allParticles
    662771  set ParticleLHEFInputArray Delphes/allParticlesLHEF
    663772  set JetInputArray JetEnergyScale/jets
    664  
     773
    665774  set DeltaR 0.5
    666775  set PartonPTMin 1.0
    … …  
    683792
    684793  # based on arXiv:1211.4462
    685  
     794
    686795  # default efficiency formula (misidentification rate)
    687796  add EfficiencyFormula {0} {0.01+0.000038*pt}
Note: See TracChangeset for help on using the changeset viewer.