Fork me on GitHub

Changeset 058006a in git


Ignore:
Timestamp:
Aug 26, 2016, 2:24:54 PM (8 years ago)
Author:
Alexandre Mertens <alexandre.mertens@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
f01b690
Parents:
ebf057e (diff), 3986893 (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 branch 'dev_01' of https://github.com/delphes/delphes into ValidationWorkWithPavel

Files:
1 added
1 deleted
2 edited

Legend:

Unmodified
Added
Removed
  • cards/delphes_card_CMS.tcl

    rebf057e r058006a  
    1515
    1616  TrackMerger
    17   Calorimeter
     17 
     18  ECal
     19  HCal
     20 
     21  TowerMerger
    1822  EFlowMerger
    1923
     
    3539 
    3640  FastJetFinder
    37   FastCaloJetFinder
    3841
    3942  JetEnergyScale
     
    193196}
    194197
     198
     199
    195200#############
    196 # Calorimeter
     201#   ECAL
    197202#############
    198203
    199 module Calorimeter Calorimeter {
     204module SimpleCalorimeter ECal {
    200205  set ParticleInputArray ParticlePropagator/stableParticles
    201206  set TrackInputArray TrackMerger/tracks
    202207
    203   set TowerOutputArray towers
    204   set PhotonOutputArray photons
    205 
     208  set TowerOutputArray ecalTowers
    206209  set EFlowTrackOutputArray eflowTracks
    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
     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
     305module 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
    215317
    216318  set SmearTowerCenter true
     
    250352
    251353  # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    252   add EnergyFraction {0} {0.0 1.0}
     354  add EnergyFraction {0} {1.0}
    253355  # energy fractions for e, gamma and pi0
    254   add EnergyFraction {11} {1.0 0.0}
    255   add EnergyFraction {22} {1.0 0.0}
    256   add EnergyFraction {111} {1.0 0.0}
     356  add EnergyFraction {11} {0.0}
     357  add EnergyFraction {22} {0.0}
     358  add EnergyFraction {111} {0.0}
    257359  # energy fractions for muon, neutrinos and neutralinos
    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}
     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}
    267369  # energy fractions for K0short and Lambda
    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)}
     370  add EnergyFraction {310} {0.7}
     371  add EnergyFraction {3122} {0.7}
    276372
    277373  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    278   set HCalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
     374  set ResolutionFormula {                      (abs(eta) <= 3.0) * sqrt(energy^2*0.050^2 + energy*1.50^2) +
    279375                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
    280 }
     376
     377}
     378
     379
     380#################
     381# Electron filter
     382#################
     383
     384module 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
     396module Merger TowerMerger {
     397# add InputArray InputArray
     398  add InputArray ECal/ecalTowers
     399  add InputArray HCal/hcalTowers
     400  set OutputArray towers
     401}
     402
     403
    281404
    282405####################
     
    286409module Merger EFlowMerger {
    287410# add InputArray InputArray
    288   add InputArray Calorimeter/eflowTracks
    289   add InputArray Calorimeter/eflowPhotons
    290   add InputArray Calorimeter/eflowNeutralHadrons
     411  add InputArray HCal/eflowTracks
     412  add InputArray ECal/eflowPhotons
     413  add InputArray HCal/eflowNeutralHadrons
    291414  set OutputArray eflow
    292415}
     
    297420
    298421module Efficiency PhotonEfficiency {
    299   set InputArray Calorimeter/eflowPhotons
     422  set InputArray ECal/eflowPhotons
    300423  set OutputArray photons
    301424
     
    326449}
    327450
    328 #################
    329 # Electron filter
    330 #################
    331 
    332 module PdgCodeFilter ElectronFilter {
    333   set InputArray Calorimeter/eflowTracks
    334   set OutputArray electrons
    335   set Invert true
    336   add PdgCode {11}
    337   add PdgCode {-11}
    338 }
    339451
    340452#####################
     
    498610}
    499611
    500 ################
    501 # caloJet finder
    502 ################
    503 
    504 module 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 
    516612##################
    517613# Jet Energy Scale
     
    617713
    618714  add Branch TrackMerger/tracks Track Track
    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
     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
    624720
    625721  add Branch GenJetFinder/jets GenJet Jet
    626722  add Branch GenMissingET/momentum GenMissingET MissingET
    627 
     723 
    628724  add Branch UniqueObjectFinder/jets Jet Jet
    629   add Branch FastCaloJetFinder/jets CaloJet Jet
    630725  add Branch UniqueObjectFinder/electrons Electron Electron
    631726  add Branch UniqueObjectFinder/photons Photon Photon
  • examples/validation.sh

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