From 6dea38ddb01a15ff048c6fd4330d78dc0fe5596b Mon Sep 17 00:00:00 2001 From: Artur Glavic Date: Tue, 10 Nov 2026 11:20:08 +0100 Subject: [PATCH] Change supermirror to m=6 with variable mSM --- eos/normalization.py | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/eos/normalization.py b/eos/normalization.py index aeb030f..abf362a 100644 --- a/eos/normalization.py +++ b/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]