void read_and_draw_vdepl_Perugia2015_fl1e15(const char* filename) {
  std::fstream file;
  file.open(filename);
  if (! file.good() ){
    std::cout << filename << " is not good: " << file.rdstate() << "\n";
    exit(1);
  }
  
  int nmax = 10;
  int i = 0;
  TString** models = (TString**) (new TString[nmax]);
  double* vdepls = new double[nmax];
  double* evdepls = new double[nmax];
  double vdepl,evdepl;
  std::string model;
  while(1) {
    file >> model >> vdepl >> evdepl;
    if ( file.eof() ) break;
    vdepls[i]=vdepl;
    evdepls[i] = evdepl;
    models[i] = new TString(model.c_str());
    std::cout << model << "\n";
    std::cout << ((models)[i])->Data() << "\n";
    std::cout << models[i] << "\n";
    i++;
  }
  std::cout << i << " lines\n";
  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[i] << "\n";
    std::cout << ((models)[i])->Data() << "\n";
    hVdepl->GetXaxis()->SetBinLabel(j+1,((models)[i])->Data());
  }
  TCanvas *c1 = new TCanvas("c1","",800,600);
  c1->cd();
  hVdepl->Draw("e");

  gStyle->SetOptStat(0);

  file.close();
}