source: ETALON/SPESO/ana2015/data_get_vals.m @ 622

function [valsused,vals]=data_get_vals(var1,var2)

%%%Each variable constant is associated to a value through this function
%%% ND 2016


global allData;
global dvar;
define_data_var();

check_data_vars_plotted();

idx=find(allData.selected==1);

nloop=0;
for iloop=1:2       
  if (nargin>=iloop)
     if (iloop==1)
       varval(iloop)=var1;
     elseif (iloop==2)
     varval(iloop)=var2;
     else
       error('Error in for loop');
     end
     if (varval(iloop)>dvar.LOG_VAR)
       log_requested(iloop)=1;
       varval(iloop)=varval(iloop)-dvar.LOG_VAR;
     else
       log_requested(iloop)=0;
     end  
     varidx(iloop)=mod(varval(iloop),1000);
     vartype(iloop)=floor(varval(iloop)/1000)*1000;
     nloop=iloop;
   end %if
end %for






if (~isempty(allData.selection))
fprintf('Selection: %s \n %d events.\n',allData.selection,length(idx));
end

if (isempty(idx))
warning_sp([ 'Plot2D: The selection eliminated all events!!!' allData.selection ])
valsused=zeros(4,1);
vals=zeros(4,1);
else

valsused=zeros(4,1);
vals=zeros(4,length(idx));
for iloop=1:nloop
switch vartype(iloop)
    case dvar.BEAMINFO
        vals(iloop,:)=allData.beamInfo(idx,varidx(iloop))'; %'
    case dvar.DATAINFO        
            switch varval(iloop)
            case dvar.DATA_VAL_LPM
            warning_sp('Using dvar.DATA_VAL_LPM is obsolete; use dvar.DATA1_VAL_LPM or dvar.DATA2_VAL_LPM instead!!!')
            valsused(iloop,1)=1;
            valsused(iloop+2,1)=1;
                vals(iloop,:)=allData.ch1signalValLPM(idx);
                vals(iloop+2,:)=allData.ch2signalValLPM(idx);
            case dvar.DATA_PEAK_LPM
            warning_sp('Using dvar.DATA_PEAK_LPM is obsolete; use dvar.DATA1_PEAK_LPM or dvar.DATA2_PEAK_LPM instead!!!')
            valsused(iloop,1)=1;
            valsused(iloop+2,1)=1;
                vals(iloop,:)=allData.ch1signalNegPeakLPM(idx);
                vals(iloop+2,:)=allData.ch2signalNegPeakLPM(idx);
            case dvar.DATA_VAL_SPM
            warning_sp('Using dvar.DATA_VAL_SPM is obsolete; use dvar.DATA1_VAL_SPM or dvar.DATA2_VAL_SPM instead!!!')
            valsused(iloop,1)=1;
            valsused(iloop+2,1)=1;
                vals(iloop,:)=allData.ch1signalValSPM(idx);
                vals(iloop+2,:)=allData.ch2signalValSPM(idx);
            case dvar.DATA_PEAK_SPM
            warning_sp('Using dvar.DATA_PEAK_SPM is obsolete; use dvar.DATA1_PEAK_SPM or dvar.DATA2_PEAK_SPM instead!!!')
            valsused(iloop,1)=1;
            valsused(iloop+2,1)=1;
                vals(iloop,:)=allData.ch1signalNegPeakSPM(idx);
                vals(iloop+2,:)=allData.ch2signalNegPeakSPM(idx);

            case dvar.DATA1_VAL_LPM
                valsused(iloop,1)=1;
                vals(iloop,:)=allData.ch1signalValLPM(idx);
            case dvar.DATA1_PEAK_LPM
                valsused(iloop,1)=1;
                vals(iloop,:)=allData.ch1signalNegPeakLPM(idx);
            case dvar.DATA1_VAL_SPM
                valsused(iloop,1)=1;
                vals(iloop,:)=allData.ch1signalValSPM(idx);
            case dvar.DATA1_PEAK_SPM
                valsused(iloop,1)=1;
                vals(iloop,:)=allData.ch1signalNegPeakSPM(idx);

            case dvar.DATA2_VAL_LPM
                valsused(iloop,1)=1;
                vals(iloop,:)=allData.ch2signalValLPM(idx);
            case dvar.DATA2_PEAK_LPM
                valsused(iloop,1)=1;
                vals(iloop,:)=allData.ch2signalNegPeakLPM(idx);
            case dvar.DATA2_VAL_SPM
                valsused(iloop,1)=1;
                vals(iloop,:)=allData.ch2signalValSPM(idx);
            case dvar.DATA2_PEAK_SPM
                valsused(iloop,1)=1;
                vals(iloop,:)=allData.ch2signalNegPeakSPM(idx);
            % Select time in minutes    
            case dvar.DATA_TIME_MINUTES
                valsused(iloop,1)=1;
                %convert time HHMMSS to minutes
                vals(iloop,:)=( fix(allData.hoursNum(idx)/10000)*60+...
                                fix(mod(allData.hoursNum(idx),10000)/100)+...
                                mod(allData.hoursNum(idx),100)/60);
            case dvar.DATA_TIME_30
                valsused(iloop,1)=1;
                %convert time HHMMSS to minutes
                vals(iloop,:)=mod(( fix(allData.hoursNum(idx)/10000)*60+...
                                fix(mod(allData.hoursNum(idx),10000)/100)+...
                                mod(allData.hoursNum(idx),100)/60),30);
            
            otherwise
            warning_sp([ 'The varidx ' num2str(varidx(iloop)) ' is unknown.' ]);
                error('end');
            end %switch
    otherwise
            warning_sp([ 'The vartype ' num2str(vartype) ' is unknown.' ]);
        error('end');
end %switch
if (log_requested(iloop)==1)
  vals(iloop,:)=log(abs(vals(iloop,:)));
 end
end %for

end %function