Fork me on GitHub

Changes in / [058006a:ebf057e] in git


Ignore:
Files:
1 added
1 deleted
2 edited

Legend:

Unmodified
Added
Removed

  • cards/delphes_card_CMS.tcl

    r058006a rebf057e  
    1515
    1616  TrackMerger
    17  
    18   ECal
    19   HCal
    20  
    21   TowerMerger
     17  Calorimeter
    2218  EFlowMerger
    2319
    … …  
    3935 
    4036  FastJetFinder
     37  FastCaloJetFinder
    4138
    4239  JetEnergyScale
    … …  
    196193}
    197194
    198 
    199 
    200195#############
    201 #   ECAL
     196# Calorimeter
    202197#############
    203198
    204 module SimpleCalorimeter ECal {
     199module Calorimeter Calorimeter {
    205200  set ParticleInputArray ParticlePropagator/stableParticles
    206201  set TrackInputArray TrackMerger/tracks
    207202
    208   set TowerOutputArray ecalTowers
     203  set TowerOutputArray towers
     204  set PhotonOutputArray photons
     205
    209206  set EFlowTrackOutputArray eflowTracks
    210   set EFlowTowerOutputArray eflowPhotons
    211 
    212   set IsEcal true
    213 
    214   set EnergyMin 0.5
    215   set EnergySignificanceMin 2.0
    216 
    217   set SmearTowerCenter true
    218 
    219   set pi [expr {acos(-1)}]
    220 
    221   # lists of the edges of each tower in eta and phi
    222   # each list starts with the lower edge of the first tower
    223   # the list ends with the higher edged of the last tower
    224 
    225   # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
    226 
    227   set PhiBins {}
    228   for {set i -180} {$i <= 180} {incr i} {
    229     add PhiBins [expr {$i * $pi/180.0}]
    230   }
    231 
    232   # 0.02 unit in eta up to eta = 1.5 (barrel)
    233   for {set i -85} {$i <= 86} {incr i} {
    234     set eta [expr {$i * 0.0174}]
    235     add EtaPhiBins $eta $PhiBins
    236   }
    237 
    238   # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
    239 
    240   set PhiBins {}
    241   for {set i -180} {$i <= 180} {incr i} {
    242     add PhiBins [expr {$i * $pi/180.0}]
    243   }
    244 
    245   # 0.02 unit in eta up to eta = 3
    246   for {set i 1} {$i <= 84} {incr i} {
    247     set eta [expr { -2.958 + $i * 0.0174}]
    248     add EtaPhiBins $eta $PhiBins
    249   }
    250 
    251   for {set i 1} {$i <= 84} {incr i} {
    252     set eta [expr { 1.4964 + $i * 0.0174}]
    253     add EtaPhiBins $eta $PhiBins
    254   }
    255 
    256   # take present CMS granularity for HF
    257 
    258   # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
    259   set PhiBins {}
    260   for {set i -18} {$i <= 18} {incr i} {
    261     add PhiBins [expr {$i * $pi/18.0}]
    262   }
    263 
    264   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} {
    265     add EtaPhiBins $eta $PhiBins
    266   }
    267 
    268 
    269   add EnergyFraction {0} {0.0}
    270   # energy fractions for e, gamma and pi0
    271   add EnergyFraction {11} {1.0}
    272   add EnergyFraction {22} {1.0}
    273   add EnergyFraction {111} {1.0}
    274   # energy fractions for muon, neutrinos and neutralinos
    275   add EnergyFraction {12} {0.0}
    276   add EnergyFraction {13} {0.0}
    277   add EnergyFraction {14} {0.0}
    278   add EnergyFraction {16} {0.0}
    279   add EnergyFraction {1000022} {0.0}
    280   add EnergyFraction {1000023} {0.0}
    281   add EnergyFraction {1000025} {0.0}
    282   add EnergyFraction {1000035} {0.0}
    283   add EnergyFraction {1000045} {0.0}
    284   # energy fractions for K0short and Lambda
    285   add EnergyFraction {310} {0.3}
    286   add EnergyFraction {3122} {0.3}
    287 
    288   # set ResolutionFormula {resolution formula as a function of eta and energy}
    289 
    290   # for the ECAL barrel (|eta| < 1.5), see hep-ex/1306.2016 and 1502.02701
    291 
    292   # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    293   # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
    294   set ResolutionFormula {                      (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
    295                              (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) +
    296                              (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    297 
    298 }
    299 
    300 
    301 #############
    302 #   HCAL
    303 #############
    304 
    305 module SimpleCalorimeter HCal {
    306   set ParticleInputArray ParticlePropagator/stableParticles
    307   set TrackInputArray ECal/eflowTracks
    308 
    309   set TowerOutputArray hcalTowers
    310   set EFlowTrackOutputArray eflowTracks
    311   set EFlowTowerOutputArray eflowNeutralHadrons
    312 
    313   set IsEcal false
    314 
    315   set EnergyMin 1.0
    316   set EnergySignificanceMin 2.0
     207  set EFlowPhotonOutputArray eflowPhotons
     208  set EFlowNeutralHadronOutputArray eflowNeutralHadrons
     209
     210  set ECalEnergyMin 0.5
     211  set HCalEnergyMin 1.0
     212
     213  set ECalEnergySignificanceMin 2.0
     214  set HCalEnergySignificanceMin 2.0
    317215
    318216  set SmearTowerCenter true
    … …  
    352250
    353251  # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    354   add EnergyFraction {0} {1.0}
     252  add EnergyFraction {0} {0.0 1.0}
    355253  # energy fractions for e, gamma and pi0
    356   add EnergyFraction {11} {0.0}
    357   add EnergyFraction {22} {0.0}
    358   add EnergyFraction {111} {0.0}
     254  add EnergyFraction {11} {1.0 0.0}
     255  add EnergyFraction {22} {1.0 0.0}
     256  add EnergyFraction {111} {1.0 0.0}
    359257  # energy fractions for muon, neutrinos and neutralinos
    360   add EnergyFraction {12} {0.0}
    361   add EnergyFraction {13} {0.0}
    362   add EnergyFraction {14} {0.0}
    363   add EnergyFraction {16} {0.0}
    364   add EnergyFraction {1000022} {0.0}
    365   add EnergyFraction {1000023} {0.0}
    366   add EnergyFraction {1000025} {0.0}
    367   add EnergyFraction {1000035} {0.0}
    368   add EnergyFraction {1000045} {0.0}
     258  add EnergyFraction {12} {0.0 0.0}
     259  add EnergyFraction {13} {0.0 0.0}
     260  add EnergyFraction {14} {0.0 0.0}
     261  add EnergyFraction {16} {0.0 0.0}
     262  add EnergyFraction {1000022} {0.0 0.0}
     263  add EnergyFraction {1000023} {0.0 0.0}
     264  add EnergyFraction {1000025} {0.0 0.0}
     265  add EnergyFraction {1000035} {0.0 0.0}
     266  add EnergyFraction {1000045} {0.0 0.0}
    369267  # energy fractions for K0short and Lambda
    370   add EnergyFraction {310} {0.7}
    371   add EnergyFraction {3122} {0.7}
     268  add EnergyFraction {310} {0.3 0.7}
     269  add EnergyFraction {3122} {0.3 0.7}
     270
     271  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
     272  # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
     273  set ECalResolutionFormula {                  (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
     274                             (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) +
     275                             (abs(eta) > 2.5 && abs(eta) <= 5.0) * sqrt(energy^2*0.107^2 + energy*2.08^2)}
    372276
    373277  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    374   set ResolutionFormula {                      (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
     278  set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
    375279                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
    376 
    377 }
    378 
    379 
    380 #################
    381 # Electron filter
    382 #################
    383 
    384 module PdgCodeFilter ElectronFilter {
    385   set InputArray HCal/eflowTracks
    386   set OutputArray electrons
    387   set Invert true
    388   add PdgCode {11}
    389   add PdgCode {-11}
    390 }
    391 
    392 ###################################################
    393 # Tower Merger (in case not using e-flow algorithm)
    394 ###################################################
    395 
    396 module Merger TowerMerger {
    397 # add InputArray InputArray
    398   add InputArray ECal/ecalTowers
    399   add InputArray HCal/hcalTowers
    400   set OutputArray towers
    401 }
    402 
    403 
     280}
    404281
    405282####################
    … …  
    409286module Merger EFlowMerger {
    410287# add InputArray InputArray
    411   add InputArray HCal/eflowTracks
    412   add InputArray ECal/eflowPhotons
    413   add InputArray HCal/eflowNeutralHadrons
     288  add InputArray Calorimeter/eflowTracks
     289  add InputArray Calorimeter/eflowPhotons
     290  add InputArray Calorimeter/eflowNeutralHadrons
    414291  set OutputArray eflow
    415292}
    … …  
    420297
    421298module Efficiency PhotonEfficiency {
    422   set InputArray ECal/eflowPhotons
     299  set InputArray Calorimeter/eflowPhotons
    423300  set OutputArray photons
    424301
    … …  
    449326}
    450327
     328#################
     329# Electron filter
     330#################
     331
     332module PdgCodeFilter ElectronFilter {
     333  set InputArray Calorimeter/eflowTracks
     334  set OutputArray electrons
     335  set Invert true
     336  add PdgCode {11}
     337  add PdgCode {-11}
     338}
    451339
    452340#####################
    … …  
    610498}
    611499
     500################
     501# caloJet finder
     502################
     503
     504module FastJetFinder FastCaloJetFinder {
     505  set InputArray Calorimeter/towers
     506
     507  set OutputArray jets
     508
     509  # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt
     510  set JetAlgorithm 6
     511  set ParameterR 0.5
     512
     513  set JetPTMin 20.0
     514}
     515
    612516##################
    613517# Jet Energy Scale
    … …  
    713617
    714618  add Branch TrackMerger/tracks Track Track
    715   add Branch TowerMerger/towers Tower Tower
    716 
    717   add Branch HCal/eflowTracks EFlowTrack Track
    718   add Branch ECal/eflowPhotons EFlowPhoton Tower
    719   add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
     619  add Branch Calorimeter/towers Tower Tower
     620
     621  add Branch Calorimeter/eflowTracks EFlowTrack Track
     622  add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
     623  add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
    720624
    721625  add Branch GenJetFinder/jets GenJet Jet
    722626  add Branch GenMissingET/momentum GenMissingET MissingET
    723  
     627
    724628  add Branch UniqueObjectFinder/jets Jet Jet
     629  add Branch FastCaloJetFinder/jets CaloJet Jet
    725630  add Branch UniqueObjectFinder/electrons Electron Electron
    726631  add Branch UniqueObjectFinder/photons Photon Photon

  • examples/validation.sh

    r058006a rebf057e  
    55  cmnd="examples/Pythia8/configParticleGun_$name.cmnd"
    66  sed '/Main:spareMode1/s/=[[:space:]]*[0-9]*/= '$pid'/' examples/Pythia8/configParticleGun.cmnd > $cmnd
    7   ./DelphesPythia8 cards/validation_card.tcl $cmnd delphes_ParticleGun_$name.root
     7  ./DelphesPythia8 cards/delphes_card_CMS.tcl $cmnd delphes_ParticleGun_$name.root
    88}
    99
Note: See TracChangeset for help on using the changeset viewer.