#!/usr/bin/env python
import sys
import os
import numpy as np
from tools.tools import load_config
from qcdlib.aux import AUX
from tools.config import conf
def get_h1(x, Q2, had='p'):
return conf['transversity'].get_C(x, Q2)
if __name__ == '__main__':
from qcdlib.pdf1 import PDF
conf['aux'] = AUX()
conf['transversity'] = PDF()
x = 0.3
Q2 = 2.0
print get_h1(x, Q2, 'p')
F = conf['aux'].p2n(F)
w_tar = conf['aux'].p2n(w_tar)
return _get_FUT(x, z, Q2, pT, tar, had, F, D, w_tar, w_had)
elif tar == 'd':
return 0.5 * (get_FUT(x, z, Q2, pT, 'p', had) +
get_FUT(x, z, Q2, pT, 'n', had))
if __name__ == '__main__':
from qcdlib.pdf1 import PDF
from qcdlib.ff0 import FF
conf['aux'] = AUX()
conf['sivers'] = PDF()
conf['ffpi'] = FF('pi')
conf['ffk'] = FF('k')
x = 0.186
z = 0.379
Q2 = 4.321
mu2 = Q2
pT = 0.339
tar = 'p'
had = 'pi+'
print get_FUT(x, z, Q2, pT, 'p', 'pi+')
print get_FUT(x, z, Q2, pT, 'p', 'pi-')
print get_FUT(x, z, Q2, pT, 'p', 'pi0')
F = conf['aux'].p2n(F)
w_tar = conf['aux'].p2n(w_tar)
return _get_FUU(x, z, Q2, pT, tar, had, F, D, w_tar, w_had)
elif tar == 'd':
return 0.5 * (get_FUU(x, z, Q2, pT, 'p', had) +
get_FUU(x, z, Q2, pT, 'n', had))
if __name__ == '__main__':
from qcdlib.pdf1 import PDF
from qcdlib.ff1 import FF
conf['aux'] = AUX()
conf['boermulders'] = PDF()
conf['collinspi'] = FF('pi')
conf['collinsk'] = FF('k')
x = 0.25
z = 0.5
Q2 = 2.4
mu2 = Q2
pT = 0.3
tar = 'p'
had = 'pi+'
print get_FUU(x, z, Q2, pT, 'p', 'pi+')
print get_FUU(x, z, Q2, pT, 'p', 'pi-')
print get_FUU(x, z, Q2, pT, 'p', 'pi0')