sr0[i] = np.empty(numSys)
sr0.fill(e[1][0])
# read moments for all ensembles and resample systematic deviations
j = 0
for e2 in e[2]:
print "preparing "+e2[0]
xnmedian = np.load("/hiskp2/oehm/averx/analysis/"+e2[0]+"/xnmedian.npy")
xnhist = np.load("/hiskp2/oehm/averx/analysis/"+e2[0]+"/xnhist0.npy")
udmedian = np.load("/hiskp2/oehm/averx/analysis/"+e2[0]+"/udmedian.npy")
udhist = np.load("/hiskp2/oehm/averx/analysis/"+e2[0]+"/udhist0.npy")
ufmedian = np.load("/hiskp2/oehm/averx/analysis/"+e2[0]+"/ufmedian.npy")
ufhist = np.load("/hiskp2/oehm/averx/analysis/"+e2[0]+"/ufhist0.npy")
# resample systematic deviations
udhistre, udxnhistre = a.resample_whists2(udhist,np.random.uniform(0,1,numSys))
ufhistre, ufxnhistre = a.resample_whists2(ufhist,np.random.uniform(0,1,numSys))
uddata[i][j,:,0] = udmedian[:,0]
uddata[i][j,:,1] = xnmedian**2
ufdata[i][j,:,0] = ufmedian[:,0]
ufdata[i][j,:,1] = xnmedian**2
suddata[i][j,:,0] = udhistre
suddata[i][j,:,1] = udxnhistre**2
sufdata[i][j,:,0] = ufhistre
sufdata[i][j,:,1] = ufxnhistre**2
cov = np.cov(uddata[i][j].T)
print e[0]+" "+str(r02M2[0])+" "+str(udNx[0])+" "+str(np.std(r02M2))+" "+str(np.std(udNx))+" "+str(r02M2[0]-r02M2low)+" "+str(r02M2up-r02M2[0])+" "+str(udNx[0]-udNxlow)+" "+str(udNxup-udNx[0]) udr02M2 = np.zeros(numBootstraps) ufr02M2 = np.zeros(numBootstraps) #np.random.seed(0) #udhistre, udxnhistre = a.resample_whists(udhist2,xnhist,numSys) #ufhistre, ufxnhistre = a.resample_whists(ufhist2,xnhist,numSys) #a.plotbstraphistarray(destdir,udxnhistre,e[0]+"_udxnhistre.pdf",50) #a.plotbstraphistarray(destdir,ufxnhistre,e[0]+"_ufxnhistre.pdf",50) #a.plotbstraphistarray(destdir,udhistre,e[0]+"_udhistre.pdf",50) #a.plotbstraphistarray(destdir,ufhistre,e[0]+"_ufhistre.pdf",50) #np.random.seed(0) udhistre, udxnhistre = a.resample_whists2(udhist,sysrandnums) ufhistre, ufxnhistre = a.resample_whists2(ufhist,sysrandnums) a.plotbstraphistarray(destdir,udxnhistre,e[0]+"_udxnhistre2.pdf",50) a.plotbstraphistarray(destdir,ufxnhistre,e[0]+"_ufxnhistre2.pdf",50) a.plotbstraphistarray(destdir,udhistre,e[0]+"_udhistre2.pdf",50) a.plotbstraphistarray(destdir,ufhistre,e[0]+"_ufhistre2.pdf",50) udr02M2 = np.zeros(numSys) ufr02M2 = np.zeros(numSys) udNx = np.zeros(numSys) ufNx = np.zeros(numSys) for b in range(0,numSys): udr02M2[b] = (e[2]*udxnhistre[b])**2
