MakeNus17PredictionSysts.C

Go to the documentation of this file.

#include "CAFAna/Analysis/Calcs.h"
#include "CAFAna/Analysis/Exposures.h"
#include "CAFAna/Analysis/Prod3Loaders.h"
#include "CAFAna/Core/Loaders.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "3FlavorAna/Cuts/NumuCuts.h"
#include "NuXAna/Cuts/NusCuts.h"
#include "NuXAna/Cuts/NusCuts17.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Cuts/TimingCuts.h"
#include "CAFAna/Decomp/ProportionalDecomp.h"
#include "CAFAna/Extrap/ExtrapSterile.h"
#include "CAFAna/Prediction/PredictionCombinePeriods.h"
#include "CAFAna/Prediction/PredictionExtrap.h"
#include "CAFAna/Prediction/PredictionGenerator.h"
#include "CAFAna/Prediction/PredictionInterp.h"
#include "NuXAna/Prediction/PredictionSterile.h"
#include "NuXAna/Systs/Nus17Systs.h"
#include "CAFAna/Vars/GenieWeights.h"
#include "CAFAna/Vars/HistAxes.h"
#include "CAFAna/Vars/PPFXWeights.h"
#include "OscLib/OscCalcSterile.h"
#include "NuXAna/macros/Nus17/NusLoadProd3.h"
#include "NuXAna/Prediction/PredictionGeneratorNuX.h"
#include "NuXAna/Vars/HistAxes.h"

#include <utility>

using namespace ana;

void MakeNus17PredictionSysts(std::string outfile)
{
  TH1::AddDirectory(0);

  const Var kReweight = kXSecCVWgt2017*kPPFXFluxCVWgt;
  
  // Set up some calculators:
  // Set up oscillation calculator that uses default 3 flavor parameters
  osc::OscCalcSterile* calc3f = DefaultSterileCalc(3);

  Loaders loaders;
  loaders.SetLoaderPath(fNDmcDECAF, caf::kNEARDET, Loaders::kMC, ana::kBeam, Loaders::kNonSwap);
  loaders.SetLoaderPath(fNDdataDECAF, caf::kNEARDET, Loaders::kData, ana::kBeam, Loaders::kNonSwap);
  loaders.SetLoaderPath(fFDmcDECAF, caf::kFARDET, Loaders::kMC, ana::kBeam, Loaders::kNonSwap);
  loaders.SetLoaderPath(fFDswapDECAF, caf::kFARDET, Loaders::kMC, ana::kBeam, Loaders::kFluxSwap);
  loaders.SetLoaderPath(fFDtauDECAF, caf::kFARDET, Loaders::kMC, ana::kBeam, Loaders::kTauSwap);
  loaders.SetSpillCut(kStandardSpillCuts);

  std::vector<const ISyst*> systs = getAllNusFDSysts();

  // Set up the sterile prediction
  const SterileGenerator genNus17(
    kNus17Axis, kNus17BinsNumuCCAxis,
    kNus17FDAlt, kNus17ND, kNumuND,
    kNoShift, kReweight
  );

  PredictionInterp predInterp(systs, calc3f, genNus17, loaders);

  loaders.Go();

  // Set up output filename
  TFile* rootF = new TFile(outfile.c_str(), "RECREATE");

  TDirectory* tmp = gDirectory;
  TDirectory* saveDir = gDirectory;

  // Save all of the objects
  saveDir = rootF->mkdir("nus17_predI");
  predInterp.SaveTo(saveDir);

  tmp->cd();
  rootF->Close(); // Close the file
}