'max': 2
}
conf['lattice tabs'] = READER(conf).load_data_sets('lattice')
#####################################################
conf['path2CJ'] = '../../external/CJLIB'
conf['path2LSS'] = '../../external/LSSLIB'
conf['path2DSS'] = '../../external/DSSLIB'
# setup for inclusive dis
conf['alphaSmode'] = 'backward'
conf['Q20'] = 1
conf['order'] = 'LO'
conf['aux'] = AUX()
conf['alphaS'] = ALPHAS(conf)
conf['pdf-NLO'] = CJ(conf)
conf['dis stfuncs'] = DIS_STFUNCS(conf)
# setup tmd sidis
conf['order'] = 'LO'
conf['_pdf'] = CJ(conf)
conf['_ppdf'] = LSS(conf)
conf['_ff'] = DSS(conf)
conf['pdf'] = PDF(conf)
conf['ppdf'] = PPDF(conf)
conf['ff'] = FF(conf)
conf['transversity'] = TRANSVERSITY(conf)
conf['sivers'] = SIVERS(conf)
conf['boermulders'] = BOERMULDERS(conf)
conf['pretzelosity'] = PRETZELOSITY(conf)
conf['wormgearg'] = WORMGEARG(conf)
return res
def tmd_int(b, kt, x, Q, iq, order):
return f_tilde(b, x, Q, iq, order) * b / (2.0 * np.pi) * bessel0(b * kt)
def tmd(kt, x, Q, iq, order):
res, err = integrate.quad(tmd_int, 0, np.inf, args = (kt, x, Q, iq, order))
return res
if __name__=='__main__':
conf={}
conf['path2CJ']='./'
conf['order']='NLO'
cjNLO=CJ(conf)
print cjNLO.get_f(0.1, 1.62)
sys.exit()
if target == 'n':
C = self.p2n(C)
return C
if __name__ == '__main__':
conf['path2CJ'] = '../external/CJLIB'
conf['path2LSS'] = '../external/LSSLIB'
conf['path2DSS'] = '../external/DSSLIB'
conf['path2CT10'] = '../external/PDF'
conf['order'] = 'LO'
conf['aux'] = AUX()
# conf['_pdf']=CT10()
conf['_pdf'] = CJ()
conf['_ppdf'] = LSS()
conf['_ff'] = DSS()
conf['hadron'] = 'pi'
conf['pdf'] = PDF()
conf['ppdf'] = PPDF()
conf['ff'] = FF()
conf['transversity'] = TRANSVERSITY()
conf['sivers'] = SIVERS()
conf['boermulders'] = BOERMULDERS()
conf['pretzelosity'] = PRETZELOSITY()
conf['wormgearg'] = WORMGEARG()
conf['wormgearh'] = WORMGEARH()
conf['collins'] = COLLINS()
conf['dcollinsdz'] = DCOLLINSDZ()
import NLO.Pg.fchn2C,NLO.Ppp.fchn2C
import NLO.Pg.fchn3C,NLO.Ppp.fchn3C
import NLO.Pg.fchn4C,NLO.Ppp.fchn4C
import NLO.Pg.fchn5C,NLO.Ppp.fchn5C
import NLO.Pg.fchn6C,NLO.Ppp.fchn6C
from external.CJLIB.CJ import CJ
from external.DSSLIB.DSS import DSS
import vegas
from numba import jit
conf['order']='NLO'
conf['Q20'] = 1.0
conf['aux']=aux.AUX()
aS=alphaS.ALPHAS()
pdf=CJ({'path2CJ':'external/CJLIB','order':'NLO'})
ff=DSS({'path2DSS':'external/DSSLIB','order':'NLO'})
zero=1e-10
alfa=1/137.036
M=0.93891897
eU= 2./3
eD=-1./3
eq3=np.array([0,eU,-eU,eD,-eD,eD,-eD, 0, 0, 0, 0])
eq4=np.array([0,eU,-eU,eD,-eD,eD,-eD,eU,-eU, 0, 0])
eq5=np.array([0,eU,-eU,eD,-eD,eD,-eD,eU,-eU,eD,-eD])
# g, u, ub, d, db, s, sb, c, cb, b, bb
qqp=np.ones((10,10))
for i in range(10): qqp[i,i]=0
for i in range(10): qqp[i,0]=0
for i in range(10): qqp[0,i]=0
# plt.show()
if __name__ == '__main__':
import qcdlib.aux
import qcdlib.alphaS
#conf={}
#cwd = os.getcwd()
conf['path2CJ'] = '../../external/CJLIB'
#conf['path2CT10'] ='../../external/PDF'
conf['path2LSS'] = '../../external/LSSLIB'
conf['path2DSS'] = '../../external/DSSLIB'
conf['order'] = 'NLO'
conf['aux'] = AUX()
conf['_pdf'] = CJ(conf)
#conf['_pdf'] =CT10(conf)
conf['_ppdf'] = LSS(conf)
conf['_ff'] = DSS(conf)
conf['gk'] = GK()
conf['pdf'] = PDF()
conf['ppdf'] = PPDF()
conf['ff'] = FF()
conf['Q20'] = conf['gk'].Q0**2
conf['alphaSmode'] = 'backward'
conf['alphaS'] = qcdlib.alphaS.ALPHAS()
stfuncs = STFUNCS()
x = 0.25
z = 0.5
+ self.Cg(z, lambda y: self.glue(self.get_PDFs(y)), x)
def get_F2(self, x, Q2, hadron, method='gauss', n=20):
if (x, Q2, hadron) not in self.storage:
self.hadron = hadron
self.Nf = conf['alphaS'].get_Nf(Q2)
alpi = conf['alphaS'].get_alphaS(Q2) / (2. * np.pi)
self.get_PDFs = lambda y: conf['pdf-NLO'].get_f(y, Q2)
LO = self.qplus(self.get_PDFs(x))
#integrand=lambda z:self.integrand(x,z,Q2)
# NLO=self.integrator(integrand,x,1,method,n=n)
self.storage[(x, Q2, hadron)] = x * (LO) # +alpi*NLO)
return self.storage[(x, Q2, hadron)]
if __name__ == '__main__':
conf['alphaSmode'] = 'backward'
conf['order'] = 'NLO'
conf['Q20'] = 1
conf['aux'] = AUX()
conf['alphaS'] = ALPHAS()
conf['path2CJ'] = '../../external/CJLIB'
conf['pdf-NLO'] = CJ()
stfuncs = STFUNCS()
x = 0.5
Q2 = 1000.
print stfuncs.get_F2(x, Q2, 'p', method='gauss')
print stfuncs.get_F2(x, Q2, 'n', method='gauss')