Changeset 798 in ETALON for BPM

Timestamp:

Sep 6, 2018, 4:40:35 PM (6 years ago)

Author:

delerue

Message:

Added Codesys directory

Location:

BPM

Files:

• print_datas.py (modified) (2 diffs)
• print_datas_lib.py (modified) (6 diffs)

BPM/print_datas.py

@@ -11 +11 @@
 basename='bpm_name_bpm_number_x_motor_step_x_motor_mm_y_motor_step_y_motor_mm_Va_Vb_Vc_Vd_Sum_x_libera_mm_y_libera_mm_'
 
+#Big dataset 1 (along X)
+#graph_fr(basedir+basename+"BPM0-E_BPM1-impr_BPM2-C_100pts-horizontal_3pts-vertical_20180801_1.txt", "x_motor_mm", "x_libera_mm",lin_reg = "yes", residu = "yes", rms = "yes",legende = "yes" )#, centrage = "yes", y_minimal = -20., y_maximal = -10.)
+#plt.show()
+#graph_fr(basedir+basename+"BPM0-E_BPM1-impr_BPM2-C_100pts-horizontal_3pts-vertical_20180801_1.txt", "x_motor_mm", "x_libera_mm",lin_reg = "yes", residu = "yes", rms = "yes",legende = "yes" , x_minimal = -50000, x_maximal = 50000) #, centrage = "yes"
+#plt.show()
+#graph_fr(basedir+basename+"BPM0-E_BPM1-impr_BPM2-C_100pts-horizontal_3pts-vertical_20180801_1.txt", "x_motor_mm", "x_libera_mm",lin_reg = "yes", residu = "yes", rms = "yes",legende = "yes" , x_minimal = -50000, x_maximal = 50000) #, centrage = "yes"
+
+
+#Decrochage
+graph_fr(basedir+basename+"BPM0-E_BPM1-impr_BPM2-C_100pts-horizontal_3pts-vertical_20180801_1.txt", "x_motor_mm", "x_libera_mm",lin_reg = "yes", residu = "yes", rms = "yes",legende = "yes", x_minimal = -200000., x_maximal = -150000. )#, centrage = "yes", y_minimal = -20., y_maximal = -10.)
+graph_fr(basedir+basename+"BPM0-E_BPM1-impr_BPM2-C_100pts-horizontal_3pts-vertical_20180801_1.txt", "x_motor_mm", "Va",lin_reg = "yes", residu = "yes", rms = "yes",legende = "yes", x_minimal = -200000., x_maximal = -150000. )#, centrage = "yes", y_minimal = -20., y_maximal = -10.)
+graph_fr(basedir+basename+"BPM0-E_BPM1-impr_BPM2-C_100pts-horizontal_3pts-vertical_20180801_1.txt", "x_motor_mm", "Vb",lin_reg = "yes", residu = "yes", rms = "yes",legende = "yes", x_minimal = -200000., x_maximal = -150000. )#, centrage = "yes", y_minimal = -20., y_maximal = -10.)
+graph_fr(basedir+basename+"BPM0-E_BPM1-impr_BPM2-C_100pts-horizontal_3pts-vertical_20180801_1.txt", "x_motor_mm", "Vc",lin_reg = "yes", residu = "yes", rms = "yes",legende = "yes", x_minimal = -200000., x_maximal = -150000. )#, centrage = "yes", y_minimal = -20., y_maximal = -10.)
+graph_fr(basedir+basename+"BPM0-E_BPM1-impr_BPM2-C_100pts-horizontal_3pts-vertical_20180801_1.txt", "x_motor_mm", "Vd",lin_reg = "yes", residu = "yes", rms = "yes",legende = "yes", x_minimal = -200000., x_maximal = -150000. )#, centrage = "yes", y_minimal = -20., y_maximal = -10.)
+plt.show()
+
+#Big data set 2 (along y)
+#graph_fr(basedir+basename+"BPM0-E_BPM1-impr_BPM2-C_3pts-horizontal_100pts-vertical_20180802_0.txt", "y_motor_step", "y_libera_mm",lin_reg = "yes", residu = "yes", rms = "yes",legende = "yes" , save_figure = 'y_20180802_0')#, centrage = "yes", y_minimal = -20., y_maximal = -10.)
+#plt.show()
+#graph_fr(basedir+basename+"BPM0-E_BPM1-impr_BPM2-C_3pts-horizontal_100pts-vertical_20180802_0.txt", "y_motor_mm", "y_libera_mm",lin_reg = "yes", residu = "yes", rms = "yes",legende = "yes" , centrage = "yes", y_minimal = -20., y_maximal = -10., save_figure = 'y_zoom_20180802_0' )
+#plt.show()
+
+
+#Ultra Big data set 3 (along x)
+#graph_fr(basedir+basename+"BPM0-E_BPM1-impr_BPM2-C_100pts-horizontal_5pts-vertical_20180803_0.txt", "x_motor_step", "x_libera_mm",lin_reg = "yes", residu = "yes", rms = "yes",legende = "yes", save_figure = 'x_20180803_0' )#, centrage = "yes", y_minimal = -20., y_maximal = -10.)
+#plt.show()
+
+
+#tests de decrochage
 #graph_fr(basedir+basename+"test_decrochage_BPM0-E_BPM1-impr_BPM2-C_High_Stat_20180803_0.txt", "y_motor_step", "y_libera_mm",lin_reg = "yes", residu = "yes", rms = "yes",legende = None )#, centrage = "yes", y_minimal = -20., y_maximal = -10.)
 #graph_fr(basedir+basename+"test_decrochage_BPM0-E_BPM1-impr_BPM2-C_High_Stat_20180803_1.txt", "y_motor_step", "y_libera_mm",lin_reg = "yes", residu = "yes", rms = "yes",legende = None )#, centrage = "yes", y_minimal = -20., y_maximal = -10.)
@@ -20 +49 @@
     print(ifile)
     graph_fr(basedir+basename+"test_decrochage_BPM0-E_BPM1-impr_BPM2-C_20180803_"+str(ifile)+".txt", "y_motor_step", "y_libera_mm",lin_reg = "yes", residu = "yes", rms = "yes",legende = None )#, centrage = "yes", y_minimal = -20., y_maximal = -10.)
-""""
+"""
 
 #plt.show(block = False)

BPM/print_datas_lib.py

@@ -8 +8 @@
 global_index_figure = 0
 
-figure_size_in_hundreds_pixels=(8,6)
+figure_size_in_hundreds_pixels=(12,8)
 
 
@@ -262 +262 @@
     if y_data_name == "x_motor_step" or y_data_name == "y_motor_step":
         Y = [i[y_data_number]/1000000. for i in data]
-    if legende is None:
+    if not(lin_reg is None):
         plt.plot(X,Y, color+".")
-        plt.errorbar(X,Y,yerr = rms, ecolor = color)
-    else:
-        plt.plot(X,Y, color+".", label = legende)
-        plt.errorbar(X,Y,yerr = rms, ecolor = color)
-    if not(lin_reg is None):
         slope, intercept, r_value, p_value, std_err = stats.linregress(X,Y) # return tuple (pente,ordonnee a l'origine, coef de correlation, p-value, erreur standard de l'estimation)
         if legende is None:
@@ -274 +269 @@
             plt.plot(X,[slope*i + intercept for i in X], color+"--")
         else:
-            plt.plot(X,[slope*i + intercept for i in X], color+"--", label = legende + " : regression lineaire, \nerreur = "+str(round(std_err,3))+",\ncoefficient de correlation  = "+str(round(r_value,3))+ "\npente de la regression lineaire = " + str(round(slope,3)))            
-            plt.legend(fontsize=15)
+            plt.plot(X,[slope*i + intercept for i in X], color+"--", label = legende + " : reg. lin, \nerr = "+str(round(std_err,3))+",\ncoef correl = "+str(round(r_value,3))+ "\npente = " + str(round(slope,3)))            
+            plt.legend(fontsize=15,bbox_to_anchor=(1.11, 1.11))            
+    else:
+        if legende is None:
+            plt.plot(X,Y, color+".")
+            plt.errorbar(X,Y,yerr = rms, ecolor = color)
+        else:
+            plt.plot(X,Y, color+".", label = legende)
+            plt.errorbar(X,Y,yerr = rms, ecolor = color)
+
     x_label = plot_legend_fr(x_data_name)
     y_label = plot_legend_fr(y_data_name)
@@ -512 +515 @@
 
 
-def graph_fr(fichier, x_data_name, y_data_name, lin_reg = None, residu = None, rms = None, x_begin_mm = None, x_end_mm = None, y_begin_mm = None, y_end_mm = None, x_minimal = None, x_maximal = None, y_minimal = None,  y_maximal = None, centrage = None, one_bpm = None, legende = "yes"):
+def graph_fr(fichier, x_data_name, y_data_name, lin_reg = None, residu = None, rms = None, x_begin_mm = None, x_end_mm = None, y_begin_mm = None, y_end_mm = None, x_minimal = None, x_maximal = None, y_minimal = None,  y_maximal = None, centrage = None, one_bpm = None, legende = "yes", save_figure = None):
     (list_data_name,data) = read_data(fichier)
     if not(x_data_name in list_data_name):
@@ -526 +529 @@
     linear_fonction_on_data(list_data_name, data, x_data_name, begin_mm = x_begin_mm, end_mm = x_end_mm)
     linear_fonction_on_data(list_data_name, data, y_data_name, begin_mm = y_begin_mm, end_mm = y_end_mm)
-    data = cut_data(list_data_name, data, x_data_name, inf = x_maximal, equal = x_maximal)
-    data = cut_data(list_data_name, data, x_data_name, sup = x_minimal, equal = x_minimal)
-    data = cut_data(list_data_name, data, y_data_name, inf = y_maximal, equal = y_maximal)
-    data = cut_data(list_data_name, data, y_data_name, sup = y_minimal, equal = y_minimal)
+    print("Data size")
+    print(len(data))
+    if not(x_maximal is None):
+        data = cut_data(list_data_name, data, x_data_name, inf = x_maximal, equal = x_maximal)
+        print("x_maximal "+str(x_maximal))
+        print(len(data))
+    if not(x_minimal is None):
+        data = cut_data(list_data_name, data, x_data_name, sup = x_minimal, equal = x_minimal)
+        print("x_minimal "+str(x_minimal))
+        print(len(data))
+    if not(y_maximal is None):
+        print("y_maximal "+str(y_maximal))
+        data = cut_data(list_data_name, data, y_data_name, inf = y_maximal, equal = y_maximal)
+        print(len(data))
+    if not(y_minimal is None):
+        print("y_minimal "+str(y_minimal))
+        data = cut_data(list_data_name, data, y_data_name, sup = y_minimal, equal = y_minimal)
+        print(len(data))
     if data == []:
         print("cut to high")
@@ -551 +568 @@
                 if not(legende is None):
                     the_legend=all_bpm_name[bpm_number]
+                    print(the_legend)
                 else:
                     the_legend=None
                 print_graph_fr(list_data_name, cut_data(list_data_name, data, "bpm_name", equal =  all_bpm_name[bpm_number]), x_data_name, y_data_name, global_index_figure, color = color_list[bpm_number], legende = the_legend, lin_reg = lin_reg, rms = rms, centrage = centrage)
+                if not(save_figure is None):
+                    plt.savefig(save_figure+'_'+str(global_index_figure)+'.png')
+                    print("Figure saved")
+                    print(save_figure+'_'+str(global_index_figure)+'.png')
                 if not(residu is None):
                     print_residu_fr(list_data_name,  cut_data(list_data_name, data, "bpm_name", equal = all_bpm_name[bpm_number]), x_data_name, y_data_name, global_index_figure + 1, color_list[bpm_number], legende = the_legend, rms = rms, centrage = centrage)
+                    if not(save_figure is None):
+                        plt.savefig(save_figure+'_residu_'+str(global_index_figure)+'.png')
             global_index_figure += 2
         else:
             all_oposit_y_value = all_value_of_data(list_data_name, data,oposit_y_data_name)
-            if not(legende is None):
-                the_legend=all_bpm_name[bpm_number]
-            else:
-                the_legend=None
             for index_value in range(len(all_oposit_y_value)):
                 data_cut =  cut_data(list_data_name, data, oposit_y_data_name, equal = all_oposit_y_value[index_value])
                 for bpm_number in range(len(all_bpm_name)):
+                    if not(legende is None):
+                        the_legend=all_bpm_name[bpm_number]
+                    else:
+                        the_legend=None
                     #print(cut_data(list_data_name, data_cut, "bpm_name", equal =  all_bpm_name[bpm_number]))
                     print_graph_fr(list_data_name, cut_data(list_data_name, data_cut, "bpm_name", equal =  all_bpm_name[bpm_number]), x_data_name, y_data_name, global_index_figure + 2*index_value, color = color_list[bpm_number], legende = the_legend, lin_reg = lin_reg, rms = rms, title = oposit_y_data_name + " = " + str(all_oposit_y_value[index_value]), centrage = centrage)
                     if not(residu is None):
-                        print_residu_fr(list_data_name, cut_data(list_data_name, data_cut, "bpm_name", equal = all_bpm_name[bpm_number]), x_data_name, y_data_name, global_index_figure + 2*index_value + 1, color_list[bpm_number], legende = the_legend, rms = rms, title = oposit_y_data_name + " = " + str(all_oposit_y_value[index_value]), centrage = centrage)
-            global_index_figure += 2*len(all_oposit_y_value)
+                        print_residu_fr(list_data_name, cut_data(list_data_name, data_cut, "bpm_name", equal = all_bpm_name[bpm_number]), x_data_name, y_data_name, global_index_figure + 2*index_value + 1, color_list[bpm_number], legende = the_legend, rms = rms, title = oposit_y_data_name + " = " + str(all_oposit_y_value[index_value]), centrage = centrage)                        
+                if not(save_figure is None):
+                    plt.figure(global_index_figure + 2*index_value)
+                    plt.savefig(save_figure+'_'+str(global_index_figure + 2*index_value)+'.png')
+                    print("Figure saved")
+                    print(save_figure+'_'+str(global_index_figure + 2*index_value)+'.png')
+                    if not(residu is None):
+                        plt.figure(global_index_figure + 2*index_value+1)
+                        plt.savefig(save_figure+'_residu_'+str(global_index_figure + 2*index_value)+'.png')
+                global_index_figure += 2*len(all_oposit_y_value)
     else: # need to be done, same as befor with "bpm_number" to replace "bpm_name" 
         print("ERROR TYPE OF DATA")