#include "TString.h"
#include "TH1.h"
#include "TAxis.h"
#include "TStyle.h"
#include "TCanvas.h"

#include <cstdlib>
#include <string>
#include <iostream>
#include <fstream>

void read_and_draw_vdepl_Perugia2015_fl1e15(const char* filename = "vdepl_Perugia2015_fl1e15.txt") {
  std::fstream file;
  file.open(filename);
  if (! file.good() ){
    std::cout << filename << " is not good: " << file.rdstate() << std::endl;
    exit(1);
  }
  
  int nmax = 10;
  TString* models = new TString[nmax];
  double* vdepls = new double[nmax];
  double* evdepls = new double[nmax];
  double vdepl,evdepl;
  std::string model;
  int i = 0;
  while(file >> model >> vdepl >> evdepl) {
    vdepls[i]=vdepl;
    evdepls[i] = evdepl;
    models[i] = model.c_str();
    //std::cout << model << std::endl;
    //std::cout << (models[i]).Data() << std::endl;
    //std::cout << models[i] << std::endl;
    i++;
    if (i >= nmax) break;
  }
  
  file.close();
  
  //std::cout << i << " lines" << std::endl;
  TH1D* hVdepl = new TH1D("hVdepl",";;V_{depl} [V]",i,0,i);
  for (int j = 0; j<i; j++) {
    hVdepl->SetBinContent(j+1,vdepls[j]);
    hVdepl->SetBinError(j+1,evdepls[j]);
    //std::cout << models[j] << std::endl;
    //std::cout << (models[j]).Data() << std::endl;
    hVdepl->GetXaxis()->SetBinLabel(j+1,(models[j]).Data());
  }
  
  gStyle->SetOptStat(0);
  
  TCanvas *c1 = new TCanvas("c1","",800,600);
  // c1->cd();
  hVdepl->Draw("e");
  
  delete[] evdepls;
  delete[] vdepls;
  delete[] models;
}