test_nuwro.C

Go to the documentation of this file.

#include "CAFAna/Core/Binning.h"
#include "CAFAna/Core/Ratio.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "3FlavorAna/Cuts/NumuCuts.h"
#include "CAFAna/Systs/NuWROSysts.h"
#include "CAFAna/Vars/Vars.h"
#include "3FlavorAna/Vars/Binnings.h"

using namespace ana;

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

void test_nuwro()
{
  SpectrumLoader loader("$NOVA_PROD/mc/FA14-12-29/genie/nd/caf/000103/10377/neardet_genie_fhc_nonswap*s9*.caf.root");

  HistAxis axis("Reconstructed energy (GeV)", kNumuEnergyBinning, kCCE);

  Spectrum sNom(loader, axis, kNumuND, kNoShift);
  Spectrum sUp(loader, axis, kNumuND, SystShifts(&kNuWROReweightEnu, +1));
  Spectrum sDn(loader, axis, kNumuND, SystShifts(&kNuWROReweightEnu, -1));

  loader.Go();

  sNom.ToTH1(18e20)->Draw("hist");

  TH1* hup = sUp.ToTH1(18e20);
  hup->SetLineColor(kRed);
  hup->Draw("hist same");

  TH1* hdn = sDn.ToTH1(18e20);
  hdn->SetLineColor(kBlue);
  hdn->Draw("hist same");

  new TCanvas;
  Ratio rUp(sUp, sNom);
  TH1* hrup = rUp.ToTH1();
  hrup->SetLineColor(kRed);
  hrup->Draw("hist");
  hrup->GetYaxis()->SetRangeUser(.8, 1.2);
  Ratio rDn(sDn, sNom);
  TH1* hrdn = rDn.ToTH1();
  hrdn->SetLineColor(kBlue);
  hrdn->Draw("hist same");
}