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")
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