evwgh::WeightCalc Class Reference

#include "WeightCalc.h"

Inheritance diagram for evwgh::WeightCalc:

Public Member Functions
virtual void	Configure (fhicl::ParameterSet const &pset)=0
virtual std::vector< std::vector< double > >	GetWeight (art::Event &e)=0
void	SetName (std::string name)
std::string	GetName ()

Static Public Member Functions
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)
	Applies Gaussian smearing to a set of data. More...
static std::vector< double >	MultiGaussianSmearing (std::vector< double > const &centralValue, TMatrixD *const &inputCovarianceMatrix, std::vector< double > rand)
	Over load of the above function that only returns a single varied parameter set. More...
static std::vector< double >	MultiGaussianSmearing (std::vector< double > const &centralValue, TMatrixD *const &LowerTriangleCovarianceMatrix, bool isDecomposed, std::vector< double > rand)

Private Attributes
std::string	fName

Detailed Description

Definition at line 18 of file WeightCalc.h.

Member Function Documentation

virtual void evwgh::WeightCalc::Configure ( fhicl::ParameterSet const &  pset )

pure virtual

Implemented in evwgh::GenieWeightCalc.

std::string evwgh::WeightCalc::GetName ( )

inline

Definition at line 24 of file WeightCalc.h.

References fName, MultiGaussianSmearing(), and lar::dump::vector().

Referenced by evwgh::GenieWeightCalc::Configure().

{return fName;}

virtual std::vector<std::vector<double> > evwgh::WeightCalc::GetWeight ( art::Event &  e )

pure virtual

Implemented in evwgh::GenieWeightCalc.

Referenced by evwgh::Weight_t::GetWeight().

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

static

Applies Gaussian smearing to a set of data.

Parameters

centralValues	the values to be smeared
inputCovarianceMatrix	covariance matrix for smearing
n_multisims	number of sets of smeared values to be produced

Returns

a set of n_multisims value sets smeared from the central value

If centralValues is of dimension N, inputCovarianceMatrix needs to be NxN, and each of the returned data sets will be also of dimension N.

Referenced by GetName().

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

static

Over load of the above function that only returns a single varied parameter set.

Definition at line 106 of file WeightCalc.cxx.

References col, and art::errors::StdException.

  {

    //compute the smeared central values
    std::vector<double> smearedCentralValues;

    //
    //perform Choleskey Decomposition
    //
    // Best description I have found 
    //     is in the PDG (Monte Carlo techniques, Algorithms, Gaussian distribution)
    //
    //  http://pdg.lbl.gov/2016/reviews/rpp2016-rev-monte-carlo-techniques.pdf (Page 5)
    //
    //
    TDecompChol dc = TDecompChol(*(inputCovarianceMatrix));
    if(!dc.Decompose())
      {
        throw art::Exception(art::errors::StdException)
          << "Cannot decompose covariance matrix to begin smearing.";
        return smearedCentralValues;
      }

    //Get upper triangular matrix. This maintains the relations in the
    //covariance matrix, but simplifies the structure.
    TMatrixD U = dc.GetU();

 
    for(unsigned int col = 0; col < centralValue.size(); ++col)
      {
        //find the weight of each col of the upper triangular cov. matrix
        double weightFromU = 0.;
        
        for(unsigned int row = 0; row < col+1; ++row)
          {
            weightFromU += U(row,col)*rand[row];
          }

        //multiply this weight by each col of the central values to obtain
        //the gaussian smeared and constrained central values
        // smearedCentralValues.push_back(weightFromU * centralValue[col]);
        smearedCentralValues.push_back(weightFromU + centralValue[col]);
      }
  return smearedCentralValues;
  } // WeightCalc::MultiGaussianSmearing() - Gives a single set of smeared parameters

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

static

Definition at line 163 of file WeightCalc.cxx.

References col, and art::errors::StdException.

  {
    
    //compute the smeared central values
    std::vector<double> smearedCentralValues;
    
    // This guards against accidentally 
    if(!isDecomposed){
      throw art::Exception(art::errors::StdException)
        << "Must supply the decomposed lower triangular covariance matrix.";
      return smearedCentralValues;      
    } 

     
    for(unsigned int col = 0; col < centralValue.size(); ++col)
      {
        //find the weight of each col of the upper triangular cov. matrix
        double weightFromU = 0.;
        
        for(unsigned int row = 0; row < col+1; ++row)
          {
            weightFromU += LowerTriangleCovarianceMatrix[0][row][col]*rand[row];
          }

        //multiply this weight by each col of the central values to obtain
        //the gaussian smeared and constrained central values
        // smearedCentralValues.push_back(weightFromU * centralValue[col]);
        smearedCentralValues.push_back(weightFromU + centralValue[col]);
      }
  return smearedCentralValues;
  } // WeightCalc::MultiGaussianSmearing() - Gives a single set of smeared parameters

void evwgh::WeightCalc::SetName ( std::string  name )

inline

Definition at line 23 of file WeightCalc.h.

References fName.

{fName=name;}

Member Data Documentation

std::string evwgh::WeightCalc::fName

private

Definition at line 54 of file WeightCalc.h.

Referenced by GetName(), and SetName().

---

The documentation for this class was generated from the following files:

• larsim/v06_53_00/source/larsim/EventWeight/Base/WeightCalc.h
• larsim/v06_53_00/source/larsim/EventWeight/Base/WeightCalc.cxx