void spectrum(int i=0) {
  gROOT->SetStyle("Plain");

  TCanvas *c1=new TCanvas("c1","Canvas",600,600);
  gStyle->SetOptStat(0);
  gStyle->SetPalette(1);
  gStyle->SetPadLeftMargin(0.15);
  gStyle->SetPadRightMargin(0.05);
  c1->UseCurrentStyle();

  if(i==0) TFile::Open("summaryDpw.root");
  if(i==1) TFile::Open("summaryNodpw.root");
  if(i==2) TFile::Open("summaryDpw05.root");
  if(i==3) TFile::Open("summaryNodpw05.root");

  hDif->SetTitle(";Fraction of charge;Epitaxial layer volume (#mum^{3})");

  hDif->GetXaxis()->SetRange(0,60);
  hDif->GetXaxis()->SetTitleOffset(1.1);
  hDif->GetYaxis()->SetTitleOffset(1.6);
  //zAmm0->GetZaxis()->SetRangeUser(0,1.0);
  if(i<2) hDif->GetYaxis()->SetRangeUser(0,8000);
  else    hDif->GetYaxis()->SetRangeUser(0,3000);
  //Trim0Sigmas2D->GetZaxis()->SetRangeUser(0,15);

  //c1->logy();
  hDif->SetLineWidth(2);
  hDif->SetLineStyle(2);
  hDif->Draw();

  TH1F hDifO(*hDif);
  hDifO.Add(hDifC,-1.0);
  hDifO.SetLineColor(4);
  hDifO.SetLineStyle(1);
  hDifO.SetLineWidth(2);
  hDifO.Draw("same");

  hDifC->SetLineColor(2);
  hDifC->SetLineWidth(2);
  hDifC->Draw("same");

  hDif->Draw("same");

  c1->Update();
  if(i==0) c1->Print("../dpwSpectrum.png");
  if(i==1) c1->Print("../nodwSpectrum.png");
  if(i==2) c1->Print("../dpw05Spectrum.png");
  if(i==3) c1->Print("../nodpw05Spectrum.png");
}