Fork me on GitHub

Changeset d091310 in git


Ignore:
Timestamp:
Aug 26, 2016, 5:47:59 PM (8 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
1d208a2
Parents:
8f0b34c
Message:

replace calo by simple calo in CMS cards

Location:
cards
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • cards/delphes_card_CMS_NoFastJet.tcl

    r8f0b34c rd091310  
    1515
    1616  TrackMerger
     17
     18  ECal
     19  HCal
     20
    1721  Calorimeter
    1822  EFlowMerger
     23
    1924}
    2025
    … …  
    167172
    168173#############
    169 # Calorimeter
     174#   ECAL
    170175#############
    171176
    172 module Calorimeter Calorimeter {
     177module SimpleCalorimeter ECal {
    173178  set ParticleInputArray ParticlePropagator/stableParticles
    174179  set TrackInputArray TrackMerger/tracks
    175180
    176   set TowerOutputArray towers
    177   set PhotonOutputArray photons
    178 
     181  set TowerOutputArray ecalTowers
    179182  set EFlowTrackOutputArray eflowTracks
    180   set EFlowPhotonOutputArray eflowPhotons
    181   set EFlowNeutralHadronOutputArray eflowNeutralHadrons
    182 
    183   set ECalEnergyMin 0.5
    184   set HCalEnergyMin 1.0
    185 
    186   set ECalEnergySignificanceMin 1.0
    187   set HCalEnergySignificanceMin 1.0
     183  set EFlowTowerOutputArray eflowPhotons
     184
     185  set IsEcal true
     186
     187  set EnergyMin 0.5
     188  set EnergySignificanceMin 2.0
    188189
    189190  set SmearTowerCenter true
    … …  
    195196  # the list ends with the higher edged of the last tower
    196197
     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
    197299  # 5 degrees towers
    198300  set PhiBins {}
    … …  
    223325
    224326  # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    225   add EnergyFraction {0} {0.0 1.0}
     327  add EnergyFraction {0} {1.0}
    226328  # energy fractions for e, gamma and pi0
    227   add EnergyFraction {11} {1.0 0.0}
    228   add EnergyFraction {22} {1.0 0.0}
    229   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}
    230332  # energy fractions for muon, neutrinos and neutralinos
    231   add EnergyFraction {12} {0.0 0.0}
    232   add EnergyFraction {13} {0.0 0.0}
    233   add EnergyFraction {14} {0.0 0.0}
    234   add EnergyFraction {16} {0.0 0.0}
    235   add EnergyFraction {1000022} {0.0 0.0}
    236   add EnergyFraction {1000023} {0.0 0.0}
    237   add EnergyFraction {1000025} {0.0 0.0}
    238   add EnergyFraction {1000035} {0.0 0.0}
    239   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}
    240342  # energy fractions for K0short and Lambda
    241   add EnergyFraction {310} {0.3 0.7}
    242   add EnergyFraction {3122} {0.3 0.7}
    243 
    244   # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    245   set ECalResolutionFormula {                  (abs(eta) <= 3.0) * sqrt(energy^2*0.007^2 + energy*0.07^2 + 0.35^2)  +
    246                              (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}
    247345
    248346  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    249   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) +
    250348                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
    251 }
     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
    252363
    253364####################
    … …  
    257368module Merger EFlowMerger {
    258369# add InputArray InputArray
    259   add InputArray Calorimeter/eflowTracks
    260   add InputArray Calorimeter/eflowPhotons
    261   add InputArray Calorimeter/eflowNeutralHadrons
     370  add InputArray HCal/eflowTracks
     371  add InputArray ECal/eflowPhotons
     372  add InputArray HCal/eflowNeutralHadrons
    262373  set OutputArray eflow
    263374}
     375

  • cards/delphes_card_CMS_PileUp.tcl

    r8f0b34c rd091310  
    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}
    … …  
    232235
    233236#############
    234 # Calorimeter
     237#   ECAL
    235238#############
    236239
    237 module Calorimeter Calorimeter {
     240module SimpleCalorimeter ECal {
    238241  set ParticleInputArray ParticlePropagator/stableParticles
    239242  set TrackInputArray TrackMerger/tracks
    240243
    241   set TowerOutputArray towers
    242   set PhotonOutputArray photons
    243 
     244  set TowerOutputArray ecalTowers
    244245  set EFlowTrackOutputArray eflowTracks
    245   set EFlowPhotonOutputArray eflowPhotons
    246   set EFlowNeutralHadronOutputArray eflowNeutralHadrons
    247 
    248   set ECalEnergyMin 0.5
    249   set HCalEnergyMin 1.0
    250 
    251   set ECalEnergySignificanceMin 1.0
    252   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
    253353
    254354  set SmearTowerCenter true
    … …  
    288388
    289389  # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    290   add EnergyFraction {0} {0.0 1.0}
     390  add EnergyFraction {0} {1.0}
    291391  # energy fractions for e, gamma and pi0
    292   add EnergyFraction {11} {1.0 0.0}
    293   add EnergyFraction {22} {1.0 0.0}
    294   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}
    295395  # energy fractions for muon, neutrinos and neutralinos
    296   add EnergyFraction {12} {0.0 0.0}
    297   add EnergyFraction {13} {0.0 0.0}
    298   add EnergyFraction {14} {0.0 0.0}
    299   add EnergyFraction {16} {0.0 0.0}
    300   add EnergyFraction {1000022} {0.0 0.0}
    301   add EnergyFraction {1000023} {0.0 0.0}
    302   add EnergyFraction {1000025} {0.0 0.0}
    303   add EnergyFraction {1000035} {0.0 0.0}
    304   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}
    305405  # energy fractions for K0short and Lambda
    306   add EnergyFraction {310} {0.3 0.7}
    307   add EnergyFraction {3122} {0.3 0.7}
    308 
    309   # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    310   # Eta shape from arXiv:1306.2016, Energy shape from arXiv:1502.02701
    311   set ECalResolutionFormula {                  (abs(eta) <= 1.5) * (1+0.64*eta^2) * sqrt(energy^2*0.008^2 + energy*0.11^2 + 0.40^2) +
    312                              (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) +
    313                              (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}
    314408
    315409  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
    316   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) +
    317411                             (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.130^2 + energy*2.70^2)}
     412
    318413}
    319414
    … …  
    323418
    324419module PdgCodeFilter ElectronFilter {
    325   set InputArray Calorimeter/eflowTracks
     420  set InputArray HCal/eflowTracks
    326421  set OutputArray electrons
    327422  set Invert true
    … …  
    330425}
    331426
     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
    332439##########################
    333440# Track pile-up subtractor
    … …  
    336443module TrackPileUpSubtractor TrackPileUpSubtractor {
    337444# add InputArray InputArray OutputArray
    338   add InputArray Calorimeter/eflowTracks eflowTracks
     445  add InputArray HCal/eflowTracks eflowTracks
    339446  add InputArray ElectronFilter/electrons electrons
    340447  add InputArray MuonMomentumSmearing/muons muons
    … …  
    346453}
    347454
     455
    348456####################
    349 # Neutral tower merger
     457# Neutral Tower merger
    350458####################
    351459
    352460module Merger NeutralTowerMerger {
    353461# add InputArray InputArray
    354   add InputArray Calorimeter/eflowPhotons
    355   add InputArray Calorimeter/eflowNeutralHadrons
    356   set OutputArray eflowTowers
    357 }
    358 
    359 ##################################
    360 # Energy flow merger (all tracks)
    361 ##################################
    362 
    363 module Merger EFlowMergerAllTracks {
    364 # add InputArray InputArray
    365   add InputArray TrackMerger/tracks
    366   add InputArray Calorimeter/eflowPhotons
    367   add InputArray Calorimeter/eflowNeutralHadrons
    368   set OutputArray eflow
     462  add InputArray ECal/eflowPhotons
     463  add InputArray HCal/eflowNeutralHadrons
     464  set OutputArray towers
    369465}
    370466
    … …  
    374470####################
    375471
     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
    376487module Merger EFlowMerger {
    377488# add InputArray InputArray
    378489  add InputArray TrackPileUpSubtractor/eflowTracks
    379   add InputArray Calorimeter/eflowPhotons
    380   add InputArray Calorimeter/eflowNeutralHadrons
     490  add InputArray ECal/eflowPhotons
     491  add InputArray HCal/eflowNeutralHadrons
    381492  set OutputArray eflow
    382493}
    383 
    384 
    385494
    386495#############
    … …  
    453562}
    454563
    455 
    456 
    457564############
    458565# Jet finder
    … …  
    481588module PileUpJetID PileUpJetID {
    482589  set JetInputArray FastJetFinder/jets
    483   set TrackInputArray Calorimeter/eflowTracks
    484   set NeutralInputArray NeutralTowerMerger/eflowTowers
     590  set TrackInputArray HCal/eflowTracks
     591  set NeutralInputArray NeutralTowerMerger/towers
    485592
    486593  set VertexInputArray PileUpMerger/vertices
    … …  
    527634
    528635module Efficiency PhotonEfficiency {
    529   set InputArray Calorimeter/eflowPhotons
     636  set InputArray ECal/eflowPhotons
    530637  set OutputArray photons
    531638
    … …  
    659766
    660767module JetFlavorAssociation JetFlavorAssociation {
    661  
     768
    662769  set PartonInputArray Delphes/partons
    663770  set ParticleInputArray Delphes/allParticles
    664771  set ParticleLHEFInputArray Delphes/allParticlesLHEF
    665772  set JetInputArray JetEnergyScale/jets
    666  
     773
    667774  set DeltaR 0.5
    668775  set PartonPTMin 1.0
    … …  
    685792
    686793  # based on arXiv:1211.4462
    687  
     794
    688795  # default efficiency formula (misidentification rate)
    689796  add EfficiencyFormula {0} {0.01+0.000038*pt}
Note: See TracChangeset for help on using the changeset viewer.