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