MESA Released!

hp4155/MESA/lib/help.py

+import matplotlib.pyplot as plt
+import numpy as np
+import time
+from datetime import datetime
+
+import tkinter as tk
+from tkinter import filedialog
+import tkinter.messagebox
+
+import pandas as pd
+
+def normalization_factor(width):
+    factor = 10**6/width
+    return factor
+
+def add_widgets_to_list(source_dictionary,target_list):
+    for widget in source_dictionary.values():
+        target_list.append(widget)
+
+def change_state(widgets_list):
+    for widget in widgets_list:
+        widget.disabled = not widget.disabled
+
+#saving functions : save parameters and dataframe to txt
+def create_file(filename):
+    root = tk.Tk()
+    root.withdraw()
+    root.lift() #show window above all other applications
+
+    root.attributes("-topmost", True)#window stays above all other applications
+
+    file = filedialog.asksaveasfilename(defaultextension=".txt", filetypes=[("Text files","*.txt")],title = "save results path",initialfile =filename)
+
+    #check if the file path is correct(.txt)
+    while file.endswith(".txt") == False:
+        #open again filedialog with error message box
+        tk.messagebox.showerror(message='invalid filename!')
+        file = filedialog.asksaveasfilename(defaultextension=".txt", filetypes=[("Text files","*.txt")],title = "save results path",initialfile =filename)
+    root.destroy()
+    return file
+
+
+def check_values(parameters):
+    valid = True
+
+    root = tk.Tk()
+    root.withdraw()
+    root.lift() #show window above all other applications
+
+    root.attributes("-topmost", True)#window stays above all other applications
+
+    #get better variable names
+    start = parameters['start']
+    stop = parameters['stop']
+    step = parameters['step']
+
+    #check values
+    if abs(step.value) > abs(stop.value-start.value) or step.value==0:#invalid parameter setting 
+        valid = False
+        tkinter.messagebox.showerror(message="Invalid parameter setting!")
+
+    if start.value<stop.value and step.value<0: #change polarity
+        step.value =(-1)*step.value
+
+    elif start.value>stop.value and step.value>0:
+        step.value = (-1)*step.value
+
+    else:
+        pass
+
+    root.destroy()
+    return valid
+
+#Setup The device for mesa measurement 
+def setup(device):
+    device.reset()
+
+    #setup sweep measurement mode
+    device.measurement_mode('SWE')
+
+    #disable all irrelevant units
+    device.disable_not_smu()
+
+    #disable smu 2 and 4
+    device.smu_disable(2)
+    device.smu_disable(4)
+
+def measure(parameters,norm,device):
+    #delete user functions
+    device.del_user_functions()
+    
+    smu1 = device.smu_dict()
+    smu1.update(vname ='VGND',iname = 'IGND',mode = 'COMM',func='CONS')
+
+    smu3 = device.smu_dict()
+    smu3.update(vname='V',iname='I',mode ='V',func='VAR1')
+
+    device.setup_smu(1,smu1)
+    device.setup_smu(3,smu3)
+
+    #define user functions
+    device.user_function('Imm','mA/mm',f'I*{norm}')
+    device.user_function('ABSImm','mA/mm','ABS(Imm)')
+
+    var1 = device.var1_dict()
+    var1.update(
+        mode = parameters['hyst'].value,
+        start = parameters['start'].value,
+        stop = parameters['stop'].value,
+        step = parameters['step'].value,
+        comp = parameters['comp'].value,
+        pcomp = 0
+    )
+
+    device.setup_var1(var1)
+    device.integration_time(parameters['integration'].value)
+
+
+    #display
+    device.display_variable('X','V')
+    device.display_variable('Y1','ABSImm')
+
+    device.axis_scale('Y1','LOG')
+
+    if parameters['start'].value < parameters['stop'].value:
+        device.display_variable_min_max('X','MIN',parameters['start'].value)
+        device.display_variable_min_max('X','MAX',parameters['stop'].value)
+
+    else:
+        device.display_variable_min_max('X','MAX',parameters['start'].value)
+        device.display_variable_min_max('X','MIN',parameters['stop'].value)
+
+    if parameters['comp'].value !=0:
+        device.display_variable_min_max('Y1','MIN',0)
+        device.display_variable_min_max('Y1','MAX',abs(parameters['comp'].value)*norm)
+
+    
+    variables_list = ['V','I','Imm','ABSImm']
+    device.variables_to_save(variables_list)
+
+    device.single_measurement()
+    while device.operation_completed() == False:
+        pass
+    
+    device.autoscaling()
+    
+    values = dict([(variable,device.return_values(variable)) for variable in variables_list])
+    return values
+
+def plot_results(x,y):
+    fig, ax1 = plt.subplots()
+
+    color = 'tab:red'
+    ax1.set_xlabel('V/V') 
+    ax1.set_ylabel('I/mA/mm', color = color)
+    ax1.set_yscale('log')
+
+    ax1.plot(x,y, color = color) 
+    ax1.tick_params(axis ='y', labelcolor = color)
+
+    fig.suptitle('MESA Plot', fontweight ="bold")
+    fig.tight_layout()
+    fig.show() #interactve Figures
+
+def save_to_file(values,sample,parameters,file):
+    with open(file,'w') as f:
+        date = str(datetime.today().replace(microsecond=0))
+        f.write(f"Mesa Measurement at {date}"+"\n")
+        f.write(f"Series:{sample['processing_number'].value}"+"\n")
+        f.write(f"Sample:{sample['sample'].value}"+"\n")
+        f.write(f"Field:{sample['field'].value}"+"\n")
+        f.write(f"Width/um:{sample['width'].value}"+"\n\n")
+
+        f.write('Parameters\n')
+        f.write(f"V from {parameters['start'].value}V to {parameters['stop'].value}V with step {parameters['step'].value}V"+"\n")
+        f.write(f"I Compliance/A:{parameters['comp'].value}"+"\n")
+        f.write(f"Integration Time:{parameters['integration'].value}"+"\n")
+        
+        f.write("\nResults\n")
+
+    #create pandas dataframe
+    df = pd.DataFrame(values)
+
+    df = df.drop(columns=['ABSImm'])
+    df = df.rename(columns={'V':'V/V','I':'I/A','Imm':'Imm/mA/mm'})
+
+    df.to_csv(file,sep= " ",mode = 'a')
\ No newline at end of file

hp4155/MESA/lib/interface.py

+import ipywidgets as widgets
+from ipywidgets import GridspecLayout,Layout
+from IPython.display import clear_output
+import sys
+import os
+
+width = "50%"
+height = 'auto'
+style = {'description_width': 'initial'}
+floatbox_width = "80%"
+
+def parameters_interface():
+    grid = GridspecLayout(5,4)
+    grid[0,:]=widgets.Label("V (SMU3)",layout=Layout(height='auto', width='auto'))
+    grid[0,:].style.font_weight='bold'
+
+
+    #first line
+    grid[1,0]=widgets.Label("Start(V)",layout=Layout(height='auto', width='auto'))
+    grid[1,1]=widgets.Label("Step(V)",layout=Layout(height='auto', width='auto'))
+    grid[1,2]=widgets.Label("Stop(V)",layout=Layout(height='auto', width='auto'))
+
+    #second line
+    grid[2,0]=widgets.BoundedFloatText(value=0,min=-100,max=100,step=1,layout=Layout(height='auto', width=floatbox_width))
+    grid[2,1]=widgets.BoundedFloatText(value=0.5,min=-200,max=200,step=1,layout=Layout(height='auto', width=floatbox_width))
+    grid[2,2]=widgets.BoundedFloatText(value=100,min=-100,max=100,step=1,layout=Layout(height='auto', width=floatbox_width))
+
+    #third line 
+    grid[3,0]=widgets.Label("Compliance(A)",layout=Layout(height='auto', width='auto'))
+    grid[3,1] =widgets.Label("Integration Time",layout=Layout(height='auto', width='auto'))#mind the gap
+    grid[3,2] =widgets.Label("Hysterisis",layout=Layout(height='auto', width='auto'))#mind the gap
+
+    #fourth line
+    grid[4,0]=widgets.BoundedFloatText(value=0.02,min=-0.1,max=0.1,step=0.01,layout=Layout(height='auto', width=floatbox_width))
+    grid[4,1]=widgets.Dropdown(options=["SHORt","MEDium","LONG"],value="MEDium",style =style,layout=Layout(height='auto', width=floatbox_width))
+    grid[4,2]=widgets.Dropdown(options=['SINGle','DOUBle'],value='SINGle',layout=Layout(height='auto', width=floatbox_width))#mind the gap
+
+    parameters = {
+        'start': grid[2,0],
+        'step': grid[2,1],
+        'stop': grid[2,2],
+        'comp': grid[4,0],
+        'hyst':grid[4,2],
+        'integration':grid[4,1]
+    }
+    
+
+    display(grid)
+    return parameters
+
+
+def sample_information_interface():
+    width = '90%'
+    sample_information=widgets.Label("Sample Information",layout=Layout(height=height, width='50%'))
+    sample_information.style.font_weight='bold'
+    information_grid=GridspecLayout(3,2)
+    
+    for i in range(3):
+        for j in range(2):
+            if i ==2 and j == 1:
+               #information_grid[i,j]=widgets.Dropdown(options=['mA/mm','uA/um'],value='mA/mm',layout=Layout(height=height, width=width))#mind the gap
+                pass
+            elif i == 2 and j == 0:                
+                information_grid[i,j]=widgets.BoundedFloatText(
+                    value=100,
+                    min=1e-3,
+                    max=sys.float_info.max,step=1,
+                    layout=Layout(height=height, width=width)
+                )
+            else:
+                information_grid[i,j]=widgets.Text(layout=Layout(height=height, width=width))
+
+    information_grid[0,0].description = "Processing-Nr:"
+    information_grid[1,0].description = "Sample:"
+    information_grid[2,0].description = "Device Width(um):"
+    information_grid[0,1].description = "Field(XYY):"
+    information_grid[1,1].description = "Device:"
+    information_grid[1,1].disabled = True
+    #information_grid[2,1].description = "Normalization:"
+
+    for i in range(3):
+        for j in range(2):
+            try:
+                information_grid[i,j].style = style
+            except:
+                pass
+
+    image_title = widgets.Label("SMU Configuration",layout=Layout(height=height, width='50%'))
+    image_title.style.font_weight='bold'
+
+    filename = os.getcwd()+r"\lib\MESA_SMU_configuration.png"
+    #print(filename)
+
+    file = open(filename, "rb")
+    image = file.read()
+    image_widget =widgets.Image(
+        value = image,
+        format='png',
+        width='auto',
+        height='auto',
+    )
+    
+    #display(widgets.HBox([sample_information,image_title]))
+    #display(widgets.HBox([information_grid,image_widget]))
+    vbox1 =widgets.VBox([sample_information,information_grid])
+    vbox2 =widgets.VBox([image_title,image_widget])
+
+    display(widgets.HBox([vbox1,vbox2]))
+
+    information = {
+        'processing_number': information_grid[0,0],
+        'sample' : information_grid[1,0],
+        'field': information_grid[0,1],
+        #'device':information_grid[1,1],    
+        'width': information_grid[2,0],
+    }
+
+    return information
+            
\ No newline at end of file

hp4155/MESA/mesa.py

+import sys
+sys.path.insert(0, './lib')
+sys.path.insert(0, '..') #append parent directory
+
+import hp4155a
+from interface import *
+from help import *
+
+sample = sample_information_interface()
+parameters =parameters_interface()
+button = widgets.Button(description ='Start Measurement')
+output = widgets.Output()
+
+display(button,output)
+
+all_widgets = [button]
+add_widgets_to_list(sample,all_widgets)
+add_widgets_to_list(parameters,all_widgets)
+
+device = hp4155a.HP4155a('GPIB0::17::INSTR')
+setup(device)
+def on_button_clicked(b):
+    with output:
+        clear_output(wait = True)
+        change_state(all_widgets)
+        norm=normalization_factor(sample['width'].value)
+        
+        valid = check_values(parameters)
+        if valid == True:
+            values = measure(parameters,norm,device)
+            plot_results(values['V'],values['ABSImm'])
+            filename = f"MESA_ISO_{sample['field'].value}.txt"
+            file = create_file(filename)
+            save_to_file(values,sample,parameters,file)
+        
+        change_state(all_widgets)
+
+button.on_click(on_button_clicked)
+    
\ No newline at end of file

hp4155/MESA/mesa_interface.ipynb

+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": 1,
+   "id": "1b607685-438c-466e-9671-140ae6c8d9f9",
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "application/vnd.jupyter.widget-view+json": {
+       "model_id": "1e8f68e7a1a44930b5f8aa7274a71bd1",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "HBox(children=(VBox(children=(Label(value='Sample Information', layout=Layout(height='auto', width='50%'), sty…"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "application/vnd.jupyter.widget-view+json": {
+       "model_id": "59dc533f8cc2493597699a3a72882bbf",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "GridspecLayout(children=(Label(value='V (SMU3)', layout=Layout(grid_area='widget001', height='auto', width='au…"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "application/vnd.jupyter.widget-view+json": {
+       "model_id": "b79166daed814a348b6538ae9102b66c",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "Button(description='Start Measurement', style=ButtonStyle())"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "application/vnd.jupyter.widget-view+json": {
+       "model_id": "68d3460573264361a55153b8c1f43020",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "Output()"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "%matplotlib widget\n",
+    "%run mesa.py"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "c62885a4-1186-4770-937f-fc8b5c484b98",
+   "metadata": {},
+   "outputs": [],
+   "source": []
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3 (ipykernel)",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.11.4"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}