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')
print get_FX(x, z, Q2, pT, 'd', 'pi+')
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)
def get_AN(xF, pT, rs, nx=10):
return get_SFP(xF, pT, rs, nx=10) + get_SGP(xF, pT, rs, nx=10)
# 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
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
xF2 = xF * xF
xT2 = xT * xT
def test():