WeightCalc.h

Go to the documentation of this file.

#ifndef _WEIGHTCALC_H_
#define _WEIGHTCALC_H_
#include "art/Framework/Principal/Event.h"
#include "fhiclcpp/ParameterSet.h"

#include "art/Framework/Services/Optional/RandomNumberGenerator.h"
#include "CLHEP/Random/RandomEngine.h"
#include "CLHEP/Random/RandGaussQ.h"

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

//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) = 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_