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

void main()
{
  float mz = 91.2;
  float Ecm_min = 150.;
  float Ecm_max = 350.;
  int num_points = 101;
  float mh;
  cout << "enter mh" << endl;
  cin >> mh;

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

}