'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')