Changeset 312 in ETALON for SPESO/ana2015/data_plot2d.m

Timestamp:

Feb 18, 2016, 4:33:49 PM (8 years ago)

Author:

delerue

Message:

Latest version of the SPESO analysis code

File:

• SPESO/ana2015/data_plot2d.m (modified) (1 diff)

SPESO/ana2015/data_plot2d.m

@@ -2 +2 @@
 
 global allData;
-global dvar;
-define_data_var();
 
-idx=find(allData.selected==1);
+[valsused,vals]=data_get_vals(var1,var2);
 
-varidx(1)=mod(var1,1000);
-vartype(1)=floor(var1/1000)*1000;
+if (isempty(vals(1,:))||isempty(vals(2,:)))
+  warning_sp([ ' Warning one of the two variable is empty: ' length(vals(1,:)) ' / '  length(vals(2,:))  ])
+else
+hold on
+if (valsused(4,1))
+   plot(vals(1,:),vals(2,:),'dm')
+   plot(vals(1,:),-10*vals(4,:),'sg')
+   legend('Det. 1','Det. 2 x-10', 'Location','Best')
+   else
+   plot(vals(1,:),vals(2,:),'ob')
+end % there is a var in pos 4
 
-varidx(2)=mod(var2,1000);
-vartype(2)=floor(var2/1000)*1000;
-
-if (~isempty(allData.selection))
-fprintf('Selection: %s \n %d events.\n',allData.selection,length(idx));
-end
-
-vals=zeros(2,length(idx));
-for iloop=1:2
-switch vartype(iloop)
-    case dvar.BEAMINFO
-        vals(iloop,:)=allData.beamInfo(idx,varidx(iloop))';
-    case dvar.DATAINFO
-        switch varidx(iloop)
-            case 1
-                vals(iloop,:)=allData.ch2signalVal(idx);
-            case 2
-                vals(iloop,:)=allData.ch2signalNegPeak(idx);
-        end
-    otherwise
-        warning_sp([ 'The vartype ' vartype ' is unknown.' ]);
-        error('end');
-end %switch
-end %for
-
-plot(vals(1,:),vals(2,:),'o')
 title({ [ get_var_name(var1) ' vs ' get_var_name(var2) ];allData.selection})
 xlabel(get_var_name(var1))
 ylabel(get_var_name(var2))
+end %not empty
 end %function