getFlatWeightSpectra.C

Go to the documentation of this file.

// I found that the energy estimator used by the nuebar oscillation analysis
// was not sufficient for this analysis due to training on a sample of 
// neutrinos with energy less than 4.5 GeV.
//
// This script marks the first piece of training a new estimator.
// Following docdb-26696, we need to normalize the true energy distribution
// to avoid biasing the energy estimation toward our peak beam energy.
// 
// Grab true energy distributions
// dddoyle@colostate.edu

#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/Var.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Cuts/TruthCuts.h"
#include "3FlavorAna/Cuts/NueCutsSecondAna.h"
#include "3FlavorAna/Vars/NueVars.h"
#include "3FlavorAna/Vars/NueVarsExtra.h"
#include "CAFAna/Vars/XsecTunes.h"
#include "CAFAna/Vars/PPFXWeights.h"

#include "NDAna/nuebarcc_inc/NuebarCCIncVars.h"
#include "NDAna/nuebarcc_inc/NuebarCCIncCuts.h"
#include "NDAna/nuebarcc_inc/NuebarCCIncExtra.h"

using namespace ana;

// selection used in the nuebar CC inclusive analysis
Var weights = kPPFXFluxCVWgt * kXSecCVWgt2018;
std::string training_data = NominalMC_energy_training_RHC_prod4;

const Binning trueE_bins = Binning::Simple(50, 0, 10);


void getFlatWeightSpectra()
{
  Cut cuts = kcPresel && kIsTrueSig; // ktIsNuebarCCST && kTrueFiducialST
  SpectrumLoaderBase* loader = new SpectrumLoader(training_data);

  Spectrum* trueE = new Spectrum("True #nu Energy [GeV]",
                                 trueE_bins,
                                 *loader,
                                 kTrueE, 
                                 cuts,
                                 kNoShift,
                                 weights);

  loader->Go();

  
  TFile* out_file = new TFile("flat_weight_energy_spectra.root", "recreate");
  trueE->SaveTo(out_file, "trueE");

}