From 56bd8e4647b687b0490f81e7e25a948dbe8c504e Mon Sep 17 00:00:00 2001
From: guillaume <guillaume@E112-PC.FEUPsig.fe.up.pt>
Date: Tue, 21 Mar 2017 15:53:52 +0000
Subject: [PATCH] update AnalyzeThis

---
 @AnalyzeThis/plot.m | 45 +++++++++++++++++++++++++++++++++++++++++++--
 1 file changed, 43 insertions(+), 2 deletions(-)

diff --git a/@AnalyzeThis/plot.m b/@AnalyzeThis/plot.m
index 205113e..ecd7c35 100644
--- a/@AnalyzeThis/plot.m
+++ b/@AnalyzeThis/plot.m
@@ -1,5 +1,5 @@
 function plot(obj)
-
+figure
 XY = obj.TraceIn.Data;
 s  = subplot(2,1, 1);
 s.Parent.Name = [obj.TraceIn.Title, ' - Results AnalyzeThis'];
@@ -29,18 +29,59 @@ function plot(obj)
 hold on
 stem(PL.Data(:,1), PL.Data(:,2), 'r')
 hold off
-    
+
 xlabel([obj.TraceIn.AxisX.Label, ' / ', obj.TraceIn.AxisX.Unit]);
 ylabel([obj.TraceIn.AxisY.Label, ' / ', obj.TraceIn.AxisY.Unit]);
 
 c = uicontextmenu;
 plotline.UIContextMenu = c;
+
+uD = plotline.UserData;
+if ~isempty(uD)
+    dothat
+end
 % Set c to be the plot line's UIContextMenu
 % Create menu items for the uicontextmenu
 uimenu(c,'Label','FiguresOfMerits' ,'Callback',@getFOM);
     function getFOM(~,~)
+        h = msgbox('Select peak start', '','modal');
+        uiwait(h)
+        axis(sb);
+        [TmS, ~] = ginput(1);
+        if isempty(TmS)
+            return
+        end
+        h = msgbox('Select peak end', '','modal');
+        uiwait(h)
         axis(sb);
+        [TmE, ~] = ginput(1);
+        if isempty(TmE)
+            return
+        end
+        XLim(1) = min(TmS, TmE);
+        XLim(2) = max(TmS, TmE);
+        IdX(1)  = find(XY(:,1) <= XLim(1), 1, 'last');
+        IdX(2)  = find(XY(:,1) >= XLim(2), 1, 'first');
         
+        XY = obj.XY;
+        uD = plotline.UserData;
+        if isempty(uD)
+            Id = 1;
+            [Imax, idImax] = max(XY(IdX(1):IdX(2), 2));
+            tmax = XY(idImax+IdX(1)-1, 1);
+            x = XY(IdX(1):IdX(2),1); y = XY(IdX(1):IdX(2),2);
+            M0 = trapz(x, y);
+            M1 = trapz(x, x.*y)/M0;
+            M2 = trapz(x, (x-M1).^2.*y)/M0;
+            T     = table;
+            T.Id  = Id;
+            T.Imax = Imax;
+            T.tmax = tmax;
+            T.M0   = M0;
+            T.M1   = M1;
+            T.M2   = M2;
+            disp(T)
+        end
     end
 
 end