Input file¶
In this page we provide a template to build an input.py
file to be read by resman.py
Headers¶
import os
conf={}
#--setup qcd evolution
conf['dglap mode']='truncated'
conf['alphaSmode']='backward'
conf['order'] = 'NLO'
conf['Q20'] = 1.27**2
#--setup posterior sampling
conf['bootstrap']=False
conf['flat par']=False
conf['ftol']=1e-8
#--setup for idis
conf['tmc'] = False
conf['ht'] = False
conf['nuc smearing']=True
conf['sidis nuc smearing']=False
conf['hq']=False
#--grids
conf['path2idistab'] = '%s/grids/grids-idis/distab/'%os.environ['FITPACK']
conf['path2dytab'] = '%s/grids/grids-dy/'%os.environ['FITPACK']
conf['path2jettab'] = '%s/grids/grids-jets/'%os.environ['FITPACK']
conf['path2pidistab'] = '%s/grids/grids-pidis/'%os.environ['FITPACK']
Datasets¶
First we create a key for all the datasets
conf['datasets']={}
Unpolarized DIS¶
#--lepton-hadron reactions
Q2cut=1.3**2
W2cut=4.0
pt_cut = 4.0
##--IDIS
conf['datasets']['idis']={}
conf['datasets']['idis']['filters']=[]
conf['datasets']['idis']['filters'].append("Q2>%f"%Q2cut)
conf['datasets']['idis']['filters'].append("W2>%f"%W2cut)
conf['datasets']['idis']['xlsx']={}
conf['datasets']['idis']['xlsx'][10010]='idis/expdata/10010.xlsx' # proton | F2 | SLAC
conf['datasets']['idis']['xlsx'][10011]='idis/expdata/10011.xlsx' # deuteron | F2 | SLAC
conf['datasets']['idis']['xlsx'][10016]='idis/expdata/10016.xlsx' # proton | F2 | BCDMS
conf['datasets']['idis']['xlsx'][10017]='idis/expdata/10017.xlsx' # deuteron | F2 | BCDMS
conf['datasets']['idis']['xlsx'][10020]='idis/expdata/10020.xlsx' # proton | F2 | NMC
conf['datasets']['idis']['xlsx'][10021]='idis/expdata/10021.xlsx' # d/p | F2d/F2p | NMC
conf['datasets']['idis']['xlsx'][10026]='idis/expdata/10026.xlsx' # proton | sigma red | HERA II NC e+ (1)
conf['datasets']['idis']['xlsx'][10027]='idis/expdata/10027.xlsx' # proton | sigma red | HERA II NC e+ (2)
conf['datasets']['idis']['xlsx'][10028]='idis/expdata/10028.xlsx' # proton | sigma red | HERA II NC e+ (3)
conf['datasets']['idis']['xlsx'][10029]='idis/expdata/10029.xlsx' # proton | sigma red | HERA II NC e+ (4)
conf['datasets']['idis']['xlsx'][10030]='idis/expdata/10030.xlsx' # proton | sigma red | HERA II NC e-
conf['datasets']['idis']['xlsx'][10031]='idis/expdata/10031.xlsx' # proton | sigma red | HERA II CC e+
conf['datasets']['idis']['xlsx'][10032]='idis/expdata/10032.xlsx' # proton | sigma red | HERA II NC e-
conf['datasets']['idis']['norm']={}
conf['datasets']['idis']['norm'][10010]={'value': 1, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['idis']['norm'][10011]={'value': 1, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['idis']['norm'][10016]={'value': 9, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['idis']['norm'][10017]={'value': 1, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['idis']['norm'][10020]={'value': 1, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
Polarized DIS¶
conf['datasets']['pidis']={}
conf['datasets']['pidis']['filters']=[]
conf['datasets']['pidis']['filters'].append("Q2>%f"%Q2cut)
conf['datasets']['pidis']['filters'].append("W2>%f"%W2cut)
conf['datasets']['pidis']['xlsx']={}
conf['datasets']['pidis']['xlsx'][10001]='pidis/expdata/10001.xlsx' # 10001 | deuteron | A1 | COMPASS |
conf['datasets']['pidis']['xlsx'][10006]='pidis/expdata/10006.xlsx' # 10006 | deuteron | Apa | HERMES |
conf['datasets']['pidis']['xlsx'][10020]='pidis/expdata/10020.xlsx' # 10020 | deuteron | Ape | SLAC(E143) |
conf['datasets']['pidis']['xlsx'][10021]='pidis/expdata/10021.xlsx' # 10021 | deuteron | Apa | SLAC(E143) |
conf['datasets']['pidis']['xlsx'][10026]='pidis/expdata/10026.xlsx' # 10026 | deuteron | Ape | SLAC(E155) |
conf['datasets']['pidis']['xlsx'][10027]='pidis/expdata/10027.xlsx' # 10027 | deuteron | Apa | SLAC(E155) |
#conf['datasets']['pidis']['xlsx'][10030]='pidis/expdata/10030.xlsx' # 10030 | deuteron | Atpe | SLAC(E155x) |
conf['datasets']['pidis']['xlsx'][10033]='pidis/expdata/10033.xlsx' # 10033 | deuteron | A1 | SMC |
conf['datasets']['pidis']['xlsx'][10034]='pidis/expdata/10034.xlsx' # 10034 | deuteron | A1 | SMC |
conf['datasets']['pidis']['xlsx'][10002]='pidis/expdata/10002.xlsx' # 10002 | proton | A1 | COMPASS |
conf['datasets']['pidis']['xlsx'][10003]='pidis/expdata/10003.xlsx' # 10003 | proton | A1 | COMPASS |
conf['datasets']['pidis']['xlsx'][10004]='pidis/expdata/10004.xlsx' # 10004 | proton | A1 | EMC |
conf['datasets']['pidis']['xlsx'][10007]='pidis/expdata/10007.xlsx' # 10007 | proton | Apa | HERMES |
conf['datasets']['pidis']['xlsx'][10008]='pidis/expdata/10008.xlsx' # 10008 | proton | A2 | HERMES |
conf['datasets']['pidis']['xlsx'][10022]='pidis/expdata/10022.xlsx' # 10022 | proton | Apa | SLAC(E143) |
conf['datasets']['pidis']['xlsx'][10023]='pidis/expdata/10023.xlsx' # 10023 | proton | Ape | SLAC(E143) |
conf['datasets']['pidis']['xlsx'][10028]='pidis/expdata/10028.xlsx' # 10028 | proton | Ape | SLAC(E155) |
conf['datasets']['pidis']['xlsx'][10029]='pidis/expdata/10029.xlsx' # 10029 | proton | Apa | SLAC(E155) |
#conf['datasets']['pidis']['xlsx'][10031]='pidis/expdata/10031.xlsx' # 10031 | proton | Atpe | SLAC(E155x) |
conf['datasets']['pidis']['xlsx'][10032]='pidis/expdata/10032.xlsx' # 10032 | proton | Apa | SLACE80E130 |
conf['datasets']['pidis']['xlsx'][10035]='pidis/expdata/10035.xlsx' # 10035 | proton | A1 | SMC |
conf['datasets']['pidis']['xlsx'][10036]='pidis/expdata/10036.xlsx' # 10036 | proton | A1 | SMC |
conf['datasets']['pidis']['xlsx'][10005]='pidis/expdata/10005.xlsx' # 10005 | neutron | A1 | HERMES |
#conf['datasets']['pidis']['xlsx'][10018]='pidis/expdata/10018.xlsx' # 10018 | helium | A1 | SLAC(E142) |
#conf['datasets']['pidis']['xlsx'][10019]='pidis/expdata/10019.xlsx' # 10019 | helium | A2 | SLAC(E142) |
#conf['datasets']['pidis']['xlsx'][10024]='pidis/expdata/10024.xlsx' # 10024 | helium | Ape | SLAC(E154) |
#conf['datasets']['pidis']['xlsx'][10025]='pidis/expdata/10025.xlsx' # 10025 | helium | Apa | SLAC(E154) |
conf['datasets']['pidis']['norm']={}
conf['datasets']['pidis']['norm'][10001]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10006]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10020]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10021]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10026]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10027]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
#conf['datasets']['pidis']['norm'][10030]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10033]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10034]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10002]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10003]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10004]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10007]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10008]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10022]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10023]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10028]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10029]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
#conf['datasets']['pidis']['norm'][10031]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10032]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10035]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10036]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
conf['datasets']['pidis']['norm'][10005]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
#conf['datasets']['pidis']['norm'][10018]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
#conf['datasets']['pidis']['norm'][10019]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
#conf['datasets']['pidis']['norm'][10024]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
#conf['datasets']['pidis']['norm'][10025]={'value': 1.00000, 'min': 8.00000e-01, 'max': 1.20000e+00, 'fixed': False}
Unpolarized SIDIS¶
conf['datasets']['sidis']={}
conf['datasets']['sidis']['filters']=[]
conf['datasets']['sidis']['filters'].append("Q2>%f"%Q2cut)
conf['datasets']['sidis']['filters'].append("W2>%f"%W2cut)
conf['datasets']['sidis']['filters'].append('Z>0.2 and Z<0.8')
conf['datasets']['sidis']['xlsx']={}
conf['datasets']['sidis']['xlsx'][1005]='sidis/expdata/1005.xlsx' # deuteron , mult , pi+ , COMPASS
conf['datasets']['sidis']['xlsx'][1006]='sidis/expdata/1006.xlsx' # deuteron , mult , pi- , COMPASS
conf['datasets']['sidis']['xlsx'][2005]='sidis/expdata/2005.xlsx' # deuteron , mult , K+ , COMPASS
conf['datasets']['sidis']['xlsx'][2006]='sidis/expdata/2006.xlsx' # deuteron , mult , K- , COMPASS
conf['datasets']['sidis']['xlsx'][1007]='sidis/expdata/1007.xlsx' # proton , mult , pi+ , HERMES, Q2-z
conf['datasets']['sidis']['xlsx'][1008]='sidis/expdata/1008.xlsx' # proton , mult , pi- , HERMES, Q2-z
conf['datasets']['sidis']['xlsx'][1009]='sidis/expdata/1009.xlsx' # deuteron , mult , pi+ , HERMES, Q2-z
conf['datasets']['sidis']['xlsx'][1010]='sidis/expdata/1010.xlsx' # deuteron , mult , pi- , HERMES, Q2-z
conf['datasets']['sidis']['xlsx'][2007]='sidis/expdata/2007.xlsx' # proton , mult , K+ , HERMES, Q2-z
conf['datasets']['sidis']['xlsx'][2008]='sidis/expdata/2008.xlsx' # proton , mult , K- , HERMES, Q2-z
conf['datasets']['sidis']['xlsx'][2009]='sidis/expdata/2009.xlsx' # deuteron , mult , K+ , HERMES, Q2-z
conf['datasets']['sidis']['xlsx'][2010]='sidis/expdata/2010.xlsx' # deuteron , mult , K- , HERMES, Q2-z
conf['datasets']['sidis']['norm']={}
conf['datasets']['sidis']['norm'][1005]={'value': 1.0000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sidis']['norm'][1006]={'value': 1.0000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sidis']['norm'][1007]={'value': 1.0000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sidis']['norm'][1008]={'value': 1.0000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sidis']['norm'][1009]={'value': 1.0000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sidis']['norm'][1010]={'value': 1.0000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sidis']['norm'][2005]={'value': 1.0000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sidis']['norm'][2006]={'value': 1.0000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sidis']['norm'][2007]={'value': 1.0000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sidis']['norm'][2008]={'value': 1.0000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sidis']['norm'][2009]={'value': 1.0000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sidis']['norm'][2010]={'value': 1.0000000,'fixed':False,'min':0.5,'max':1.5}
Polarized SIDIS¶
conf['datasets']['psidis']={}
conf['datasets']['psidis']['filters']=[]
conf['datasets']['psidis']['filters'].append("Q2>%f"%Q2cut)
conf['datasets']['psidis']['filters'].append("W2>%f"%W2cut)
conf['datasets']['psidis']['filters'].append('Z>0.2 and Z<0.8')
conf['datasets']['psidis']['xlsx']={}
conf['datasets']['psidis']['xlsx'][20004]='psidis/expdata/20004.xlsx' # 20004 | proton | A1pi+ | HERMES
conf['datasets']['psidis']['xlsx'][20005]='psidis/expdata/20005.xlsx' # 20005 | proton | A1pi- | HERMES
conf['datasets']['psidis']['xlsx'][20008]='psidis/expdata/20008.xlsx' # 20008 | deuteron | A1pi+ | HERMES
conf['datasets']['psidis']['xlsx'][20009]='psidis/expdata/20009.xlsx' # 20009 | deuteron | A1pi- | HERMES
conf['datasets']['psidis']['xlsx'][20012]='psidis/expdata/20012.xlsx' # 20012 | deuteron | A1K+ | HERMES
conf['datasets']['psidis']['xlsx'][20013]='psidis/expdata/20013.xlsx' # 20013 | deuteron | A1K- | HERMES
conf['datasets']['psidis']['xlsx'][20014]='psidis/expdata/20014.xlsx' # 20014 | deuteron | A1Ksum | HERMES
conf['datasets']['psidis']['xlsx'][20017]='psidis/expdata/20017.xlsx' # 20017 | proton | A1pi+ | COMPASS
conf['datasets']['psidis']['xlsx'][20018]='psidis/expdata/20018.xlsx' # 20018 | proton | A1pi- | COMPASS
conf['datasets']['psidis']['xlsx'][20019]='psidis/expdata/20019.xlsx' # 20019 | proton | A1K+ | COMPASS
conf['datasets']['psidis']['xlsx'][20020]='psidis/expdata/20020.xlsx' # 20020 | proton | A1K- | COMPASS
conf['datasets']['psidis']['xlsx'][20021]='psidis/expdata/20021.xlsx' # 20021 | deuteron | A1pi+ | COMPASS
conf['datasets']['psidis']['xlsx'][20022]='psidis/expdata/20022.xlsx' # 20022 | deuteron | A1pi- | COMPASS
conf['datasets']['psidis']['xlsx'][20025]='psidis/expdata/20025.xlsx' # 20025 | deuteron | A1K+ | COMPASS
conf['datasets']['psidis']['xlsx'][20026]='psidis/expdata/20026.xlsx' # 20026 | deuteron | A1K- | COMPASS
conf['datasets']['psidis']['norm']={}
conf['datasets']['psidis']['norm'][20004]={'value': 1.000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['psidis']['norm'][20005]={'value': 1.000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['psidis']['norm'][20008]={'value': 1.000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['psidis']['norm'][20009]={'value': 1.000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['psidis']['norm'][20012]={'value': 1.000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['psidis']['norm'][20013]={'value': 1.000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['psidis']['norm'][20014]={'value': 1.000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['psidis']['norm'][20017]={'value': 1.000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['psidis']['norm'][20018]={'value': 1.000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['psidis']['norm'][20019]={'value': 1.000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['psidis']['norm'][20020]={'value': 1.000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['psidis']['norm'][20021]={'value': 1.000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['psidis']['norm'][20022]={'value': 1.000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['psidis']['norm'][20025]={'value': 1.000000,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['psidis']['norm'][20026]={'value': 1.000000,'fixed':False,'min':0.5,'max':1.5}
Drell-Yan¶
conf['datasets']['dy']={}
conf['datasets']['dy']['filters']=[]
conf['datasets']['dy']['filters'].append("Q2>36")
conf['datasets']['dy']['xlsx']={}
conf['datasets']['dy']['xlsx'][10001]='dy/expdata/10001.xlsx'
conf['datasets']['dy']['xlsx'][10002]='dy/expdata/10002.xlsx'
conf['datasets']['dy']['norm']={}
conf['datasets']['dy']['norm'][10001]={'value': 1,'fixed':False,'min': 0.8,'max': 1.2}
conf['datasets']['dy']['norm'][10002]={'value': 1,'fixed':False,'min': 0.8,'max': 1.2}
SIA pion¶
conf['datasets']['sia']={}
conf['datasets']['sia']['filters']=[]
conf['datasets']['sia']['filters'].append('z>0.2 and z<0.9')
conf['datasets']['sia']['xlsx']={}
conf['datasets']['sia']['xlsx'][1001]='sia/expdata/1001.xlsx' # hadron: pion exp: TASSO
conf['datasets']['sia']['xlsx'][1002]='sia/expdata/1002.xlsx' # hadron: pion exp: TASSO
conf['datasets']['sia']['xlsx'][1003]='sia/expdata/1003.xlsx' # hadron: pion exp: TASSO
conf['datasets']['sia']['xlsx'][1004]='sia/expdata/1004.xlsx' # hadron: pion exp: TASSO
conf['datasets']['sia']['xlsx'][1005]='sia/expdata/1005.xlsx' # hadron: pion exp: TASSO
conf['datasets']['sia']['xlsx'][1006]='sia/expdata/1006.xlsx' # hadron: pion exp: TASSO
conf['datasets']['sia']['xlsx'][1007]='sia/expdata/1007.xlsx' # hadron: pion exp: TPC
conf['datasets']['sia']['xlsx'][1008]='sia/expdata/1008.xlsx' # hadron: pion exp: TPC
conf['datasets']['sia']['xlsx'][1012]='sia/expdata/1012.xlsx' # hadron: pion exp: HRS
conf['datasets']['sia']['xlsx'][1013]='sia/expdata/1013.xlsx' # hadron: pion exp: TOPAZ
conf['datasets']['sia']['xlsx'][1014]='sia/expdata/1014.xlsx' # hadron: pion exp: SLD
conf['datasets']['sia']['xlsx'][1018]='sia/expdata/1018.xlsx' # hadron: pion exp: ALEPH
conf['datasets']['sia']['xlsx'][1019]='sia/expdata/1019.xlsx' # hadron: pion exp: OPAL
conf['datasets']['sia']['xlsx'][1025]='sia/expdata/1025.xlsx' # hadron: pion exp: DELPHI
conf['datasets']['sia']['xlsx'][1028]='sia/expdata/1028.xlsx' # hadron: pion exp: BABAR
conf['datasets']['sia']['xlsx'][1029]='sia/expdata/1029.xlsx' # hadron: pion exp: BELL
conf['datasets']['sia']['xlsx'][1030]='sia/expdata/1030.xlsx' # hadron: pion exp: ARGUS
conf['datasets']['sia']['norm']={}
conf['datasets']['sia']['norm'][1001]={'value': 1.10478e+00,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][1002]={'value': 9.82581e-01,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][1003]={'value': 1.03054e+00,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][1005]={'value': 1.03419e+00,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][1006]={'value': 9.79162e-01,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][1014]={'value': 9.97770e-01,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][1018]={'value': 1.02378e+00,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][1028]={'value': 9.76001e-01,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][1029]={'value': 8.68358e-01,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][1030]={'value': 1.01509e+00,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][1007]={'value': 1.06806e+00,'fixed':False,'min':0.5,'max':1.5,'dN':0.1}
conf['datasets']['sia']['norm'][1008]={'value': 1.06603e+00,'fixed':False,'min':0.5,'max':1.5,'dN':0.1}
conf['datasets']['sia']['norm'][1019]={'value': 1.03100e+00,'fixed':False,'min':0.5,'max':1.5,'dN':0.1}
##--SIA pion HQ
conf['datasets']['sia']['xlsx'][1010]='sia/expdata/1010.xlsx' # hadron: pion exp: TPC(c)
conf['datasets']['sia']['xlsx'][1011]='sia/expdata/1011.xlsx' # hadron: pion exp: TPC(b)
conf['datasets']['sia']['xlsx'][1016]='sia/expdata/1016.xlsx' # hadron: pion exp: SLD(c)
conf['datasets']['sia']['xlsx'][1017]='sia/expdata/1017.xlsx' # hadron: pion exp: SLD(b)
conf['datasets']['sia']['xlsx'][1023]='sia/expdata/1023.xlsx' # hadron: pion exp: OPAL(c)
conf['datasets']['sia']['xlsx'][1024]='sia/expdata/1024.xlsx' # hadron: pion exp: OPAL(b)
conf['datasets']['sia']['xlsx'][1027]='sia/expdata/1027.xlsx' # hadron: pion exp: DELPHI(b)
conf['datasets']['sia']['norm'][1016]={'value': 1.18920e+00,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][1017]={'value': 1.00345e+00,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][1023]={'value': 1.33434e+00,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][1024]={'value': 1.19151e+00,'fixed':False,'min':0.5,'max':1.5}
SIA Kaon¶
conf['datasets']['sia']['xlsx'][2030]='sia/expdata/2030.xlsx' # hadron: kaon exp: ARGUS
conf['datasets']['sia']['xlsx'][2028]='sia/expdata/2028.xlsx' # hadron: kaon exp: BABAR
conf['datasets']['sia']['xlsx'][2029]='sia/expdata/2029.xlsx' # hadron: kaon exp: BELL
conf['datasets']['sia']['xlsx'][2001]='sia/expdata/2001.xlsx' # hadron: kaon exp: TASSO
conf['datasets']['sia']['xlsx'][2002]='sia/expdata/2002.xlsx' # hadron: kaon exp: TASSO
conf['datasets']['sia']['xlsx'][2003]='sia/expdata/2003.xlsx' # hadron: kaon exp: TASSO
conf['datasets']['sia']['xlsx'][2004]='sia/expdata/2004.xlsx' # hadron: kaon exp: TASSO
conf['datasets']['sia']['xlsx'][2005]='sia/expdata/2005.xlsx' # hadron: kaon exp: TASSO
conf['datasets']['sia']['xlsx'][2006]='sia/expdata/2006.xlsx' # hadron: kaon exp: TASSO
conf['datasets']['sia']['xlsx'][2007]='sia/expdata/2007.xlsx' # hadron: kaon exp: TPC
conf['datasets']['sia']['xlsx'][2008]='sia/expdata/2008.xlsx' # hadron: kaon exp: TPC
conf['datasets']['sia']['xlsx'][2012]='sia/expdata/2012.xlsx' # hadron: kaon exp: HRS
conf['datasets']['sia']['xlsx'][2013]='sia/expdata/2013.xlsx' # hadron: kaon exp: TOPAZ
conf['datasets']['sia']['xlsx'][2014]='sia/expdata/2014.xlsx' # hadron: kaon exp: SLD
conf['datasets']['sia']['xlsx'][2018]='sia/expdata/2018.xlsx' # hadron: kaon exp: ALEPH
conf['datasets']['sia']['xlsx'][2019]='sia/expdata/2019.xlsx' # hadron: kaon exp: OPAL
conf['datasets']['sia']['xlsx'][2025]='sia/expdata/2025.xlsx' # hadron: kaon exp: DELPHI
conf['datasets']['sia']['xlsx'][2031]='sia/expdata/2031.xlsx' # hadron: kaon exp: DELPHI
conf['datasets']['sia']['norm'][2030]={'value': 1.00482e+00,'fixed':False,'min':0.5,'max':1.5,'dN':0.1}
conf['datasets']['sia']['norm'][2028]={'value': 9.97435e-01,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][2029]={'value': 1.00000e+00,'fixed':True,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][2002]={'value': 9.83394e-01,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][2003]={'value': 9.94421e-01,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][2005]={'value': 9.92876e-01,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][2014]={'value': 9.12186e-01,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][2018]={'value': 8.88453e-01,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][2019]={'value': 8.40305e-01,'fixed':False,'min':0.5,'max':1.5,'dN':0.1}
conf['datasets']['sia']['norm'][2007]={'value': 1.20488e+00,'fixed':False,'min':0.5,'max':1.5,'dN':0.1}
conf['datasets']['sia']['norm'][2008]={'value': 1.10314e+00,'fixed':False,'min':0.5,'max':1.5,'dN':0.1}
##--SIA kaon HQ
conf['datasets']['sia']['xlsx'][2016]='sia/expdata/2016.xlsx' # hadron: kaon exp: SLD(c)
conf['datasets']['sia']['xlsx'][2017]='sia/expdata/2017.xlsx' # hadron: kaon exp: SLD(b)
conf['datasets']['sia']['xlsx'][2023]='sia/expdata/2023.xlsx' # hadron: kaon exp: OPAL(c)
conf['datasets']['sia']['xlsx'][2024]='sia/expdata/2024.xlsx' # hadron: kaon exp: OPAL(b)
conf['datasets']['sia']['xlsx'][2027]='sia/expdata/2027.xlsx' # hadron: kaon exp: DELPHI(b)
conf['datasets']['sia']['norm'][2016]={'value': 1.03996e+00,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][2017]={'value': 1.00015e+00,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][2023]={'value': 1.35922e+00,'fixed':False,'min':0.5,'max':1.5}
conf['datasets']['sia']['norm'][2024]={'value': 1.38163e+00,'fixed':False,'min':0.5,'max':1.5}
Unpolarized JETs¶
conf['datasets']['jet'] = {}
conf['datasets']['jet']['xlsx'] = {}
conf['datasets']['jet']['xlsx'][10001] = 'jets/expdata/10001.xlsx' ## D0 dataset
conf['datasets']['jet']['xlsx'][10002] = 'jets/expdata/10002.xlsx' ## CDF dataset
conf['datasets']['jet']['xlsx'][10003] = 'jets/expdata/10003.xlsx' ## STAR MB dataset
conf['datasets']['jet']['xlsx'][10004] = 'jets/expdata/10004.xlsx' ## STAR HT dataset
conf['datasets']['jet']['norm'] = {}
conf['datasets']['jet']['norm'][10001] = {'value': 1, 'fixed': True, 'min': 0.8, 'max': 1.2}
conf['datasets']['jet']['norm'][10002] = {'value': 1, 'fixed': False, 'min': 0.8, 'max': 1.2}
conf['datasets']['jet']['norm'][10003] = {'value': 1, 'fixed': True, 'min': 0.8, 'max': 1.2}
conf['datasets']['jet']['norm'][10004] = {'value': 1, 'fixed': True, 'min': 0.8, 'max': 1.2}
conf['channels'] = {'q_qp': 0, 'q_qbp': 1, 'q_q': 2, 'q_qb': 3, 'q_g': 4, 'g_g': 5}
conf['qr'] = {}
conf['qr']['parameter'] = {'epsilon': 1e-11, 'block_size': 10, 'power_scheme': 1}
conf['fit'] = {'method': 'fixed', 'f_scale': 0.25, 'r_scale': 0.25}
Polarized JETs¶
conf['datasets']['pjet'] = {}
conf['datasets']['pjet']['filters'] = []
conf['datasets']['pjet']['filters'].append("pT>%f" % pt_cut)
conf['datasets']['pjet']['xlsx'] = {}
conf['datasets']['pjet']['xlsx'][20001] = 'pjets/expdata/20001.xlsx' ## STAR 2006 paper on 2003 and 2004 data
conf['datasets']['pjet']['xlsx'][20002] = 'pjets/expdata/20002.xlsx' ## STAR 2012 paper on 2005 data
conf['datasets']['pjet']['xlsx'][20003] = 'pjets/expdata/20003.xlsx' ## STAR 2012 paper on 2006 data
conf['datasets']['pjet']['xlsx'][20004] = 'pjets/expdata/20004.xlsx' ## STAR 2015 paper on 2009 data
conf['datasets']['pjet']['xlsx'][20005] = 'pjets/expdata/20005.xlsx' ## PHENIX 2011 paper on 2005 data
conf['datasets']['pjet']['xlsx'][99999] = 'pjets/expdata/99999.xlsx' ## STAR 2019 paper on 2012 data
conf['datasets']['pjet']['norm'] = {}
conf['datasets']['pjet']['norm'][20001] = {'value': 1, 'fixed': True, 'min': 0.8, 'max': 1.2}
conf['datasets']['pjet']['norm'][20002] = {'value': 1, 'fixed': False, 'min': 0.8, 'max': 1.2}
conf['datasets']['pjet']['norm'][20003] = {'value': 1, 'fixed': False, 'min': 0.8, 'max': 1.2}
conf['datasets']['pjet']['norm'][20004] = {'value': 1, 'fixed': False, 'min': 0.8, 'max': 1.2}
# conf['datasets']['pjet']['norm'][20005] = {'value': 1, 'fixed': False, 'min': 0.8, 'max': 1.2}
conf['datasets']['pjet']['norm'][20006] = {'value': 1, 'fixed': False, 'min': 0.8, 'max': 1.2}
conf['channels'] = {'q_qp': 0, 'q_qbp': 1, 'q_q': 2, 'q_qb': 3, 'q_g': 4, 'g_g': 5}
conf['qr'] = {}
conf['qr']['parameter'] = {'epsilon': 1e-11, 'block_size': 10, 'power_scheme': 1}
conf['fit'] = {'method': 'fixed', 'f_scale': 0.25, 'r_scale': 0.25}
Parameters¶
#--parameters
conf['params'] = {}
pdfs¶
conf['pdf parametrization'] = 2
conf['params']['pdf'] = {}
conf['params']['pdf']['g1 N'] ={'value': 3.87592e-01, 'min': None, 'max': None, 'fixed': True }
conf['params']['pdf']['g1 a'] ={'value': -6.60217e-01, 'min': -1.9, 'max': 1, 'fixed': False}
conf['params']['pdf']['g1 b'] ={'value': 8.12091e+00, 'min': 0, 'max': 10, 'fixed': False}
conf['params']['pdf']['uv1 N'] ={'value': 3.47549e-01, 'min': None, 'max': None, 'fixed': True }
conf['params']['pdf']['uv1 a'] ={'value': -1.14055e-01, 'min': -0.5, 'max': 1, 'fixed': False}
conf['params']['pdf']['uv1 b'] ={'value': 3.21230e+00, 'min': 0, 'max': 10, 'fixed': False}
conf['params']['pdf']['dv1 N'] ={'value': 1.52089e-01, 'min': None, 'max': None, 'fixed': True }
conf['params']['pdf']['dv1 a'] ={'value': -3.80099e-02, 'min': -0.5, 'max': 1, 'fixed': False}
conf['params']['pdf']['dv1 b'] ={'value': 4.36319e+00, 'min': 0, 'max': 10, 'fixed': False}
conf['params']['pdf']['db1 N'] ={'value': 2.08640e-02, 'min': 0, 'max': 1, 'fixed': False}
conf['params']['pdf']['db1 a'] ={'value': -4.05372e-01, 'min': -1, 'max': 1, 'fixed': False}
conf['params']['pdf']['db1 b'] ={'value': 1.00000e+01, 'min': 0, 'max': 10, 'fixed': False}
conf['params']['pdf']['ub1 N'] ={'value': 1.64626e-02, 'min': 0, 'max': 1, 'fixed': False}
conf['params']['pdf']['ub1 a'] ={'value': 0.000000000, 'min': -1, 'max': 1, 'fixed': False}
conf['params']['pdf']['ub1 b'] ={'value': 1.00000e+01, 'min': 0, 'max': 10, 'fixed': False}
conf['params']['pdf']['s1 N'] ={'value': 0.00000e+00, 'min': 0, 'max': 1, 'fixed': True }
conf['params']['pdf']['s1 a'] ={'value': -0.19319e+00, 'min': -0.5, 'max': 1, 'fixed': False}
conf['params']['pdf']['s1 b'] ={'value': 4.99999e+00, 'min': 0, 'max': 10, 'fixed': False}
conf['params']['pdf']['sb1 N'] ={'value': 0.00000e+00, 'min': 0, 'max': 1, 'fixed': False}
conf['params']['pdf']['sb1 a'] ={'value': -0.19317e+00, 'min': -0.5, 'max': 1, 'fixed': False}
conf['params']['pdf']['sb1 b'] ={'value': 4.99988e+00, 'min': 0, 'max': 10, 'fixed': False}
conf['params']['pdf']['sea1 N'] ={'value': 2.08640e-02, 'min': 0, 'max': 1, 'fixed': False}
conf['params']['pdf']['sea1 a'] ={'value': -1.500000000, 'min': -1.9, 'max': -1, 'fixed': False}
conf['params']['pdf']['sea1 b'] ={'value': 1.00000e+01, 'min': 0, 'max': 10, 'fixed': False}
conf['params']['pdf']['sea2 N'] ={'value': 2.08640e-02, 'min': 0, 'max': 1, 'fixed': 'sea1 N'}
conf['params']['pdf']['sea2 a'] ={'value': -1.500000000, 'min': -1.9, 'max': -1, 'fixed': 'sea1 a'}
conf['params']['pdf']['sea2 b'] ={'value': 1.00000e+01, 'min': 0, 'max': 10, 'fixed': 'sea1 b'}
spin pdfs¶
conf['ppdf parametrization']=0
conf['su2+su3']=False
conf['params']['ppdf']={}
conf['params']['ppdf']['g1 N'] ={'value': 3.87592e-01, 'min': -10, 'max': 10, 'fixed': False}
conf['params']['ppdf']['g1 a'] ={'value': -6.60217e-01, 'min': -1, 'max': 2, 'fixed': False}
conf['params']['ppdf']['g1 b'] ={'value': 8.12091e+00, 'min': 0, 'max': 10, 'fixed': False}
conf['params']['ppdf']['g1 d'] ={'value': 0.21230e+00, 'min': -10, 'max': 10, 'fixed': False}
conf['params']['ppdf']['up1 N'] ={'value': 3.47549e-01, 'min': -10, 'max': 10, 'fixed': False}
conf['params']['ppdf']['up1 a'] ={'value': -1.14055e-01, 'min': -1, 'max': 2, 'fixed': False}
conf['params']['ppdf']['up1 b'] ={'value': 3.21230e+00, 'min': 0, 'max': 10, 'fixed': False}
conf['params']['ppdf']['up1 d'] ={'value': 0.21230e+00, 'min': -10, 'max': 10, 'fixed': False}
conf['params']['ppdf']['dp1 N'] ={'value': 1.52089e-01, 'min': -10, 'max': 10, 'fixed': False}
conf['params']['ppdf']['dp1 a'] ={'value': -3.80099e-02, 'min': -1, 'max': 2, 'fixed': False}
conf['params']['ppdf']['dp1 b'] ={'value': 4.36319e+00, 'min': 0, 'max': 10, 'fixed': False}
conf['params']['ppdf']['dp1 d'] ={'value': 0.36319e+00, 'min': -10, 'max': 10, 'fixed': False}
conf['params']['ppdf']['sp1 N'] ={'value': 1.52089e-01, 'min': -10, 'max': 10, 'fixed': False}
conf['params']['ppdf']['sp1 a'] ={'value': -3.80099e-02, 'min': -1, 'max': 2, 'fixed': False}
conf['params']['ppdf']['sp1 b'] ={'value': 4.36319e+00, 'min': 0, 'max': 10, 'fixed': False}
conf['params']['ppdf']['sp1 d'] ={'value': 0.36319e+00, 'min': -10, 'max': 10, 'fixed': False}
conf['params']['ppdf']['um1 N'] ={'value': 3.47549e-01, 'min': -10, 'max': 10, 'fixed': False}
conf['params']['ppdf']['um1 a'] ={'value': -1.14055e-01, 'min': -1, 'max': 2, 'fixed': False}
conf['params']['ppdf']['um1 b'] ={'value': 3.21230e+00, 'min': 0, 'max': 10, 'fixed': False}
conf['params']['ppdf']['um1 d'] ={'value': 0.21230e+00, 'min': -10, 'max': 10, 'fixed': False}
conf['params']['ppdf']['dm1 N'] ={'value': 1.52089e-01, 'min': -10, 'max': 10, 'fixed': False}
conf['params']['ppdf']['dm1 a'] ={'value': -3.80099e-02, 'min': -1, 'max': 2, 'fixed': False}
conf['params']['ppdf']['dm1 b'] ={'value': 4.36319e+00, 'min': 0, 'max': 10, 'fixed': False}
conf['params']['ppdf']['dm1 d'] ={'value': 0.36319e+00, 'min': -10, 'max': 10, 'fixed': False}
conf['params']['ppdf']['sm1 N'] ={'value': 1.52089e-01, 'min': -10, 'max': 10, 'fixed': False}
conf['params']['ppdf']['sm1 a'] ={'value': -3.80099e-02, 'min': -1, 'max': 2, 'fixed': False}
conf['params']['ppdf']['sm1 b'] ={'value': 4.36319e+00, 'min': 0, 'max': 10, 'fixed': False}
conf['params']['ppdf']['sm1 d'] ={'value': 0.36319e+00, 'min': -10, 'max': 10, 'fixed': False}
pion fragamentation functions¶
conf['ffpion parametrization']=0
conf['params']['ffpion']={}
conf['params']['ffpion']['g1 N'] ={'value': 2.95437e-01, 'min': 0.0 , 'max': 1,'fixed':False}
conf['params']['ffpion']['g1 a'] ={'value': 1.00469e+00, 'min': -1.8 , 'max': 2,'fixed':False}
conf['params']['ffpion']['g1 b'] ={'value': 6.85766e+00, 'min': 0 , 'max': 10,'fixed':False}
conf['params']['ffpion']['u1 N'] ={'value': 2.67821e-02, 'min': 0.0 , 'max': 1,'fixed':False}
conf['params']['ffpion']['u1 a'] ={'value': 1.76877e-01, 'min': -1.8 , 'max': 2,'fixed':False}
conf['params']['ffpion']['u1 b'] ={'value': 4.81521e+00, 'min': 0 , 'max': 10,'fixed':False}
conf['params']['ffpion']['d1 N'] ={'value': 2.99974e-01, 'min': 0.0 , 'max': 1,'fixed':False}
conf['params']['ffpion']['d1 a'] ={'value': -6.89477e-01, 'min': -1.8 , 'max': 2,'fixed':False}
conf['params']['ffpion']['d1 b'] ={'value': 4.79992e+00, 'min': 0 , 'max': 10,'fixed':False}
conf['params']['ffpion']['s1 N'] ={'value': 1.54863e-01, 'min': 0.0 , 'max': 1,'fixed':'d1 N'}
conf['params']['ffpion']['s1 a'] ={'value': 3.00305e-01, 'min': -1.8 , 'max': 2,'fixed':'d1 a'}
conf['params']['ffpion']['s1 b'] ={'value': 1.83178e+00, 'min': 0 , 'max': 10,'fixed':'d1 b'}
conf['params']['ffpion']['c1 N'] ={'value': 1.84550e-01, 'min': 0.0 , 'max': 1,'fixed':False}
conf['params']['ffpion']['c1 a'] ={'value': -5.05798e-02, 'min': -1.8 , 'max': 2,'fixed':False}
conf['params']['ffpion']['c1 b'] ={'value': 3.19952e+00, 'min': 0 , 'max': 10,'fixed':False}
conf['params']['ffpion']['b1 N'] ={'value': 3.74125e-01, 'min': 0.0 , 'max': 1,'fixed':False}
conf['params']['ffpion']['b1 a'] ={'value': -1.59541e+00, 'min': -1.8 , 'max': 2,'fixed':False}
conf['params']['ffpion']['b1 b'] ={'value': 4.50102e+00, 'min': 0 , 'max': 10,'fixed':False}
conf['params']['ffpion']['ub1 N'] ={'value': 2.99974e-01, 'min': 0.0 , 'max': 1,'fixed':'d1 N'}
conf['params']['ffpion']['ub1 a'] ={'value': -6.89477e-01, 'min': -1.8 , 'max': 2,'fixed':'d1 a'}
conf['params']['ffpion']['ub1 b'] ={'value': 4.79992e+00, 'min': 0 , 'max': 10,'fixed':'d1 b'}
conf['params']['ffpion']['db1 N'] ={'value': 2.67821e-02, 'min': 0.0 , 'max': 1,'fixed':'u1 N'}
conf['params']['ffpion']['db1 a'] ={'value': 1.76877e-01, 'min': -1.8 , 'max': 2,'fixed':'u1 a'}
conf['params']['ffpion']['db1 b'] ={'value': 4.81521e+00, 'min': 0 , 'max': 10,'fixed':'u1 b'}
conf['params']['ffpion']['sb1 N'] ={'value': 1.54863e-01, 'min': 0.0 , 'max': 1,'fixed':'d1 N'}
conf['params']['ffpion']['sb1 a'] ={'value': 3.00305e-01, 'min': -1.8 , 'max': 2,'fixed':'d1 a'}
conf['params']['ffpion']['sb1 b'] ={'value': 1.83178e+00, 'min': 0 , 'max': 10,'fixed':'d1 b'}
conf['params']['ffpion']['cb1 N'] ={'value': 1.84550e-01, 'min': 0.0 , 'max': 1,'fixed':'c1 N'}
conf['params']['ffpion']['cb1 a'] ={'value': -5.05798e-02, 'min': -1.8 , 'max': 2,'fixed':'c1 a'}
conf['params']['ffpion']['cb1 b'] ={'value': 3.19952e+00, 'min': 0 , 'max': 10,'fixed':'c1 b'}
conf['params']['ffpion']['bb1 N'] ={'value': 3.74125e-01, 'min': 0.0 , 'max': 1,'fixed':'b1 N'}
conf['params']['ffpion']['bb1 a'] ={'value': -1.59541e+00, 'min': -1.8 , 'max': 2,'fixed':'b1 a'}
conf['params']['ffpion']['bb1 b'] ={'value': 4.50102e+00, 'min': 0 , 'max': 10,'fixed':'b1 b'}
kaon fragamentation functions¶
conf['params']['ffkaon']={}
conf['ffkaon parametrization']=0
conf['params']['ffkaon']['g1 N'] ={'value': 2.33320e-01, 'min': 0.0 , 'max': 1,'fixed':False}
conf['params']['ffkaon']['g1 a'] ={'value': 1.48737e+00, 'min': -1.8 , 'max': 2,'fixed':False}
conf['params']['ffkaon']['g1 b'] ={'value': 9.62755e+00, 'min': 0 , 'max': 10,'fixed':False}
conf['params']['ffkaon']['u1 N'] ={'value': 4.03672e-02, 'min': 0.0 , 'max': 1,'fixed':False}
conf['params']['ffkaon']['u1 a'] ={'value': 1.26356e+00, 'min': -1.8 , 'max': 2,'fixed':False}
conf['params']['ffkaon']['u1 b'] ={'value': 1.62596e+00, 'min': 0 , 'max': 10,'fixed':False}
conf['params']['ffkaon']['d1 N'] ={'value': 2.51671e-01, 'min': 0.0 , 'max': 1,'fixed':False}
conf['params']['ffkaon']['d1 a'] ={'value': -1.43444e+00, 'min': -1.8 , 'max': 2,'fixed':False}
conf['params']['ffkaon']['d1 b'] ={'value': 1.65143e+00, 'min': 0 , 'max': 10,'fixed':False}
conf['params']['ffkaon']['s1 N'] ={'value': 2.51671e-01, 'min': 0.0 , 'max': 1,'fixed':'d1 N'}
conf['params']['ffkaon']['s1 a'] ={'value': -1.43444e+00, 'min': -1.8 , 'max': 2,'fixed':'d1 a'}
conf['params']['ffkaon']['s1 b'] ={'value': 1.65143e+00, 'min': 0 , 'max': 10,'fixed':'d1 b'}
conf['params']['ffkaon']['c1 N'] ={'value': 7.76923e-01, 'min': 0.0 , 'max': 1,'fixed':False}
conf['params']['ffkaon']['c1 a'] ={'value': -1.80000e+00, 'min': -1.8 , 'max': 2,'fixed':False}
conf['params']['ffkaon']['c1 b'] ={'value': 2.50452e+00, 'min': 0 , 'max': 10,'fixed':False}
conf['params']['ffkaon']['b1 N'] ={'value': 5.66971e-01, 'min': 0.0 , 'max': 1,'fixed':False}
conf['params']['ffkaon']['b1 a'] ={'value': -1.80000e+00, 'min': -1.8 , 'max': 2,'fixed':False}
conf['params']['ffkaon']['b1 b'] ={'value': 3.41727e+00, 'min': 0 , 'max': 10,'fixed':False}
conf['params']['ffkaon']['ub1 N'] ={'value': 4.03672e-02, 'min': 0.0 , 'max': 1,'fixed':'d1 N'}
conf['params']['ffkaon']['ub1 a'] ={'value': 2.26356e+00, 'min': -1.8 , 'max': 2,'fixed':'d1 a'}
conf['params']['ffkaon']['ub1 b'] ={'value': 1.62596e+00, 'min': 0 , 'max': 10,'fixed':'d1 b'}
conf['params']['ffkaon']['db1 N'] ={'value': 2.51671e-01, 'min': 0.0 , 'max': 1,'fixed':'d1 N'}
conf['params']['ffkaon']['db1 a'] ={'value': -1.43444e+00, 'min': -1.8 , 'max': 2,'fixed':'d1 a'}
conf['params']['ffkaon']['db1 b'] ={'value': 1.65143e+00, 'min': 0 , 'max': 10,'fixed':'d1 b'}
conf['params']['ffkaon']['sb1 N'] ={'value': 2.51671e-01, 'min': 0.0 , 'max': 1,'fixed':False}
conf['params']['ffkaon']['sb1 a'] ={'value': -1.43444e+00, 'min': -1.8 , 'max': 2,'fixed':False}
conf['params']['ffkaon']['sb1 b'] ={'value': 1.65143e+00, 'min': 0 , 'max': 10,'fixed':False}
conf['params']['ffkaon']['cb1 N'] ={'value': 7.76923e-01, 'min': 0.0 , 'max': 1,'fixed':'c1 N'}
conf['params']['ffkaon']['cb1 a'] ={'value': -1.80000e+00, 'min': -1.8 , 'max': 2,'fixed':'c1 a'}
conf['params']['ffkaon']['cb1 b'] ={'value': 2.50452e+00, 'min': 0 , 'max': 10,'fixed':'c1 b'}
conf['params']['ffkaon']['bb1 N'] ={'value': 5.66971e-01, 'min': 0.0 , 'max': 1,'fixed':'b1 N'}
conf['params']['ffkaon']['bb1 a'] ={'value': -1.80000e+00, 'min': -1.8 , 'max': 2,'fixed':'b1 a'}
conf['params']['ffkaon']['bb1 b'] ={'value': 3.41727e+00, 'min': 0 , 'max': 10,'fixed':'b1 b'}
Steps¶
We provide here examples to define a step in the multistep strategy.
The steps are currentlty interpretable by the script maxlike.py.
Some explanation are in order:
dep: a list of the steps that are used as priors.
active distributions: a list of distributions like pdfs or ffs that are needed
passive distributions: remove from the active set the distributions whose parameters are fixed according from the priors specified at
dep
#--steps
conf['steps'] = {}
#--idis and no hera
conf['steps'][1]={}
conf['steps'][1]['dep']=[]
conf['steps'][1]['active distributions']=['pdf']
conf['steps'][1]['datasets']={}
conf['steps'][1]['datasets']['idis']=[]
conf['steps'][1]['datasets']['idis'].append(10010) # proton | F2 | SLAC
conf['steps'][1]['datasets']['idis'].append(10011) # deuteron | F2 | SLAC
conf['steps'][1]['datasets']['idis'].append(10016) # proton | F2 | BCDMS
conf['steps'][1]['datasets']['idis'].append(10017) # deuteron | F2 | BCDMS
conf['steps'][1]['datasets']['idis'].append(10020) # proton | F2 | NMC
conf['steps'][1]['datasets']['idis'].append(10021) # d/p | F2d/F2p | NMC
#--idis + hera
conf['steps'][2]={}
conf['steps'][2]['dep']=[1]
conf['steps'][2]['active distributions']=['pdf']
conf['steps'][2]['datasets']={}
conf['steps'][2]['datasets']['idis']=[]
conf['steps'][2]['datasets']['idis'].append(10010) # proton | F2 | SLAC
conf['steps'][2]['datasets']['idis'].append(10011) # deuteron | F2 | SLAC
conf['steps'][2]['datasets']['idis'].append(10016) # proton | F2 | BCDMS
conf['steps'][2]['datasets']['idis'].append(10017) # deuteron | F2 | BCDMS
conf['steps'][2]['datasets']['idis'].append(10020) # proton | F2 | NMC
conf['steps'][2]['datasets']['idis'].append(10021) # d/p | F2d/F2p | NMC
conf['steps'][2]['datasets']['idis'].append(10026) # proton | sigma red | HERA II NC e+ (1)
conf['steps'][2]['datasets']['idis'].append(10027) # proton | sigma red | HERA II NC e+ (2)
conf['steps'][2]['datasets']['idis'].append(10028) # proton | sigma red | HERA II NC e+ (3)
conf['steps'][2]['datasets']['idis'].append(10029) # proton | sigma red | HERA II NC e+ (4)
conf['steps'][2]['datasets']['idis'].append(10030) # proton | sigma red | HERA II NC e-
conf['steps'][2]['datasets']['idis'].append(10031) # proton | sigma red | HERA II CC e+
conf['steps'][2]['datasets']['idis'].append(10032) # proton | sigma red | HERA II NC e-
#--idis + dy
conf['steps'][3]={}
conf['steps'][3]['dep']=[2]
conf['steps'][3]['active distributions']=['pdf']
conf['steps'][3]['datasets']={}
conf['steps'][3]['datasets']['idis']=[]
conf['steps'][3]['datasets']['idis'].append(10010) # proton | F2 | SLAC
conf['steps'][3]['datasets']['idis'].append(10011) # deuteron | F2 | SLAC
conf['steps'][3]['datasets']['idis'].append(10016) # proton | F2 | BCDMS
conf['steps'][3]['datasets']['idis'].append(10017) # deuteron | F2 | BCDMS
conf['steps'][3]['datasets']['idis'].append(10020) # proton | F2 | NMC
conf['steps'][3]['datasets']['idis'].append(10021) # d/p | F2d/F2p | NMC
conf['steps'][3]['datasets']['idis'].append(10026) # proton | sigma red | HERA II NC e+ (1)
conf['steps'][3]['datasets']['idis'].append(10027) # proton | sigma red | HERA II NC e+ (2)
conf['steps'][3]['datasets']['idis'].append(10028) # proton | sigma red | HERA II NC e+ (3)
conf['steps'][3]['datasets']['idis'].append(10029) # proton | sigma red | HERA II NC e+ (4)
conf['steps'][3]['datasets']['idis'].append(10030) # proton | sigma red | HERA II NC e-
conf['steps'][3]['datasets']['idis'].append(10031) # proton | sigma red | HERA II CC e+
conf['steps'][3]['datasets']['idis'].append(10032) # proton | sigma red | HERA II NC e-
conf['steps'][3]['datasets']['dy']=[]
conf['steps'][3]['datasets']['dy'].append(10001)
conf['steps'][3]['datasets']['dy'].append(10002)