- Timestamp:
- Apr 19, 2016, 4:17:00 PM (9 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 3110d6a
- Parents:
- 730bc64
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
cards/CMS_PhaseII_200PU.tcl
r730bc64 r16e0e46 341 341 342 342 set PhiBins {} 343 for {set i -1 57} {$i <= 157} {incr i} {344 add PhiBins [expr {$i * $pi/1 57.0}]343 for {set i -180} {$i <= 180} {incr i} { 344 add PhiBins [expr {$i * $pi/180.0}] 345 345 } 346 346 347 347 # 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.0 2}]348 for {set i -85} {$i <= 86} {incr i} { 349 set eta [expr {$i * 0.0174}] 350 350 add EtaPhiBins $eta $PhiBins 351 351 } 352 352 353 # assume 0.0 5 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) 354 354 355 355 set PhiBins {} 356 for {set i - 65} {$i <= 65} {incr i} {357 add PhiBins [expr {$i * $pi/ 65.0}]358 } 359 360 # 0.0 5unit in eta up to eta = 3361 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}] 363 363 add EtaPhiBins $eta $PhiBins 364 364 } 365 365 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}] 368 368 add EtaPhiBins $eta $PhiBins 369 369 } … … 377 377 } 378 378 379 foreach eta {-5 -4.7 -4.525 -4.35 -4.175 -4 -3.825 -3.65 -3.475 -3.3 -3.125 - 33.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} { 380 380 add EtaPhiBins $eta $PhiBins 381 381 } … … 406 406 # for the endcaps (1.5 < |eta| < 3.0), we take HGCAL see LHCC-P-008, Fig. 3.39, p.117 407 407 408 set ResolutionFormula { (abs(eta) <= 1.50) * sqrt(energy^2*0.00 3^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) + \ 409 409 (abs(eta) > 1.50 && abs(eta) <= 1.75) * sqrt(energy^2*0.006^2 + energy*0.20^2) + \ 410 410 (abs(eta) > 1.75 && abs(eta) <= 2.15) * sqrt(energy^2*0.007^2 + energy*0.21^2) + \ 411 411 (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.9 7^2)}412 (abs(eta) > 3.0 && abs(eta) <= 5.0) * sqrt(energy^2*0.08^2 + energy*1.98^2)} 413 413 414 414 } … … 500 500 501 501 # 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)} 505 505 506 506 }
Note:
See TracChangeset
for help on using the changeset viewer.