def TDEVTest():
# dataFile=unicode('C:\Users\Levit\Experiment Data\德令哈测试\\20171226\零基线实验\\20171227015305-tdc2_5_filterN_coindence_filtered_250-350s_residual-1-0.1s-ps.txt','utf8')
dataFile = unicode(
'C:\Users\Levit\Experiment Data\双站数据\\20180121\\result\\synCoincidence-124-216--17-1-Coin-紫台WGS84-atm-factor-haiji_laser改正_filtered_residual-2-0.1s-ps.txt',
'utf8')
dataList=fileToList.fileToList(dataFile)
saveFile = dataFile[:-4] + '_TDEV.txt'
tdev=TDEV(dataList,0.1)
fileToList.listToFileFloat(tdev, saveFile)
def twoLocalTDCDataProcess(dataFile1,dataFile2,fineTimeFile1,fineTimeFile2,order,tau):
channel1=int(fineTimeFile1[-8])
channel2=int(fineTimeFile2[-8])
saveFile1 = dataFile1[:-4] + '_channel%s.txt'%channel1
dataList1 = TDCDataConvert.TDCDataParse(dataFile1, fineTimeFile1, 8, str(40+channel1-1))
saveFile2 = dataFile2[:-4] + '_channel%s.txt' % channel2
dataList2 = TDCDataConvert.TDCDataParse(dataFile2, fineTimeFile2, 8, str(40+channel2-1))
fileToList.listToFile(dataList1, saveFile1)
fileToList.listToFile(dataList2, saveFile2)
num1=len(dataList1)
num2=len(dataList2)
if num1!=num2:
print 'data lenght is not equal!'
return 0
xa=[]
xb=[]
ya=[]
for i in range(num1):
xa.append(dataList1[i][0])
xb.append([dataList1[i][0]])
ya.append(dataList1[i][0]-dataList2[i][0])
xa, ya, fitList, residual =fitting.polyFitSegment(xa, ya, order, 10000)
del xa[:]
del fitList[:]
xOneSec, residualOneSec = filter.normalByTime(xb, residual, 1000000000000)
fileToList.listToFileLong(residualOneSec, saveFile1[:-4] + '-%s_residual-%s-1s-ps.txt'%(channel2,order))
xTau,residualTau=filter.normalByTime(xb, residual, tau)
residualFile=saveFile1[:-4] + '-%s_residual-%s-%.3fs-ps.txt' % (channel2, order,tau/1000000000000.0)
fileToList.listToFileLong(residualTau, residualFile)
fig = plt.figure()
ax1 = fig.add_subplot(111)
for i,item in enumerate(xOneSec):
xOneSec[i]=item/1000000000000.0
residualOneSec[i][0]=residualOneSec[i][0]*1000000000000
ax1.plot(xOneSec,residualOneSec, color='g', linestyle='-', marker='*')
ax1.xaxis.grid(True, which='major') # x坐标轴的网格使用主刻度
ax1.yaxis.grid(True, which='major') # y坐标轴的网格使用次刻度show()
plt.show()
tdevFile = residualFile[:-4] + '_TDEV.txt'
tdev = varianceStatistics.TDEV(residualTau, tau/1000000000000.0)
fileToList.listToFileFloat(tdev, tdevFile)
print 'TDEV calculation finished!'
fig = plt.figure()
dataPlot.logPlotAx(tdev, fig, 'r','--','s', '11.25 tdev')
dataFile1 = unicode(
'C:\Users\Levit\Experiment Data\Rakon晶振测试数据\两TDC测试\\20171114200130-tdc2-13-2k-500s-3_residual-1-0.01s-ps_TDEV.txt',
'utf8')
List1 = fileToList.fileToList(dataFile1)
dataPlot.logPlotAx(List1, fig, 'g', '--', 'o', '11.14')
plt.show()
def fitting_segment(startTime,endTime,coindenceFile):
coindenceList=fileToList.fileToList(coindenceFile)
second=1000000000000.0
tau=100000000000
order=2
num = len(coindenceList)
xa = []
xb = []
ya = []
for i in range(num):
if coindenceList[i][0]/second>=startTime and coindenceList[i][0]/second<endTime:
xa.append(coindenceList[i][0])
ya.append(coindenceList[i][2])
# xa, ya, coindenceList, fitList, residual = filter.fitFilter(coindenceList, 2000 / 1000000000000.0, 2, 1)
xa, ya, fitList, residual = fitting.polyFitSegment(xa, ya, order, 10000)
for i in range(len(xa)):
xb.append([xa[i]])
xOneSec, residualOneSec = filter.normalByTime(xb, residual, second)
del xa[:]
del fitList[:]
del coindenceList[:]
fileToList.listToFileLong(residualOneSec, coindenceFile[:-4] + '_%s-%s-1s-ps.txt' % (startTime, endTime))
xTau, residualTau = filter.normalByTime(xb, residual, tau)
residualFile = coindenceFile[:-4] + '_%s-%s-%.3fs-ps.txt' % (startTime, endTime, tau / 1000000000000.0)
fileToList.listToFileLong(residualTau, residualFile)
fig = plt.figure()
ax1 = fig.add_subplot(111)
for i, item in enumerate(xOneSec):
xOneSec[i] = item / 1000000000000.0
residualOneSec[i][0] = residualOneSec[i][0] * 1000000000000
ax1.plot(xOneSec, residualOneSec, color='g', linestyle='-', marker='*')
ax1.xaxis.grid(True, which='major') # x坐标轴的网格使用主刻度
ax1.yaxis.grid(True, which='major') # y坐标轴的网格使用次刻度show()
plt.show()
tdevFile = residualFile[:-4] + '_TDEV.txt'
tdev = varianceStatistics.TDEV(residualTau, tau / 1000000000000.0)
fileToList.listToFileFloat(tdev, tdevFile)
print 'TDEV calculation finished!'
tdevName='TDEV %s-%s second'%(startTime,endTime)
fig = plt.figure()
dataPlot.logPlotAx(tdev, fig, 'r', '--', 's', tdevName)
plt.show()
def MDEVTest():
dataFile=unicode('C:\Users\Levit\Experiment Data\Rakon晶振测试数据\两TDC测试\\20171121\\20171121230538-tdc2-13-2k-500s-dating-1_residual-2-0.01s-ps.txt','utf8')
dataList=fileToList.fileToList(dataFile)
saveFile = dataFile[:-4] + '_MDEV.txt'
mdev=MDEV(dataList,0.01)
fileToList.listToFileFloat(mdev, saveFile)
def ADEVTest():
dataFile=unicode('C:\Users\Levit\Experiment Data\Rakon晶振测试数据\两TDC测试\\20171114200130-tdc2-13-2k-500s-3_residual-1-0.01s-ps.txt','utf8')
saveFile=dataFile[:-4]+'_ADEV.txt'
dataList=fileToList.fileToList(dataFile)
adev=ADEV(dataList,0.01)
fileToList.listToFileFloat(adev,saveFile)
def twoLightTDCDataProcess(channel1,channel2,delay,window,dataFile1,dataFile2,fineTimeFile1,fineTimeFile2,order,tau,tdevName,tdevComp,tdevCompFile):
saveFile1 = dataFile1[:-4] + '_channel_%s.txt'%channel1
dataList1 = TDCDataConvert.TDCDataParse(dataFile1, fineTimeFile1, 8, str(40+channel1-1))
fileToList.listToFile(dataList1, saveFile1)
dataList1_filter=filter.freqFilter(dataList1,10000200,6,300000)
dataList1_filterN=filter.reflectNoiseFilter(dataList1_filter,1000000,0)
fileToList.listToFile(dataList1_filterN, saveFile1[:-4]+'_filterN.txt')
saveFile2 = dataFile2[:-4] + '_channel_%s.txt' % channel2
dataList2 = TDCDataConvert.TDCDataParse(dataFile2, fineTimeFile2, 8, str(40+channel2 - 1))
fileToList.listToFile(dataList2, saveFile2)
dataList2_filter = filter.freqFilter(dataList2, 10000200, 6, 300000)
dataList2_filterN = filter.reflectNoiseFilter(dataList2_filter, 1000000, 0)
TDCTest.countBySec(dataList2)
TDCTest.countBySec(dataList1)
TDCTest.countBySec(dataList2_filterN)
TDCTest.countBySec(dataList1_filterN)
fileToList.listToFile(dataList2_filterN, saveFile2[:-4] + '_filterN.txt')
saveCoinFile = dataFile1[:-4] + '_%s_filterN_coindence.txt'%channel1
coindenceList,averSecCount=TDCTest.coindenceTest(dataList1_filterN, dataList2_filterN, delay, window, saveCoinFile)
tdevName=tdevName+' %s'%averSecCount
del dataList1,dataList1_filter,dataList1_filterN
del dataList2,dataList2_filter,dataList2_filterN
num=len(coindenceList)
xa = []
xb = []
ya = []
for i in range(num):
xa.append(coindenceList[i][0])
ya.append(coindenceList[i][2])
xa, ya, coindenceList, fitList, residual = filter.fitFilter(coindenceList, 2000 / 1000000000000.0, 2, 1)
xa, ya, fitList, residual = fitting.polyFitSegment(xa, ya, order, 10000)
fileToList.listToFile(coindenceList,saveCoinFile[:-4]+'_filtered.txt')
for i in range(len(xa)):
xb.append([xa[i]])
xOneSec, residualOneSec = filter.normalByTime(xb, residual, 1000000000000)
del xa[:]
del fitList[:]
del coindenceList[:]
fileToList.listToFileLong(residualOneSec, saveFile1[:-4] + '-%s_residual-%s-1s-ps.txt' % (channel2, order))
xTau, residualTau = filter.normalByTime(xb, residual, tau)
residualFile = saveFile1[:-4] + '-%s_residual-%s-%.3fs-ps.txt' % (channel2, order, tau / 1000000000000.0)
fileToList.listToFileLong(residualTau, residualFile)
fig = plt.figure()
ax1 = fig.add_subplot(111)
for i, item in enumerate(xOneSec):
xOneSec[i] = item / 1000000000000.0
residualOneSec[i][0] = residualOneSec[i][0] * 1000000000000
ax1.plot(xOneSec, residualOneSec, color='g', linestyle='-', marker='*')
ax1.xaxis.grid(True, which='major') # x坐标轴的网格使用主刻度
ax1.yaxis.grid(True, which='major') # y坐标轴的网格使用次刻度show()
plt.show()
tdevFile = residualFile[:-4] + '_TDEV.txt'
tdev = varianceStatistics.TDEV(residualTau, tau / 1000000000000.0)
fileToList.listToFileFloat(tdev, tdevFile)
print 'TDEV calculation finished!'
fig = plt.figure()
dataPlot.logPlotAx(tdev, fig, 'r', '--', 's', tdevName)
List1 = fileToList.fileToList(tdevCompFile)
dataPlot.logPlotAx(List1, fig, 'g', '--', 'o', tdevComp)
plt.show()