void gausOut(int iShow){

  TF1 *f0 = new TF1("f0","gaus");

  int nBins = 100;
  TH1D *hGaus = new TH1D("hGaus","",nBins,-3.0,3.0);
  TH1D *hGausSqrt = new TH1D("hGausSqrt","",nBins,-3.0,3.0);

  f0->SetParameter(0,1.0);
  f0->SetParameter(1,0.0);
  f0->SetParameter(2,1.0);

  hGaus->Eval(f0);

  double intVal = hGaus->Integral();

  cout<<"intVal = "<<intVal<<endl;

  hGaus->Scale(1.0/intVal);

  cout<<"intVal New = "<<hGaus->Integral()<<endl;

  hGaus->Draw();

  
  for (int iBin = 1; iBin <= nBins; iBin++) {
    cout<<hGaus->GetBinContent(iBin)<<", ";
    if (iBin%10 == 0) cout<<endl;
    hGausSqrt->SetBinContent(iBin,sqrt(hGaus->GetBinContent(iBin)));
  }
  cout<<"Centers ->"<<endl;
  for (int iBin = 1; iBin <= nBins; iBin++) {
    cout<<hGaus->GetBinCenter(iBin)<<", ";
    if (iBin%10 == 0) cout<<endl;
  }
  hGaus->Draw();

}