import os import traceback import math print "\nStarting delay scan - Parameters: ", print START,END,SIZE,ENERGIES,RUNTYPE ############################################################################### # Preparing for the Scan ############################################################################### '''if ID1 == 1: caput ("X11PHS-E:OPT","PGM+ID1") elif ID2 == 1: caput ("X11PHS-E:OPT","PGM+ID2") else: caput ("X11PHS-E:OPT","PGM+ID1+ID2")''' file_prefix = time.strftime("%y%m%d") input_path = "/sls/X11MA/Data1/public/X11MA/temp/"+file_prefix+"/" output_path = input_path x=ENERGIES y = [] x=x.split(',') for i in x: y.append(float(i)) def switchpol(activeID, runtype): global pol_str if activeID == 1: caput(OTF_OFF1,OFFSET1) caput(OTF_OFF2,OFFSET2-40) #detune ID2 if runtype in ["+/-", "+"]: pol_str = "circ +" elif runtype in ["LH/LV", "LH"]: pol_str = "Lin. Horizontal" elif activeID == 2: caput(OTF_OFF1,OFFSET1-40) #detune ID1 caput(OTF_OFF2,OFFSET2) if runtype in ["+/-", "-"]: pol_str = "circ -" elif runtype in ["LH/LV", "LV"]: pol_str = "Lin. Vertical" else: raise Exception("Invalid parameter") LPP_DELAY_NEW = 0 LPP_FINE_DELAY_NEW = 0 def Convert_LPP_Delays(NewDelay): global LPP_DELAY_NEW, LPP_FINE_DELAY_NEW LPP_DELAY_NEW = math.trunc(((NewDelay*100.0+DELAY*BucketSize.read()*100+FINEDELAY)/100.0)/BucketSize.read()) LPP_FINE_DELAY_NEW = math.trunc(NewDelay*100.0+DELAY*BucketSize.read()*100+FINEDELAY-LPP_DELAY_NEW*BucketSize.read()*100.0) pol_str = None polswitch = 1 number_of_scans = 1 if RUNTYPE in ["+/-", "+", "-"]: caput(OTF_MODE1,1) # circ + in ID1 caput(OTF_MODE2,2) # circ - in ID2 wait_channel(OTF_DONE, 1, type = 'i') if RUNTYPE == "+/-": number_of_scans = 2 * ROUNDS else: number_of_scans = ROUNDS elif RUNTYPE in ["LH/LV", "LH", "LV"]: caput(OTF_MODE1,0) caput(OTF_MODE2,0) wait_channel(OTF_DONE, 1, type = 'i') caput(OTF_ALPHA1, 0.0) # LH in ID1 caput(OTF_ALPHA2, 90.0) # LV in ID2 wait_channel(OTF_DONE, 1, type = 'i') if RUNTYPE == "LH/LV": number_of_scans = 2 * ROUNDS else: number_of_scans = ROUNDS else: raise Exception("Invalid run type: " + RUNTYPE) if RUNTYPE in ["-", "LV"]: switchpol(2, RUNTYPE) # tune ID2 --> polarization: C- or LV polswitch = 0 elif RUNTYPE in ["+/-", "+", "LH/LV", "LH"]: switchpol(1, RUNTYPE) # tune ID1 --> polarization: C+ or LH time.sleep(1.0) caput (TIME_DELAY_SET, START) caput (TIME_DELAY_START, 1) if TIME_DELAY_SET == 0: caput (LPP_DELAY, DELAY) caput (LPP_FINE_DELAY, FINE_DELAY) else: Convert_LPP_Delays(caget(TIME_DELAY_SET, type='d')) caput (LPP_DELAY, LPP_DELAY_NEW) caput (LPP_FINE_DELAY, LPP_FINE_DELAY_NEW) wait_channel(TIME_DELAY_COMPLETE, 'SCAN Complete', timeout= abs((caget(TIME_DELAY_SET, type='d')-caget(TIME_DELAY_VAL, type='d')))*10) def Timing(Step_size): New_Time = caget(TIME_DELAY_VAL, type='d') + Step_Size caput (TIME_DELAY_SET, New_Time) caput (TIME_DELAY_START, 1) Convert_LPP_Delays(New_Time) caput (LPP_DELAY, LPP_DELAY_NEW) caput (LPP_FINE_DELAY, LPP_FINE_DELAY_NEW) wait_channel(TIME_DELAY_COMPLETE, 'SCAN Complete', timeout= abs((caget(TIME_DELAY_SET, type='d')-caget(TIME_DELAY_VAL, type='d')))*10) ############################################################################### # Scanning ############################################################################### '''task = None running = False''' MCP1 = [] MCP2 = [] sep = "\t" line_sep = "\r\n" '''def _startPlot(type): global running sep = "\t" line_sep = "\r\n" print "Starting plot: type " + str(type) running = True p = plot(None,name="Energy")[0] s = p.getSeries(0) cur = 0 time.sleep(3.0) while running: try: if otf_start.read() == 0: break''' Delay_Time = START '''for Delay_Time in xrange(START, END+0.1, STEP): e = energy.read() if (MCP_1 == 1 or MCP_2 == 1): aK1=keithley_1a.read() aP=FillingPattern.read() aPr = [] for z in aP: aPr.append(z) SumaP=0 CamaP=0 for z in xrange(0,480): SumaP = SumaP + aPr[z] for z in xrange(460,470): CamaP = CamaP + aPr[z] aK1norm = aK1 * CamaP / SumaP if MCP_1 == 1: a1=MCPArray1.read() a1r = [] for i in a1: a1r.append(i) Sum1a1 = 0 Cam1a1 = 0 Sum2a1 = 0 Cam2a1 = 0 MCP1.append(line_sep) MCP1.append(e) for i in xrange(0,480): Sum1a1 = Sum1a1 + a1r[i] for i in xrange(Cam_start,Cam_end+1): Cam1a1 = Cam1a1 + a1r[i] MCP1.append(Sum1a1) MCP1.append(Cam1a1) MCP1.append(aK1norm) if NrCounters.read() > 460: for i in xrange(480,960): Sum2a1 = Sum2a1 + a1r[i] for i in xrange(Cam_start+480,Cam_end+481): Cam2a1 = Cam2a1 + a1r[i] MCP1.append(Sum2a1) MCP1.append(Cam2a1) if Save_array_data == 1: MCP1.append(a1r) if MCP_2 == 1: a2=MCPArray2.read() a2r = [] for j in a2: a2r.append(j) Sum1a2 = 0 Cam1a2 = 0 Sum2a2 = 0 Cam2a2 = 0 MCP2.append(line_sep) MCP2.append(e) for j in xrange(0,480): Sum1a2 = Sum1a2 + a2r[j] for j in xrange(Cam_start,Cam_end+1): Cam1a2 = Cam1a2 + a2r[j] MCP2.append(Sum1a2) MCP2.append(Cam1a2) MCP2.append(aK1norm) if NrCounters.read() > 460: for j in xrange(480,960): Sum2a2 = Sum2a2 + a2r[j] for j in xrange(Cam_start+480,Cam_end+481): Cam2a2 = Cam2a2 + a2r[j] MCP2.append(Sum2a2) MCP2.append(Cam2a2) if Save_array_data == 1: MCP2.append(a2r) time.sleep(SamplingTime.read()*0.001) if (abs(e-cur)) > 0.1: v = abs((keithley_2a.read() / ((keithley_1a if (type==1) else keithley_3a).read() ))) s.appendData(e,v) cur = e # time.sleep(0.2) except: pass print "Done Plotting" if MCP_1 == 1: output_file_MCP1 = output_path+"MCP_1_"+file_prefix+"_" + suffix + ".dat" MCP1out = open(output_file_MCP1, "a+") sMCP1 = sep.join(str(x) for x in MCP1) + line_sep # MCP1.write("%s" + sep % i) sMCP1out = sep + "rbkenergy" + sep + "Sum1" + sep + "Cam1" + sep + "Keithley1_norm" if NrCounters.read() > 460: sMCP1out = sMCP1out + sep + "Sum2" + sep + "Cam2" if Save_array_data == 1: sMCP1out = sMCP1out + sep + "Array" sMCP1out = sMCP1out + line_sep + sMCP1 MCP1out.write(sMCP1out) MCP1out.close() print "Saved MCP signal data" if MCP_2 == 1: output_file_MCP2 = output_path+"MCP_2_"+file_prefix+"_" + suffix + ".dat" MCP2out = open(output_file_MCP2, "a+") sMCP2 = sep.join(str(y) for y in MCP2) + line_sep sMCP2out = sep + "rbkenergy" + sep + "Sum1" + sep + "Cam1" + sep + "Keithley1_norm" if NrCounters.read() > 460: sMCP2out = sMCP2out + sep + "Sum2" + sep + "Cam2" if Save_array_data == 1: sMCP2out = sMCP2out + sep + "Array" sMCP2out = sMCP2out + line_sep + sMCP2 MCP2out.write(sMCP2out) MCP2out.close() print "Saved MCP I_0 data" ''' ''' def startPlot(type = 1): global task task = fork((_startPlot,(type,)),) def stopPlot(): global task, running running = False ret = join(task) ''' print "Finished Delay scan" print("Success")