#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 ( )

Definition at line 5 of file PDFWeightsHelper.cc.

                                    {
 
 }

Member Function Documentation

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

Definition at line 35 of file PDFWeightsHelper.cc.

References neigenvectors(), and transformation_.

Referenced by PDFWeightsTest::analyze().

                                                                                                        {
   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;
   }  
     
 }

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

Definition at line 9 of file PDFWeightsHelper.cc.

References Exception, edm::FileInPath::fullPath(), geometryCSVtoXML::line, neigenvectors(), edm::FileInPath::relativePath(), AlCaHLTBitMon_QueryRunRegistry::string, transformation_, and heppy_batch::val.

Referenced by PDFWeightsTest::PDFWeightsTest().

                                                                                                           {  
   
   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);
     }
     
   }
   
 }