NueSignalOnlyExtrap.cxx

Go to the documentation of this file.

#include "3FlavorAna/Extrap/NueSignalOnlyExtrap.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"

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

  NueSignalOnlyExtrap::NueSignalOnlyExtrap(SpectrumLoaderBase& nearMC,
                                           SpectrumLoaderBase& farMCswap,
                                           SpectrumLoaderBase& farMCnonswap,
                                           SpectrumLoaderBase& farMCtauswap,
                                           const IDecomp& numuDecomp,
                                           const HistAxis& axis,
                                           const HistAxis& axisNumuND,
                                           const Cut& fdcut,
                                           const Cut& numuNDcut,
                                           const SystShifts& shiftMC,
                                           const Var& weight)
      : ModularExtrap(farMCswap,
                      farMCnonswap,
                      farMCtauswap,
                      axis,

                      // Need background extrap comps to be empty
                      // First set all comps to be empty
                      // Then reset signal comps below
                      !kNoCut, //fdcut,
                      shiftMC,
                      weight)
  {
    // mu -> e ----
    fMEextrap = std::unique_ptr<ModularExtrapComponent>(
        new TruthReweight(
            nearMC, axis, axisNumuND, fdcut, shiftMC, weight,
            "mu -> e", "#nu_{#mu} #rightarrow #nu_{e}",
            numuNDcut, numuDecomp,                          // numu selection in ND
            DecompResult::numu, kIsNumuCC && !kIsAntiNu,    // numu truth in ND
            farMCswap, kIsSig && !kIsAntiNu                 // mu->e in FD
        )
    );
    // mubar -> ebar
    fMEAntiextrap = std::unique_ptr<ModularExtrapComponent>(
        new TruthReweight(
            nearMC, axis, axisNumuND, fdcut, shiftMC, weight,
            "mubar -> ebar", "#bar{#nu}_{#mu} #rightarrow #bar{#nu}_{e}",
            numuNDcut, numuDecomp,                          // numu selection in ND
            DecompResult::numubar, kIsNumuCC && kIsAntiNu,  // numubar truth in ND
            farMCswap, kIsSig && kIsAntiNu                  // mubar->ebar in FD
        )
    );
  }

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

  NueSignalOnlyExtrap::NueSignalOnlyExtrap(Loaders& loaders,
                                           const IDecomp& numuDecomp,
                                           const HistAxis& axis,
                                           const HistAxis& axisNumuND,
                                           const Cut& fdcut,
                                           const Cut& numuNDcut,
                                           const SystShifts& shiftMC,
                                           const Var& weight)
      : NueSignalOnlyExtrap(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),
                            numuDecomp,
                            axis,
                            axisNumuND,
                            fdcut,
                            numuNDcut,
                            shiftMC,
                            weight)
  {}

}