if __name__ == '__main__':
conf['datasets'] = {}
conf['datasets']['lattice'] = {}
conf['datasets']['lattice']['xlsx'] = {}
conf['datasets']['lattice']['norm'] = {}
# lattice gT(u-d)
conf['datasets']['lattice']['xlsx'][
1000] = '../../database/lattice/1000.xlsx'
conf['datasets']['lattice']['norm'][1000] = {
'value': 1,
'fixed': True,
'min': 0,
'max': 2
}
conf['lattice tabs'] = READER(conf).load_data_sets('lattice')
#####################################################
conf['path2CJ'] = '../../external/CJLIB'
conf['path2LSS'] = '../../external/LSSLIB'
conf['path2DSS'] = '../../external/DSSLIB'
# setup for inclusive dis
conf['alphaSmode'] = 'backward'
conf['Q20'] = 1
conf['order'] = 'LO'
conf['aux'] = AUX()
conf['alphaS'] = ALPHAS(conf)
conf['pdf-NLO'] = CJ(conf)
conf['dis stfuncs'] = DIS_STFUNCS(conf)
4006] = 'sidis/expdata/4006.xlsx' # | deuteron | k+ | AUTcollins | compass | x
conf['datasets']['sidis']['xlsx'][
4008] = 'sidis/expdata/4008.xlsx' # | deuteron | k+ | AUTcollins | compass | z
conf['datasets']['sidis']['xlsx'][
4007] = 'sidis/expdata/4007.xlsx' # | deuteron | k+ | AUTcollins | compass | pt
conf['datasets']['sidis']['xlsx'][
4009] = 'sidis/expdata/4009.xlsx' # | deuteron | k- | AUTcollins | compass | x
conf['datasets']['sidis']['xlsx'][
4011] = 'sidis/expdata/4011.xlsx' # | deuteron | k- | AUTcollins | compass | z
conf['datasets']['sidis']['xlsx'][
4010] = 'sidis/expdata/4010.xlsx' # | deuteron | k- | AUTcollins | compass | pt
#conf['datasets']['sidis']['xlsx'][3001]='sidis/expdata/3001.xlsx' # | neutron | pi- | AUTcollins | jlab | x
#conf['datasets']['sidis']['xlsx'][3002]='sidis/expdata/3002.xlsx' # | neutron | pi+ | AUTcollins | jlab | x
conf['datasets']['sidis']['norm'] = {}
for k in conf['datasets']['sidis']['xlsx']:
conf['datasets']['sidis']['norm'][k] = {
'value': 1,
'fixed': True,
'min': 0,
'max': 1
}
conf['datasets']['sidis']['filters'] = {}
conf['sidis tabs'] = READER().load_data_sets('sidis')
conf['residuals'] = RESIDUALS()
print conf['residuals'].get_residuals()
#conf['residuals'].gen_report(verb=1, level=1)
conf['datasets'] = {}
conf['datasets']['sia'] = {}
conf['datasets']['sia']['xlsx'] = {}
# babar | pi,pi | AUL-0 | 9 | z1,z2,pT0 |
conf['datasets']['sia']['xlsx'][
1000] = '../../database/sia/expdata/1000.xlsx'
# babar | pi,pi | AUC-0 | 9 | z1,z2,pT0 |
conf['datasets']['sia']['xlsx'][
1001] = '../../database/sia/expdata/1001.xlsx'
# conf['datasets']['sia']['xlsx'][1002]='../database/sia/expdata/1002.xlsx' # babar | pi,pi | AUC-0 | 36 | z1,z2 |
# conf['datasets']['sia']['xlsx'][1003]='../database/sia/expdata/1003.xlsx' # babar | pi,pi | AUL-0 | 36 | z1,z2 |
# conf['datasets']['sia']['xlsx'][1004]='../database/sia/expdata/1004.xlsx' # belle | pi,pi | AUT-0-CCP | 16 | z1,z2,qT |
# conf['datasets']['sia']['xlsx'][1005]='../database/sia/expdata/1005.xlsx' # belle | pi,pi | AUT-0 | 16 | z1,z2,qT |
conf['datasets']['sia']['norm'] = {}
conf['datasets']['sia']['norm'][1000] = {'value': 1, 'fixed': False}
conf['datasets']['sia']['norm'][1001] = {'value': 1, 'fixed': False}
# conf['datasets']['sia']['norm'][1002]={'value':1,'fixed':False}
# conf['datasets']['sia']['norm'][1003]={'value':1,'fixed':False}
# conf['datasets']['sia']['norm'][1004]={'value':1,'fixed':False}
# conf['datasets']['sia']['norm'][1005]={'value':1,'fixed':False}
# conf['datasets']['sia']['filters']=[]
# conf['datasets']['sia']['filters'].append("z<0.6")
# conf['datasets']['sia']['filters'].append("Q2>1.69")
#conf['datasets']['sia']['filters'].append("pT>0.2 and pT<0.8")
conf['sia tabs'] = READER().load_data_sets('sia')
conf['residuals'] = RESIDUALS()
conf['residuals'].get_residuals()
conf['residuals'].gen_report(verb=1, level=1)
from reader import READER
conf['aux'] = AUX()
conf['pdf'] = pdf0.PDF('p')
conf['pdfpi-'] = pdf0.PDF('pi-')
conf['sivers'] = pdf1.PDF()
conf['datasets']={}
conf['datasets']['dy']={}
conf['datasets']['dy']['xlsx']={}
# COMPASS Sivers
conf['datasets']['dy']['xlsx'][1000]='dy/expdata/1000.xlsx'
conf['datasets']['dy']['xlsx'][1001]='dy/expdata/1001.xlsx'
conf['datasets']['dy']['xlsx'][1002]='dy/expdata/1002.xlsx'
conf['datasets']['dy']['xlsx'][1003]='dy/expdata/1003.xlsx'
conf['datasets']['dy']['norm']={}
for k in conf['datasets']['dy']['xlsx']: conf['datasets']['dy']['norm'][k]={'value':1,'fixed':True,'min':0,'max':1}
conf['datasets']['dy']['filters']={}
conf['dy tabs'] = READER().load_data_sets('dy')
conf['residuals']= RESIDUALS()
print(conf['residuals'].get_residuals())
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')
conf['residuals']= RESIDUALS()
print(conf['residuals'].get_residuals())
#conf['residuals'].gen_report(verb=1, level=1)
conf['aux'] = AUX()
conf['pdf'] = pdf0.PDF('p')
conf['sivers'] = pdf1.PDF()
conf['datasets'] = {}
conf['datasets']['wz'] = {}
conf['datasets']['wz']['xlsx'] = {}
# STAR Sivers
conf['datasets']['wz']['xlsx'][2000] = 'wz/expdata/2000.xlsx'
conf['datasets']['wz']['xlsx'][2001] = 'wz/expdata/2001.xlsx'
conf['datasets']['wz']['xlsx'][2001] = 'wz/expdata/2002.xlsx'
conf['datasets']['wz']['norm'] = {}
for k in conf['datasets']['wz']['xlsx']:
conf['datasets']['wz']['norm'][k] = {
'value': 1,
'fixed': True,
'min': 0,
'max': 1
}
conf['datasets']['wz']['filters'] = {}
conf['wz tabs'] = READER().load_data_sets('wz')
conf['residuals'] = RESIDUALS()
print(conf['residuals'].get_residuals())