heppy::PdfWeightProducerTool Class Reference

#include <PdfWeightProducerTool.h>

Public Member Functions
void	addPdfSet (const std::string &name)
void	beginJob ()
const std::vector< double > &	getWeights (const std::string &name) const
	PdfWeightProducerTool ()
void	processEvent (const GenEventInfoProduct &pdfstuff)

Private Attributes
std::vector< std::string >	pdfs_
std::map< std::string, std::vector< double > >	weights_

Detailed Description

Definition at line 16 of file PdfWeightProducerTool.h.

Constructor & Destructor Documentation

PdfWeightProducerTool()

heppy::PdfWeightProducerTool::PdfWeightProducerTool ( )

inline

Definition at line 18 of file PdfWeightProducerTool.h.

{}

Member Function Documentation

addPdfSet()

void heppy::PdfWeightProducerTool::addPdfSet

(

const std::string &

name

)

Definition at line 18 of file PdfWeightProducerTool.cc.

References Skims_PA_cff::name, pdfs_, and weights_.

                                                             {
    pdfs_.push_back(name);
    weights_[name] = std::vector<double>();
  }

beginJob()

void heppy::PdfWeightProducerTool::beginJob

(

void

)

Definition at line 23 of file PdfWeightProducerTool.cc.

References mps_fire::i, LHAPDF::initPDFSet(), dqmiodumpmetadata::n, and pdfs_.

                                       {
    for (unsigned int i = 0, n = pdfs_.size(); i < n; ++i) {
      LHAPDF::initPDFSet(i + 1, pdfs_[i]);
    }
  }

getWeights()

const std::vector< double > & heppy::PdfWeightProducerTool::getWeights

(

const std::string &

name

)

const

Definition at line 56 of file PdfWeightProducerTool.cc.

References cms::cuda::assert(), match(), Skims_PA_cff::name, and weights_.

                                                                                        {
    std::map<std::string, std::vector<double> >::const_iterator match = weights_.find(name);
    assert(match != weights_.end());
    return match->second;
  }

processEvent()

void heppy::PdfWeightProducerTool::processEvent

(

const GenEventInfoProduct &

pdfstuff

)

Definition at line 29 of file PdfWeightProducerTool.cc.

References mps_fire::i, gen::PdfInfo::id, globals_cff::id1, globals_cff::id2, dqmiolumiharvest::j, dqmiodumpmetadata::n, LHAPDF::numberPDF(), GenEventInfoProduct::pdf(), pdfs_, gen::PdfInfo::scalePDF, LHAPDF::usePDFMember(), HLT_2022v12_cff::weights, weights_, gen::PdfInfo::x, testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, LHAPDF::xfx(), and gen::PdfInfo::xPDF.

                                                                              {
    float Q = pdfstuff.pdf()->scalePDF;

    int id1 = pdfstuff.pdf()->id.first;
    double x1 = pdfstuff.pdf()->x.first;
    double pdf1 = pdfstuff.pdf()->xPDF.first;

    int id2 = pdfstuff.pdf()->id.second;
    double x2 = pdfstuff.pdf()->x.second;
    double pdf2 = pdfstuff.pdf()->xPDF.second;

    for (unsigned int i = 0, n = pdfs_.size(); i < n; ++i) {
      std::vector<double> &weights = weights_[pdfs_[i]];
      unsigned int nweights = 1;
      if (LHAPDF::numberPDF(i + 1) > 1)
        nweights += LHAPDF::numberPDF(i + 1);
      weights.resize(nweights);

      for (unsigned int j = 0; j < nweights; ++j) {
        LHAPDF::usePDFMember(i + 1, j);
        double newpdf1 = LHAPDF::xfx(i + 1, x1, Q, id1) / x1;
        double newpdf2 = LHAPDF::xfx(i + 1, x2, Q, id2) / x2;
        weights[j] = newpdf1 / pdf1 * newpdf2 / pdf2;
      }
    }
  }

Member Data Documentation

pdfs_

std::vector<std::string> heppy::PdfWeightProducerTool::pdfs_

private

Definition at line 25 of file PdfWeightProducerTool.h.

Referenced by addPdfSet(), beginJob(), and processEvent().

weights_

std::map<std::string, std::vector<double> > heppy::PdfWeightProducerTool::weights_

private

Definition at line 26 of file PdfWeightProducerTool.h.

Referenced by addPdfSet(), getWeights(), and processEvent().