WeightCalc.h

Go to the documentation of this file.

#ifndef _WEIGHTCALC_H_
#define _WEIGHTCALC_H_

#include "art/Framework/Principal/fwd.h"
namespace fhicl {
  class ParameterSet;
}
namespace CLHEP {
  class HepRandomEngine;
  class RandGaussQ;
}

#include "TMatrixD.h"
#include <map>
#include <string>

//weight calc base
namespace evwgh {
  typedef std::map<std::string, std::vector<double>> WeightMap_t;

  class WeightCalc {
  public:
    virtual void Configure(fhicl::ParameterSet const& pset, CLHEP::HepRandomEngine&) = 0;
    virtual std::vector<std::vector<double>> GetWeight(art::Event& e) = 0;
    void SetName(std::string name) { fName = name; }
    std::string GetName() { return fName; }

    static std::vector<std::vector<double>> MultiGaussianSmearing(
      std::vector<double> const& centralValues,
      std::vector<std::vector<double>> const& inputCovarianceMatrix,
      int n_multisims,
      CLHEP::RandGaussQ& GaussRandom);

    static std::vector<double> MultiGaussianSmearing(std::vector<double> const& centralValue,
                                                     TMatrixD* const& inputCovarianceMatrix,
                                                     std::vector<double> rand);

    static std::vector<double> MultiGaussianSmearing(std::vector<double> const& centralValue,
                                                     TMatrixD* const& LowerTriangleCovarianceMatrix,
                                                     bool isDecomposed,
                                                     std::vector<double> rand);

  private:
    std::string fName;
  };

}

#endif // _WEIGHTCALC_H_