test_prediction_interp.C

Go to the documentation of this file.

#include "CAFAna/Analysis/Calcs.h"
#include "CAFAna/Core/Loaders.h"
#include "CAFAna/Core/LoadFromFile.h"
#include "3FlavorAna/Cuts/NumuCuts.h"
#include "CAFAna/Prediction/PredictionInterp.h"
#include "3FlavorAna/Systs/EnergySysts.h"
using namespace ana;

#include "OscLib/IOscCalc.h"

#include "TCanvas.h"
#include "TGraph.h"
#include "TH2.h"
#include "3FlavorAna/Vars/NumuVars.h"

void test_prediction_interp()
{
  Loaders loaders;
  loaders.SetLoaderPath("defname: prod_caf_S15-05-22_fd_genie_fhc_nonswap with stride 100",
                        caf::kFARDET, Loaders::kMC, kBeam, Loaders::kNonSwap);

  osc::IOscCalc* calc = DefaultOscCalc();

  const Binning bins = Binning::Simple(10, 0, 5);
  NoExtrapGenerator gen(HistAxis("CC E", bins, kCCE), kNumuCC);
  PredictionInterp pred({&kEnergyScaleSyst}, calc, gen, loaders);

  loaders.Go();

  pred.DebugPlots(calc);

  // Make sure we get the same answer after going to a file and back
  SaveToFile(pred, "pred.root", "pred");
  PredictionInterp* pred2 = LoadFromFile<PredictionInterp>("pred.root", "pred").release();
  pred2->DebugPlots(calc);
}