NCDecomp.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Decomp/IDecomp.h"

#include "CAFAna/Core/Loaders.h"

namespace ana
{
  /// Uses MC for \f$ \nu_\mu \f$ CC and \f$ \nu_e \f$ CC components, assigns remainder of data to NC
  class NCDecomp: public IDecomp
  {
  public:
    NCDecomp(SpectrumLoaderBase& loaderMC,
               SpectrumLoaderBase& loaderData,
               const std::string& label,
               const Binning& bins,
               const Var& var,
               const Cut& cut,
               const SystShifts& shiftMC = kNoShift,
               const SystShifts& shiftData = kNoShift,
               const Var& wei = kUnweighted);

    NCDecomp(SpectrumLoaderBase& loaderMC,
               SpectrumLoaderBase& loaderData,
               const HistAxis& axis,
               const Cut& cut,
               const SystShifts& shiftMC = kNoShift,
               const SystShifts& shiftData = kNoShift,
               const Var& wei = kUnweighted);

    NCDecomp(Loaders& loaders,
               const HistAxis& axis,
               const Cut& cut,
               const SystShifts& shiftMC = kNoShift,
               const SystShifts& shiftData = kNoShift,
               const Var& wei = kUnweighted);

    Spectrum NumuComponent()     const override {return fNumu;}
    Spectrum AntiNumuComponent() const override {return fAntiNumu;}
    Spectrum NCTotalComponent()  const override;
    Spectrum NCComponent()  const override;
    Spectrum NCAntiComponent()  const override;
    Spectrum NueComponent()      const override {return fNue;}
    Spectrum AntiNueComponent()  const override {return fAntiNue;}
    Spectrum Data_Component()    const override {return fData;}

    void SaveTo(TDirectory* dir, const std::string& name) const override;
    static std::unique_ptr<NCDecomp> LoadFrom(TDirectory* dir, const std::string& name);

  protected:
    NCDecomp()
      : fData(Spectrum::Uninitialized()),
        fNCTot   (Spectrum::Uninitialized()),
        fNue (Spectrum::Uninitialized()),
        fAntiNue(Spectrum::Uninitialized()),
        fNumu(Spectrum::Uninitialized()),
        fAntiNumu(Spectrum::Uninitialized()),
        fNotNC(Spectrum::Uninitialized())
    {};

    Spectrum fData;
    Spectrum fNCTot;
    Spectrum fNue;
    Spectrum fAntiNue;

    Spectrum fNumu;
    Spectrum fAntiNumu;

    Spectrum fNotNC;
  };
}