def setup_core(self):
conf['aux'] = qcdlib.aux.AUX()
conf['mellin'] = qcdlib.mellin.MELLIN(npts=4)
conf['alphaS'] = qcdlib.alphaS.ALPHAS()
if 'pdf parametrization' in conf:
if conf['pdf parametrization'] == 0: conf['pdf'] = pdf0.PDF()
if conf['pdf parametrization'] == 1: conf['pdf'] = pdf1.PDF()
if conf['pdf parametrization'] == 2: conf['pdf'] = pdf2.PDF()
if conf['pdf parametrization'] == 3: conf['pdf'] = pdf3.PDF()
if 'pdf' in conf['params']: conf['pdf'] = pdf0.PDF()
if 'pdfpi-' in conf['params']: conf['pdfpi-'] = pdf0.PDF('pi-')
if 'transversity' in conf['params']: conf['transversity'] = pdf1.PDF()
if 'sivers' in conf['params']: conf['sivers'] = pdf1.PDF()
if 'boermulders' in conf['params']: conf['boermulders'] = pdf1.PDF()
if 'ffpi' in conf['params']: conf['ffpi'] = ff0.FF('pi')
if 'ffk' in conf['params']: conf['ffk'] = ff0.FF('k')
if 'collinspi' in conf['params']: conf['collinspi'] = ff1.FF('pi')
if 'collinsk' in conf['params']: conf['collinsk'] = ff1.FF('k')
if 'Htildepi' in conf['params']: conf['Htildepi'] = ff1.FF('pi')
if 'Htildek' in conf['params']: conf['Htildek'] = ff1.FF('k')
if 'transversity+' in conf['params']:
conf['transversity'] = pdf2.PDF('h1') # Transversity
if 'collinspi+' in conf['params']:
conf['collinspi'] = ff2.FF('Col') # Collins
if 'sivers+' in conf['params']:
conf['sivers'] = pdf2.PDF('Siv') # Sivers
def setup_core(self):
conf['aux'] = qcdlib.aux.AUX()
conf['mellin'] = qcdlib.mellin.MELLIN(npts=4)
conf['alphaS'] = qcdlib.alphaS.ALPHAS()
if 'pdf parametrization' in conf:
if conf['pdf parametrization'] == 0: conf['pdf'] = pdf0.PDF()
if conf['pdf parametrization'] == 1: conf['pdf'] = pdf1.PDF()
if conf['pdf parametrization'] == 2: conf['pdf'] = pdf2.PDF()
if conf['pdf parametrization'] == 3: conf['pdf'] = pdf3.PDF()
if 'version' in conf: version = conf['version']
else: version = 0 #--for back compatibility
if 'pdf' in conf['params']: conf['pdf'] = pdf0.PDF()
if 'pdfpi-' in conf['params']: conf['pdfpi-'] = pdf0.PDF('pi-')
if 'transversity' in conf['params']:
if version == 0: conf['transversity'] = pdf1.PDF()
if version == 'JAM20+': conf['transversity'] = pdf2.PDF('h1')
if 'sivers' in conf['params']:
if version == 0:
conf['sivers'] = pdf1.PDF()
conf['dsivers'] = pdf1.PDF('deriv')
if version == 'JAM20+':
conf['sivers'] = pdf2.PDF('Siv')
conf['dsivers'] = pdf2.PDF('Siv', 'deriv')
if 'boermulders' in conf['params']: conf['boermulders'] = pdf1.PDF()
if 'ffpi' in conf['params']: conf['ffpi'] = ff0.FF('pi')
if 'ffk' in conf['params']: conf['ffk'] = ff0.FF('k')
if 'collinspi' in conf['params']:
if version == 0:
conf['collinspi'] = ff1.FF('pi')
conf['dcollinspi'] = ff1.FF('pi', 'deriv')
if version == 'JAM20+':
conf['collinspi'] = ff2.FF('Col', 'pi')
conf['dcollinspi'] = ff2.FF('Col', 'pi', 'deriv')
if 'collinsk' in conf['params']:
if version == 0:
conf['collinsk'] = ff1.FF('k')
conf['dcollinsk'] = ff1.FF('k', 'deriv')
if version == 'JAM20+':
conf['collinsk'] = ff2.FF('Col', 'k')
conf['dcollinsk'] = ff2.FF('Col', 'k', 'deriv')
if 'Htildepi' in conf['params']:
if version == 0:
conf['Htildepi'] = ff1.FF('pi')
if version == 'JAM20+':
conf['Htildepi'] = ff2.FF(
'Col',
'pi') # Htilde (using same splitting functions as Collins)
if 'Htildek' in conf['params']:
if version == 0:
conf['Htildek'] = ff1.FF('k')
if version == 'JAM20+':
conf['Htildek'] = ff2.FF(
'Col',
'k') # Htilde (using same splitting functions as Collins)
for l in L:
print l
return L
if __name__ == '__main__':
from qcdlib.interpolator import INTERPOLATOR
from qcdlib import pdf0, ff0, pdf1, ff1
from qcdlib.aux import AUX
from reader import READER
conf['aux'] = AUX()
conf['pdf'] = pdf0.PDF()
conf['transversity'] = pdf1.PDF()
conf['sivers'] = pdf1.PDF()
conf['ffpi'] = ff0.FF('pi')
conf['ffk'] = ff0.FF('k')
conf['collinspi'] = ff1.FF('pi')
conf['collinsk'] = ff1.FF('k')
conf['datasets'] = {}
conf['datasets']['sidis'] = {}
conf['datasets']['sidis']['xlsx'] = {}
# upol
conf['datasets']['sidis']['xlsx'][
1000] = 'sidis/expdata/1000.xlsx' # | proton | pi+ | M_Hermes | hermes
conf['datasets']['sidis']['xlsx'][
for l in L:
print l
return L
if __name__ == '__main__':
from qcdlib.interpolator import INTERPOLATOR
from qcdlib import pdf0, ff0, pdf1, ff1
from qcdlib.aux import AUX
from reader import READER
conf['aux'] = AUX()
conf['pdf'] = pdf0.PDF()
conf['transversity'] = pdf1.PDF()
conf['sivers'] = pdf1.PDF()
conf['boermulders'] = pdf1.PDF()
conf['ffpi'] = ff0.FF('pi')
conf['ffk'] = ff0.FF('k')
conf['collinspi'] = ff1.FF('pi')
conf['collinsk'] = ff1.FF('k')
conf['datasets'] = {}
conf['datasets']['sidis'] = {}
conf['datasets']['sidis']['xlsx'] = {}
''' For AUUcos2 testing purposes, most entries have been commented out.
Uncomment for full program functionality. '''
'''
# upol
Q2 = conf['aux'].mZ2
xA = np.sqrt(Q2/energy) * np.exp(+y) # A is polarized
xB = np.sqrt(Q2/energy) * np.exp(-y) # B is unpolarized
# build structure function
return get_FUTWZ(xA,xB,qT,hadronA,hadronB,boson,TransversePolarizationA,TransversePolarizationB)
if __name__ == '__main__':
from qcdlib import pdf0
from qcdlib import pdf1
conf['aux']= AUX()
conf['sivers']=pdf1.PDF()
conf['pdf']=pdf0.PDF('p')
conf['pdfpi-']=pdf0.PDF('pi-')
xA = 0.25
xB = 0.5
Q2 = 16
qT = 0.3
hadronA = 'p'
hadronB = 'p'
TransversePolarizationA = True
TransversePolarizationB = False
print(get_FUTWZ(xA,xB,qT,hadronA,hadronB,'W+',TransversePolarizationA,TransversePolarizationB))
return L
elif verb == 1:
for l in L:
print(l)
return L
if __name__ == '__main__':
from qcdlib.interpolator import INTERPOLATOR
from qcdlib import pdf1
from qcdlib.aux import AUX
from reader import READER
conf['aux'] = AUX()
conf['transversity'] = pdf1.PDF()
conf['datasets']={}
conf['datasets']['SB']={}
conf['datasets']['SB']['xlsx']={}
# SB
conf['datasets']['SB']['xlsx'][1000]='Soffer_Bound/expdata/1000.xlsx' # | u | SB
conf['datasets']['SB']['xlsx'][2000]='Soffer_Bound/expdata/2000.xlsx' # | d | SB
conf['datasets']['SB']['norm']={}
for k in conf['datasets']['SB']['xlsx']: conf['datasets']['SB']['norm'][k]={'value':1,'fixed':True,'min':0,'max':1}
conf['datasets']['SB']['filters']={}
conf['SB tabs'] = READER().load_data_sets('SB')
