Changeset 787 in ETALON

Timestamp:

Jul 27, 2018, 11:43:24 AM (6 years ago)

Author:

moutardier

Message:

huge change in acquisition program

File:

• BPM/initialise_motors.py (modified) (7 diffs)

BPM/initialise_motors.py

@@ -5 +5 @@
 __author__ = 'delerue'
 
-
 from bpm_constants import *
-from scope_function import *
+from read_adc import *
 import numpy as np
 import matplotlib.pyplot as plt
-
+from datetime import datetime
 
 VERBOSE_DEFAULT=0
@@ -17 +16 @@
 write_to_motor_dated_log = data_dir + './write_to_motor_log'
 
-from writeclient import *
+#from writeclient import *
 import time
 
@@ -59 +58 @@
             print('Be carful with disable vertical motor')
         else :
-            move_motor_absolute(2,10000,-10000)
-            while motor_get(2, 'status') != 'READY':
-                time.sleep(1)
-            print('Say ok when motor can be put in position 0')
-            x = raw_input()
-            while not (x == 'OK' or x == 'ok' or x == 'Ok'):
-                print('Say ok when motor can be put in position 0')
-                x = raw_input()
-            move_motor_absolute(2,5000,0)
+            move_motor_absolute(2,2000,0)
             while motor_get(2, 'status') != 'READY':
                 time.sleep(1)
@@ -212 +203 @@
         time.sleep(1)
 
+
 def motor_move_absolute(motor, speed, steps):
     "Moves a motor with the given parameters"
@@ -226 +218 @@
         time.sleep(1)
 
+
 def motor_wait_for_ready(motor,verbose=VERBOSE_DEFAULT):
     motor_status=motor_get(motor,'status')
@@ -257 +250 @@
     if (x == 'NO'):
         print('Please don\'t forbid to move BPM at maxima right when you look from motor side befor to shut down the motor')
-        exit()
-    
+        exit()    
     else:
         motor_enable(1)
@@ -277 +269 @@
     if not os.path.exists(data_name_and_path) :
         os.makedirs(data_name_and_path)
-    data_name_and_path = data_name_and_path+"position_vs_tension_"+bpm_name+"_"+start_time_str+"_"+str(idata)+".txt"
-    while os.path.exists(data_name_and_path) :
+    data_name_and_path = data_name_and_path+"BPM-number_x-motor-step_x-motor-mm_y-motor-step_y-motor_mm_Va_Vb_Vc_Vd_Sum_X_Y_"+bpm_name+"_"+start_time_str+"_"
+    while os.path.exists(data_name_and_path+str(idata)+".txt") :
         idata+=1
-        data_name_and_path = data_dir+"data/position_vs_tension_"+bpm_name+"_"+start_time_str+"_"+str(idata)+".txt"
-    return(data_name_and_path)
+    return(data_name_and_path+str(idata)+".txt")
+
+def write_first_line(data_name_and_path, list_of_data):
+    length = len(list_of_data)
+    fichier = open(data_name_and_path, "w")
+    for index in range(length):
+        fichier.write(list_of_data[index])
+        if index != length - 1:
+            fichier.write(" ")
+    fichier.write("\n")
+    fichier.close()
     
+def horizontal_acquisition(begin, end, number_of_point, # cordinate of the the first and the last point of data taken horizontally and total number of point 
+                           title_bpm_name, # a string that will appear in the name of the datafile 
+                           statistic_number = 1, # nuber of repetition of taking data whith 1 set of parameter, as default no repetition
+                           list_of_data = ["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"], # list of the name of the data that the program will read, as default it read all
+                           bpm_name1 = "BPM_E",bpm_name2 = "BPM_impr",bpm_name3 = "BPM_C",bpm_name4 = None): #name of the BPM used (as default, it the configuration of the BPM when this code was writen 
+
+    data_name_and_path = path(title_bpm_name) #create the name of the datafile
+    write_first_line(data_name_and_path, list_of_data) # write first line of the datafile with the name of the data
+    bpm_list = initialise_bpm(bpm_name1,bpm_name2,bpm_name3,bpm_name4) # initialise BPM (program on read_adc.py)
+    y_position = 0 # set an abirtary position for y since I can't get motor step
+    motor_enable(1)
+    for x_position in np.linspace(begin,end,number_of_point):
+        move_motor_absolute(1,5000,x_position)
+        print('x_motor going to ', x_position)
+        while (motor_get(1, 'status') != 'READY') or (motor_get(2, 'status') != 'READY'):
+            time.sleep(1)
+        time.sleep(1)
+        acquisition_libera(x_position, y_position, data_name_and_path, statistic_number, list_of_data, bpm_list) # read_the data (program on read_adc.py)
+    motor_disable(1)
+
+
+
+def vertical_acquisition(begin, end, number_of_point, # cordinate of the the first and the last point of data taken vertically and total number of point 
+                           title_bpm_name, # a string that will appear in the name of the datafile 
+                           statistic_number = 1, # nuber of repetition of taking data whith 1 set of parameter, as default no repetition
+                           list_of_data = ["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"], # list of the name of the data that the program will read, as default it read all
+                           bpm_name1 = "BPM_E",bpm_name2 = "BPM_impr",bpm_name3 = "BPM_C",bpm_name4 = None): #name of the BPM used (as default, it the configuration of the BPM when this code was writen
     
-def horizontal_acquisition(begin,end,pas,bpm_name):
-    if pas == 0: #position test
-        motor_enable(1)
-        move_motor_absolute(1,10000,begin)
-        while motor_get(1, 'status') != 'READY':
-                time.sleep(1)
-        print("position ok ?")
-        time.sleep(15)
-        move_motor_absolute(1,10000,end)
-        while motor_get(1, 'status') != 'READY':
-            time.sleep(1)
-        print("position ok ?")
-        motor_disable(1)
-    else :
-        if (end - begin)*pas < 0 :
-            pas = -pas
-        data_name_and_path = path(bpm_name)
-        motor_enable(1)
-        L = [[],[],[],[],[]]
-        for i in range(begin, end,pas):
-            move_motor_absolute(1,5000,i)
-            while motor_get(1, 'status') != 'READY':
-                time.sleep(1)
-            time.sleep(1)
-            l = read_scope_and_write(i, data_name_and_path)
-            for j in range(5):
-                L[j].append(l[j])
-        time.sleep(1)
-        move_motor_absolute(1,5000,end) # pour assuré le dernier pas, en aller ou retour
-        while motor_get(1, 'status') != 'READY':
-            time.sleep(1)
-        time.sleep(1)
-        l = read_scope_and_write(end, data_name_and_path)
-        for j in range(5):
-            L[j].append(l[j])
-        motor_disable(1)
-        for i in range(4):
-            plt.plot(L[0],L[i+1])
-        plt.show()
-
-
-
-
-def vertical_acquisition(begin,end,pas,bpm_name):
-    if pas == 0: #position test
-        motor_enable(2)
-        move_motor_absolute(2,10000,begin)
-        while motor_get(2, 'status') != 'READY':
-            time.sleep(1)
-        print("position ok ?")
-        time.sleep(15)
-        move_motor_absolute(2,10000,end)
-        while motor_get(2, 'status') != 'READY':
-                time.sleep(1)
-        print("position ok ?")
-        motor_disable(2)
-    else:
-        if (end - begin)*pas < 0 :
-            pas = -pas
-        data_name_and_path = path(bpm_name)
-        motor_enable(2)
-        L = [[],[],[],[],[]]
-        for i in range(begin,end,pas):
-            move_motor_absolute(2,5000,i)
-            while motor_get(2, 'status') != 'READY':
+    data_name_and_path = path(title_bpm_name) #create the name of the datafile
+    write_first_line(data_name_and_path, list_of_data)  # write first line of the datafile with the name of the data
+    bpm_list = initialise_bpm(bpm_name1,bpm_name2,bpm_name3,bpm_name4) # initialise BPM (program on read_adc.py)
+    x_position = 0 # set an abirtary position for x since I can't get motor step
+    motor_enable(2)
+    for y_position in np.linspace(begin, end, number_of_point):            
+        motor_move_absolute(2,2000,y_position)
+        print('y_moteur going to', y_position)
+        while (motor_get(1, 'status') != 'READY') or (motor_get(2, 'status') != 'READY'):
+            time.sleep(1)
+        time.sleep(1)
+        acquisition_libera(x_position, y_position, data_name_and_path, statistic_number, list_of_data, bpm_list) # read_the data (program on read_adc.py)
+    motor_disable(2)
+
+
+def acquisition_2D(begin1,end1,number_of_point1, # cordinate of the the first and the last point of data taken horizontally and total number of point 
+                   begin2,end2,number_of_point2, # cordinate of the the first and the last point of data taken vertically and total number of point 
+                   title_bpm_name, # a string that will appear in the name of the datafile
+                   statistic_number = 1,# nuber of repetition of taking data whith 1 set of parameter, as default no repetition
+                   list_of_data = ["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"], # list of the name of the data that the program will read, as default it read all
+                   bpm_name1 = "BPM_E",bpm_name2 = "BPM_impr",bpm_name3 = "BPM_C",bpm_name4 = None):#name of the BPM used (as default, it the configuration of the BPM when this code was writen
+    
+    data_name_and_path = path(title_bpm_name)#create the name of the datafile
+    write_first_line(data_name_and_path, list_of_data)# write first line of the datafile with the name of the data
+    bpm_list = initialise_bpm(bpm_name1,bpm_name2,bpm_name3,bpm_name4)# initialise BPM (program on read_adc.py)
+    motor_enable(1)
+    motor_enable(2)
+    for y_position in np.linspace(begin2, end2, number_of_point2):            
+         motor_move_absolute(2,2000,y_position)
+        print('moteur going to', y_position)
+        for x_position in np.linspace(begin1,end1,number_of_point1):
+            move_motor_absolute(1,5000,x_position)
+            print('x_motor going to ', x_position)
+            while (motor_get(1, 'status') != 'READY') or (motor_get(2, 'status') != 'READY'): #waits will the motor move
                 time.sleep(1)
             time.sleep(1)
-            l = read_scope_and_write(i, data_name_and_path)
-            for j in range(5):
-                L[j].append(l[j])
-        time.sleep(1)
-        move_motor_absolute(2,5000,end) # pour assuré le dernier pas, en aller ou retour
-        while motor_get(2, 'status') != 'READY':
-            time.sleep(1)
-        time.sleep(1)
-        l = read_scope_and_write(end, data_name_and_path)
-        for j in range(5):
-            L[j].append(l[j])
-        time.sleep(1)
-        for i in range(4):
-            plt.plot(L[0],L[i+1])
-        plt.show()
-        motor_disable(2)
-
-def acqui_column(begin, end, pas, data_path, other_position): # use in aquisition to read a column
-    if (end - begin)*pas < 0 :
-        pas = -pas
-    for i in range(begin,end,pas):
-        move_motor_absolute(2,5000,i)
-        while motor_get(2, 'status') != 'READY':
-            time.sleep(1)
-        time.sleep(1)
-        l = read_scope_and_write_2D(other_position, i, data_path)
-    time.sleep(1)
-    move_motor_absolute(2,5000,end) # pour assuré le dernier pas, en aller ou retour
-    while motor_get(2, 'status') != 'READY':
-        time.sleep(1)
-    time.sleep(1)
-    read_scope_and_write_2D(other_position, end, data_path)    
-
-def acquisition(begin1,end1,pas1,begin2,end2,pas2,bpm_name):
-    if pas1 == 0 or pas2 == 0: #position test
-        motor_enable(1)
-        motor_enable(2)
-        move_motor_absolute(1,10000,begin1)
-        move_motor_absolute(2,10000,begin2)
-        while (motor_get(1, 'status') != 'READY') or (motor_get(2, 'status') != 'READY'):
-            time.sleep(1)
-        print("position ok ?")
-        time.sleep(15)
-        move_motor_absolute(1,10000,end1)
-        move_motor_absolute(2,10000,end2)
-        while (motor_get(1, 'status') != 'READY') or (motor_get(2, 'status') != 'READY'):
-            time.sleep(1)
-        print("position ok ?")
-        motor_disable(1)
-        motor_disable(2)
-    else:
-        if (end1 - begin1)*pas1 < 0 :
-            pas1 = -pas1
-        if (end2 - begin2)*pas2 < 0 :
-            pas2 = -pas2
-        data_name_and_path = path(bpm_name)
-        motor_enable(1)
-        motor_enable(2)
-        move_motor_absolute(1,10000,begin1)
-        move_motor_absolute(2,10000,begin2)
-        while (motor_get(1, 'status') != 'READY') or (motor_get(2, 'status') != 'READY'):
-            time.sleep(1)
-        print("bagin position ok ?")
-        time.sleep(15)
-#        L = [[],[],[],[],[],[]]
-        for i in range(begin1, end1, pas1):            
-            move_motor_absolute(1,5000,i)
-            while motor_get(1, 'status') != 'READY':
-                time.sleep(1)
-            acqui_column(begin2, end2, pas2, data_name_and_path, i)
-            begin2,end2 = end2,begin2
-        move_motor_absolute(1,5000,end1) # pour assuré le dernier pas, en aller ou retour
-        while motor_get(1, 'status') != 'READY': 
-            time.sleep(1)
-        time.sleep(1)
-        acqui_column(begin2, end2, pas2, data_name_and_path, end1)
-        print("end position ok ?")
-        time.sleep(15)
-        motor_disable(1)
-        motor_disable(2)
-#        for j in range(6):
-#L[j].append(l[j])
-#        for i in range(4):
-#            plt.plot(L[0],L[i+2])
-#        plt.show()
-
-
-
+            acquisition_libera(x_position, y_position, data_name_and_path, statistic_number, list_of_data, bpm_list) # read_the data (program on read_adc.py)
+        begin1,end1 = end1,begin1 #to 
+    motor_disable(1)
+    motor_disable(2)
     
-#controller_hello()
-#move_motor_relative(1,1000,-1000)
-#time.sleep(1)
-#move_motor_relative(1,1000,1000)
+#write_first_line("test.txt", ["test","de","si","ca","marche","bien","ou","pas"])
+
+#acquisition_2D(0,1,3,-3.1415,3.1415,2,"test_acquisition_2D")
+#horizontal_acquisition(0,1,3,"hori")
+#vertical_acquisition(-3.1415,3.1415,2,"test_verti")