#include <iostream.h>
#include <math.h>
#include "sigma_zz.h"

void main()
{

  float mz = 91.2;
  float Ecm_min = 150.;
  float Ecm_max = 350.;
  int num_points = 101;

  for (int i=0; i<num_points; i++){
    float Ecm = (Ecm_max - Ecm_min)*float(i)/float(num_points-1) + Ecm_min;
    float sigma = sigma_zz(Ecm, mz);
    cout << Ecm << "  "  << sigma << endl;
  }

}