Fork me on GitHub

Changeset e0a1a44 in git for cards/delphes_card_ATLAS.tcl


Ignore:
Timestamp:
Sep 30, 2016, 2:41:40 PM (8 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
4d128b9
Parents:
b4ec6ac
Message:

Separate ECAL and HCAL also for ATLAS card

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/delphes_card_ATLAS.tcl

    rb4ec6ac re0a1a44  
    44
    55set ExecutionPath {
    6   ParticlePropagator
     6 ParticlePropagator
    77
    88  ChargedHadronTrackingEfficiency
    … …  
    1515
    1616  TrackMerger
     17
     18  ECal
     19  HCal
     20
    1721  Calorimeter
    1822  EFlowMerger
    … …  
    3337  GenJetFinder
    3438  GenMissingET
    35  
     39
    3640  FastJetFinder
    3741
    3842  JetEnergyScale
    39  
     43
    4044  JetFlavorAssociation
    4145
    … …  
    189193}
    190194
     195
    191196#############
    192 # Calorimeter
     197#   ECAL
    193198#############
    194199
    195 module Calorimeter Calorimeter {
     200module SimpleCalorimeter ECal {
    196201  set ParticleInputArray ParticlePropagator/stableParticles
    197202  set TrackInputArray TrackMerger/tracks
    198203
    199   set TowerOutputArray towers
    200   set PhotonOutputArray photons
    201 
     204  set TowerOutputArray ecalTowers
    202205  set EFlowTrackOutputArray eflowTracks
    203   set EFlowPhotonOutputArray eflowPhotons
    204   set EFlowNeutralHadronOutputArray eflowNeutralHadrons
    205 
    206   set ECalEnergyMin 0.5
    207   set HCalEnergyMin 1.0
    208 
    209   set ECalEnergySignificanceMin 1.0
    210   set HCalEnergySignificanceMin 1.0
     206  set EFlowTowerOutputArray eflowPhotons
     207
     208  set IsEcal true
     209
     210  set EnergyMin 0.5
     211  set EnergySignificanceMin 2.0
    211212
    212213  set SmearTowerCenter true
    213214
    214215  set pi [expr {acos(-1)}]
     216
     217  # lists of the edges of each tower in eta and phi
     218  # each list starts with the lower edge of the first tower
     219  # the list ends with the higher edged of the last tower
     220
     221  # assume 0.02 x 0.02 resolution in eta,phi in the barrel |eta| < 1.5
     222
     223  set PhiBins {}
     224  for {set i -180} {$i <= 180} {incr i} {
     225    add PhiBins [expr {$i * $pi/180.0}]
     226  }
     227
     228  # 0.02 unit in eta up to eta = 1.5 (barrel)
     229  for {set i -85} {$i <= 86} {incr i} {
     230    set eta [expr {$i * 0.0174}]
     231    add EtaPhiBins $eta $PhiBins
     232  }
     233
     234  # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0
     235  set PhiBins {}
     236  for {set i -180} {$i <= 180} {incr i} {
     237    add PhiBins [expr {$i * $pi/180.0}]
     238  }
     239
     240  # 0.02 unit in eta up to eta = 3
     241  for {set i 1} {$i <= 84} {incr i} {
     242    set eta [expr { -2.958 + $i * 0.0174}]
     243    add EtaPhiBins $eta $PhiBins
     244  }
     245
     246  for {set i 1} {$i <= 84} {incr i} {
     247    set eta [expr { 1.4964 + $i * 0.0174}]
     248    add EtaPhiBins $eta $PhiBins
     249  }
     250
     251  # take present CMS granularity for HF
     252
     253  # 0.175 x (0.175 - 0.35) resolution in eta,phi in the HF 3.0 < |eta| < 5.0
     254  set PhiBins {}
     255  for {set i -18} {$i <= 18} {incr i} {
     256    add PhiBins [expr {$i * $pi/18.0}]
     257  }
     258
     259  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} {
     260    add EtaPhiBins $eta $PhiBins
     261  }
     262
     263
     264  add EnergyFraction {0} {0.0}
     265  # energy fractions for e, gamma and pi0
     266  add EnergyFraction {11} {1.0}
     267  add EnergyFraction {22} {1.0}
     268  add EnergyFraction {111} {1.0}
     269  # energy fractions for muon, neutrinos and neutralinos
     270  add EnergyFraction {12} {0.0}
     271  add EnergyFraction {13} {0.0}
     272  add EnergyFraction {14} {0.0}
     273  add EnergyFraction {16} {0.0}
     274  add EnergyFraction {1000022} {0.0}
     275  add EnergyFraction {1000023} {0.0}
     276  add EnergyFraction {1000025} {0.0}
     277  add EnergyFraction {1000035} {0.0}
     278  add EnergyFraction {1000045} {0.0}
     279  # energy fractions for K0short and Lambda
     280  add EnergyFraction {310} {0.3}
     281  add EnergyFraction {3122} {0.3}
     282
     283  # set ResolutionFormula {resolution formula as a function of eta and energy}
     284
     285  # set ECalResolutionFormula {resolution formula as a function of eta and energy}
     286  # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
     287  # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
     288  # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
     289  set ResolutionFormula {                      (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) +
     290                             (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
     291
     292
     293}
     294
     295
     296
     297#############
     298#   HCAL
     299#############
     300
     301module SimpleCalorimeter HCal {
     302  set ParticleInputArray ParticlePropagator/stableParticles
     303  set TrackInputArray ECal/eflowTracks
     304
     305  set TowerOutputArray hcalTowers
     306  set EFlowTrackOutputArray eflowTracks
     307  set EFlowTowerOutputArray eflowNeutralHadrons
     308
     309  set IsEcal false
     310
     311  set EnergyMin 1.0
     312  set EnergySignificanceMin 2.0
     313
     314  set SmearTowerCenter true
     315
     316 set pi [expr {acos(-1)}]
    215317
    216318  # lists of the edges of each tower in eta and phi
    … …  
    237339
    238340  # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    239   add EnergyFraction {0} {0.0 1.0}
     341  add EnergyFraction {0} {1.0}
    240342  # energy fractions for e, gamma and pi0
    241   add EnergyFraction {11} {1.0 0.0}
    242   add EnergyFraction {22} {1.0 0.0}
    243   add EnergyFraction {111} {1.0 0.0}
     343  add EnergyFraction {11} {0.0}
     344  add EnergyFraction {22} {0.0}
     345  add EnergyFraction {111} {0.0}
    244346  # energy fractions for muon, neutrinos and neutralinos
    245   add EnergyFraction {12} {0.0 0.0}
    246   add EnergyFraction {13} {0.0 0.0}
    247   add EnergyFraction {14} {0.0 0.0}
    248   add EnergyFraction {16} {0.0 0.0}
    249   add EnergyFraction {1000022} {0.0 0.0}
    250   add EnergyFraction {1000023} {0.0 0.0}
    251   add EnergyFraction {1000025} {0.0 0.0}
    252   add EnergyFraction {1000035} {0.0 0.0}
    253   add EnergyFraction {1000045} {0.0 0.0}
     347  add EnergyFraction {12} {0.0}
     348  add EnergyFraction {13} {0.0}
     349  add EnergyFraction {14} {0.0}
     350  add EnergyFraction {16} {0.0}
     351  add EnergyFraction {1000022} {0.0}
     352  add EnergyFraction {1000023} {0.0}
     353  add EnergyFraction {1000025} {0.0}
     354  add EnergyFraction {1000035} {0.0}
     355  add EnergyFraction {1000045} {0.0}
    254356  # energy fractions for K0short and Lambda
    255   add EnergyFraction {310} {0.3 0.7}
    256   add EnergyFraction {3122} {0.3 0.7}
    257 
    258   # set ECalResolutionFormula {resolution formula as a function of eta and energy}
    259   # http://arxiv.org/pdf/physics/0608012v1 jinst8_08_s08003
    260   # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
    261   # http://www.physics.utoronto.ca/~krieger/procs/ComoProceedings.pdf
    262   set ECalResolutionFormula {                  (abs(eta) <= 3.2) * sqrt(energy^2*0.0017^2 + energy*0.101^2) +
    263                              (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.0350^2 + energy*0.285^2)}
    264 
    265   # set HCalResolutionFormula {resolution formula as a function of eta and energy}
     357  add EnergyFraction {310} {0.7}
     358  add EnergyFraction {3122} {0.7}
     359
    266360  # http://arxiv.org/pdf/hep-ex/0004009v1
    267361  # http://villaolmo.mib.infn.it/ICATPP9th_2005/Calorimetry/Schram.p.pdf
    268   set HCalResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
     362  # set HCalResolutionFormula {resolution formula as a function of eta and energy}
     363  set ResolutionFormula {                      (abs(eta) <= 1.7) * sqrt(energy^2*0.0302^2 + energy*0.5205^2 + 1.59^2) +
    269364                             (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.0500^2 + energy*0.706^2) +
    270365                             (abs(eta) > 3.2 && abs(eta) <= 4.9) * sqrt(energy^2*0.09420^2 + energy*1.00^2)}
    271366}
     367
     368
     369#################
     370# Electron filter
     371#################
     372
     373module PdgCodeFilter ElectronFilter {
     374  set InputArray HCal/eflowTracks
     375  set OutputArray electrons
     376  set Invert true
     377  add PdgCode {11}
     378  add PdgCode {-11}
     379}
     380
     381
     382###################################################
     383# Tower Merger (in case not using e-flow algorithm)
     384###################################################
     385
     386module Merger Calorimeter {
     387# add InputArray InputArray
     388  add InputArray ECal/ecalTowers
     389  add InputArray HCal/hcalTowers
     390  set OutputArray towers
     391}
     392
     393
    272394
    273395####################
    … …  
    277399module Merger EFlowMerger {
    278400# add InputArray InputArray
    279   add InputArray Calorimeter/eflowTracks
    280   add InputArray Calorimeter/eflowPhotons
    281   add InputArray Calorimeter/eflowNeutralHadrons
     401  add InputArray HCal/eflowTracks
     402  add InputArray ECal/eflowPhotons
     403  add InputArray HCal/eflowNeutralHadrons
    282404  set OutputArray eflow
    283405}
    … …  
    288410
    289411module Efficiency PhotonEfficiency {
    290   set InputArray Calorimeter/eflowPhotons
     412  set InputArray ECal/eflowPhotons
    291413  set OutputArray photons
    292414
    … …  
    322444
    323445module PdgCodeFilter ElectronFilter {
    324   set InputArray Calorimeter/eflowTracks
     446  set InputArray HCal/eflowTracks
    325447  set OutputArray electrons
    326448  set Invert true
    … …  
    403525module Merger MissingET {
    404526# add InputArray InputArray
    405   add InputArray EFlowMerger/eflow
     527  add InputArray Calorimeter/towers
    406528  set MomentumOutputArray momentum
    407529}
    … …  
    503625
    504626module JetFlavorAssociation JetFlavorAssociation {
    505  
     627
    506628  set PartonInputArray Delphes/partons
    507629  set ParticleInputArray Delphes/allParticles
    508630  set ParticleLHEFInputArray Delphes/allParticlesLHEF
    509631  set JetInputArray JetEnergyScale/jets
    510  
     632
    511633  set DeltaR 0.5
    512634  set PartonPTMin 1.0
    … …  
    545667
    546668module TrackCountingTauTagging TauTagging {
    547  
     669
    548670  set ParticleInputArray Delphes/allParticles
    549671  set PartonInputArray Delphes/partons
    … …  
    555677
    556678  set TrackPTMin 1.0
    557  
     679
    558680  set TauPTMin 1.0
    559681  set TauEtaMax 2.5
    560682
    561   # instructions: {n-prongs} {eff} 
    562  
     683  # instructions: {n-prongs} {eff}
     684
    563685  # 1 - one prong efficiency
    564686  # 2 - two or more efficiency
    565687  # -1 - one prong mistag rate
    566688  # -2 - two or more mistag rate
    567  
     689
    568690  set BitNumber 0
    569  
     691
    570692  # taken from ATL-PHYS-PUB-2015-045 (medium working point)
    571693  add EfficiencyFormula {1} {0.70}
    … …  
    604726  add Branch Calorimeter/towers Tower Tower
    605727
    606   add Branch Calorimeter/eflowTracks EFlowTrack Track
    607   add Branch Calorimeter/eflowPhotons EFlowPhoton Tower
    608   add Branch Calorimeter/eflowNeutralHadrons EFlowNeutralHadron Tower
     728  add Branch HCal/eflowTracks EFlowTrack Track
     729  add Branch ECal/eflowPhotons EFlowPhoton Tower
     730  add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
    609731
    610732  add Branch GenJetFinder/jets GenJet Jet
Note: See TracChangeset for help on using the changeset viewer.