parentdir = os.path.dirname(currentdir)
sys.path.insert(0,parentdir)
import a
import numpy as np
numBoots = 200
#O0123_1 = a.readcor("/hiskp2/oehm/correlators/momentum3/A80/O0123/sigma1/O0123_sigma1",48)
#O0123_2 = a.readcor("/hiskp2/oehm/correlators/momentum3/A80/O0123/sigma2/O0123_sigma2",48)
#O0123_3 = a.readcor("/hiskp2/oehm/correlators/momentum3/A80/O0123/sigma3/O0123_sigma3",48)
#O0123_4 = a.readcor("/hiskp2/oehm/correlators/momentum3/A80/O0123/sigma4/O0123_sigma4",48)
#O0123_5 = a.readcor("/hiskp2/oehm/correlators/momentum3/A80/O0123/sigma5/O0123_sigma5",48)
O0123_5t8 = a.readcor("/hiskp2/oehm/correlators/momentum3/A80/O0123/sigma5-8/O0123_sigma5-8",48)
#O0123_1_bstraps = a.boot(O0123_1,numBoots)
#O0123_2_bstraps = a.boot(O0123_2,numBoots)
#O0123_3_bstraps = a.boot(O0123_3,numBoots)
#O0123_4_bstraps = a.boot(O0123_4,numBoots)
#O0123_5_bstraps = a.boot(O0123_5,numBoots)
O0123_5t8_bstraps = a.boot(O0123_5t8,numBoots)
#a.nsettofile2(O0123_1_bstraps[0,:],np.std(O0123_1_bstraps,axis=0),"/hiskp2/oehm/correlators/momentum3/A80/O0123/analysis/O0123_1")
#a.nsettofile2(O0123_2_bstraps[0,:],np.std(O0123_2_bstraps,axis=0),"/hiskp2/oehm/correlators/momentum3/A80/O0123/analysis/O0123_2")
#a.nsettofile2(O0123_3_bstraps[0,:],np.std(O0123_3_bstraps,axis=0),"/hiskp2/oehm/correlators/momentum3/A80/O0123/analysis/O0123_3")
#a.nsettofile2(O0123_4_bstraps[0,:],np.std(O0123_4_bstraps,axis=0),"/hiskp2/oehm/correlators/momentum3/A80/O0123/analysis/O0123_4")
#a.nsettofile2(O0123_5_bstraps[0,:],np.std(O0123_5_bstraps,axis=0),"/hiskp2/oehm/correlators/momentum3/A80/O0123/analysis/O0123_5")
a.nsettofile2(O0123_5t8_bstraps[0,:],np.std(O0123_5t8_bstraps,axis=0),"/hiskp2/oehm/correlators/momentum3/A80/O0123/analysis/O0123_5t8")
if correlator == "corr2": for t in range(24,48): x_bstraps[:,t] = -x_bstraps[:,t] _effmass, effmass_error, effmass_nans = a.effmass(x_bstraps) _effmass2, effmass_error2, effmass_nans2 = a.effmass2(x_bstraps) _effmasswerr = np.zeros((2,_effmass.shape[0])) _effmasswerr[0] = _effmass _effmasswerr[1] = effmass_error np.save(os.path.join(destdir,correlator+"_effmass_50.npy"),_effmasswerr) a.nsettofile2(x_bstraps[0],np.std(x_bstraps,axis=0),os.path.join(destdir,correlator+"_mean_50")) a.nsettofile2(_effmass,effmass_error,os.path.join(destdir,correlator+"_effmass_50")) # plot corr log path = os.path.join(destdir,correlator+"_log_50.pdf") a.ensure_dir(path) plotpdf = PdfPages(path) cut = 0.158655254 orig = np.log(x_bstraps[0,:]) up = np.percentile(np.log(x_bstraps),(1-cut)*100,axis=0)-orig
errors = a.autocorr_bstrap(data,int(taskr.rparam("start")),int(taskr.rparam("end")),numBstraps)
a.plotautocorr_bstrap(errors,os.path.join(destdir,"autocorr_bstrap.pdf"))
elif task == "mwerr":
data = a.readcor(runr.getopcf(op,sigma),latT)
if runr.istrue("wpe"):
data = a.sym_wope(data)
data = data[:,1:len(data[0])-1]
else:
data = a.sym(data)
databstraps = a.boot(data,numBstraps)
a.nsettofile2(databstraps[0],np.std(databstraps,axis=0),os.path.join(destdir,op+"_mean_werr"))
elif task == "fitplatval2":
doasym = 0
if taskr.istrue("asym"):
doasym = 1
wpe = 0
if runr.istrue("wpe"):
wpe = 1
print "wpe: "+str(wpe)
a.fitplatval2(runr.getopcf(op,sigma),destdir,latT,numBstraps,doasym,wpe)
elif task == "plotbstraphists":
doasym = 0
if taskr.istrue("asym"):
