F = conf['aux'].p2n(F)
return _get_FX(x, z, Q2, pT, target, hadron, F, D)
elif target == 'd':
return 0.5 * (get_FX(x, z, Q2, pT, 'p', hadron) +
get_FX(x, z, Q2, pT, 'n', hadron))
if __name__ == '__main__':
from qcdlib.ff1 import FF as FF1
from qcdlib.pdf1 import PDF as PDF1
conf['aux'] = AUX()
conf['Htildepi'] = FF1('pi')
conf['Htildek'] = FF1('k')
conf['transversity'] = PDF1()
x = 0.25
z = 0.5
Q2 = 2.4
pT = 0.3
print get_FX(x, z, Q2, pT, 'p', 'pi+')
print get_FX(x, z, Q2, pT, 'p', 'pi-')
print get_FX(x, z, Q2, pT, 'p', 'pi0')
print get_FX(x, z, Q2, pT, 'n', 'pi+')
print get_FX(x, z, Q2, pT, 'n', 'pi-')
print get_FX(x, z, Q2, pT, 'n', 'pi0')
sig = dblquad(lambda x, z: get_dsigST(x, z, xF, pT, rs, tar, had),
zmin, 1., xmin, lambda x: 1.)[0]
#print had,sig
return sig
if __name__ == '__main__':
from qcdlib.ff0 import FF as FF0
from qcdlib.ff1 import FF as FF1
from qcdlib.pdf0 import PDF as PDF0
from qcdlib.pdf1 import PDF as PDF1
conf['aux'] = AUX()
conf['pdf'] = PDF0()
conf['collinspi'] = FF1('pi')
conf['collinsk'] = FF1('k')
conf['Htildepi'] = FF1('pi')
conf['Htildek'] = FF1('k')
conf['transversity'] = PDF1()
conf['sivers'] = PDF1()
conf['ffpi'] = FF0('pi')
conf['ffk'] = FF0('k')
rs = 200.
tar = 'p'
#had = 'pi+'
had = 'jet'
pT = 3.0
xF = 0.2375
xT = 2. * pT / rs
def get_ANQS(xF, pT, rs, tar, had):
numQS = get_numQS(xF, pT, rs, tar, had)
denom = get_denom(xF, pT, rs, tar, had)
return numQS / denom
if __name__ == '__main__':
from qcdlib.ff0 import FF as FF0
from qcdlib.ff1 import FF as FF1
from qcdlib.pdf0 import PDF as PDF0
from qcdlib.pdf1 import PDF as PDF1
conf['aux'] = AUX()
conf['pdf'] = PDF0()
conf['collinspi'] = FF1('pi')
conf['dcollinspi'] = FF1('pi', 'deriv')
conf['dcollinsk'] = FF1('k', 'deriv')
conf['Htildepi'] = FF1('pi')
conf['transversity'] = PDF1()
conf['sivers'] = PDF1()
conf['dsivers'] = PDF1('deriv')
conf['ffpi'] = FF0('pi')
rs = 7.25
tar = 'p'
had = 'pi+'
pT = 1.0
xF = 0.3
Q2 = 2.0
if __name__ == '__main__':
from qcdlib.interpolator import INTERPOLATOR
from qcdlib.ff0 import FF as FF0
from qcdlib.ff1 import FF as FF1
conf['aux'] = AUX()
conf['cpipff'] = INTERPOLATOR('dsspipNLO_0000')
conf['cpimff'] = INTERPOLATOR('dsspimNLO_0000')
conf['cKpff'] = INTERPOLATOR('dssKpNLO_0000')
conf['cKmff'] = INTERPOLATOR('dssKmNLO_0000')
conf['lam2'] = 0.4
conf['Q02'] = 1.0
conf['ffpi'] = FF0('pi')
conf['ffpi'] = FF0('k')
conf['collinspi'] = FF1('pi')
conf['collinsk'] = FF1('k')
z1 = 0.5
z2 = 0.3
Q2 = 2.5
pT = 1.0
hadron1 = 'pi+'
hadron2 = 'pi-'
#for i in [1,2]:
print 1, ZX(1, z1, z2, Q2, pT, hadron1, hadron2)
print 2, ZX(2, z1, z2, Q2, pT, hadron1, hadron2)
