[chronojump] Showing corrected speed in allSet graph

[Xavier Padullés <xpadulles src gnome org>]

commit 4c77a1d74be156f7c1f24ca46662c7e7448fbadd
Author: Xavier Padullés <x padulles gmail com>
Date:   Tue Jul 18 13:36:35 2017 +0200

    Showing corrected speed in allSet graph

 encoder/graph.R |   30 +++++++++++++++++++++++++++---
 1 files changed, 27 insertions(+), 3 deletions(-)
---
diff --git a/encoder/graph.R b/encoder/graph.R
index ef1252a..9db2d0b 100644
--- a/encoder/graph.R
+++ b/encoder/graph.R
@@ -2902,11 +2902,35 @@ doProcess <- function(options)
                        }
 
                        #showSpeed only on gravitatory until speed is fixed on this experimental graph
-                       if (showSpeed && ! isInertial(op$EncoderConfigurationName)) {
+                       #if (showSpeed && ! isInertial(op$EncoderConfigurationName)) {
+                       if (showSpeed ) {
                                par(new=T)
-                               ylimHeight = max(abs(range(speed$y)))
+                               
+                               correctedSpeed = speed$y[0:(curves[1,1] -1)]    #The first part of the graph 
isn't corrected
+                               #The speed at the excentric phase should be negative
+                               #Seeing if the speed at the excentric phase is positive or negative
+                               if(speed$y[curves[1,1] + 300] > 0)
+                               {
+                                       sign = -1                               #If the speed is positive it 
must be changeed the sign of it
+                               } else 
+                               {
+                                       sign = 1                                #If the speed is negative it 
must to be manteined
+                               }
+                               
+                               
+                               #adding the consecutive repetitions to the correctedSpeed
+                               for (n in 2:length(curves[,1]))
+                               {
+                                       correctedSpeed = c(correctedSpeed, sign * (speed$y[curves[n -1, 
1]:(curves[n,1] -1 )]))
+                                       sign = sign * (-1)
+                               }
+                               
+                               correctedSpeed = c(correctedSpeed, sign 
*speed$y[curves[length(curves[,1]),1]: length(speed$y)])
+                               print(paste("CorrectedSpeed length =", length(correctedSpeed)))
+                               ylimHeight = max(abs(range(correctedSpeed)))
                                ylim=c(- 1.05 * ylimHeight, 1.05 * ylimHeight)  #put 0 in the middle, and 
have 5% margin at each side
-                               plot(speed$y, col=cols[1], ylim=ylim, type="l", xlab="",ylab="",axes=F)
+                               print(paste("ylim =",ylim))
+                               plot(correctedSpeed, col=cols[1], ylim=ylim, type="l", xlab="",ylab="",axes=F)
                                axis(4, col=cols[1], lty=lty[1], line=axisLineRight, lwd=1, padj=-.5)
                                axisLineRight = axisLineRight +2
                        }