|
|
|
#include <HFDarkening.h>
Public Member Functions | |
| float | degradation (float mrad) |
| float | dose (int layer, float radius) |
| HFDarkening () | |
| float | int_lumi (float intlumi) |
| ~HFDarkening () | |
Private Attributes | |
| float | dose_layer_depth [10][10] |
| float | radius [10] |
Definition at line 4 of file HFDarkening.h.
| HFDarkening::HFDarkening | ( | ) |
Definition at line 12 of file HFDarkening.cc.
References dose_layer_depth, i, and radius.
{
//HF area of consideration is 1100cm from interaction point to 1300cm in z-axis
//Radius (cm) - 12.5 cm from Beam pipe up to 170cm to top of PMT robox.
//Dose in MRad
for(int i = 0; i < 10; ++i) radius[i] = 12.5+i*15.75;
//Depth 0 - 1100cm-1120 - calculated @ 1110cm
dose_layer_depth[0][0] = 45.8; dose_layer_depth[0][1] = 39.5;
dose_layer_depth[0][2] = 9.97; dose_layer_depth[0][3] = 6.46;
dose_layer_depth[0][4] = 3.32; dose_layer_depth[0][5] = 2.21;
dose_layer_depth[0][6] = 0.737; dose_layer_depth[0][7] = 0.0557;
dose_layer_depth[0][8] = 0.00734; dose_layer_depth[0][9] = 0.00212;
//Depth 1 - 1120-1140 - @1130
dose_layer_depth[1][0] = 45.8; dose_layer_depth[1][1] = 39.5;
dose_layer_depth[1][2] = 9.97; dose_layer_depth[1][3] = 6.46;
dose_layer_depth[1][4] = 3.32; dose_layer_depth[1][5] = 2.21;
dose_layer_depth[1][6] = 0.737; dose_layer_depth[1][7] = 0.0557;
dose_layer_depth[1][8] = 0.00734; dose_layer_depth[1][9] = 0.00212;
//Depth 2 - 1140-1160 - @1150
dose_layer_depth[2][0] = 74.8; dose_layer_depth[2][1] = 42.3;
dose_layer_depth[2][2] = 8.17; dose_layer_depth[2][3] = 5.12;
dose_layer_depth[2][4] = 2.21; dose_layer_depth[2][5] = 1.48;
dose_layer_depth[2][6] = 0.400; dose_layer_depth[2][7] = 0.0388;
dose_layer_depth[2][8] = 0.0050; dose_layer_depth[2][9] = 0.00202;
//Depth 3 - 1160-1180 - @1170
dose_layer_depth[3][0] = 73.9; dose_layer_depth[3][1] = 18.1;
dose_layer_depth[3][2] = 3.81; dose_layer_depth[3][3] = 1.95;
dose_layer_depth[3][4] = 0.839; dose_layer_depth[3][5] = 0.559;
dose_layer_depth[3][6] = 0.200; dose_layer_depth[3][7] = 0.0386;
dose_layer_depth[3][8] = 0.00582; dose_layer_depth[3][9] = 0.00216;
//Depth 4 - 1180-1200 - @1190
dose_layer_depth[4][0] = 68.0; dose_layer_depth[4][1] = 9.75;
dose_layer_depth[4][2] = 2.19; dose_layer_depth[4][3] = 0.974;
dose_layer_depth[4][4] = 0.426; dose_layer_depth[4][5] = 0.265;
dose_layer_depth[4][6] = 0.124; dose_layer_depth[4][7] = 0.0258;
dose_layer_depth[4][8] = 0.00565; dose_layer_depth[4][9] = 0.00206;
//Depth 5 - 1200-1220 - @1210
dose_layer_depth[5][0] = 73.1; dose_layer_depth[5][1] = 5.42;
dose_layer_depth[5][2] = 1.24; dose_layer_depth[5][3] = 0.507;
dose_layer_depth[5][4] = 0.21; dose_layer_depth[5][5] = 0.129;
dose_layer_depth[5][6] = 0.0625; dose_layer_depth[5][7] = 0.0220;
dose_layer_depth[5][8] = 0.00404; dose_layer_depth[5][9] = 0.00156;
//Depth 6 - 1220-1240 - @1230
dose_layer_depth[6][0] = 66.1; dose_layer_depth[6][1] = 2.79;
dose_layer_depth[6][2] = 0.650; dose_layer_depth[6][3] = 0.257;
dose_layer_depth[6][4] = 0.111; dose_layer_depth[6][5] = 0.0624;
dose_layer_depth[6][6] = 0.0335; dose_layer_depth[6][7] = 0.0134;
dose_layer_depth[6][8] = 0.00280; dose_layer_depth[6][9] = 0.00127;
//Depth 7 - 1240-1260 - @1250
dose_layer_depth[7][0] = 68.7; dose_layer_depth[7][1] = 1.69;
dose_layer_depth[7][2] = 0.350; dose_layer_depth[7][3] = 0.137;
dose_layer_depth[7][4] = 0.0614; dose_layer_depth[7][5] = 0.0321;
dose_layer_depth[7][6] = 0.0167; dose_layer_depth[7][7] = 0.00621;
dose_layer_depth[7][8] = 0.00187; dose_layer_depth[7][9] = 0.000867;
//Depth 8 - 1260-1280 - @1270
dose_layer_depth[8][0] = 62.5; dose_layer_depth[8][1] = 2.39;
dose_layer_depth[8][2] = 0.252; dose_layer_depth[8][3] = 0.135;
dose_layer_depth[8][4] = 0.0505; dose_layer_depth[8][5] = 0.0295;
dose_layer_depth[8][6] = 0.0159; dose_layer_depth[8][7] = 0.0131;
dose_layer_depth[8][8] = 0.00803; dose_layer_depth[8][9] = 0.00161;
//Depth 9 - 1280-1300 - @1290
dose_layer_depth[9][0] = 48.6; dose_layer_depth[9][1] = 204.0;
dose_layer_depth[9][2] = 0.260; dose_layer_depth[9][3] = 2.19;
dose_layer_depth[9][4] = 0.528; dose_layer_depth[9][5] = 0.386;
dose_layer_depth[9][6] = 0.231; dose_layer_depth[9][7] = 0.19;
dose_layer_depth[9][8] = 0.773; dose_layer_depth[9][9] = 0.00502;
}
| HFDarkening::~HFDarkening | ( | ) |
Definition at line 90 of file HFDarkening.cc.
{ }
| float HFDarkening::degradation | ( | float | mrad | ) |
Definition at line 143 of file HFDarkening.cc.
References funct::exp(), and funct::pow().
Referenced by HCalSD::ProcessHits().
| float HFDarkening::dose | ( | int | layer, |
| float | radius | ||
| ) |
Definition at line 92 of file HFDarkening.cc.
References dose_layer_depth, and radius.
Referenced by HCalSD::ProcessHits().
{
// float dose_acquired = 0.;
int radius = (int) floor((Radius-12.5)/10.0+0.5);
if (layer == 15) {
return 0.;
}
return dose_layer_depth[layer][radius];
/*
for (int i = 0; i<= layer; ++i) {
dose_acquired+=dose_layer_depth[i][radius];
}
if((layer<-1)||(layer>17)) return (0.);
if((layer==-1)||(layer==0)) {
if(ind<0) return (dose_layer_depth[0][0]); if(ind>=9) return (dose_layer_depth[0][10]);
if((ind>=0)&&(ind<=9)) {float delta = (dose_layer_depth[0][ind+1]-dose_layer_depth[0][ind])/10; return (dose_layer_depth[0][ind]+delta*(Radius-radius[ind]));}
}
if((layer>=1)&&(layer<=3)) {
if(ind<0) return (dose_layer_depth[5][0]); if(ind>=9) return (dose_layer_depth[5][10]);
if((ind>=0)&&(ind<=9)) {float delta = (dose_layer_depth[5][ind+1]-dose_layer_depth[5][ind])/10; return (dose_layer_depth[5][ind]+delta*(Radius-radius[ind]));}
}
if((layer==4)||(layer==5)) {
if(ind<0) return (dose_layer_depth[10][0]); if(ind>=9) return (dose_layer_depth[10][10]);
if((ind>=0)&&(ind<=9)) {float delta = (dose_layer_depth[10][ind+1]-dose_layer_depth[10][ind])/10; return (dose_layer_depth[10][ind]+delta*(Radius-radius[ind]));}
}
if((layer>=6)&&(layer<=8)) {
if(ind<0) return (dose_l6_l8[0]); if(ind>=23) return (dose_l6_l8[23]);
if((ind>=0)&&(ind<=23)) {float delta = (dose_l6_l8[ind+1]-dose_l6_l8[ind])/10; return (dose_l6_l8[ind]+delta*(Radius-radius[ind]));}
}
if((layer==9)||(layer==10)) {
if(ind<0) return (dose_l9_l10[0]); if(ind>=23) return (dose_l9_l10[23]);
if((ind>=0)&&(ind<=23)) {float delta = (dose_l9_l10[ind+1]-dose_l9_l10[ind])/10; return (dose_l9_l10[ind]+delta*(Radius-radius[ind]));}
}
if((layer>=11)&&(layer<=13)) {
if(ind<0) return (dose_l11_l13[0]); if(ind>=23) return (dose_l11_l13[23]);
if((ind>=0)&&(ind<=23)) {float delta = (dose_l11_l13[ind+1]-dose_l11_l13[ind])/10; return (dose_l11_l13[ind]+delta*(Radius-radius[ind]));}
}
if((layer==14)||(layer==15)) {
if(ind<0) return (dose_l14_l15[0]); if(ind>=23) return (dose_l14_l15[23]);
if((ind>=0)&&(ind<=23)) {float delta = (dose_l14_l15[ind+1]-dose_l14_l15[ind])/10; return (dose_l14_l15[ind]+delta*(Radius-radius[ind]));}
}
if((layer==16)||(layer==17)) {
if(ind<0) return (dose_l16_l17[0]); if(ind>=23) return (dose_l16_l17[23]);
if((ind>=0)&&(ind<=23)) {float delta = (dose_l16_l17[ind+1]-dose_l16_l17[ind])/10; return (dose_l16_l17[ind]+delta*(Radius-radius[ind]));}
}
return 0.;
*/
}
| float HFDarkening::int_lumi | ( | float | intlumi | ) |
Definition at line 147 of file HFDarkening.cc.
Referenced by HCalSD::ProcessHits().
{
return (intlumi/500.);
}
float HFDarkening::dose_layer_depth[10][10] [private] |
Definition at line 16 of file HFDarkening.h.
Referenced by dose(), and HFDarkening().
float HFDarkening::radius[10] [private] |
Definition at line 15 of file HFDarkening.h.
Referenced by dose(), and HFDarkening().
1.7.3