import numpy as np
def c8_anapole(mchi,sigma_ref):
      hbarc2=0.389e-27
      gp=5.585694713
      gn=-3.82608545
      e_tau=np.array([1,1])
      m_N=0.931
      m_chi_N=mchi*m_N/(mchi+m_N)
      eg_over_lambda2=np.sqrt(4*np.pi*sigma_ref/(2*m_chi_N**2*hbarc2))
      return 2*eg_over_lambda2*e_tau


def c9_anapole(mchi,sigma_ref):
      hbarc2=0.389e-27
      gp=5.585694713
      gn=-3.82608545
      g_tau=np.array([gp+gn,gp-gn])
      m_N=0.931
      m_chi_N=mchi*m_N/(mchi+m_N)
      eg_over_lambda2=np.sqrt(4*np.pi*sigma_ref/(2*m_chi_N**2*hbarc2))
      return -eg_over_lambda2*g_tau