std::complex<double> raremodel::vectorNonres()
{
  double mmu = pdg::muon_M;
  double mb = pdg::bquark_mass;
  double ms = pdg::squark_mass;
  double mK = pdg::kaon_M;
  double mB = pdg::Bplus_M;


  double betaplus = sqrt(1-4*(mmu*mmu)/((*m)*(*m)));
  double betaplus_function = sqrt(betaplus*(1-(1.0/3.0)*betaplus*betaplus));

  double triangle_d = triangle();
  triangle_d=triangle_d*triangle_d*triangle_d;
  triangle_d=sqrt(triangle_d);
  triangle_d=sqrt(triangle_d);

  double phasespace = triangle_d*sqrt(2*(*m));


  double C9component = fplusFactorLattice()*(*C9eff);

  double C7component = 2*((mb+ms)/(mK+mB))*ftransverseFactorLattice()*(*C7eff);

  std::complex<double> nonresamp = std::complex<double> (phasespace*betaplus_function*(C9component+C7component), 0. );
  return nonresamp;