#include <iostream>
#include <sstream>
#include <iomanip>
#include <vector>

void wager(std::string hname)
{

  std::vector<double> V, eV;

  if (hname=="j_bjet_Pt_binned") 
    {
      std::cout<<" pippo"<<std::endl;
      V.push_back(0.974);  eV.push_back(0.004);
      V.push_back(0.974);  eV.push_back(0.004);
      V.push_back(0.974);  eV.push_back(0.005);
      V.push_back(0.972);  eV.push_back(0.008);
      V.push_back(0.975);  eV.push_back(0.012);
      V.push_back(0.970);  eV.push_back(0.033);
    }
  else if (hname=="j_bjet_Y_binned") 
    {
      V.push_back(0.974);  eV.push_back(0.007);
      V.push_back(0.975);  eV.push_back(0.007);
      V.push_back(0.974);  eV.push_back(0.007);
      V.push_back(0.975);  eV.push_back(0.007);
      V.push_back(0.974);  eV.push_back(0.005);
      V.push_back(0.974);  eV.push_back(0.006);
      V.push_back(0.973);  eV.push_back(0.007);
      V.push_back(0.972);  eV.push_back(0.009);
    }
  else if (hname=="j_bjet_ZPt_binned") 
    {
      V.push_back(0.975);  eV.push_back(0.007);
      V.push_back(0.978);  eV.push_back(0.006);
      V.push_back(0.975);  eV.push_back(0.006);
      V.push_back(0.975);  eV.push_back(0.004);
      V.push_back(0.970);  eV.push_back(0.006);
      V.push_back(0.967);  eV.push_back(0.007);
      V.push_back(0.972);  eV.push_back(0.009);
      V.push_back(0.965);  eV.push_back(0.022);
    }
  else if (hname=="j_bjet_dPhiF20_binned") 
    {
      V.push_back(0.974);  eV.push_back(0.020);
      V.push_back(0.969);  eV.push_back(0.019);
      V.push_back(0.978);  eV.push_back(0.015);
      V.push_back(0.973);  eV.push_back(0.012);
      V.push_back(0.979);  eV.push_back(0.010);
      V.push_back(0.976);  eV.push_back(0.007);
      V.push_back(0.973);  eV.push_back(0.007);
      V.push_back(0.971);  eV.push_back(0.008);
    }
  else if (hname=="j_bjet_dR620_binned") 
    {
      V.push_back(0.973);  eV.push_back(0.015);
      V.push_back(0.974);  eV.push_back(0.012);
      V.push_back(0.970);  eV.push_back(0.009);
      V.push_back(0.975);  eV.push_back(0.007);
      V.push_back(0.976);  eV.push_back(0.004);
      V.push_back(0.972);  eV.push_back(0.004);
      V.push_back(0.973);  eV.push_back(0.010);
      V.push_back(0.982);  eV.push_back(0.022);
      V.push_back(0.980);  eV.push_back(0.051);
    }		 
  else  
    {
      std::cout<<" error "<<std::endl;
      return;
    }


       std::cout<<" here we are "<<std::endl;
      double average = 0.;
      double eaverage = 0.;
      double waverage = 0.;
      double ewaverage = 0.;
      double sumw = 0.;
      unsigned int np = 0;
      
      std::vector<double>::const_iterator it =V.begin();
      std::vector<double>::const_iterator ite=V.end();
      std::vector<double>::const_iterator itErr=eV.begin();
      std::cout<<" vsize = "<<V.size()<<std::endl;
      
      for (;it!=ite; ++it)
	{
	  ++np;
	  average += (*it);
	  double w = (*itErr);
	  //average = V[np];
	  //double w = eV[np];
	  w = 1./(w*w);

	  waverage += w*(*it);
	  sumw += w;

	  std::cout<<" component n. "<<np<<" --- adding "<<(*it)<<"+/- "<<(*itErr)<<" ====> average  sum = "<< average<<"  average = "<< average/double(np)<<std::endl;
	  std::cout<<" component n. "<<np<<" --- adding "<<(*it)<<"+/- "<<(*itErr)<<" ====> waverage sum = "<<waverage<<" waverage = "<<waverage/sumw<<std::endl;
	  ++itErr;
	}

      std::cout<<" *** "<<std::endl;
      average  = average/double(np);
      waverage = waverage/sumw ;
      std::cout<<" Final   average = "<< average<<std::endl;
      std::cout<<" Final  waverage = "<<waverage<<" +\- "<<1./sqrt(sumw)<<std::endl;
      std::cout<<" *** "<<std::endl;

      double residual=0.; 
      double wresidual=0.;
      std::vector<double>::const_iterator itn =V.begin();
      std::vector<double>::const_iterator itne =V.end();
      std::vector<double>::const_iterator itErrn=eV.begin();
      np=0;
      for (;itn!=itne; ++itn)
	{
	  ++np;
	  residual  = ((*itn)-average)/(*itErrn);
	  wresidual = ((*itn)-waverage)/(*itErrn);

	  std::cout<<" component n. "<<np<<" is "<<(*itn)<<" --- residual "<< residual<<std::endl;
	  std::cout<<" component n. "<<np<<" is "<<(*itn)<<" --- residual "<<wresidual<<std::endl;
	  ++itErrn;
	}
}
////// results 
////// bjet pT : Final  waverage = 0.973852 +- 0.00230334
////// bjet |y|: Final  waverage = 0.973852 +- 0.00230334
////// Z pT    : Final  waverage = 0.973625 +- 0.00223016
////// dPhi    : Final  waverage = 0.974247 +- 0.00346896
////// dR      : Final  waverage = 0.973865 +- 0.00234666