import sys
import numpy as np
import os
import math
import os,inspect
currentdir = os.path.dirname(os.path.abspath(inspect.getfile(inspect.currentframe())))
parentdir = os.path.dirname(currentdir)
sys.path.insert(0,parentdir)
import a
sigma=5
print "sigma: "+str(sigma)
averx2val = lambda x,y,z: -x/y*2*(np.sqrt(3))**3/((2*math.pi/24*(sigma)/5)**3)
a.calcplatval("/hiskp2/oehm/correlators/momentum3/A80/O0123/sigma5-8/O0123_sigma5-8","/hiskp2/oehm/correlators/momentum3/A80/x/sigma"+str(sigma)+"/x"+str(sigma),"/hiskp2/oehm/correlators/momentum3/A80/O0123/analysis/sigma5-8/",48,1500,averx2val)
a.fitplatval3(runr.getopcf(op,sigma),runr.getopcf("x",sigma),runr.getopcf("x0",sigma),destdir,latT,numBstraps,valfunc,doasym = doasym)
elif task == "calcplatval":
sigmaterm = runr.getsigma(sigma)
if op == "O44":
valfunc = lambda x,y,z: -x/y*z/(z**2+((2*math.pi/24*sigmaterm)**2)/3)
elif op == "O012":
valfunc = lambda x,y,z: -x/y/((2*math.pi/24*sigmaterm)**2)
elif op == "O0ii":
valfunc = lambda x,y,z: -x/y*2/((2*math.pi/24*sigmaterm)**2)
print runr.getopcf(op,sigma)+","+runr.getopcf("x",sigma)+","+destdir+","+str(latT)+","+str(numBstraps)+","+str(doasym)
a.calcplatval(runr.getopcf(op,sigma),runr.getopcf("x",sigma),destdir,latT,numBstraps,valfunc,wpe)
elif task == "autocorr_bstrap":
data = a.readcor(runr.getopcf(op,sigma),latT)
if runr.istrue("wpe"):
data = a.sym_wope(data)
else:
data = a.sym(data)
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"):
