Merge branch 'main' into off_specular

eos/normalization.py

@@ -21,6 +21,7 @@ class LZNormalisation:
     norm: np.ndarray
 
     def __init__(self, reference:EventDatasetProtocol, normalisationMethod: NormalisationMethod, grid: LZGrid):
+        mSM = 6
         self.angle = reference.geometry.nu-reference.geometry.mu
         lamda_e = reference.data.events.lamda
         detZ_e = reference.data.events.detZ
@@ -36,10 +37,10 @@ class LZNormalisation:
             theta_lz = grid.lz()*theta_z
             qz_lz = 4.0*np.pi*np.sin(np.deg2rad(theta_lz))/lamda_lz
             # TODO: introduce variable for `m` and propably for the slope
-            # Correct reflectivity of m=5 supermirror
+            # Correct reflectivity of m=mSM supermirror
             Rsm_lz = np.ones(np.shape(qz_lz))
             Rsm_lz = np.where(qz_lz>0.0217, 1-(qz_lz-0.0217)*(0.0625/0.0217), Rsm_lz)
-            Rsm_lz = np.where(qz_lz>0.0217*5, np.nan, Rsm_lz)
+            Rsm_lz = np.where(qz_lz>0.0217*mSM, np.nan, Rsm_lz)
             self.norm = norm_lz/Rsm_lz
         self.file_list = [os.path.basename(entry) for entry in reference.file_list]