Changeset 777 in ETALON

Timestamp:

Jun 14, 2018, 11:25:20 AM (6 years ago)

Author:

moutardier

Message:

add program to initilase_motor.py to generat data_path and add exeception to acquisition's program to test position

Location:

BPM

Files:

• initialise_motors.py (modified) (3 diffs)
• motor_control.py (modified) (2 diffs)
• print_datas.py (modified) (4 diffs)
• scope_function.py (modified) (2 diffs)

BPM/initialise_motors.py

@@ -58 +58 @@
         if (x == 'no' or x == 'n' or  x == 'NO' or x == 'N'):
             print('Be carful with disable vertical motor')
-            exit()
-    msg = str(motor) + 'CD'
-    motor_send(msg, verbose)
+        else :
+            msg = str(motor) + 'CD'
+            motor_send(msg, verbose)
+    else :
+        msg = str(motor) + 'CD'
+        motor_send(msg, verbose)
 
 
@@ -254 +257 @@
     motor_disable(1)
 
-def horizontal_aquisition(begin,end,pas,bpm_name):
+def path(bpm_name): # create the path of the file what will contain data
     idata = 0
-    start_time_str=datetime.now().strftime("%Y%m%d_%H%M%S")
-    data_name_and_path = data_dir+"data/position_vs_tension_"+bpm_name+"_"+start_time_str+"_"+str(idata)+".txt"
+    start_time_str=datetime.now().strftime("%Y%m%d")
+    data_name_and_path = data_dir+"data/"
+    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) :
         idata+=1
         data_name_and_path = data_dir+"position_vs_tension_"+bpm_name+start_time_str+"_"+str(idata)+".txt"
-    fichier  = open(data_name_and_path, "w")
-    motor_enable(1)
-    L = [[],[],[],[],[]]
-    for i in range(begin, end,pas):
-        move_motor_absolute(1,5000,i)
+    return(data_name_and_path)
+    
+    
+def horizontal_aquisition(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)
-        time.sleep(1)
-        l = read_scope_and_write(i, fichier)
+        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])
-    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, fichier)
-    for j in range(5):
-        L[j].append(l[j])
-    time.sleep(1)
-    motor_disable(1)
-    fichier.close()
-#    for i in range(4):
-#         plt.plot(L[0],L[i+1])
-#    plt.show()
+        motor_disable(1)
+        for i in range(4):
+            plt.plot(L[0],L[i+1])
+        plt.show()
 
 
@@ -292 +314 @@
 
 def vertical_aquisition(begin,end,pas,bpm_name):
-    idata = 0
-    start_time_str=datetime.now().strftime("%Y%m%d_%H%M%S")
-    data_name_and_path = data_dir+"data/position_vs_tension_"+bpm_name+"_"+start_time_str+"_"+str(idata)+".txt"
-    while os.path.exists(data_name_and_path) :
-        idata+=1
-        data_name_and_path = data_dir+"position_vs_tension_"+bpm_name+start_time_str+"_"+str(idata)+".txt"
-    fichier  = open(data_name_and_path, "w")
-    motor_enable(2)
-    L = [[],[],[],[],[]]
-    for i in range(begin, end,pas):
-        move_motor_absolute(1,5000,i)
+    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)
-        time.sleep(1)
-        l = read_scope_and_write(i, fichier)
+        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':
+                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)
-    
-    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, fichier)
-    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)
-    fichier.close()
-
-
-"""
-def aquisition(begin1,end1,begin2,end2,pas1,pas2,bpm_name):
-    idata = 0
-    start_time_str=datetime.now().strftime("%Y%m%d_%H%M%S")
-    data_name_and_path = data_dir+"data/position_vs_tension_"+bpm_name+"_"+start_time_str+"_"+str(idata)+".txt"
-    while os.path.exists(data_name_and_path) :
-        idata+=1
-        data_name_and_path = data_dir+"position_vs_tension_"+bpm_name+start_time_str+"_"+str(idata)+".txt"
-    fichier  = open(data_name_and_path, "w")
-    motor_enable(1)
-    motor_enable(2)
-    L = [[],[],[],[],[]]
-    for i in range(begin1, end1,pas1):
-        move_motor_absolute(1,5000,i)
-        while motor_get(1, 'status') != 'READY':
+        time.sleep(1)
+        for i in range(4):
+            plt.plot(L[0],L[i+1])
+        plt.show()
+        motor_disable(2)
+
+
+
+
+def aquisition(begin1,end1,pas2,begin2,end2,pas1,bpm_name):
+    if pas == 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') and (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') and (motor_get(2, 'status') != 'READY'):
+            time.sleep(1)
+        print("position ok ?")
+        motor_disable(1)
+        motor_disable(2)
+    else:
+        if (end - begin)*pas < 0 :
+            pas = -pas
+        data_name_and_path = path(bpm_name)
+        motor_enable(1)
+        motor_enable(2)
+        L = [[],[],[],[],[]]
+        for i in range(begin1, end1,pas1):
+            move_motor_absolute(1,5000,i)
+            while motor_get(1, 'status') != 'READY':
                 time.sleep(1)
-        time.sleep(1)
-        for j in range(begin2, end2,pas2):
-            move_motor_absolute(2,5000,i)
-            while motor_get(2, 'status') != 'READY':
-                time.sleep(1)
-            time.sleep(1)
-            l = read_scope_and_write(i, fichier)
-            for k in range(5):
-                L[k].append(l[k])
-            time.sleep(1)
-    move_motor_absolute(1,5000,end1)    
-    move_motor_absolute(2,5000,end2) # pour assuré le dernier pas, en aller ou retour
-    while (motor_get(1, 'status') != 'READY') & (motor_get(2, 'status') != 'READY'): 
-        time.sleep(1)
-    time.sleep(1)
-    l = read_scope_and_write(end, fichier)
-    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()
-    fichier.close()
-"""
+            time.sleep(1)
+            for j in range(begin2, end2,pas2):
+                move_motor_absolute(2,5000,i)
+                while motor_get(2, 'status') != 'READY':
+                    time.sleep(1)
+                time.sleep(1)
+                l = read_scope_and_write_2D(i, data_name_and_path)
+                for k in range(5):
+                    L[k].append(l[k])
+                time.sleep(1)
+        move_motor_absolute(1,5000,end1)    
+        move_motor_absolute(2,5000,end2) # pour assuré le dernier pas, en aller ou retour
+        while (motor_get(1, 'status') != 'READY') & (motor_get(2, 'status') != 'READY'): 
+            time.sleep(1)
+        time.sleep(1)
+        l = read_scope_and_write_2D(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()
+
 
 

BPM/motor_control.py

@@ -15 +15 @@
     time.sleep(1)
 time.sleep(5)
-    
+""" 
 move_motor_absolute(2,5000,0)
 while motor_get(1, 'status') != 'READY':
     time.sleep(1)
-
+"""
 #controller_hello()
 #move_motor_relative(1,1000,-1000)
@@ -34 +34 @@
 
 
-turn_on()
+#turn_on()
 
 vertical_aquisition(-110000,90000,4000,"ch1-3_bpm_ref_ch2-4_bpm_impr_vertical_acquisition")
 
-vertical_aquisition(90000,-110000,-4000,"ch1-3_bpm_ref_ch2-4_bpm_impr_vertical_acquisition")
+#vertical_aquisition(90000,-110000,4000,"ch1-3_bpm_ref_ch2-4_bpm_impr_vertical_acquisition")
 
-#vertical_aquisition(0,4000,2000,"test")
+#vertical_aquisition(0,-4000,2000,"test")
+#vertical_aquisition(4000,0,2000,"test2")
 
 #horizontal_aquisition(400000,0,-4000,"ch1-3_bpm_ref_ch2-4_bpm_impr")    
 
-turn_off()
+#turn_off()
 
 #fichier = open("24-05-2018_aquisition_from_0_to_400000_step_of_1000.txt", "r")

BPM/print_datas.py

@@ -8 +8 @@
         L.append(float(i))
     return(L)
+
+def product_list(L,x):
+    Lx = []
+    for i in L:
+        Lx.append(i*x)
+    return(Lx)
+
 
 def print_data(fichier):
@@ -42 +49 @@
     for i in range(5):
         L[i] = float_list(L[i])
+    L[0] = product_list(L[0],10**-6)
+    for i in range(1,5):
+        L[i] = product_list(L[i],1000)
     print_graph(L)
 
@@ -66 +76 @@
     
     plt.legend()
-    plt.xlabel("length in motor step")
-    plt.ylabel("tension in V")
+    plt.xlabel("length in million motor step")
+    plt.ylabel("tension in mV")
     plt.show()
 
@@ -80 +90 @@
 
 
-print_data("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_20180611_160153_0.txt")
-print_data("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_20180611_163036_0.txt")
+#print_data("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_20180611_160153_0.txt")
+#print_data("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_20180611_163036_0.txt")
+
+#print_data("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_vertical_acquisition_20180612_0.txt")
+print_data("data/position_vs_tension_ch1-3_bpm_ref_ch2-4_bpm_impr_vertical_acquisition_20180612_0.txt")

BPM/scope_function.py

@@ -68 +68 @@
     return(L)
     
-def read_scope_and_write(position, fichier):
+def read_scope_and_write(position, path):
+    fichier  = open(path, "a")
     fichier.write(str(position)+' ')
     L = [position]
@@ -78 +79 @@
         fichier.write(valeur+' ')
     fichier.write('\n')
+    fichier.close()
     return(L)