JointAna2018Systs.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 EAnaType2018 {kNueAna2018, kNumuAna2018, kJointAna2018};
 enum BeamType2018 {kFHC, kRHC, kBoth};

  std::string AnaTypeToString(const EAnaType2018 ana);

  void AddJointAna2018XSecSysts(std::vector<const ISyst*> &systs,
        const EAnaType2018 ana = kJointAna2018,
        bool smallgenie = true);

  void AddJointAna2018BeamSysts (std::vector<const ISyst*> &systs,
         const EAnaType2018 ana = kJointAna2018);

  void AddJointAna2018FileSysts (std::vector<const ISyst*> &systs,
         const EAnaType2018 ana = kJointAna2018, const BeamType2018 beam = kBoth);

  void AddJointAna2018NormSysts(std::vector<const ISyst*> &systs,
         const EAnaType2018 ana = kJointAna2018);

  void AddJointAna2018OtherSysts(std::vector<const ISyst*> &systs,
         const EAnaType2018 ana = kJointAna2018, const BeamType2018 beam = kBoth);

  void AddJointAna2018NotCorrelSysts(std::vector<const ISyst*> &systs,
         const EAnaType2018 ana = kJointAna2018, const BeamType2018 beam = kBoth);

  void AddNonLoadable2018Systs(std::vector<const ISyst*> &systs,
         const EAnaType2018 ana = kJointAna2018);

  void AddWrongSignSysts(std::vector<const ISyst*> &systs,
                         const EAnaType2018 ana = kJointAna2018);

  std::vector<const ISyst *> getAna2018LargeXsecSysts(
    const EAnaType2018 ana = kJointAna2018);
  std::vector<const ISyst *> getAna2018SmallXsecSysts(
    const EAnaType2018 ana = kJointAna2018);
  const SummedSyst * getAna2018SummedSmallXsecSysts(
    const EAnaType2018 ana = kJointAna2018);
  
//////////////////////////////////////////////////////////////////////

  std::vector<const ISyst*> getAllAna2018Systs(
    const EAnaType2018 ana = kJointAna2018, const bool smallgenie = true, const BeamType2018 beam = kBoth, bool isFit = false); //fit or syst loading? some systs are not in files (like rock)


  std::vector<const ISyst*> getAllCPTSysts
    (const EAnaType2018 ana = kJointAna2018, const bool smallgenie = true, const BeamType2018 beam = kBoth, bool isFit = false);

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

  std::vector<const ISyst *> getAna2018NueOnly( const BeamType2018 beam);
  std::vector<const ISyst *> getAna2018NumuOnly( const BeamType2018 beam);

 
} // namespace