PredictionSystNumu2017.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Prediction/PredictionInterp.h"

#include "CAFAna/Core/ISyst.h"

namespace ana
{
  /// Doesn't actually shift anything. Contains info about what files to load
  /// shifted predictions from.
  class DummyNumu2017Syst: public ISyst
  {
  public:
    /// \a shortName must match exactly the component in the filename
    DummyNumu2017Syst(const std::string& shortName,
                     const std::string& longName,
                     const std::vector<int>& shifts, bool isAltSample)
      : ISyst(shortName, longName), fShifts(shifts),
        fIsAltSample(isAltSample) {}
    virtual void Shift(double, caf::SRProxy*, double&) const {}

    std::vector<int> Shifts() const {return fShifts;}
    bool IsAltSample() const {return fIsAltSample;}

  protected:
    std::vector<int> fShifts;
    bool fIsAltSample;
  };

  // TODO: sort out what systs we'll need
  extern const DummyNumu2017Syst kNumu2017MECq0ShapeSyst;
  extern const DummyNumu2017Syst kNumu2017MaCCQE_reducedSyst;
  extern const DummyNumu2017Syst kNumu2017RPAShapeRESSyst;
  extern const DummyNumu2017Syst kNumu2017RPAShapeEnhSyst;
  extern const DummyNumu2017Syst kNumu2017MaCCRESSyst;
  extern const DummyNumu2017Syst kNumu2017MaNCRESSyst;
  extern const DummyNumu2017Syst kNumu2017MvCCRESSyst;
  extern const DummyNumu2017Syst kNumu2017DISvnCC1piSyst;
  extern const DummyNumu2017Syst kNumu2017AllSmallSysts;
  extern const DummyNumu2017Syst kNumu2017NumuSmallSysts;

  extern const DummyNumu2017Syst kNumu2017AllBeamSysts;

  extern const DummyNumu2017Syst kNumu2017RelHadEScale;
  extern const DummyNumu2017Syst kNumu2017AbsHadEScale;
  extern const DummyNumu2017Syst kNumu2017RelMuEScale;
  extern const DummyNumu2017Syst kNumu2017AbsMuEScale;

  extern const DummyNumu2017Syst kNumu2017AbsCalibSyst;
  extern const DummyNumu2017Syst kNumu2017RelCalibSyst;
  extern const DummyNumu2017Syst kNumu2017CalibShapeSyst;
  extern const DummyNumu2017Syst kNumu2017LightLevelSyst;
  extern const DummyNumu2017Syst kNumu2017CherenkovSyst;

  extern const DummyNumu2017Syst kNumu2017PPFXPC00Syst;
  extern const DummyNumu2017Syst kNumu2017PPFXPC01Syst;
  extern const DummyNumu2017Syst kNumu2017PPFXPC02Syst;
  extern const DummyNumu2017Syst kNumu2017PPFXPC03Syst;
  extern const DummyNumu2017Syst kNumu2017PPFXPC04Syst;

  //extern const DummyNumu2017Syst kNumu2017AbsCalibSyst2Sig;

  extern const DummyNumu2017Syst kNumu2017NormSystBoth;
  extern const DummyNumu2017Syst kNumu2017NormSystSig;
  extern const DummyNumu2017Syst kNumu2017NormSystBkg;


  /// These are the systs that \ref PredictionSystNumu2017 will be able to
  /// handle.
  const std::vector<const DummyNumu2017Syst*> kNumu2017Systs =
    {
      &kNumu2017MECq0ShapeSyst,
      &kNumu2017MaCCQE_reducedSyst,
      &kNumu2017RPAShapeRESSyst,
      &kNumu2017RPAShapeEnhSyst,
      &kNumu2017MaCCRESSyst,
      &kNumu2017MaNCRESSyst,
      &kNumu2017MvCCRESSyst,
      &kNumu2017DISvnCC1piSyst,
      &kNumu2017AllSmallSysts,
      &kNumu2017NumuSmallSysts,

      &kNumu2017AllBeamSysts,

      &kNumu2017RelHadEScale,
      &kNumu2017AbsHadEScale,
      &kNumu2017RelMuEScale,
      &kNumu2017AbsMuEScale,

      &kNumu2017AbsCalibSyst,
      &kNumu2017RelCalibSyst,
      &kNumu2017CalibShapeSyst,
      &kNumu2017LightLevelSyst,
      &kNumu2017CherenkovSyst,

      &kNumu2017PPFXPC00Syst,
      &kNumu2017PPFXPC01Syst,
      &kNumu2017PPFXPC02Syst,
      &kNumu2017PPFXPC03Syst,
      &kNumu2017PPFXPC04Syst,

//      &kNumu2017AbsCalibSyst2Sig,

      &kNumu2017NormSystBoth
//      &kNumu2017NormSystSig,
//      &kNumu2017NormSystBkg
    };

  enum class ENumu2017ExtrapType{
    kNumuNoExtrap,
    kNuMu
  };

  /// Loads shifted spectra from files
  class PredictionSystNumu2017: public PredictionInterp
  {
  public:
    PredictionSystNumu2017(ENumu2017ExtrapType extrap,
                           osc::IOscCalc* osc,
                           int WhichQuant = 1,
                           const std::vector<const DummyNumu2017Syst*>& systs = kNumu2017Systs,
                           const std::string fname = "/nova/ana/nu_mu_ana/Ana2017/Predictions/provisional/pred_numu_reduced_v5.root");

    virtual ~PredictionSystNumu2017();

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

  protected:
    PredictionSystNumu2017() {}

    void AddNormSysts(const std::vector<const DummyNumu2017Syst*>& systs);

  };
}