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')