#include <HEDarkening.h>
Public Member Functions | |
float | degradation (float intlumi, int ieta, int lay) |
HEDarkening () | |
~HEDarkening () | |
Private Attributes | |
int | ieta_shift |
float | lumiscale [maxEta][maxLay] |
Definition at line 15 of file HEDarkening.h.
HEDarkening::HEDarkening | ( | ) |
Definition at line 13 of file HEDarkening.cc.
References i, ieta_shift, j, lumiscale, maxEta, and maxLay.
{ //HE starts at tower 16 ieta_shift = 16; //scale parameter in fb-1 for exponential darkening from radiation damage for each tile //fits from laser data for L1 and L7 //L0,Lm1 = L1, L2-L6 interpolated, L8-L17 extrapolated //from Vladimir Epshteyn float _lumiscale[maxEta][maxLay] = { {1194.9, 1194.9, 1194.9, 1651.5, 2282.7, 3155.2, 4361.0, 6027.8, 8331.5, 11515.7, 15916.8, 22000.0, 30408.2, 42029.8, 58093.1, 80295.6, 110983.5, 153400.1, 212027.7}, {952.8, 952.8, 952.8, 1293.9, 1757.1, 2386.1, 3240.3, 4400.3, 5975.4, 8114.5, 11019.3, 14963.9, 20320.6, 27594.9, 37473.2, 50887.7, 69104.3, 93841.9, 127435.0}, {759.8, 759.8, 759.8, 1013.8, 1352.5, 1804.5, 2407.6, 3212.2, 4285.7, 5717.9, 7628.7, 10178.1, 13579.5, 18117.6, 24172.3, 32250.4, 43028.0, 57407.4, 76592.2}, {605.9, 605.9, 605.9, 794.2, 1041.1, 1364.7, 1788.9, 2344.9, 3073.7, 4029.1, 5281.4, 6922.9, 9074.7, 11895.2, 15592.4, 20438.8, 26791.5, 35118.8, 46034.2}, {483.2, 483.2, 483.2, 622.3, 801.4, 1032.1, 1329.2, 1711.8, 2204.5, 2839.1, 3656.3, 4708.8, 6064.3, 7809.9, 10058.0, 12953.2, 16681.8, 21483.8, 27667.9}, {385.3, 385.3, 385.3, 487.5, 616.9, 780.5, 987.6, 1249.6, 1581.1, 2000.6, 2531.3, 3202.8, 4052.5, 5127.6, 6487.9, 8209.2, 10387.0, 13142.6, 16629.3}, {307.3, 307.3, 307.3, 382.0, 474.8, 590.3, 733.8, 912.2, 1134.0, 1409.7, 1752.4, 2178.5, 2708.1, 3366.6, 4185.1, 5202.6, 6467.5, 8039.9, 9994.7}, {245.0, 245.0, 245.0, 299.3, 365.5, 446.4, 545.2, 665.9, 813.3, 993.3, 1213.2, 1481.8, 1809.7, 2210.3, 2699.6, 3297.2, 4027.0, 4918.4, 6007.1}, {195.4, 195.4, 195.4, 234.5, 281.3, 337.6, 405.1, 486.1, 583.3, 700.0, 839.9, 1007.9, 1209.4, 1451.2, 1741.4, 2089.6, 2507.4, 3008.8, 3610.5}, {155.8, 155.8, 155.8, 183.7, 216.6, 255.3, 301.0, 354.9, 418.4, 493.2, 581.5, 685.5, 808.2, 952.8, 1123.3, 1324.3, 1561.3, 1840.6, 2170.0}, {124.3, 124.3, 124.3, 143.9, 166.7, 193.1, 223.6, 259.0, 300.1, 347.5, 402.6, 466.3, 540.1, 625.6, 724.6, 839.3, 972.1, 1126.0, 1304.2}, {99.1, 99.1, 99.1, 112.8, 128.3, 146.0, 166.2, 189.1, 215.2, 244.9, 278.7, 317.2, 360.9, 410.7, 467.4, 531.9, 605.3, 688.8, 783.9}, {79.0, 79.0, 79.0, 88.3, 98.8, 110.4, 123.5, 138.0, 154.3, 172.6, 192.9, 215.7, 241.2, 269.7, 301.5, 337.1, 376.9, 421.4, 471.1}, {63.0, 63.0, 63.0, 69.2, 76.0, 83.5, 91.7, 100.8, 110.7, 121.6, 133.6, 146.7, 161.2, 177.0, 194.5, 213.6, 234.7, 257.8, 283.2} }; //store array for(int j = 0; j < maxEta; j++){ for(int i = 0; i < maxLay; i++){ lumiscale[j][i] = _lumiscale[j][i]; } } }
HEDarkening::~HEDarkening | ( | ) |
Definition at line 47 of file HEDarkening.cc.
{ }
float HEDarkening::degradation | ( | float | intlumi, |
int | ieta, | ||
int | lay | ||
) |
Definition at line 49 of file HEDarkening.cc.
References create_public_lumi_plots::exp, ieta_shift, lumiscale, maxEta, and maxLay.
Referenced by HCalSD::getEnergyDeposit(), and HERecalibration::initialize().
{ //no lumi, no darkening if(intlumi <= 0) return 1.; //shift ieta tower index ieta -= ieta_shift; //if outside eta range, no darkening if(ieta < 0 || ieta >= maxEta) return 1.; //shift layer index by 1 to act as array index lay += 1; //if outside layer range, no darkening if(lay < 0 || lay >= maxLay) return 1.; //return darkening factor return (exp(-intlumi/lumiscale[ieta][lay])); }
int HEDarkening::ieta_shift [private] |
Definition at line 24 of file HEDarkening.h.
Referenced by degradation(), and HEDarkening().
float HEDarkening::lumiscale[maxEta][maxLay] [private] |
Definition at line 25 of file HEDarkening.h.
Referenced by degradation(), and HEDarkening().