diff --git a/sim/detectors/rich/CbmRich.cxx b/sim/detectors/rich/CbmRich.cxx index 99d2d269a070ecd2cbfcd62e738b96b07dc301ac..b1c785cdc4941546cb433b7d70f547b5b9add841 100644 --- a/sim/detectors/rich/CbmRich.cxx +++ b/sim/detectors/rich/CbmRich.cxx @@ -314,22 +314,22 @@ void CbmRich::SetRichGlassPropertiesForGeant4() } std::vector<Double_t> energyAr = { - 6.1997, 5.9045, 5.6361, 5.391, 5.1664, 4.9598, 4.769, 4.5924, 4.4284, 4.2757, 4.1331, 3.9998, 3.8748, - 3.7574, 3.6469, 3.5427, 3.4443, 3.3512, 3.263, 3.1793, 3.0998, 3.0242, 2.9522, 2.8836, 2.818, 2.7554, - 2.6955, 2.6382, 2.5832, 2.5305, 2.4799, 2.4313, 2.3845, 2.3395, 2.2962, 2.2544, 2.2142, 2.1753, 2.1378, - 2.1016, 2.0666, 2.0327, 1.9999, 1.9682, 1.9374, 1.9076, 1.8787, 1.8507, 1.8234, 1.797, 1.7713, 1.7464, - 1.7221, 1.6985, 1.6756, 1.6533, 1.6315, 1.6103, 1.5897, 1.5695, 1.5499}; + 1.5499, 1.5695, 1.5897, 1.6103, 1.6315, 1.6533, 1.6756, 1.6985, 1.7221, 1.7464, 1.7713, 1.7970, 1.8234, + 1.8507, 1.8787, 1.9076, 1.9374, 1.9682, 1.9999, 2.0327, 2.0666, 2.1016, 2.1378, 2.1753, 2.2142, 2.2544, + 2.2962, 2.3395, 2.3845, 2.4313, 2.4799, 2.5305, 2.5832, 2.6382, 2.6955, 2.7554, 2.8180, 2.8836, 2.9522, + 3.0242, 3.0998, 3.1793, 3.2630, 3.3512, 3.4443, 3.5427, 3.6469, 3.7574, 3.8748, 3.9998, 4.1331, 4.2757, + 4.4284, 4.5924, 4.7690, 4.9598, 5.1664, 5.3910, 5.6361, 5.9045, 6.1997}; // Transform to GeV for (Double_t& energy : energyAr) { energy *= 1.e-9; } std::vector<Double_t> reflectivityAr = { - 0.22529, 0.1862, 0.15901, 0.14713, 0.13963, 0.13898, 0.13762, 0.13622, 0.13868, 0.13951, 0.14031, 0.14073, 0.1409, - 0.13977, 0.14205, 0.14072, 0.1396, 0.13826, 0.13665, 0.13513, 0.13463, 0.13287, 0.13182, 0.13084, 0.12824, 0.12601, - 0.12622, 0.12681, 0.12193, 0.12011, 0.12109, 0.11908, 0.11526, 0.11364, 0.11385, 0.12015, 0.11935, 0.11712, 0.1208, - 0.12021, 0.11909, 0.11783, 0.12257, 0.1215, 0.12199, 0.12494, 0.13101, 0.12089, 0.12284, 0.12569, 0.13136, 0.13307, - 0.13705, 0.13844, 0.13753, 0.14416, 0.14761, 0.14953, 0.15218, 0.15315, 0.15719}; + 0.15719, 0.15315, 0.15218, 0.14953, 0.14761, 0.14416, 0.13753, 0.13844, 0.13705, 0.13307, 0.13136, 0.12569, 0.12284, + 0.12089, 0.13101, 0.12494, 0.12199, 0.12150, 0.12257, 0.11783, 0.11909, 0.12021, 0.12080, 0.11712, 0.11935, 0.12015, + 0.11385, 0.11364, 0.11526, 0.11908, 0.12109, 0.12011, 0.12193, 0.12681, 0.12622, 0.12601, 0.12824, 0.13084, 0.13182, + 0.13287, 0.13463, 0.13513, 0.13665, 0.13826, 0.13960, 0.14072, 0.14205, 0.13977, 0.14090, 0.14073, 0.14031, 0.13951, + 0.13868, 0.13622, 0.13762, 0.13898, 0.13963, 0.14713, 0.15901, 0.18620, 0.22529}; for (Double_t& reflectivity : reflectivityAr) { reflectivity = 1. - reflectivity;