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

heppy::PdfWeightProducerTool::PdfWeightProducerTool ( )

inline

Definition at line 18 of file PdfWeightProducerTool.h.

References addPdfSet(), beginJob(), getWeights(), dataset::name, processEvent(), and AlCaHLTBitMon_QueryRunRegistry::string.

{}

Member Function Documentation

void heppy::PdfWeightProducerTool::addPdfSet

(

const std::string &

name

)

Definition at line 18 of file PdfWeightProducerTool.cc.

References dataset::name.

Referenced by PdfWeightProducerTool().

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

void heppy::PdfWeightProducerTool::beginJob

(

void

)

Definition at line 23 of file PdfWeightProducerTool.cc.

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

Referenced by PdfWeightProducerTool().

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

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

(

const std::string &

name

)

const

Definition at line 55 of file PdfWeightProducerTool.cc.

References match().

Referenced by PdfWeightProducerTool().

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

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, gen::n, LHAPDF::numberPDF(), GenEventInfoProduct::pdf(), class-composition::Q, gen::PdfInfo::scalePDF, LHAPDF::usePDFMember(), create_public_lumi_plots::weights, gen::PdfInfo::x, globals_cff::x1, globals_cff::x2, LHAPDF::xfx(), and gen::PdfInfo::xPDF.

Referenced by PdfWeightProducerTool().

                                                                             {
    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

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

private

Definition at line 24 of file PdfWeightProducerTool.h.

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

private

Definition at line 25 of file PdfWeightProducerTool.h.