SetupControl.m

classdef SetupControl < handle
    properties
        equipment
    end
    
    methods
        %Init Functions
        function obj = SetupControl
            obj.equipment.LI=[]; %Lock-Ins
            obj.equipment.SM=[]; %Source Meters
            obj.equipment.EM=[]; %Electrometers
            obj.equipment.ITC503=[];%TemperatureController
            obj.equipment.spec=[];%spectrometer
            obj.equipment.WPControl=[]; %WaveplateStageControl
            obj.equipment.WPFig=[]; %WaveplateStageControlFigure
        end
        function InitComms(obj)
            instrreset;
            %pyversion C:\Python27\python.exe
            try obj.equipment.LI = gpib('ni',0,8); fopen(obj.equipment.LI); disp("LI 1 FOUND") %Lock-in SR830
            catch; obj.equipment.LI = []; disp("LI 1 ERROR")
            end
            if length(obj.equipment.LI)==1
                try obj.equipment.LI(2) = gpib('ni',0,9); fopen(obj.equipment.LI(2)); disp("LI 2 FOUND") %Lock-in SR830
                catch; obj.equipment.LI (2) = []; disp("LI 2 ERROR")
                end
            end
            try obj.equipment.SM = gpib('ni',0,28); fopen(obj.equipment.SM(1)); disp("SM 1 FOUND") %Source meter
            catch; obj.equipment.SM = []; disp("SM 1 ERROR")
            end
            if length(obj.equipment.SM)==1
                try obj.equipment.SM(2) = gpib('ni',0,26); fopen(obj.equipment.SM(2)); disp("SM 2 FOUND") %Source meter
                catch; obj.equipment.SM (2) = []; disp("SM 2 ERROR")
                end
            end
            try obj.equipment.EM = gpib('ni',0,18); fopen(obj.equipment.EM(1)); disp("EM 1 FOUND") %Electrometer
            catch; obj.equipment.EM = []; disp("EM 1 ERROR")
            end
            if length(obj.equipment.EM)==1
                try obj.equipment.EM(2) = gpib('ni',0,14); fopen(obj.equipment.EM(2)); disp("EM 2 FOUND") %Electrometer
                catch; obj.equipment.EM (2) = []; disp("EM 2 ERROR")
                end
            end
            try obj.equipment.ITC503 = gpib('ni',0,13); fopen(obj.equipment.ITC503);obj.equipment.ITC503.EOSMode='read&write';obj.equipment.ITC503.EOSCharCode='CR'; disp("ITC503 1 FOUND") %TemperatureController
            catch; obj.equipment.ITC503 = []; disp("ITC503 ERROR")
            end
            try obj.equipment.spec = py.stellarnet_driver.array_get_spec(0);disp("Spectrometer FOUND")% calling spectrometer device id
            catch; obj.equipment.spec = []; disp("Spectrometer ERROR")
            end
            
            try 
                [obj.equipment.WPControl,obj.equipment.WPFig] = obj.WaveplateInit(); disp("Waveplate FOUND")% calling waveplate device id
                obj.WaveplateHome;
            catch; disp("Waveplate ERROR")
            end
        end
        function Name = IDN(obj,instr,ind)
            switch instr
                case 'LI'
                    if ind<=length(obj.equipment.LI)
                        Name=query(obj.equipment.LI(ind),'*IDN?');
                    else
                        Name='';
                    end
                case 'SM'
                    if ind<=length(obj.equipment.SM)
                        Name=query(obj.equipment.SM(ind),'*IDN?');
                    else
                        Name='';
                    end
                case 'ITC503'
                    if ind<=length(obj.equipment.ITC503)
                        Name=query(obj.equipment.ITC503,'V');
                    else
                        Name='';
                    end
            end
        end
        
        %SM Functions
        function r = SM_ReadI(obj,ind,channel)
            if ind > length(obj.equipment.SM)
                r= 'NaN';
            else
                switch channel % select channel
                    case 1
                        r=query(obj.equipment.SM(ind),'print(smua.measure.i())');
                    case 2
                        r=query(obj.equipment.SM(ind),'print(smub.measure.i())');
                end
            end
        end
        function r = SM_ReadV(obj,ind,channel)
            if ind > length(obj.equipment.SM)
                r= 'NaN';
            else
                switch channel % select channel
                    case 1
                        r=query(obj.equipment.SM(ind),'print(smua.measure.v())');
                    case 2
                        r=query(obj.equipment.SM(ind),'print(smub.measure.v())');
                end
            end
        end
        function [] = SM_RampV(obj,ind,channel,Vstart,Vend,stepsize,delay)
            if Vstart>Vend
                stepsize=-abs(stepsize);
            elseif Vstart<Vend
                stepsize=abs(stepsize);
            end
            Vstart
            stepsize
            Vend
            V=[Vstart:stepsize:Vend Vend] %define ramp
            disp(['ramping channel ',num2str(channel),' to V = ',num2str(Vend)])
            switch channel % select channel
                case 1
                    command='smua.source.levelv = ';
                case 2
                    command='smub.source.levelv = ';
            end
            for n=1:length(V) % do ramp
                message=strcat(command,num2str(V(n)));
                fprintf(obj.equipment.SM(ind),message);
                pause(delay)
            end
            disp('ramp done')
        end
        function [] = SM_RampI(obj,ind,channel,Istart,Iend,stepsize,delay)
            if Istart>Iend
                stepsize=-abs(stepsize);
            elseif Istart<Iend
                stepsize=abs(stepsize);
            end
            I=[Istart:stepsize:Iend Iend]; %define ramp
            disp(['ramping channel ',num2str(channel),' to V = ',num2str(Iend)])
            switch channel % select channel
                case 1
                    command='smua.source.leveli = ';
                case 2
                    command='smub.source.leveli = ';
            end
            for n=1:length(I) % do ramp
                message=strcat(command,num2str(I(n)));
                fprintf(obj.equipment.SM(ind),message);
                pause(delay)
            end
            disp('ramp done')
        end
        function [] = SM_SetV(obj,ind,channel,V)
            switch channel % select channel
                case 1
                    command='smua.source.levelv = ';
                case 2
                    command='smub.source.levelv = ';
            end
            message=strcat(command,num2str(V));
            fprintf(obj.equipment.SM(ind),message);
        end
        function [] = SM_SetI(obj,ind,channel,I)
            switch channel % select channel
                case 1
                    command='smua.source.leveli = ';
                case 2
                    command='smub.source.leveli = ';
            end
            message=strcat(command,num2str(I));
            fprintf(obj.equipment.SM(ind),message);
        end
        
        %ITC503 Functions
        function x=queryITC(obj,command)
            clrdevice(obj.equipment.ITC503)
            x=query(obj.equipment.ITC503,command);
        end
        function T = ITC503_ReadT(obj)
            queryITC(obj,'C3');
            x=queryITC(obj,'R1');
            T=str2double(extractAfter(x,1));% read temperature
        end
        function St = ITC503_ReadSetT(obj)
            queryITC(obj,'C3');
            x=queryITC(obj,'R0');
            St=str2double(extractAfter(x,1));% read set temperature
        end
        function stabilizationT = ITC503_SetT(obj,SetT,Tol,Time)
            queryITC(obj,'C3');%Remote Mode
            queryITC(obj,['T' num2str(SetT)]); %set temperature
            queryITC(obj,'A1');
            queryITC(obj,'L1');%Auto-PID
            Setpoint=str2double(extractAfter(queryITC(obj,'R0'),1));
            check=0;
            while check<Time
                T=str2double(extractAfter(queryITC(obj,'R1'),1))% read temperature
                stabilization=0;
                if Setpoint<T+Tol && Setpoint>T-Tol
                    check=check+1;
                else
                    check=0;
                end
                pause(1)
            end
            stabilizationT=1;
        end
        
        %LI Functions
        function [xr1,yt2] = LI_Read(obj,mode)
            xr1 = [];
            yt2 = [];
            % Mode should be either xy or rt
            if ~(strcmp(mode,'xy') || strcmp(mode,'rt') || strcmp(mode,'CH'))
                error('GPIBManager.m: Mode should be either "xy" or "rt".')
                return
            end
            for h = obj.equipment.LI
                if strcmp(mode,'xy')
                    value = query(h, 'SNAP?1,2');
                elseif strcmp(mode,'rt')
                    value = query(h, 'SNAP?3,4');
                elseif strcmp(mode,'CH')
                    value = query(h, 'SNAP?10,11');
                end
                
                % Convert string to double 2x1 array
                value = str2num(value);
                xr1(end+1) = value(1);
                yt2(end+1) = value(2);
            end
        end
        function [FREQ] = LI_FreqRead(obj)
            FREQ=[];
            for h = obj.equipment.LI
                value=query(h, 'FREQ?');
                % Convert string to double 2x1 array
                value = str2num(value);
                FREQ(end+1) = value(1);
            end
        end
        function [] = LI_FreqSet(obj,ind,FREQ)
            h = obj.equipment.LI(ind);
            message=strcat('FREQ ',num2str(FREQ));
            fprintf(h, message);
            % Convert string to double 2x1 array   
        end
        function [SENS] = LI_SensRead(obj)
            senses=["2nV","5nV","10nV","20nV","50nV","100nV","200nV",...
                "500nV","1muV","2muV","5muV","10muV","20muV","50muV",...
                "100muV","200muV","500muV","1mV","2mV","5mV","10mV",...
                "20mV","50mV","100mV","200mV","500mV","1V"];
            for h = obj.equipment.LI
                value=query(h, 'SENS?');
                % Convert string to double 2x1 array
                value = str2num(value);
            end
            SENS=senses(value+1);
        end
        function []=LI_SensSet(obj,ind,SENS)
            h = obj.equipment.LI(ind);
            message=strcat('SENS',num2str(SENS));
            fprintf(h, message);
        end
        %EM Functions
        function EM_Init(obj,ind,mode)
            %fprintf(obj.equipment.EM(ind),'*RST')
            fprintf(obj.equipment.EM(ind),'VOLT:GUAR OFF')% funcion interna?
            fprintf(obj.equipment.EM(ind),'SENS:FUNC ''VOLT''')% mode
            fprintf(obj.equipment.EM(ind),'VOLT:RANG:AUTO OFF')% auto range
            fprintf(obj.equipment.EM(ind),'SYST:ZCOR OFF')% offset corr off
            fprintf(obj.equipment.EM(ind),'SYST:ZCH OFF')% internal zero off
        end
        function EM_Start(obj,ind,mode)
            %fprintf(obj.equipment.EM(ind),'VOLT:GUAR ON')% funcion interna?
            %fprintf(obj.equipment.EM(ind),'SENS:FUNC ''VOLT''')% mode
            %fprintf(obj.equipment.EM(ind),'VOLT:RANG:AUTO ON')% auto range
            %fprintf(obj.equipment.EM(ind),'SYST:ZCOR ON')% offset corr on
            %fprintf(obj.equipment.EM(ind),'SYST:ZCH ON')% internal zero on
            %fprintf(obj.equipment.EM(ind),'SYST:ZCOR OFF')% offset corr off
            %fprintf(obj.equipment.EM(ind),'SYST:ZCH OFF')% internal zero off
        end
        function r = EM_Read(obj,ind)
            EM_Start(obj,ind,1);
            r=query(obj.equipment.EM(ind),'READ?'); %read signal
        end
        
        %MS257 Functions
        function r = MS257move (~,WL)
            [exec_state,output]=system(['.resources\MS257com\MS257com.exe -m ',num2str(WL)]);
            if exec_state==0
                disp(['MS257 Moving wl to ',num2str(WL),' nm'])
                r=str2double(output(17:end));
            else
                r=[];
                error(output)
            end
        end
        function MS257scan (~,startwl,endwl)
            [exec_state,output]=system(['.resources\MS257com\MS257com.exe -s ',num2str(startwl),' ',num2str(endwl)]);
            if exec_state==0
                disp(['MS257 Starting scan from ',num2str(startwl),' to ',num2str(endwl),'.'])
            else
                error(output)
            end
        end
        function MS257abort (~)
            [exec_state,output]=system('.resources\MS257com\MS257com.exe -a');
            if exec_state==0
                disp('MS257 Aborting scan.')
            else
                error(output)
            end
        end
        function MS257status (~)
            [exec_state,output]=system('.resources\MS257com\MS257com.exe -st');
            if exec_state==0
                disp('MS257 Status.')
            else
                error(output)
            end
        end
        
        %LA Functions
        function r = LAread(~)
            [exec_state,output]=system('.resources\NKTcom\NKTcom.exe -r');
            if exec_state==0
                r = output;
            else
                r=999;
                error(output)
            end
        end
        function LAgo(~,PW)
            system(['.resources\NKTcom\NKTcom.exe -p ',num2str(PW)]);
        end
        function LAon(~)
            system('.resources\NKTcom\NKTcom.exe -on');
        end
        function LAoff(~)
            system('.resources\NKTcom\NKTcom.exe -off ');
        end
        
        %Spectrometer Functions
        function [maxWL,counts]=SPread(obj)
            % setting parameters
            inttime =  int64(100);
            xtiming = int64(1);
            scansavg = int64(1);
            smoothing = int64(1);
            param = py.stellarnet_driver.setparam(obj.equipment.spec,inttime,xtiming,scansavg,smoothing); % calling lib to set parameter
            wav = py.stellarnet_driver.array_get_wav(obj.equipment.spec); %getting wavelengths
            spectrum = py.stellarnet_driver.array_spectrum(obj.equipment.spec); % getting spectrum
            data = double(py.array.array('d',py.numpy.nditer(spectrum))); %d is for double, coverting spectrum to matlab type
            x = double(py.array.array('d',py.numpy.nditer(wav))); %d is for double, coverting wavelengths to matlab type
            LF=@(I,x0,gamma,x) I./(1+((x-x0)/gamma).^2);
            xmax=mean(x(data==max(data)));
            datanorm=data./max(data);
            options = optimoptions('lsqcurvefit','MaxFunEvals',10000,'MaxIter',2000,'Display','off');
            Fit=@(p,x) LF(p(1),p(2),p(3),x)+p(4);
            guess=[1 xmax 7 10];lb=[0 400 0 0];ub=[1 1200 20 inf];
            p=lsqcurvefit(Fit,guess,x,datanorm,lb,ub,options);
            maxWL=p(2);
            counts=max(data);
            %figure;hold on;p([1 4])=p([1 4]).*max(data);plot(x,data);plot(x,Fit(p,x));hold off
        end
        
        %Waveplate Functions
        function [] = WaveplateHome(obj)
            timeout = 20; % timeout for waiting the move to be completed
            %h.MoveJog(0,1); % Jog
            
            % Set target position and move
            obj.equipment.WPControl.MoveHome(0,1==1);
        end
        function [WPControl,WPFig] = WaveplateInit(obj)
            % Create Matlab Figure Container
            fpos    = get(0,'DefaultFigurePosition'); % figure default position
            fpos(3) = 650; % figure window size;Width
            fpos(4) = 450; % Height
            
            WPFig = figure('Position', fpos,...
                'Menu','None',...
                'Name','APT GUI',...
                'Visible','on');
            
            % Create ActiveX Controller
            WPControl = actxcontrol('MGMOTOR.MGMotorCtrl.1',[20 20 600 400 ], WPFig);
            
            % Initialize
            % Start Control
            WPControl.StartCtrl;
            
            % Set the Serial Number
            SN = 27253957; % put in the serial number of the hardware
            set(WPControl,'HWSerialNum', SN);
            
            % Indentify the device
            WPControl.Identify;
            pause(5); % waiting for the GUI to load up;
        end
        function [] = WaveplateMove(obj,position)
            timeout = 10; % timeout for waiting the move to be completed
            %h.MoveJog(0,1); % Jog
            
            % Set target position and move
            obj.equipment.WPControl.SetAbsMovePos(0,position);
            obj.equipment.WPControl.MoveAbsolute(0,1==1);
        end
    end
end