Fork me on GitHub

Changeset 16e0e46 in git for cards


Ignore:
Timestamp:
Apr 19, 2016, 4:17:00 PM (9 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
3110d6a
Parents:
730bc64
Message:

updated calo parameters

File:
1 edited

Legend:

Unmodified
Added
Removed
  • cards/CMS_PhaseII_200PU.tcl

    r730bc64 r16e0e46  
    341341 
    342342  set PhiBins {}
    343   for {set i -157} {$i <= 157} {incr i} {
    344     add PhiBins [expr {$i * $pi/157.0}]
     343  for {set i -180} {$i <= 180} {incr i} {
     344    add PhiBins [expr {$i * $pi/180.0}]
    345345  }
    346346
    347347  # 0.02 unit in eta up to eta = 1.5 (barrel)
    348   for {set i -74} {$i <= 75} {incr i} {
    349     set eta [expr {$i * 0.02}]
     348  for {set i -85} {$i <= 86} {incr i} {
     349    set eta [expr {$i * 0.0174}]
    350350    add EtaPhiBins $eta $PhiBins
    351351  }
    352352
    353   # assume 0.05 x 0.05 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL)
     353  # assume 0.02 x 0.02 resolution in eta,phi in the endcaps 1.5 < |eta| < 3.0 (HGCAL- ECAL)
    354354 
    355355  set PhiBins {}
    356   for {set i -65} {$i <= 65} {incr i} {
    357     add PhiBins [expr {$i * $pi/65.0}]
    358   }
    359 
    360   # 0.05 unit in eta up to eta = 3
    361   for {set i 1} {$i <= 30} {incr i} {
    362     set eta [expr { -3.00 + $i * 0.05}]
     356  for {set i -180} {$i <= 180} {incr i} {
     357    add PhiBins [expr {$i * $pi/180.0}]
     358  }
     359
     360  # 0.02 unit in eta up to eta = 3
     361  for {set i 1} {$i <= 84} {incr i} {
     362    set eta [expr { -2.958 + $i * 0.0174}]
    363363    add EtaPhiBins $eta $PhiBins
    364364  }
    365365
    366   for {set i 1} {$i <= 30} {incr i} {
    367     set eta [expr { 1.5 + $i * 0.05}]
     366  for {set i 1} {$i <= 84} {incr i} {
     367    set eta [expr { 1.4964 + $i * 0.0174}]
    368368    add EtaPhiBins $eta $PhiBins
    369369  }
     
    377377  }
    378378 
    379   foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 -3 3.125 3.3 3.475 3.65 3.825 4 4.175 4.35 4.525 4.7 5} {
     379  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} {
    380380    add EtaPhiBins $eta $PhiBins
    381381  }
     
    406406  # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL  see LHCC-P-008, Fig. 3.39, p.117
    407407
    408   set ResolutionFormula {                     (abs(eta) <= 1.50) * sqrt(energy^2*0.003^2 + energy*0.028^2 + 0.12^2) + \
     408  set ResolutionFormula {                     (abs(eta) <= 1.50) * sqrt(energy^2*0.005^2 + energy*0.027^2 + 0.15^2) + \
    409409                           (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \
    410410                           (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \
    411411                           (abs(eta) > 2.15 && abs(eta) <= 3.00) * sqrt(energy^2*0.008^2 + energy*0.24^2) + \
    412                            (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.97^2)}
     412                           (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)}
    413413
    414414}
     
    500500
    501501# set ResolutionFormula {resolution formula as a function of eta and energy}
    502   set ResolutionFormula {                  (abs(eta) <= 1.7) * sqrt(energy^2*0.03^2 + energy*0.52^2 + 1.59^2) + \
    503                                                    (abs(eta) > 1.7 && abs(eta) <= 3.2) * sqrt(energy^2*0.05^2 + energy*0.706^2) + \
    504                                                    (abs(eta) > 3.0 && abs(eta) <= 4.9) * sqrt(energy^2*0.05^2 + energy*1.00^2)}
     502  set ResolutionFormula {                    (abs(eta) <= 1.5) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
     503                                                   (abs(eta) > 1.5 && abs(eta) <= 3.0) * sqrt(energy^2*0.05^2 + energy*1.00^2) + \
     504                                                   (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.11^2 + energy*2.80^2)}
    505505
    506506}
Note: See TracChangeset for help on using the changeset viewer.