GeniePCASyst.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"

class TH1D;

namespace ana
{
  class GeniePCASyst : public ISyst
  {
  public:
    GeniePCASyst(const std::string& fname,
                 const std::string& shortname,
                 const std::string& latexname);
    
    ~GeniePCASyst();
    
    void Shift(double sigma,
               caf::SRProxy* sr,
               double& weight) const override;

    TH1D* getHist(int f,int d,int s, int fl,int p);

  protected:
    enum{kFHC, kRHC, kNumFluxType};
    enum{kND, kFD, kNumDets};
    enum{kNueSel, kNumuSel, kNumSels};
    enum{kSigQE, kSigNonQE, kBkgFlav, kNC, kNumFlavors};
    enum ESign{kMinus, kPlus, kNumSigns};
    
    void InitializeHistograms() const;
    
    double GetWeight(const caf::SRProxy* sr, ESign sign) const;
    
    std::string fFileName;
    std::string fHistName;
    
    mutable TH1D* fHistos[kNumFluxType][kNumDets][kNumSels][kNumFlavors][kNumSigns];
    
  };
  
  // GeniePCASyst* GetGeniePrincipals2020(int PCIdx);
  GeniePCASyst* GetGeniePrincipals2020Small(int PCIdx);
  
  GeniePCASyst* GetGeniePrincipals2018(int PCIdx);
  GeniePCASyst* GetGeniePrincipals2018Small(int PCIdx);
}