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