demo_flat.C

Go to the documentation of this file.

#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Analysis/Calcs.h"
#include "CAFAna/Analysis/Prod4Loaders.h"
#include "CAFAna/Prediction/PredictionNoExtrap.h"

#include "3FlavorAna/Vars/HistAxes.h"
#include "3FlavorAna/Vars/NumuVars.h"
#include "3FlavorAna/Cuts/NumuCuts2020.h"

using namespace ana;

#include "OscLib/IOscCalc.h"

#include "TH1.h"
#include "TPad.h"

#include <iostream>

// "" = classic nested CAFs. Other possible strings are "multitree" and
// "singletree"
void demo_flat(std::string flat = "")
{
  const std::string prefix = flat.empty() ? "" : "bckhouse_"+flat+"_flatcaf_";

  SpectrumLoader loader_nonswap (prefix+"prod_sumdecaf_R17-11-14-prod4reco.neutron-respin.c_fd_genie_nonswap_fhc_nova_v08_period5_v1_numu2018");
  SpectrumLoader loader_fluxswap(prefix+"prod_sumdecaf_R17-11-14-prod4reco.neutron-respin.c_fd_genie_fluxswap_fhc_nova_v08_period5_v1_numu2018");

  // All the interesting parts of the cut are apparently failed by these prod4
  // files (they want prod5).
  PredictionNoExtrap pred(loader_nonswap, loader_fluxswap, kNullLoader, kNumuCCOptimisedAxis, kNumuQuality && kNumuContainFD2020); // && kNumu2020PID && kNumu2020CosRej && k3flavor2020FDVeto);

  loader_nonswap.Go();
  loader_fluxswap.Go();

  pred.Predict(DefaultOscCalc()).ToTH1(9e20, kBlack, kSolid, kPOT, kBinDensity)->Draw("hist");
  std::cout << pred.Predict(DefaultOscCalc()).Integral(9e20) << " events" << std::endl;
  gPad->Print("demo_flat_spect.pdf");
}