NueAcceptSysts.h

Go to the documentation of this file.

#pragma once

#include "CAFAna/Core/ISyst.h"
#include "CAFAna/Core/Utilities.h"
#include "CAFAna/Core/Var.h"
#include "3FlavorAna/Cuts/QuantileCuts2020.h"

class TH1;

namespace ana
{
  /// FHC systs
  
  class NueAcceptSystBkg2018FHC: public ISyst
  {
  public:
    NueAcceptSystBkg2018FHC() : ISyst("accept_bkg_FHC", "Acceptance Bkg FHC") {}
    
    void Shift(double sigma, caf::SRProxy* sr, double& weight) const override;
  };
  
  extern const NueAcceptSystBkg2018FHC kNueAcceptSystBkg2018FHC;
  
  class NueAcceptSystSignalKin2018FHC: public ISyst
  {
  public:
    NueAcceptSystSignalKin2018FHC() : ISyst("accept_signalkin_FHC_2018", "Acceptance ND to FD Kinematics Signal FHC 2018"), fWeightHistFHC() {}
    
    void Shift(double sigma, caf::SRProxy* sr, double& weight) const override;
    
  protected:
    mutable TH1* fWeightHistFHC;    
  };
  
  extern const NueAcceptSystSignalKin2018FHC kNueAcceptSystSignalKin2018FHC;

  class NueAcceptSystSignalKin2020FHC: public ISyst
  {
  public:
  NueAcceptSystSignalKin2020FHC(int extrap = ana::kExtrapEmpty, std::string shortName = "accept_signalkin_FHC_2020", std::string latexName = "Acceptance ND to FD Kinematics Signal FHC 2020") : ISyst(shortName, latexName), extrapVar(extrap), fWeightHistFHC() {}
    
    void Shift(double sigma, caf::SRProxy* sr, double& weight) const override;

  protected:
    int extrapVar;
    mutable TH1* fWeightHistFHC;    
  };
  
  extern const NueAcceptSystSignalKin2020FHC kNueAcceptSystSignalKin2020FHC;  
  extern const NueAcceptSystSignalKin2020FHC kNueAcceptSystSignalKinPtExtrap2020FHC;
  
  
  /// RHC systs
  class NueAcceptSystBkg2018RHC: public ISyst
  {
  public:
    NueAcceptSystBkg2018RHC() : ISyst("accept_bkg_RHC", "Acceptance Bkg RHC") {}
    
    void Shift(double sigma, caf::SRProxy* sr, double& weight) const override;
  };
  
  extern const NueAcceptSystBkg2018RHC kNueAcceptSystBkg2018RHC;
  
  class NueAcceptSystSignalKin2018RHC: public ISyst
  {
  public:
    NueAcceptSystSignalKin2018RHC() : ISyst("accept_signalkin_RHC_2018", "Acceptance ND to FD Kinematics Signal RHC 2018"), fWeightHistRHC() {}
    
    void Shift(double sigma, caf::SRProxy* sr, double& weight) const override;
    
  protected:
    mutable TH1* fWeightHistRHC;
  };
  
  extern const NueAcceptSystSignalKin2018RHC kNueAcceptSystSignalKin2018RHC;

  class NueAcceptSystSignalKin2020RHC: public ISyst
  {
  public:
    NueAcceptSystSignalKin2020RHC(int extrap = ana::kExtrapEmpty, std::string shortName = "accept_signalkin_RHC_2020", std::string latexName = "Acceptance ND to FD Kinematics Signal RHC 2020") : ISyst(shortName, latexName), extrapVar(extrap), fWeightHistRHC() {}

    void Shift(double sigma, caf::SRProxy* sr, double& weight) const override;
    
  protected:
    int extrapVar;
    mutable TH1* fWeightHistRHC;
  };
  
  extern const NueAcceptSystSignalKin2020RHC kNueAcceptSystSignalKin2020RHC;
  extern const NueAcceptSystSignalKin2020RHC kNueAcceptSystSignalKinPtExtrap2020RHC;
}