Fork me on GitHub

Changeset a2983ec in git for cards/delphes_card_LHCb.tcl


Ignore:
Timestamp:
Dec 22, 2014, 11:38:33 AM (10 years ago)
Author:
Pavel Demin <pavel.demin@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
8c72f09
Parents:
4e09c3a
Message:

add parameter SmearTowerCenter to all cards

File:
1 edited

Legend:

Unmodified
Added
Removed

  • cards/delphes_card_LHCb.tcl

    r4e09c3a ra2983ec  
    88
    99set ExecutionPath {
    10  
     10
    1111  ParticlePropagator
    12  
     12
    1313  ChargedHadronMomentumSmearing
    1414  ElectronEnergySmearing
    … …  
    1717  TrackMerger
    1818  ImpactParameterSmearing
    19  
     19
    2020  IdentificationMap
    21  
     21
    2222  ECal
    2323  HCal
    … …  
    4141
    4242  # radius of the magnetic field coverage, in m
    43  
     43
    4444  set Radius 3.31
    45  
    46  
     45
     46
    4747  # half-length of the magnetic field coverage, in m
    4848  set HalfLength 12.0
    … …  
    5353  # Need to veto anything with theta > 0.269 rad  -> eta = 2
    5454  #                            theta < 0.0135 rad -> eta = 5
    55  
    56   # tracker and calos are at approx 0.269 rad, R = 12*tan(0.269) 
    57 
    58 }
    59 
    60  
     55
     56  # tracker and calos are at approx 0.269 rad, R = 12*tan(0.269)
     57
     58}
     59
     60
    6161########################################
    6262# Momentum resolution for charged tracks
    … …  
    8484
    8585  # resolution formula for electrons
    86   set ResolutionFormula { (eta > 2.0  && eta <= 5.0) * (energy > 0.1   && energy <= 8.0) * (energy*0.05) + \
     86  set ResolutionFormula { (eta > 2.0  && eta <= 5.0) * (energy > 0.1   && energy <= 8.0) * (energy*0.05) +
    8787                          (eta > 2.0  && eta <= 5.0) * (energy > 8.0)                    *  sqrt(energy^2*0.015^2 + energy*0.10^2)}
    8888  }
    … …  
    139139  set InputArray ImpactParameterSmearing/tracks
    140140  set OutputArray tracks
    141  
     141
    142142  # {PID in} {PID out} {formula}
    143143  # make sure "PID in" and "PID out" have the same charge (e.g {-13} {211} or {-321} {211})
    144144  # {211} {-13} is equivalent to {-211} {13} (and needs to be written once only...)
    145  
    146  
    147  
    148  
    149  
    150  
    151   # --- pions --- 
    152  
    153   add EfficiencyFormula {211} {211} {      (eta <= 2.0)                                  * (0.00) + \
    154                                            (eta > 2.0  && eta <= 5.0) *       (pt < 0.8) * (0.00) + \
    155                                            (eta > 2.0  && eta <= 5.0) *       (pt >= 0.8)* (0.95) + \
     145
     146
     147
     148
     149
     150
     151  # --- pions ---
     152
     153  add EfficiencyFormula {211} {211} {      (eta <= 2.0)                                  * (0.00) +
     154                                           (eta > 2.0  && eta <= 5.0) *       (pt < 0.8) * (0.00) +
     155                                           (eta > 2.0  && eta <= 5.0) *       (pt >= 0.8)* (0.95) +
    156156                                           (eta > 5.0)                                   * (0.00)}
    157  
    158   add EfficiencyFormula {211} {-13} {      (eta <= 2.0)                                 * (0.00) + \
    159                                            (eta > 2.0  && eta <= 5.0) *      (pt < 0.8) * (0.00) + \
    160                                            (eta > 2.0  && eta <= 5.0) *      (pt >= 0.8)* (0.005 + 0.0663*exp(-0.13*pt*cosh(eta))) + \
     157
     158  add EfficiencyFormula {211} {-13} {      (eta <= 2.0)                                 * (0.00) +
     159                                           (eta > 2.0  && eta <= 5.0) *      (pt < 0.8) * (0.00) +
     160                                           (eta > 2.0  && eta <= 5.0) *      (pt >= 0.8)* (0.005 + 0.0663*exp(-0.13*pt*cosh(eta))) +
    161161                                           (eta > 5.0)                                  * (0.00)}
    162162
    163  
     163
    164164 # --- kaons ---
    165  
    166  
    167   add EfficiencyFormula {321} {321} {      (eta <= 2.0)                                  * (0.00) + \
    168                                            (eta > 2.0  && eta <= 5.0) *       (pt < 0.8) * (0.00) + \
    169                                            (eta > 2.0  && eta <= 5.0) *       (pt >= 0.8)* (0.95) + \
     165
     166
     167  add EfficiencyFormula {321} {321} {      (eta <= 2.0)                                  * (0.00) +
     168                                           (eta > 2.0  && eta <= 5.0) *       (pt < 0.8) * (0.00) +
     169                                           (eta > 2.0  && eta <= 5.0) *       (pt >= 0.8)* (0.95) +
    170170                                           (eta > 5.0)                                   * (0.00)}
    171                                            
    172   add EfficiencyFormula {321} {-13} {      (eta <= 2.0)                                 * (0.00) + \
    173                                            (eta > 2.0  && eta <= 5.0) *      (pt < 0.8) * (0.00) + \
    174                                            (eta > 2.0  && eta <= 5.0) *      (pt >= 0.8)* (0.005 + 0.086*exp(-0.11*pt*cosh(eta))) + \
     171
     172  add EfficiencyFormula {321} {-13} {      (eta <= 2.0)                                 * (0.00) +
     173                                           (eta > 2.0  && eta <= 5.0) *      (pt < 0.8) * (0.00) +
     174                                           (eta > 2.0  && eta <= 5.0) *      (pt >= 0.8)* (0.005 + 0.086*exp(-0.11*pt*cosh(eta))) +
    175175                                           (eta > 5.0)                                  * (0.00)}
    176176
    177177
    178178 # --- protons ---
    179  
    180  
    181   add EfficiencyFormula {2212} {2212} {    (eta <= 2.0)                                  * (0.00) + \
    182                                            (eta > 2.0  && eta <= 5.0) *       (pt < 0.8) * (0.00) + \
    183                                            (eta > 2.0  && eta <= 5.0) *       (pt >= 0.8)* (0.95) + \
     179
     180
     181  add EfficiencyFormula {2212} {2212} {    (eta <= 2.0)                                  * (0.00) +
     182                                           (eta > 2.0  && eta <= 5.0) *       (pt < 0.8) * (0.00) +
     183                                           (eta > 2.0  && eta <= 5.0) *       (pt >= 0.8)* (0.95) +
    184184                                           (eta > 5.0)                                   * (0.00)}
    185                                            
    186   add EfficiencyFormula {2212} {-13} {     (eta <= 2.0)                                 * (0.00) + \
    187                                            (eta > 2.0  && eta <= 5.0) *      (pt < 0.8) * (0.00) + \
    188                                            (eta > 2.0  && eta <= 5.0) *      (pt >= 0.8)* (0.002) + \
     185
     186  add EfficiencyFormula {2212} {-13} {     (eta <= 2.0)                                 * (0.00) +
     187                                           (eta > 2.0  && eta <= 5.0) *      (pt < 0.8) * (0.00) +
     188                                           (eta > 2.0  && eta <= 5.0) *      (pt >= 0.8)* (0.002) +
    189189                                           (eta > 5.0)                                  * (0.00)}
    190190
    … …  
    195195
    196196
    197   add EfficiencyFormula {-13} {-13} {      (eta <= 2.0)                                * (0.00) + \
    198                                            (eta > 2.0  && eta <= 5.0) *      (pt < 0.8)* (0.00) + \
    199                                            (eta > 2.0  && eta <= 5.0) *     (pt >= 0.8)* (0.97) + \
     197  add EfficiencyFormula {-13} {-13} {      (eta <= 2.0)                                * (0.00) +
     198                                           (eta > 2.0  && eta <= 5.0) *      (pt < 0.8)* (0.00) +
     199                                           (eta > 2.0  && eta <= 5.0) *     (pt >= 0.8)* (0.97) +
    200200                                           (eta > 5.0)                                 * (0.00)}
    201201
    202  
     202
    203203 # efficiency for other charged particles (should be always {0} {0} {formula})
    204  
    205   add EfficiencyFormula {0} {0}     {      (eta <= 2.0)                                * (0.00) + \
    206                                            (eta > 2.0  && eta <= 5.0) *     (pt < 0.8) * (0.00) + \
    207                                            (eta > 2.0  && eta <= 5.0) *     (pt > 0.8) * (0.95) + \
     204
     205  add EfficiencyFormula {0} {0}     {      (eta <= 2.0)                                * (0.00) +
     206                                           (eta > 2.0  && eta <= 5.0) *     (pt < 0.8) * (0.00) +
     207                                           (eta > 2.0  && eta <= 5.0) *     (pt > 0.8) * (0.95) +
    208208                                           (eta > 5.0)                                 * (0.00)}
    209209
    … …  
    223223  set TowerOutputArray ecalTowers
    224224  set EFlowTowerOutputArray eflowPhotons
    225  
     225
     226  set EnergyMin 0.0
     227  set EnergySignificanceMin 0.0
     228
     229  set SmearTowerCenter true
     230
    226231  set pi [expr {acos(-1)}]
    227232
    … …  
    241246    add EtaPhiBins $eta $PhiBins
    242247  }
    243  
     248
    244249    # 1.25 degree towers
    245250  set PhiBins {}
    … …  
    285290  add EnergyFraction {3122} {0.3}
    286291
    287    set ResolutionFormula {(eta <= 5.0 && eta > 2.0) * sqrt(energy^2*0.015^2 + energy*0.10^2)}
     292  set ResolutionFormula {(eta <= 5.0 && eta > 2.0) * sqrt(energy^2*0.015^2 + energy*0.10^2)}
    288293}
    289294
    … …  
    298303  set TowerOutputArray hcalTowers
    299304  set EFlowTowerOutputArray eflowNeutralHadrons
    300  
     305
     306  set EnergyMin 0.0
     307  set EnergySignificanceMin 0.0
     308
     309  set SmearTowerCenter true
     310
    301311  set pi [expr {acos(-1)}]
    302312
    … …  
    330340
    331341
    332    
     342
    333343  # default energy fractions {abs(PDG code)} {Fecal Fhcal}
    334344  add EnergyFraction {0} {1.0}
    … …  
    352362
    353363  set ResolutionFormula { (eta <= 5.0 && eta > 2.0) * sqrt(energy^2*0.05^2 + energy*0.80^2)}
    354   }
     364}
    355365
    356366
    … …  
    360370
    361371# tracks, towers and eflow objects are not stored by default in the output.
    362 # if needed (for jet constituent or other studies), uncomment the relevant 
    363 # "add Branch ..." lines. 
     372# if needed (for jet constituent or other studies), uncomment the relevant
     373# "add Branch ..." lines.
    364374
    365375module TreeWriter TreeWriter {
    366376# add Branch InputArray BranchName BranchClass
    367  
     377
    368378  add Branch Delphes/allParticles Particle GenParticle
    369379
    … …  
    371381  add Branch HCal/eflowNeutralHadrons NeutralHadron Tower
    372382  add Branch ECal/eflowPhotons Photon Photon
    373  
    374 }
    375 
     383
     384}
     385
Note: See TracChangeset for help on using the changeset viewer.