3FlavorAna2020Systs.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Core/ISyst.h"
#include "CAFAna/Systs/SummedSysts.h"

#include <vector>

namespace ana
{

  enum class EAnaType2020 {kNueAna, kNumuAna, k3FlavorAna};
  enum class BeamType2020 {kFHC, kRHC, kBoth};

  std::string AnaTypeToString(const EAnaType2020 ana);

  void Add3FlavorAna2020XSecSysts(std::vector<const ISyst*> &systs,
        const EAnaType2020 ana = EAnaType2020::k3FlavorAna,
        bool smallgenie = true);

  void Add3FlavorAna2020BeamSysts (std::vector<const ISyst*> &systs,
         const EAnaType2020 ana = EAnaType2020::k3FlavorAna);

  void Add3FlavorAna2020CalibSysts( std::vector<const ISyst*>& systs );
  void Add3FlavorAna2020LightSysts( std::vector<const ISyst*>& systs );
  void Add3FlavorAna2020FileSysts ( std::vector<const ISyst*>& systs );

  void Add3FlavorAna2020MuEnergySysts(
         std::vector<const ISyst*> &systs,
         const EAnaType2020 ana = EAnaType2020::k3FlavorAna );

  void Add3FlavorAna2020LeptonAngleSysts(
         std::vector<const ISyst*> &systs,
         const bool         ptExtrap = true );

  void Add3FlavorAna2020NeutronSysts( std::vector<const ISyst*> &systs );

  void Add3FlavorAna2020NueAcceptSysts(
         std::vector<const ISyst*> &systs,
         const EAnaType2020 ana      = EAnaType2020::k3FlavorAna,
         const BeamType2020 beam     = BeamType2020::kBoth,
         const bool         ptExtrap = true );

  void Add3FlavorAna2020MichelTagSysts(
         std::vector<const ISyst*> &systs,
         const EAnaType2020 ana  = EAnaType2020::k3FlavorAna,
         const BeamType2020 beam = BeamType2020::kBoth );

  void Add3FlavorAna2020NotCorrelSysts(
         std::vector<const ISyst*> &systs,
         const EAnaType2020 ana      = EAnaType2020::k3FlavorAna,
         const BeamType2020 beam     = BeamType2020::kBoth,
         const bool         ptExtrap = true ); // for nue acceptance sytstematic

  void Add3FlavorAna2020NormSysts(
         std::vector<const ISyst*> &systs,
         const BeamType2020 beam = BeamType2020::kBoth );

  void AddNonLoadable2020Systs(
         std::vector<const ISyst*> &systs,
         const EAnaType2020 ana = EAnaType2020::k3FlavorAna);

  void Add3FlavorAna2020OtherSysts(
         std::vector<const ISyst*> &systs,
         const EAnaType2020 ana      = EAnaType2020::k3FlavorAna,
         const BeamType2020 beam     = BeamType2020::kBoth,
         const bool         ptExtrap = true ); // for nue acceptance sytstematic

  std::vector<const ISyst *> get3FlavorAna2020LargeXsecSysts(
    const EAnaType2020 ana = EAnaType2020::k3FlavorAna);
  std::vector<const ISyst *> get3FlavorAna2020SmallXsecSysts(
    const EAnaType2020 ana = EAnaType2020::k3FlavorAna);
  const SummedSyst * get3FlavorAna2020SummedSmallXsecSysts(
    const EAnaType2020 ana = EAnaType2020::k3FlavorAna);

//////////////////////////////////////////////////////////////////////

  std::vector<const ISyst*> get3FlavorAna2020AllSysts(
    const EAnaType2020 ana        = EAnaType2020::k3FlavorAna,
    const bool         smallgenie = true,
    const BeamType2020 beam       = BeamType2020::kBoth,
    const bool         isFit      = false, //fit or syst loading? some systs are not in files (like rock)
    const bool         ptExtrap   = true ); // for nue acceptance sytstematic


  std::vector<const ISyst*> get3FlavorAna2020AllCPTSysts(
    const EAnaType2020 ana        = EAnaType2020::k3FlavorAna,
    const bool         smallgenie = true,
    const BeamType2020 beam       = BeamType2020::kBoth,
    const bool         isFit      = false, //fit or syst loading? some systs are not in files (like rock)
    const bool         ptExtrap   = true ); // for nue acceptance sytstematic

//////////////////////////////////////////////////////////////////////

  std::vector<const ISyst *> getAna2020NueOnly( const BeamType2020 beam, const bool ptExtrap );
  std::vector<const ISyst *> getAna2020NumuOnly( const BeamType2020 beam );
} // namespace