//----------------------------------------------------------------------
// Conversion Rate as a function of Eta - overall!
void z5() {
  std::cout << "Conversion rate as a function of Eta, summed over all R" << std::endl;
  // Canvas
  c5 = new TCanvas("c5","Conversion rate as a function of Eta, summed over all R",800,600);
  c5->SetBorderMode(0);

  TH1F *photCnvRate;
  ntuple->Draw("abs(TruthPhotonEta) >> ph1a(25,0,4)","abs(TruthPhotonEta)<10");
  TH1F *ph1a = (TH1F*)gDirectory->Get("ph1a");
  ntuple->Draw("abs(eConvVerEta) >> ph1b(25,0,4)"," NConversions>0");
  TH1F *ph1b = (TH1F*)gDirectory->Get("ph1b");

  c5->Clear();
  c5->cd();
  photCnvRate = new TH1F("photCnvRate","Photon Conversion Rate",25,0,4);

  for(int i=1; i<25; i++){
    float totalentries = ph1a->GetBinContent(i);
    float conventries = ph1b->GetBinContent(i);
    if(totalentries!=0)
      photCnvRate->SetBinContent(i,conventries/totalentries);
    else
      photCnvRate->SetBinContent(i,0);
  }

  photCnvRate->SetTitle("Conv. rate");

  //    photCnvRate->GetXaxis()->SetRangeUser(0,3);
  //    photCnvRate->GetYaxis()->SetRangeUser(.4,.8);
  photCnvRate->GetXaxis()->SetTitle("|#eta|");
  photCnvRate->GetYaxis()->SetTitle("rate");
  photCnvRate->SetStats(kFALSE);
  photCnvRate->SetFillColor(40);
  photCnvRate->Draw();
    
  c5->Modified();
  c5->Update();
  printHisto(c5,"cnvRate_all_R");
  return;
}
//----------------------------------------------------------------------