#include <iostream>
#include <fstream>
#include <cmath>

const int NPoints = 6;//16+8;


struct TDataPoint{
  int index;
  float data;
  TDataPoint(): index(0), data(0) {};
  TDataPoint( TDataPoint &a){ index = a.index; data = a.data; };
  TDataPoint(int _index, float _data): index(_index), data(_data) {};
};


// TDataPoint data[NPoints] =
// {
//   TDataPoint(1, ),
//   TDataPoint(2, ),
//   TDataPoint(3, ),
//   TDataPoint(4, ),
//   TDataPoint(5, ),
//   TDataPoint(6, ),
//   TDataPoint(7, ),
//   TDataPoint(8, ),
//   TDataPoint(9, ),
//   TDataPoint(10, ),
//   TDataPoint(11, ),
//   TDataPoint(12, ),
//   TDataPoint(13, ),
//   TDataPoint(14, ),
//   TDataPoint(15, ),
//   TDataPoint(16, )
// };

TDataPoint data[NPoints] = // 1 iter
{

  TDataPoint(1, .016),
  TDataPoint(2, .030),
  TDataPoint(4, .064),
  TDataPoint(5, .081),
  TDataPoint(10, .184),
  TDataPoint(20, .430)
};

// TDataPoint data[NPoints] = // 3 iter
// {
// 
//   TDataPoint(1, .018),
//   TDataPoint(1, .048),
//   TDataPoint(2, .089),
//   TDataPoint(3, .109),
//   TDataPoint(4, .130),
//   TDataPoint(5, .150),
//   TDataPoint(6, .189),
//   TDataPoint(7, .210),
//   TDataPoint(8, .234),
//   TDataPoint(9, .258),
//   TDataPoint(10, .299),
//   TDataPoint(11, .347),
//   TDataPoint(12, .377),
//   TDataPoint(13, .398),
//   TDataPoint(14, .411),
//   TDataPoint(15, .467),
//   TDataPoint(16, .482),
//   TDataPoint(17, .505),
//   TDataPoint(18, .542)/*,
//   TDataPoint(19, ),
//   TDataPoint(20, )*/
// };

void MakeHisto()
{
  int maxIndex = data[NPoints-1].index;
  TH1F *h = new TH1F("Histo","Histo",maxIndex+1,-0.5,maxIndex+0.5);

  for (int i=0; i<NPoints; i+=1){
    h->SetBinContent(data[i].index+1,1./data[i].data*data[i].index);
    h->SetBinError(data[i].index+1,1.); // need for draw markers!
  }
  h->GetXaxis()->SetTitle("Number of packed mbias events");
  h->GetYaxis()->SetTitle("MBias Events / s");
  h->GetYaxis()->SetTitleOffset(1.2);
//   h->GetYaxis()->SetRangeUser(0,250);
  h->GetYaxis()->SetTitleFont(22);
  h->GetXaxis()->SetTitleFont(22);
  h->GetYaxis()->SetLabelFont(22);
  h->GetXaxis()->SetLabelFont(22);
  h->SetMarkerColor(2);
  h->SetLineColor(0);
  h->SetLineWidth(0);
  h->SetMarkerStyle(21);
  h->SetMarkerSize(1.2);
  h->SetMarkerColor(2);

  gStyle->SetCanvasColor(10);
  gStyle->SetFrameFillColor(10);
  gStyle->SetTitleColor(0);
  gStyle->SetHistFillColor(0);
  gStyle->SetOptStat(00000000);
  gStyle->SetPadGridX(1);
  gStyle->SetPadGridY(1);
  gStyle->SetOptTitle(0);
  gStyle->SetTextFont(22);

  TCanvas *c1 = new TCanvas("ITBB", "ITBB", 700, 700);
//   c1->Divide(2);
  c1->Draw();
//   c1->cd(1);
//   hc->Draw();
//   c1->cd(2);
  h->Draw();
//   TFile file(outputFile,"RECREATE");
//  c1->Write();
//   h->Write();
}