Defined silica_RIND_UI values explicitly

MonteCarlo/src/Materials.cc

@@ -121,10 +121,18 @@ void Materials::DefineOpticalProperties(void){
 
   // Refractive index and absorption length given by
   // https://www.heraeus.com/media/media/hca/doc_hca/products_and_solutions_8/optics/Data_and_Properties_Optics_fused_silica_EN.pdf
-  G4double silica_RIND_UI[50], silica_ABSL_UI[50];
-   for(int i = 0; i < nEntries_UI; i++){
-       silica_RIND_UI[i] = 1.44221 + 5.9e3*photonEnergy_UI[i] + 6.71e2*pow(photonEnergy_UI[i],2) + 2.3*pow(photonEnergy_UI[i],3) + 2.745*pow(photonEnergy_UI[i],4); //Refractive Index
-       silica_ABSL_UI[i] = 300.00*cm + i*20*cm; //Attenuation length
+
+  G4double silica_RIND_UI[50] = { 1.45740, 1.45770, 1.45800, 1.45830, 1.45860, 1.45889, 1.45919, 1.45949,
+                                  1.45979, 1.46009, 1.46039, 1.46069, 1.46103, 1.46137, 1.46171, 1.46203,
+                                  1.46235, 1.46267, 1.46299, 1.46333, 1.46370, 1.46406, 1.46442, 1.46479,
+                                  1.46515, 1.46552, 1.46588, 1.46623, 1.46663, 1.46703, 1.46743, 1.46782,
+                                  1.46822, 1.46862, 1.46901, 1.46941, 1.46982, 1.47025, 1.47069, 1.47113,
+                                  1.47157, 1.47201, 1.47245, 1.47291, 1.47338, 1.47384, 1.47430, 1.47477,
+                                  1.47523, 1.47573};
+
+  G4double silica_ABSL_UI[50];
+  for(int i = 0; i < nEntries_UI; i++){
+    silica_ABSL_UI[i] = 300.00*cm + i*20*cm; //Attenuation length
   }