demo2p5b.C

Go to the documentation of this file.

// Example of how to read spectra saved to disk in and make
// pretty plots...
//
// This is the second half of the demo2p5a script...

#include "TCanvas.h"
#include "TFile.h"  
#include "TH1.h"    

#include "CAFAna/Core/Spectrum.h"



using namespace ana;



void demo2p5b(){

  TFile inFile("save_your_spectra_to_disk.root");

  // Load the spectrum...
  std::unique_ptr<Spectrum> nhit = Spectrum::LoadFrom(inFile.GetDirectory("dir_nhit_spectra"));



  //
  // Plot the histo...
  //  
  TCanvas *cNHit = new TCanvas("cNHit","Slice NHit",1000,800);
  cNHit->cd();
  TH1D *hNHit = nhit->ToTH1(nhit->POT());
  hNHit->SetLineWidth(2);
  hNHit->SetLineColor(kGreen+2);
  hNHit->Draw("hist");

}