for m in ml:
gv_mp = data_gv[x_list['mp_%s' %str(m)]]
gv_pp = data_gv[x_list['pp_%s' %str(m)]]
fit = decay.fit_mres_bs(psql,decay_params,m,gv_mp,gv_pp)['mres_fit']
result['mres_%s' %str(m)] = fit
for s in ms:
gv_mp = data_gv[x_list['mp_%s' %str(s)]]
gv_pp = data_gv[x_list['pp_%s' %str(s)]]
fit = decay.fit_mres_bs(psql,decay_params,s,gv_mp,gv_pp)['mres_fit']
result['mres_%s' %str(s)] = fit
# mesons
for m in ml:
decay_params['decay_ward_fit']['ml'] = m
gv_SS = data_gv[x_list['ss_%s_%s' %(str(m),str(m))]]
gv_PS = data_gv[x_list['ps_%s_%s' %(str(m),str(m))]]
fit = decay.fit_decay_bs(psql,decay_params,'pion', gv_SS, gv_PS)['meson_fit']
result['meson_%s_%s' %(str(m),str(m))] = fit
for s in ms:
decay_params['decay_ward_fit']['ms'] = s
gv_SS = data_gv[x_list['ss_%s_%s' %(str(m),str(s))]]
gv_PS = data_gv[x_list['ps_%s_%s' %(str(m),str(s))]]
fit = decay.fit_decay_bs(psql,decay_params,'kaon', gv_SS, gv_PS)['meson_fit']
result['meson_%s_%s' %(str(m),str(s))] = fit
# make decay constants
decay = dict()
for m in ml:
Z0_p = result['meson_%s_%s' %(str(m),str(m))].p['Z0_p']
E0 = result['meson_%s_%s' %(str(m),str(m))].p['E0']
mres = result['mres_%s' %str(m)].p['mres']
decay['%s_%s' %(str(m),str(m))] = Z0_p*np.sqrt(2.)*(2.*float(m)+2.*mres)/E0**(3./2.)
for s in ms:
gvss_pion = gvboot0[x_ss_pion]
gvps_pion = gvboot0[x_ps_pion]
gvss_kaon = gvboot0[x_ss_kaon]
gvps_kaon = gvboot0[x_ps_kaon]
# fit #
# mres
mresl = decay.fit_mres_bs(psql, decay_params,
decay_params['decay_ward_fit']['ml'], gvmpl,
gvppl)
mress = decay.fit_mres_bs(psql, decay_params,
decay_params['decay_ward_fit']['ms'], gvmps,
gvpps)
# meson
pionfit = decay.fit_decay_bs(psql,
decay_params,
'pion',
gvss_pion,
gvps_pion,
flow=True)
kaonfit = decay.fit_decay_bs(psql,
decay_params,
'kaon',
gvss_kaon,
gvps_kaon,
flow=True)
# make 5D decay constant
fkaon = decay.decay_constant(decay_params, kaonfit['meson_fit'].p['Z0_p'],
kaonfit['meson_fit'].p['E0'],
mresl['mres_fit'].p['mres'],
mress['mres_fit'].p['mres']) / np.sqrt(2)
fpion = decay.decay_constant(decay_params, pionfit['meson_fit'].p['Z0_p'],
pionfit['meson_fit'].p['E0'],
print "len gAboot0:", np.shape(gAboot0) # concatenate boot0 = np.concatenate((mp,pp,ss_pion,ps_pion,gAboot0),axis=1) # correlated covariance gvboot0 = c51.make_gvars(boot0) # split to different project subsets gvmp = gvboot0[x_mp] gvpp = gvboot0[x_pp] gvss_pion = gvboot0[x_ss_pion] gvps_pion = gvboot0[x_ps_pion] gvgAboot0 = gvboot0[x_gAboot0] # fit # # mres mresl = decay.fit_mres_bs(psql,decay_params,decay_params['decay_ward_fit']['ml'],gvmp,gvpp) # meson pifit = decay.fit_decay_bs(psql,decay_params,'pion',gvss_pion,gvps_pion) # mN and gA #gAfit = gA.fit_gA(psql,gA_params,gvgAboot0) gAfit = gA.fit_proton(psql,gA_params,gvgAboot0) # chipt parameters priors = gv.BufferDict() priors['mpi'] = pifit['meson_fit'].p['E0'] priors['fpi'] = decay.decay_constant(decay_params,pifit['meson_fit'].p['Z0_p'],pifit['meson_fit'].p['E0'],mresl['mres_fit'].p['mres']) #priors['mN'] = gAfit['gA_fit'].p['E0'] priors['mN'] = gAfit['nucleon_fit'].p['E0'] print priors #print gv.evalcorr([pifit['meson_fit'].p['E0'],decay.decay_constant(decay_params,pifit['meson_fit'].p['Z0_p'],pifit['meson_fit'].p['E0'],mresl['mres_fit'].p['mres']),gAfit['gA_fit'].p['E0'],gAfit['gA_fit'].p['gA00']]) print priors['mN']/priors['fpi'] data = gv.BufferDict() data['gA'] = gAfit['gA_fit'].p['gA00'] # write output
gvss_pion = gvboot0[x_ss_pion]
gvps_pion = gvboot0[x_ps_pion]
gvss_kaon = gvboot0[x_ss_kaon]
gvps_kaon = gvboot0[x_ps_kaon]
gvss_etas = gvboot0[x_ss_etas]
gvps_etas = gvboot0[x_ps_etas]
gvphi_jj_5 = gvboot0[x_phi_jj_5]
gvphi_jr_5 = gvboot0[x_phi_jr_5]
gvphi_rr_5 = gvboot0[x_phi_rr_5]
gvphi_ju = gvboot0[x_phi_ju]
gvphi_js = gvboot0[x_phi_js]
gvphi_ru = gvboot0[x_phi_ru]
gvphi_rs = gvboot0[x_phi_rs]
# fit #
# meson
pifit = decay.fit_decay_bs(psql, decay_params, 'pion', gvss_pion,
gvps_pion)
kafit = decay.fit_decay_bs(psql, decay_params, 'kaon', gvss_kaon,
gvps_kaon)
etfit = decay.fit_decay_bs(psql, decay_params, 'etas', gvss_etas,
gvps_etas)
# hisq
jjfit = hisq.fit_hisq_meson(psql, hisq_params, 'phi_jj_5', gvphi_jj_5)
jrfit = hisq.fit_hisq_meson(psql, hisq_params, 'phi_jr_5', gvphi_jr_5)
rrfit = hisq.fit_hisq_meson(psql, hisq_params, 'phi_rr_5', gvphi_rr_5)
# mixed
jufit = mix.fit_mixed_meson(psql, mixed_params, 'phi_ju', gvphi_ju)
jsfit = mix.fit_mixed_meson(psql, mixed_params, 'phi_js', gvphi_js)
rufit = mix.fit_mixed_meson(psql, mixed_params, 'phi_ru', gvphi_ru)
rsfit = mix.fit_mixed_meson(psql, mixed_params, 'phi_rs', gvphi_rs)
# collect output
result = gv.BufferDict()
# concatenate
boot0 = np.concatenate((mpl, ppl, ss_pion, ps_pion), axis=1)
# correlated covariance
gvboot0 = c51.make_gvars(boot0)
# split for projects
gvmpl = gvboot0[x_mpl]
gvppl = gvboot0[x_ppl]
gvss_pion = gvboot0[x_ss_pion]
gvps_pion = gvboot0[x_ps_pion]
# fit #
# mres
mresl = decay.fit_mres_bs(psql, decay_params,
decay_params['decay_ward_fit']['ml'], gvmpl,
gvppl)
# meson
pionfit = decay.fit_decay_bs(psql, decay_params, 'pion', gvss_pion,
gvps_pion)
# make 5D decay constant
fpion = decay.decay_constant(decay_params, pionfit['meson_fit'].p['Z0_p'],
pionfit['meson_fit'].p['E0'],
mresl['mres_fit'].p['mres'])
mpion = pionfit['meson_fit'].p['E0']
# result
#print mresl['mres_fit'].p
#print mress['mres_fit'].p
#print pionfit['meson_fit'].p
#print kaonfit['meson_fit'].p
#print gv.evalcorr([mresl['mres_fit'].p['mres'], mress['mres_fit'].p['mres'], pionfit['meson_fit'].p['E0'], kaonfit['meson_fit'].p['E0']])
print "fpi:", fpion
#print gv.evalcorr([fkaon,fpion])
# chipt parameters
priors = gv.BufferDict()
for m in ml:
gv_mp = data_gv[x_list['mp_%s' %str(m)]]
gv_pp = data_gv[x_list['pp_%s' %str(m)]]
fit = decay.fit_mres_bs(psql,decay_params,m,gv_mp,gv_pp)['mres_fit']
result['mres_%s' %str(m)] = fit
for s in ms:
gv_mp = data_gv[x_list['mp_%s' %str(s)]]
gv_pp = data_gv[x_list['pp_%s' %str(s)]]
fit = decay.fit_mres_bs(psql,decay_params,s,gv_mp,gv_pp)['mres_fit']
result['mres_%s' %str(s)] = fit
# mesons
for m in ml:
decay_params['decay_ward_fit']['ml'] = m
gv_SS = data_gv[x_list['ss_%s_%s' %(str(m),str(m))]]
gv_PS = data_gv[x_list['ps_%s_%s' %(str(m),str(m))]]
fit = decay.fit_decay_bs(psql,decay_params,'pion', gv_SS, gv_PS)['meson_fit']
result['meson_%s_%s' %(str(m),str(m))] = fit
for s in ms:
decay_params['decay_ward_fit']['ms'] = s
gv_SS = data_gv[x_list['ss_%s_%s' %(str(m),str(s))]]
gv_PS = data_gv[x_list['ps_%s_%s' %(str(m),str(s))]]
fit = decay.fit_decay_bs(psql,decay_params,'kaon', gv_SS, gv_PS)['meson_fit']
result['meson_%s_%s' %(str(m),str(s))] = fit
# make decay constants
decay = dict()
for m in ml:
Z0_p = result['meson_%s_%s' %(str(m),str(m))].p['Z0_p']
E0 = result['meson_%s_%s' %(str(m),str(m))].p['E0']
mres = result['mres_%s' %str(m)].p['mres']
decay['%s_%s' %(str(m),str(m))] = Z0_p*np.sqrt(2.)*(2.*float(m)+2.*mres)/E0**(3./2.)
for s in ms: