if i==0: print "first point" m=control.find_laser() #m=(0,0,0,0) else: print "in the quartz" m=control.find_laser([steps[-1]],move_route[0]) #m=(0,0,0,0) pos=((m[0]+m[1])/2,(m[2]+m[3])/2) control.move_pd(pos[0],pos[1]) mpd=[] rpd=[] for kkk in xrange(10): r=control.read_time_pico_dmm() mpd.append(float(r[1])) rpd.append(float(r[2])) print "%s %s %s %s %s %s %s %s %s" % (r[0],r[1],r[2],pos[0],pos[1],control.global_lr[0],control.global_lr[1],control.global_pd[0],control.global_pd[1]) f.write("%s %s %s %s %s %s %s %s %s\n" % (r[0],r[1],r[2],pos[0],pos[1],control.global_lr[0],control.global_lr[1],control.global_pd[0],control.global_pd[1])) f.flush() control.move_pd(-pos[0],-pos[1]) rat=mean(mpd)/mean(rpd) if rat0<1e-6: print "Wrong rat0!" sys.exit(1) ir=math.pow(rat/rat0/RR/math.pow(bulk,L),1.0/bounce) corr0=Rs1*Rs1*math.pow(bulk,2*L)*math.pow(ir,2*bounce)/bounce print "Internal Reflectivity: ",(ir-corr0)*100,"%" reflect.append(ir) move_route.append(m)
#!/usr/bin/python import control print " ".join(control.read_time_pico_dmm())