PredictionGeneratorNuX.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Prediction/PredictionGenerator.h"

namespace ana
{

  /// Generates extrapolated NC predictions using ProportionalDecomp
  class SterileGenerator: public IPredictionGenerator
  {
    public:
      SterileGenerator(
          const HistAxis axis,
          const HistAxis axisNDNumu,
          const Cut cutFD,
          const Cut cutNDNC,
          const Cut cutNDNumu,
          const SystShifts shiftDataND = kNoShift,
          const Var wei = kUnweighted );
      std::unique_ptr<IPrediction> Generate(
          Loaders& loaders,
          const SystShifts& shiftMC = kNoShift ) const override;
      std::unique_ptr<IPrediction> Generate(
          Loaders& loaders_nd,
          Loaders& loaders_fd,
          const SystShifts& shiftMC = kNoShift ) const;
    private:
      const HistAxis fAxis;
      const HistAxis fAxisNDNumu;
      const Cut fCutFD;
      const Cut fCutNDNC;
      const Cut fCutNDNumu;
      const SystShifts fShiftDataND;
      const Var fWei;
  };

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

  // const SterileGenerator kSterileGenerator(
  //   kNCAxis, kNCBinsNumuCCAxis, kNusFD, kNusND, kNumuND, kNoShift, kTuftsWeightCC );

  //............................................................................
  /// Generates Near Detector predictions
  class NDPredictionGenerator: public IPredictionGenerator
  {
    public:
      NDPredictionGenerator(
          const HistAxis axis,
          const Cut cutND,
          const SystShifts shiftND = kNoShift,
          const Var wei = kUnweighted );
      std::unique_ptr<IPrediction> Generate(
          Loaders& loaders,
          const SystShifts& shiftMC = kNoShift ) const override;
    private:
      const HistAxis fAxis;
      const Cut fCutND;
      const SystShifts fShiftND;
      const Var fWei;
  };

  class FDPredictionGenerator: public IPredictionGenerator
  {
    public:
      FDPredictionGenerator(
          const HistAxis axis,
          const Cut cutFD,
          const SystShifts shiftFD = kNoShift,
          const Var wei = kUnweighted );
      std::unique_ptr<IPrediction> Generate(
          Loaders& loaders,
          const SystShifts& shiftMC = kNoShift ) const override;
    private:
      const HistAxis fAxis;
      const Cut fCutFD;
      const SystShifts fShiftFD;
      const Var fWei;
  };
}