tprak

Theoretical Physics Practical Training
git clone git://popovic.xyz/tprak.git
Log | Files | Refs

fit_single.py (742B)

      1 #!/usr/bin/env python3.9
      2 
      3 import numpy as np
      4 import matplotlib.pyplot as plt
      5 
      6 sig_p = lambda x: np.sqrt(1 - 4*m_p**2/x)
      7 g_s = lambda s, m_q, g_q: g_q*s/m_q**2 * (sig_p(s)/sig_p(m_q**2))**2 * np.heaviside(s, 4*m_p**2)
      8 
      9 def model(s, m_q, g_q, m_w, g_w, e_w, a, b):
     10     part1 = (m_q)**4/((m_q**2 - s)** + m_q**2*g_s(s, m_q, g_q)**2)
     11     part2 = 1 + (e_w * 2*s * (m_w**2 - s))/((m_w**2 - s)**2 + m_w**2*g_w**2)
     12     part3 = (1 + a*s + b*s)**2
     13     return part1 * part2 * part3
     14 
     15 def main():
     16     text = open('./data/KLOE-VFF.txt').read()
     17     data = np.loadtxt('./data/KLOE-VFF.txt')
     18     cov_syst = np.loadtxt('./data/KLOE-RelSystCov.txt')
     19     cov_stat = np.loadtxt('./data/KLOE-StatCov.txt')
     20     s = data[:,0]
     21     F2 = data[:, 1]
     22 
     23     m_p = 0.13957
     24