# def get_vars(rs, n):
# xF_over_pT = 2 * np.sinh(n) / rs
# return xF_over_pT
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['sivers'] = PDF1()
conf['dsivers'] = PDF1('deriv')
conf['ffpi'] = FF0('pi')
conf['ffk'] = FF0('k')
rs = 200.
tar = 'p'
pT = 2.
xF = 0.3
C_F = 4.0 / 3.0
N_C = 3.0
def test():
AN = get_AN(xF, pT, rs, nx=10)
print AN
test()
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
print get_AN(xF, pT, rs, tar, had)
print get_ANfrag(xF, pT, rs, tar, had)
print get_ANQS(xF, pT, rs, tar, had)