make_extrap_figure_hists.C

Go to the documentation of this file.

// This is the macro to make the inputs for make_extrap_figure.C
//
// It sets up the very simplest numu extrapolation and dumps the debug
// histograms.

#include "CAFAna/Core/Loaders.h"
#include "CAFAna/Core/OscCurve.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Extrap/ModularExtrap.h"
#include "CAFAna/Prediction/PredictionExtrap.h"
#include "3FlavorAna/Prediction/PredictionGenerator3Flavor.h"
#include "3FlavorAna/Vars/Binnings.h"
#include "3FlavorAna/Vars/NumuVars.h"
#include "CAFAna/Vars/PPFXWeights.h"
#include "CAFAna/Analysis/Calcs.h"
using namespace ana;

#include "OscLib/IOscCalc.h"

#include "TFile.h"

#include <iostream>

const double pot = 9e20;

void make_extrap_figure_hists()
{
  const std::string fname("prod_sumdecaf_R17-03-01-prod3reco.k_fd_genie_nonswap_fhc_nova_v08_full_nue_or_numu_or_nus_contain_v1");
  const std::string fnameSwap("prod_sumdecaf_R17-03-01-prod3reco.k_fd_genie_fluxswap_fhc_nova_v08_full_nue_or_numu_or_nus_contain_v1");
  const std::string fnameTau("prod_sumdecaf_R17-03-01-prod3reco.k_fd_genie_tau_fhc_nova_v08_full_nue_or_numu_or_nus_contain_v1");

  const std::string fnameNDData("prod_sumdecaf_R17-03-01-prod3reco.h_nd_numi_fhc_full_nue_or_numu_or_nus_contain_v1_goodruns");
  const std::string fnameNDMC("prod_sumdecaf_R17-03-01-prod3reco.d_nd_genie_nonswap_fhc_nova_v08_full_nue_or_numu_or_nus_contain_v1");

  Loaders loaders;

  loaders.SetLoaderPath(fname, caf::kFARDET, Loaders::kMC,
                        ana::kBeam, Loaders::kNonSwap);
  loaders.SetLoaderPath(fnameSwap, caf::kFARDET, Loaders::kMC,
                        ana::kBeam, Loaders::kFluxSwap);
  loaders.SetLoaderPath(fnameTau, caf::kFARDET, Loaders::kMC,
                        ana::kBeam, Loaders::kTauSwap);
  loaders.SetLoaderPath(fnameNDData, caf::kNEARDET, Loaders::kData);
  loaders.SetLoaderPath(fnameNDMC, caf::kNEARDET, Loaders::kMC);

  loaders.SetSpillCut(kStandardSpillCuts);

  const HistAxis axis("Reconstructed Energy (GeV)", kNumuEnergyBinning, kCCE);

  const NumuExtrapGenerator gen(axis, kNumuFD, kNumuND, kNoShift, kPPFXFluxCVWgt * kXSecCVWgt2017);

  IPrediction* pred = gen.Generate(loaders).release();


  loaders.Go();


  TFile* fout = new TFile("extrap_hists.root", "RECREATE");

  ModularExtrap* extrap = (ModularExtrap*)((PredictionExtrap*)pred)->GetExtrap();

  extrap->SavePlotsNumu(fout, pot);

  osc::IOscCalc* calc = DefaultOscCalc();

  OscCurve curve(calc, 14, 14);
  curve.ToTH1()->Write("osc");

  std::cout << "Wrote extrap_hists.root" << std::endl;
}