FluxWeightSystematics.hxx

#ifndef FluxWeightSystematics_h
#define FluxWeightSystematics_h

#include "EventWeightBase.hxx"
#include "BinnedParams.hxx"
#include "FluxWeighting.hxx"

/// This is a normalization systematic. It takes into account the uncertainty on the neutrino flux

class FluxWeightSystematics: public EventWeightBase, public BinnedParams {
  public:

    FluxWeightSystematics(const std::string& name);  
    virtual ~FluxWeightSystematics() {}

    /// Apply this systematic
    using EventWeightBase::ComputeWeight;
    Weight_h ComputeWeight(const ToyExperiment& toy, const AnaEventC& event, const ToyBoxB& box);
    virtual bool IsCorrectRunPeriod(int runPeriod)=0; ///Bit of a hack right now; hopefully future will have the beam polarity as a variable

protected:
  /// Access to the flux weighting.
  FluxWeighting* _flux;


};

//We don't have 3D systematics, so fake it with two 2D class; also allows neutrino or anti-neutrino 
//only analysis

class FluxWeightSystematicsNeutrino: public FluxWeightSystematics {
  public:

    FluxWeightSystematicsNeutrino();  
    virtual ~FluxWeightSystematicsNeutrino() {}
    bool IsCorrectRunPeriod(int runPeriod);

};

class FluxWeightSystematicsAntiNeutrino: public FluxWeightSystematics {
  public:

    FluxWeightSystematicsAntiNeutrino();  
    virtual ~FluxWeightSystematicsAntiNeutrino() {}
    bool IsCorrectRunPeriod(int runPeriod);


};



#endif