Two-step process (phi in medium -> K K -> K K_scatter)

double f_my_function(double *x, double * par)
{
    if (x[0] > 0.8) return 1;
    return 0;
}

fermi_rescattering() {