Nus17Systs.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Core/ISyst.h"
#include "CAFAna/Core/Utilities.h"

#include "NuXAna/Systs/NusSysts.h" // for NusChannel

#include <vector>

class TH1D;

#include "StandardRecord/Proxy/FwdDeclare.h"

namespace ana
{
  class Nus17SystFromHist: public ISyst
  {
  public:
    Nus17SystFromHist(const std::string &fname, const std::string &det,
                    const std::string &shortname, const std::string &latexname);
    ~Nus17SystFromHist();

    void Shift(double sigma,
               caf::SRProxy* sr, double &weight) const override;

    double WeightFor(NusChannel chan, double sigma, double nueenergy) const;

  private:
    std::string fFileName;
    std::string fDet;

    void LoadHists() const;

    // get oscillation channel for current event
    NusChannel GetNusChannel(caf::SRProxy* sr) const;
    // Store hist for sig and bkg for each syst, store a pair of a
    // sigma shift value and a histogram evaluated at that shift
    mutable std::vector< std::vector< std::pair<int, TH1D*> > > fHists;
  };

  //----------------------------------------------------------------------
  /// Nus group systematic applied as a flat shift only dependent on NC or CC
  class Nus17FlatSyst: public ISyst
  {
  public:
    Nus17FlatSyst(const std::string &shortname, const std::string &latexname,
                const double &ncWei, const double &bgWei);
    ~Nus17FlatSyst();

    void Shift(double sigma,
               caf::SRProxy* sr, double &weight) const override;

  private:
    double fNCWei;
    double fBGWei;
  };

  //----------------------------------------------------------------------
//  const std::string kNusAna01SystDir = "$NOVA_ANA/steriles/Ana01/Systematics";
  const std::string kNusAna17SystDir = FindCAFAnaDir() + "/data/nus17";
  //const std::string kNusAna17_SystFile = kNusAna17SystDir + "/NusEXISysts17.root";
  const std::string kNusAna17SystFile = kNusAna17SystDir + "/NusSystMaker.root";
  const std::string kNusAna17SystFileND = kNusAna17SystDir + "/NusNDISysts17.root";
  const std::string kNusAna17SystFAlt = kNusAna17SystDir + "/SystsOscParamAna.root";

  // Systematics based on docdb-15330

  // Systematics handled using shape information
//   extern const Nus17SystFromHist kNus17OscParamSysts;
/*   extern const Nus17SystFromHist kNus17BeamSysts;
   extern const Nus17SystFromHist kNus17CherenkovSyst;
   extern const Nus17SystFromHist kNus17CalibFlatSyst;
   extern const Nus17SystFromHist kNus17CalibRelNDSyst;
   extern const Nus17SystFromHist kNus17CalibShape;
   extern const Nus17SystFromHist kNus17Cherenkov_mSyst;
   extern const Nus17SystFromHist kNus17GENIESmallSysts;
   extern const Nus17SystFromHist kNus17NueCCSyst;
   extern const Nus17SystFromHist kNus17NumuCCSyst;
   extern const Nus17SystFromHist kNus17NCSyst;
   extern const Nus17SystFromHist kNus17LightLSyst;
   extern const Nus17SystFromHist kNus17AccSyst;
   extern const Nus17SystFromHist kNus17CalibRelFDSyst;  
*/

  /* extern const Nus17SystFromHist kNus17NDBeamSysts;
  extern const Nus17SystFromHist kNus17NDCherenkovSyst;
   extern const Nus17SystFromHist kNus17NDCalibFlatSyst;
   extern const Nus17SystFromHist kNus17NDCalibShapeSyst;
  extern const Nus17SystFromHist kNus17NDGENIESmallSysts;
   extern const Nus17SystFromHist kNus17NDCherenkov_mSyst;
   extern const Nus17SystFromHist kNus17NDNueCCSyst;
   extern const Nus17SystFromHist kNus17NDNumuCCSyst;
   extern const Nus17SystFromHist kNus17NDNCSyst;
   extern const Nus17SystFromHist kNus17NDLightLSyst;
   extern const Nus17SystFromHist kNus17NDCalibRelSyst;
*/


  // extern const Nus17SystFromHist kNus17FDBeamSysts;
  // extern const Nus17SystFromHist kNus17FDCherenkovSyst;
  // extern const Nus17SystFromHist kNus17FDCalibFlatSyst;
  // extern const Nus17SystFromHist kNus17FDCalibShapeSyst;
  // extern const Nus17SystFromHist kNus17FDGENIESmallSysts;
  // extern const Nus17SystFromHist kNus17FDCherenkov_mSyst;
  // extern const Nus17SystFromHist kNus17FDLightLSyst;
  // extern const Nus17SystFromHist kNus17FDCalibRelSyst;

  // Systematics applied as flat shifts
  // extern const NusFlatSyst kNusMCStatsSyst;
  // extern const Nus17FlatSyst kNus17NormSyst;
  // extern const Nus17FlatSyst kNus17OscParamSysts;
  // extern const Nus17FlatSyst kNusFDPPFXSyst;
  //----------------------------------------------------------------------
  /// Get a vector of all the nus group systs
  // std::vector<const ISyst*> getAllNusFDSysts();
  // std::vector<const ISyst*> getAllNusNDSysts();
  // std::vector<const ISyst*> getAllNusFDRSysts();
} // namespace