ufaverx = np.zeros((numBootstraps,21*N2,2))

for i in range(0,21):
	for j in range(0,N2):
		W = (1-2*abs(xn_result_mean[i][3]-0.5))**2*(1-2*abs(O012_result_mean[j][2]-0.5))**2/(xn_result_std[i][1]**2)/(O012_result_std[j][0]**2)
		for b in range(0,numBootstraps):
			udaverx[b][i*N2+j][0] = -O012_result[j][b][0]/xn_bstraps[b][24]/((2*math.pi/24*sigma/5)**2)
			ufaverx[b][i*N2+j][0] = -O012_result[j][b][0]/xnfunc([xn_result[i][b][0],xn_result[i][b][1]],24)/((2*math.pi/24*sigma/5)**2)
			udaverx[b][i*N2+j][1] = W
			ufaverx[b][i*N2+j][1] = W

#a.tofile(udaverx,path = "/hiskp2/oehm/correlators/momentum/A40.24/O012/analysis/sigma"+str(sigma)+"/udaverx2")

cut = 0.158655254

udmedian = a.histpercentile(udaverx,0.5)
udlow = a.histpercentile(udaverx,cut)
udup = a.histpercentile(udaverx,1-cut)

file = open("/hiskp2/oehm/correlators/momentum/A40.24/O012/analysis/sigma"+str(sigma)+"/result","w")
file.write("ud: "+str(udmedian[0])+"+-"+str(np.std(udmedian))+"-"+str(udmedian[0]-udlow[0])+"+"+str(udup[0]-udmedian[0])+"\n")

ufmedian = a.histpercentile(ufaverx,0.5)
uflow = a.histpercentile(ufaverx,cut)
ufup = a.histpercentile(ufaverx,1-cut)

file.write("uf: "+str(ufmedian[0])+"+-"+str(np.std(ufmedian))+"-"+str(ufmedian[0]-uflow[0])+"+"+str(ufup[0]-ufmedian[0])+"\n")

file.close()

file = open("/hiskp2/oehm/correlators/momentum/A40.24/O012/analysis/sigma"+str(sigma)+"/result2","w")
예제 #2
0
				O44_result[N2][b][2] = pvalue[b]
			N2 += 1

O44_result_mean = np.mean(O44_result,axis=1)
O44_result_std = np.std(O44_result,axis=1)

xn_hist = np.zeros((numBootstraps,21,2))

for i in range(0,21):
	W = (1-2*abs(xn_result_mean[i][3]-0.5))**2/(xn_result_std[i][1]**2)
	for b in range(0,numBootstraps):
		xn_hist[b][i][0] = xn_result[i][b][1]
		xn_hist[b][i][1] = W

cut = 0.158655254
xnmedian = a.histpercentile(xn_hist,0.5)
xnlow = a.histpercentile(xn_hist,cut)
xnup = a.histpercentile(xn_hist,1-cut)

udaverx = np.zeros((numBootstraps,(T/2+1-4)*N2,3))
ufaverx = np.zeros((numBootstraps,(T/2+1-4)*N2,3))

for i in range(0,T/2+1-4):
	for j in range(0,N2):
		W = (1-2*abs(xn_result_mean[i][3]-0.5))**2*(1-2*abs(O44_result_mean[j][2]-0.5))**2/(xn_result_std[i][1]**2)/(O44_result_std[j][0]**2)
		for b in range(0,numBootstraps):
			udaverx[b][i*N2+j][0] = -O44_result[j][b][0]/xn_result[i][b][1]/xn_bstraps[b][T/2]
			ufaverx[b][i*N2+j][0] = -O44_result[j][b][0]/xn_result[i][b][1]/xnfunc([xn_result[i][b][0],xn_result[i][b][1]],T/2)
			udaverx[b][i*N2+j][1] = xn_result[i][b][1]
			ufaverx[b][i*N2+j][1] = xn_result[i][b][1]
			udaverx[b][i*N2+j][2] = W