const double GF=1.16639e-5;
for(int i=0; i<90; i++) {
//------------------------------------------------------
  double mH =100.+10.*i;
  double tau = mH*mH/(ecm*ecm);
  double c1 = GF*pow(alphaS(mH),2)/(288.*1.414214*pi);
  double crossGeVInv2 = c1*I2(mtop/mH)*tau*partonLuminosity(tau);
  double crossSection = crossGeVInv2*197*197*1.e-6*1.e-2*1.e12;//GeV,barn,pico factor
//   ****** THIS PART MUST BE CODED ******
  

//------------------------------------------------------
    std::cout <<mH<<"     "<<crossSection<<std::endl;
}