Skip to content
Snippets Groups Projects
Commit f8a8906f authored by Simon Neuhaus's avatar Simon Neuhaus Committed by Volker Friese
Browse files

Fix for Error with Geant4 optical properties when using Geant4 11.2

parent 9adde218
No related branches found
No related tags found
1 merge request!1764Fix for Error with Geant4 optical properties when using Geant4 11.2
Pipeline #28936 passed
......@@ -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;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment