majoplot 0.1.10__tar.gz → 0.1.12__tar.gz

{majoplot-0.1.10 → majoplot-0.1.12}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: majoplot
-Version: 0.1.10
+Version: 0.1.12
 Summary: Automates scenario-specific data preprocessing and plotting workflows for condensed-matter physics labs (OriginLab COM backend + Matplotlib for preview).
 Project-URL: Homepage, https://github.com/ponyofshadows/majoplot
 Project-URL: Source, https://github.com/ponyofshadows/majoplot

{majoplot-0.1.10 → majoplot-0.1.12}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "majoplot"
-version = "0.1.10"
+version = "0.1.12"
 description = "Automates scenario-specific data preprocessing and plotting workflows for condensed-matter physics labs (OriginLab COM backend + Matplotlib for preview)."
 readme = "README.md"
 requires-python = ">=3.11"

{majoplot-0.1.10 → majoplot-0.1.12}/src/majoplot/domain/importers/PPMS_Resistivity.py

@@ -6,10 +6,16 @@ from ..utils import skip_lines_then_readline
 
 ROW_FILEOPENTIME  = 3
 ROW_SAMPLE1_NAME  = 6
+ROW_SAMPLE1_CROSS_SECTION = 7
+ROW_SAMPLE1_LENGTH = 8
 ROW_SAMPLE1_UNITS = 9
 ROW_SAMPLE2_NAME  = 10
+ROW_SAMPLE2_CROSS_SECTION = 11
+ROW_SAMPLE2_LENGTH = 12
 ROW_SAMPLE2_UNITS = 13
 ROW_SAMPLE3_NAME  = 14
+ROW_SAMPLE3_CROSS_SECTION = 15
+ROW_SAMPLE3_LENGTH = 16
 ROW_SAMPLE3_UNITS = 17
 ROW_HEADERS       = 32
 
@@ -52,8 +58,36 @@ class PPMS_Resistivity:
             return fail_signal
         labels["sample1_name"] = LabelValue(s1_name)
 
+        # S1 Cross Section
+        line = skip_lines_then_readline(raw_data_file, ROW_SAMPLE1_CROSS_SECTION - ROW_SAMPLE1_NAME -1)
+        info = line.split(",")
+        try:
+            if info[-1].strip() != "Sample1 Cross Section":
+                return fail_signal
+            s1_cross_section = info[1]
+        except IndexError:
+            return fail_signal
+        try:
+            labels["sample1_cross_section"] = LabelValue(float(s1_cross_section),unit="mm²")
+        except ValueError:
+            return fail_signal
+
+        # S1 Length
+        line = skip_lines_then_readline(raw_data_file, ROW_SAMPLE1_LENGTH - ROW_SAMPLE1_CROSS_SECTION -1)
+        info = line.split(",")
+        try:
+            if info[-1].strip() != "Sample1 Length":
+                return fail_signal
+            s1_length = info[1]
+        except IndexError:
+            return fail_signal
+        try:
+            labels["sample1_length"] = LabelValue(float(s1_length),unit="mm")
+        except ValueError:
+            return fail_signal
+
         # S1 Units
-        line = skip_lines_then_readline(raw_data_file, ROW_SAMPLE1_UNITS - ROW_SAMPLE1_NAME -1)
+        line = skip_lines_then_readline(raw_data_file, ROW_SAMPLE1_UNITS - ROW_SAMPLE1_LENGTH -1)
         info = line.split(",")
         try:
             if info[-1].strip() != "Sample1 Units":
@@ -74,8 +108,37 @@ class PPMS_Resistivity:
             return fail_signal
         labels["sample2_name"] = LabelValue(s2_name)
 
+        # S2 Cross Section
+        line = skip_lines_then_readline(raw_data_file, ROW_SAMPLE2_CROSS_SECTION - ROW_SAMPLE2_NAME -1)
+        info = line.split(",")
+        try:
+            if info[-1].strip() != "Sample2 Cross Section":
+                return fail_signal
+            s2_cross_section = info[1]
+        except IndexError:
+            return fail_signal
+        try:
+            labels["sample2_cross_section"] = LabelValue(float(s2_cross_section),unit="mm²")
+        except ValueError:
+            return fail_signal
+
+        # S1 Length
+        line = skip_lines_then_readline(raw_data_file, ROW_SAMPLE2_LENGTH - ROW_SAMPLE2_CROSS_SECTION -1)
+        info = line.split(",")
+        try:
+            if info[-1].strip() != "Sample2 Length":
+                return fail_signal
+            s2_length = info[1]
+        except IndexError:
+            return fail_signal
+        try:
+            labels["sample2_length"] = LabelValue(float(s2_length),unit="mm")
+        except ValueError:
+            return fail_signal
+
+
         # S2 Units
-        line = skip_lines_then_readline(raw_data_file, ROW_SAMPLE2_UNITS - ROW_SAMPLE2_NAME -1)
+        line = skip_lines_then_readline(raw_data_file, ROW_SAMPLE2_UNITS - ROW_SAMPLE2_LENGTH -1)
         info = line.split(",")
         try:
             if info[-1].strip() != "Sample2 Units":
@@ -96,8 +159,36 @@ class PPMS_Resistivity:
             return fail_signal
         labels["sample3_name"] = LabelValue(s3_name)
 
+        # S3 Cross Section
+        line = skip_lines_then_readline(raw_data_file, ROW_SAMPLE3_CROSS_SECTION - ROW_SAMPLE3_NAME -1)
+        info = line.split(",")
+        try:
+            if info[-1].strip() != "Sample3 Cross Section":
+                return fail_signal
+            s3_cross_section = info[1]
+        except IndexError:
+            return fail_signal
+        try:
+            labels["sample3_cross_section"] = LabelValue(float(s3_cross_section),unit="mm²")
+        except ValueError:
+            return fail_signal
+
+        # S1 Length
+        line = skip_lines_then_readline(raw_data_file, ROW_SAMPLE3_LENGTH - ROW_SAMPLE3_CROSS_SECTION -1)
+        info = line.split(",")
+        try:
+            if info[-1].strip() != "Sample3 Length":
+                return fail_signal
+            s3_length = info[1]
+        except IndexError:
+            return fail_signal
+        try:
+            labels["sample3_length"] = LabelValue(float(s3_length),unit="mm")
+        except ValueError:
+            return fail_signal
+
         # S3 Units
-        line = skip_lines_then_readline(raw_data_file, ROW_SAMPLE3_UNITS - ROW_SAMPLE3_NAME -1)
+        line = skip_lines_then_readline(raw_data_file, ROW_SAMPLE3_UNITS - ROW_SAMPLE3_LENGTH -1)
         info = line.split(",")
         try:
             if info[-1].strip() != "Sample3 Units":

{majoplot-0.1.10 → majoplot-0.1.12}/src/majoplot/domain/scenarios/PPMS_Resistivity/RT.py

@@ -57,8 +57,8 @@ class RT:
                     # clear null R points
                     s_points = np.array([point for point in s_points if point[1]])
                     # record
-                    Imin = np.min(s_points[:,1])
-                    Imax = np.max(s_points[:,1])
+                    Imin = np.min(s_points[:,2])
+                    Imax = np.max(s_points[:,2])
                     if (Imax - Imin) / Imax < 0.03:
                         Irange = f"{np.mean(s_points[:,1]):.1e}"
                     else:
@@ -92,6 +92,8 @@ class RT:
             major_grid=None,
             major_tick=TickSpec(),
             legend=LegendSpec(fontsize=5),
+            linewidth=1,
+            marker_size=2,
         )
             
 
@@ -104,13 +106,15 @@ class RT:
             title=None,
             figsize=FIGSIZE,
             linestyle_cycle= ("-",),
+
             linecolor_cycle = (
-                "#2d0b59", "#3b0f6f", "#4a136e", "#5a176e", "#6a1c6e",
-                "#7a216f", "#8b2770", "#9b2d71", "#ac3372", "#bd3973",
-                "#ce4074", "#df4775", "#f04f76", "#f86a5a", "#fb8c3c",
-                "#fdbb2d", "#fcfdbf",
+                "#515151", "#F14040", "#1A6FDF", "#37AD6B", "#B177DE",
+                "#CC9900", "#00CBCC", "#7D4E4E", "#8E8E00", "#FB6501",
+                "#6699CC", "#6FB802", "#f04f76", "#f86a5a", "#fb8c3c",
+                "#fdbb2d", "#fcfdbf","#2B2E83", "#E6007A", "#005F5F", 
+                "#1F3D2B",  "#E8C6E8","#FAB3d1",
                 ),
-            linemarker_cycle = ("o","s","^","v","d","*","x","+"),
+            linemarker_cycle = ("o",),
             alpa_cycle = (1.0,),
         )
     

majoplot-0.1.12/src/majoplot/domain/scenarios/PPMS_Resistivity/RT_Resistivity.py

@@ -0,0 +1,136 @@
+from __future__ import annotations
+
+from ...base import *
+
+
+FIGSIZE = (8, 6)
+
+T = "Temperature (K)"
+H = "Magnetic Field (Oe)"
+R = "Resistivity (Ω·m)"
+RI ={
+    1: {"R":"Bridge 1 Resistance (Ohms)", "I":"Bridge 1 Excitation (uA)"},
+    2: {"R":"Bridge 2 Resistance (Ohms)", "I":"Bridge 2 Excitation (uA)"},
+    3: {"R":"Bridge 3 Resistance (Ohms)", "I":"Bridge 3 Excitation (uA)"},
+}
+
+ResistivityI ={
+    1: {"R":"Bridge 1 Resistivity (Ω·m)", "I":"Bridge 1 Excitation (uA)"},
+    2: {"R":"Bridge 2 Resistivity (Ω·m)", "I":"Bridge 2 Excitation (uA)"},
+    3: {"R":"Bridge 3 Resistivity (Ω·m)", "I":"Bridge 3 Excitation (uA)"},
+}
+
+class RT_Resistivity:
+    data_summary_label_names = ["H"]
+    axes_label_names = ("date", "raw_data", "bridge", "sample_name")
+    figure_label_names = ("date", "raw_data", "bridge", "sample_name")
+    figure_summary_label_names = ("raw_data","bridge","sample_name")
+    max_axes_in_one_figure = 1
+    project_to_child_folder_label_names = {"date":"sample_name"}
+    parent_folder_name = "RT"
+
+    @classmethod
+    def preprocess(cls, raw_datas:list[Data])->list[Data]:
+        datas = []
+        for raw_data in raw_datas:
+            raw_labels = raw_data.labels
+            headers = raw_data.headers
+            raw_points = raw_data.points
+            
+            # Split by H
+            split_datas = []
+
+            last_H_stage = round(raw_points[0][headers[H]])
+            current_points = [ last_H_stage, [ raw_points[0] ] ]
+
+            for point in raw_points[1:]:
+                cur_H_stage = round(point[headers[H]])
+                if cur_H_stage != last_H_stage:
+                    split_datas.append(current_points)
+                    last_H_stage = cur_H_stage
+                    current_points = [ last_H_stage, [ point ] ]
+                else:
+                    current_points[1].append(point)
+            else:
+                split_datas.append(current_points)
+
+            # 3 bridges
+            for H_stage, points in split_datas:
+                for i in range(1,4):
+                    _headerTRI = (T,RI[i]["R"],RI[i]["I"])
+                    _headers = (T,ResistivityI[i]["R"])
+                    s_points = [ [point[headers[x]] for x in _headerTRI] for point in points]
+                    # clear null Resistance points
+                    s_points = np.array([point for point in s_points if point[1]])
+                    # calculate resistivity
+                    cross_section = raw_labels[f"sample{i}_cross_section"].value * 1e-6
+                    if cross_section <= 0:
+                        cross_section = np.nan
+                    length = raw_labels[f"sample{i}_length"].value * 1e-3
+                    if length <= 0:
+                        length = np.nan
+                    r_points = np.column_stack([s_points[:,0], s_points[:,1] * cross_section / length])
+                    # record
+                    Imin = np.min(s_points[:,2])
+                    Imax = np.max(s_points[:,2])
+                    if (Imax - Imin) / Imax < 0.03:
+                        Irange = f"{np.mean(s_points[:,1]):.1e}"
+                    else:
+                        Irange = f"{Imin:.1e}~{Imax:.1e}"
+                    labels = LabelDict()
+                    labels["instrument"] = raw_labels["instrument"]
+                    labels["raw_data"] = raw_labels["raw_data"]
+                    labels["date"] = raw_labels["date"]
+                    labels["bridge"] = LabelValue(i,unit="Bridge",unit_as_postfix=False)
+                    labels["sample_name"] = raw_labels[f"sample{i}_name"]
+                    labels["sample_units"] = raw_labels[f"sample{i}_units"]
+                    labels["H"] = LabelValue(H_stage, unit="Oe")
+                    labels["I_range"] = LabelValue(Irange,unit="μA")
+                    labels.summary_names = cls.data_summary_label_names
+                    datas.append(Data(
+                        labels=labels,
+                        _headers=_headers,
+                        points=r_points,
+                        ignore_outliers=IgnoreOutlierSpec(min_gap_base=1e-8,min_gap_multiple=10),
+                    ))
+            
+        return datas
+            
+            
+
+    @classmethod
+    def make_axes_spec(cls, axes_labels, data_pool)->AxesSpec:
+        return AxesSpec(
+            x_axis_title=T,
+            y_axis_title=R,
+            major_grid=None,
+            major_tick=TickSpec(),
+            legend=LegendSpec(fontsize=5),
+            linewidth=1,
+            marker_size=2,
+        )
+            
+
+    @classmethod
+    def make_figure_spec(cls,figure_labels, axes_pool:Iterable[Axes])->FigureSpec:
+        figure_name = figure_labels.brief_summary
+
+        return FigureSpec(
+            name=figure_name,
+            title=None,
+            figsize=FIGSIZE,
+            linestyle_cycle= ("-",),
+            linecolor_cycle = (
+                "#515151", "#F14040", "#1A6FDF", "#37AD6B", "#B177DE",
+                "#CC9900", "#00CBCC", "#7D4E4E", "#8E8E00", "#FB6501",
+                "#6699CC", "#6FB802", "#f04f76", "#f86a5a", "#fb8c3c",
+                "#fdbb2d", "#fcfdbf","#2B2E83", "#E6007A", "#005F5F", 
+                "#1F3D2B",  "#E8C6E8","#FAB3d1",
+                ),
+            linemarker_cycle = ("o","s","^","v","d","*","x","+"),
+            alpa_cycle = (1.0,),
+        )
+    
+    @classmethod
+    def make_muti_axes_spec(cls, axes_pool:list[Axes])->MutiAxesSpec|FAIL|None:
+        return None

majoplot-0.1.12/src/majoplot/domain/scenarios/PPMS_Resistivity/RT_normalized.py

@@ -0,0 +1,142 @@
+from __future__ import annotations
+
+from ...base import *
+
+
+FIGSIZE = (8, 6)
+
+T = "Temperature (K)"
+H = "Magnetic Field (Oe)"
+R = "Resistance (Ohms)"
+RI ={
+    1: {"R":"Bridge 1 Resistance (Ohms)", "I":"Bridge 1 Excitation (uA)"},
+    2: {"R":"Bridge 2 Resistance (Ohms)", "I":"Bridge 2 Excitation (uA)"},
+    3: {"R":"Bridge 3 Resistance (Ohms)", "I":"Bridge 3 Excitation (uA)"},
+}
+RI_normalized ={
+    1: {"R":"Bridge 1 Normalized Resistance", "I":"Bridge 1 Excitation (uA)"},
+    2: {"R":"Bridge 2 Normalied Resistance", "I":"Bridge 2 Excitation (uA)"},
+    3: {"R":"Bridge 3 Normalized Resistance", "I":"Bridge 3 Excitation (uA)"},
+}
+
+class RT_normalized:
+    data_summary_label_names = ["H"]
+    axes_label_names = ("date", "raw_data", "bridge", "sample_name")
+    figure_label_names = ("date", "raw_data", "bridge", "sample_name")
+    figure_summary_label_names = ("raw_data","bridge","sample_name")
+    max_axes_in_one_figure = 1
+    project_to_child_folder_label_names = {"date":"sample_name"}
+    parent_folder_name = "RT"
+
+    @classmethod
+    def preprocess(cls, raw_datas:list[Data])->list[Data]:
+        datas = []
+        for raw_data in raw_datas:
+            raw_labels = raw_data.labels
+            headers = raw_data.headers
+            raw_points = raw_data.points
+            
+            # Split by H
+            split_datas = []
+
+            last_H_stage = round(raw_points[0][headers[H]])
+            current_points = [ last_H_stage, [ raw_points[0] ] ]
+
+            for point in raw_points[1:]:
+                cur_H_stage = round(point[headers[H]])
+                if cur_H_stage != last_H_stage:
+                    split_datas.append(current_points)
+                    last_H_stage = cur_H_stage
+                    current_points = [ last_H_stage, [ point ] ]
+                else:
+                    current_points[1].append(point)
+            else:
+                split_datas.append(current_points)
+
+            # 3 bridges
+            for H_stage, points in split_datas:
+                for i in range(1,4):
+                    _headerTRI = (T,RI[i]["R"],RI[i]["I"])
+                    _headers = (T,RI_normalized[i]["R"])
+                    s_points = [ [point[headers[x]] for x in _headerTRI] for point in points]
+                    # clear null R points
+                    s_points = np.array([point for point in s_points if point[1]])
+                    # record
+                    Imin = np.min(s_points[:,2])
+                    Imax = np.max(s_points[:,2])
+                    if (Imax - Imin) / Imax < 0.03:
+                        Irange = f"{np.mean(s_points[:,1]):.1e}"
+                    else:
+                        Irange = f"{Imin:.1e}~{Imax:.1e}"
+                    
+                    labels = LabelDict()
+                    labels["instrument"] = raw_labels["instrument"]
+                    labels["raw_data"] = raw_labels["raw_data"]
+                    labels["date"] = raw_labels["date"]
+                    labels["bridge"] = LabelValue(i,unit="Bridge",unit_as_postfix=False)
+                    labels["sample_name"] = raw_labels[f"sample{i}_name"]
+                    labels["sample_units"] = raw_labels[f"sample{i}_units"]
+                    labels["H"] = LabelValue(H_stage, unit="Oe")
+                    labels["I_range"] = LabelValue(Irange,unit="μA")
+                    labels.summary_names = cls.data_summary_label_names
+                    datas.append(Data(
+                        labels=labels,
+                        _headers=_headers,
+                        points=s_points[:,0:2],
+                        ignore_outliers=IgnoreOutlierSpec(min_gap_base=1e-4,min_gap_multiple=5),
+                    ))
+        
+        # normalize
+        for data in datas:
+            data: Data
+            ymax = np.max(np.abs(data.ylim))
+            data.points = np.column_stack( [data.points[:,0], np.asarray(data.points[:,1] / ymax)])
+            data.ignore_outliers = None
+            data._ignore_outliers_spec_cache = None
+            data._points_for_plot = None
+            data._x_for_plot = None
+            data._y_for_plot = None
+            data._xlim = None
+            data._ylim = None
+
+        return datas
+            
+            
+
+    @classmethod
+    def make_axes_spec(cls, axes_labels, data_pool)->AxesSpec:
+        return AxesSpec(
+            x_axis_title=T,
+            y_axis_title=R,
+            major_grid=None,
+            major_tick=TickSpec(),
+            legend=LegendSpec(fontsize=5),
+            linewidth=1,
+            marker_size=2,
+        )
+            
+
+    @classmethod
+    def make_figure_spec(cls,figure_labels, axes_pool:Iterable[Axes])->FigureSpec:
+        figure_name = figure_labels.brief_summary
+
+        return FigureSpec(
+            name=figure_name,
+            title=None,
+            figsize=FIGSIZE,
+            linestyle_cycle= ("-",),
+
+            linecolor_cycle = (
+                "#515151", "#F14040", "#1A6FDF", "#37AD6B", "#B177DE",
+                "#CC9900", "#00CBCC", "#7D4E4E", "#8E8E00", "#FB6501",
+                "#6699CC", "#6FB802", "#f04f76", "#f86a5a", "#fb8c3c",
+                "#fdbb2d", "#fcfdbf","#2B2E83", "#E6007A", "#005F5F", 
+                "#1F3D2B",  "#E8C6E8","#FAB3d1",
+                ),
+            linemarker_cycle = ("o",),
+            alpa_cycle = (1.0,),
+        )
+    
+    @classmethod
+    def make_muti_axes_spec(cls, axes_pool:list[Axes])->MutiAxesSpec|FAIL|None:
+        return None

majoplot-0.1.10/src/majoplot/domain/scenarios/VSM/MT_origin.py → majoplot-0.1.12/src/majoplot/domain/scenarios/VSM/MT_original.py

@@ -11,7 +11,7 @@ M = "DC Moment Free Ctr (emu)"
 H = "Magnetic Field (Oe)"
 
 
-class MT_origin:
+class MT_original:
     data_summary_label_names = ["mass","H","cooling_type"]
     axes_label_names = ("material","date","raw_data", "H")
     figure_label_names = ("material","date", "raw_data","H")

{majoplot-0.1.10 → majoplot-0.1.12}/uv.lock

@@ -252,7 +252,7 @@ wheels = [
 
 [[package]]
 name = "majoplot"
-version = "0.1.9"
+version = "0.1.11"
 source = { editable = "." }
 dependencies = [
     { name = "matplotlib" },