#include <PDFWeightsHelper.h>

Public Member Functions
void	DoMC2Hessian (double nomweight, const double inweights, double outweights) const

void	Init (unsigned int nreplicas, unsigned int neigenvectors, const edm::FileInPath &incsv)

unsigned int	neigenvectors () const

	PDFWeightsHelper ()

Protected Attributes
Eigen::MatrixXd	transformation_

Detailed Description

Definition at line 10 of file PDFWeightsHelper.h.

Constructor & Destructor Documentation

◆ PDFWeightsHelper()

PDFWeightsHelper::PDFWeightsHelper ( )

Definition at line 5 of file PDFWeightsHelper.cc.

5 {}

Member Function Documentation

◆ DoMC2Hessian()

void PDFWeightsHelper::DoMC2Hessian	(	double	nomweight,
		const double *	inweights,
		double *	outweights
	)		const

Definition at line 28 of file PDFWeightsHelper.cc.

                                                                                                        {
   const unsigned int nreplicas = transformation_.rows();
   const unsigned int neigenvectors = transformation_.cols();
  
   Eigen::VectorXd inweightv(nreplicas);
   for (unsigned int irep = 0; irep < nreplicas; ++irep) {
     inweightv[irep] = inweights[irep] - nomweight;
   }
  
   Eigen::VectorXd outweightv = transformation_.transpose() * inweightv;
  
   for (unsigned int ieig = 0; ieig < neigenvectors; ++ieig) {
     outweights[ieig] = outweightv[ieig] + nomweight;
   }
 }

References neigenvectors(), and transformation_.

Referenced by PDFWeightsTest::analyze().

◆ Init()

void PDFWeightsHelper::Init	(	unsigned int	nreplicas,
		unsigned int	neigenvectors,
		const edm::FileInPath &	incsv
	)

Definition at line 7 of file PDFWeightsHelper.cc.

                                                                                                           {
   transformation_.resize(nreplicas, neigenvectors);
  
   std::ifstream instream(incsv.fullPath());
   if (!instream.is_open()) {
     throw cms::Exception("PDFWeightsHelper") << "Could not open csv file" << incsv.relativePath();
   }
  
   for (unsigned int ireplica = 0; ireplica < nreplicas; ++ireplica) {
     std::string linestr;
     getline(instream, linestr);
     std::istringstream line(linestr);
     for (unsigned int ieigen = 0; ieigen < neigenvectors; ++ieigen) {
       std::string valstr;
       getline(line, valstr, ',');
       std::istringstream val(valstr);
       val >> transformation_(ireplica, ieigen);
     }
   }
 }