NueFluxExtrap.cxx

Go to the documentation of this file.

#include "3FlavorAna/Extrap/NueFluxExtrap.h"

#include "CAFAna/Core/Loaders.h"
#include "CAFAna/Core/SpectrumLoaderBase.h"
#include "CAFAna/Core/SystShifts.h"
#include "CAFAna/Cuts/TruthCuts.h"
#include "CAFAna/Decomp/IDecomp.h"
#include "3FlavorAna/Decomp/FluxDecomp.h"
#include "3FlavorAna/Extrap/FluxReweight.h"

namespace ana
{
  //---------------------------------------------------------------------------

  NueFluxExtrap::NueFluxExtrap(SpectrumLoaderBase& nearMC,
                               SpectrumLoaderBase& farMCswap,
                               SpectrumLoaderBase& farMCnonswap,
                               SpectrumLoaderBase& farMCtauswap,
                               const FluxDecomp& nueDecomp,
                               const IDecomp& numuDecomp,
                               const HistAxis& axis,
                               const HistAxis& axisNumuND,
                               const Cut& fdcut,
                               const Cut& nueNDcut,
                               const Cut& numuNDcut,
                               const SystShifts& shiftMC,
                               const Var& weight)
      : NueRHCExtrap(nearMC,
                     farMCswap,
                     farMCnonswap,
                     farMCtauswap,
                     nueDecomp,
                     numuDecomp,
                     axis,
                     axisNumuND,
                     fdcut,
                     nueNDcut,
                     numuNDcut,
                     shiftMC,
                     weight)
  {
    fEEAntiextrap = std::unique_ptr<ModularExtrapComponent>(
        new FluxReweight(
            nearMC, axis, fdcut, shiftMC, weight,
            "ebar -> ebar", "#bar{#nu}_{e} #rightarrow #bar{#nu}_{e}",
            nueNDcut, nueDecomp,                           // nue selection in ND
            farMCnonswap, kIsBeamNue && kIsAntiNu         // e->e in FD
        )
    );

  }

 //---------------------------------------------------------------------------

  NueFluxExtrap::NueFluxExtrap(Loaders& loaders,
                               const FluxDecomp& nueDecomp,
                               const IDecomp& numuDecomp,
                               const HistAxis& axis,
                               const HistAxis& axisNumuND,
                               const Cut& fdcut,
                               const Cut& nueNDcut,
                               const Cut& numuNDcut,
                               const SystShifts& shiftMC,
                               const Var& weight)
      : NueFluxExtrap(loaders.GetLoader(caf::kNEARDET, Loaders::kMC, ana::kBeam, Loaders::kNonSwap),
                      loaders.GetLoader(caf::kFARDET,  Loaders::kMC, ana::kBeam, Loaders::kFluxSwap),
                      loaders.GetLoader(caf::kFARDET,  Loaders::kMC, ana::kBeam, Loaders::kNonSwap),
                      loaders.GetLoader(caf::kFARDET,  Loaders::kMC, ana::kBeam, Loaders::kTauSwap),
                      nueDecomp,
                      numuDecomp,
                      axis,
                      axisNumuND,
                      fdcut,
                      nueNDcut,
                      numuNDcut,
                      shiftMC,
                      weight)
  {}

}