Inheritance diagram for PdfWeightProducer:

Public Member Functions
	PdfWeightProducer (const edm::ParameterSet &)

	~PdfWeightProducer () override

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDProducer () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

std::vector< bool > const &	recordProvenanceList () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	beginJob () override

void	endJob () override

void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
std::string	fixPOWHEG_

edm::InputTag	genTag_

edm::EDGetTokenT< reco::GenParticleCollection >	genToken_

edm::InputTag	pdfInfoTag_

edm::EDGetTokenT< GenEventInfoProduct >	pdfInfoToken_

std::vector< std::string >	pdfSetNames_

std::vector< std::string >	pdfShortNames_

bool	useFirstAsDefault_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex > >

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 21 of file PdfWeightProducer.cc.

Constructor & Destructor Documentation

◆ PdfWeightProducer()

PdfWeightProducer::PdfWeightProducer ( const edm::ParameterSet & pset )

explicit

Definition at line 54 of file PdfWeightProducer.cc.

     : fixPOWHEG_(pset.getUntrackedParameter<std::string>("FixPOWHEG", "")),
       useFirstAsDefault_(pset.getUntrackedParameter<bool>("useFirstAsDefault", false)),
       genTag_(pset.getUntrackedParameter<edm::InputTag>("GenTag", edm::InputTag("genParticles"))),
       genToken_(mayConsume<reco::GenParticleCollection>(genTag_)),
       pdfInfoTag_(pset.getUntrackedParameter<edm::InputTag>("PdfInfoTag", edm::InputTag("generator"))),
       pdfInfoToken_(consumes<GenEventInfoProduct>(pdfInfoTag_)),
       pdfSetNames_(pset.getUntrackedParameter<std::vector<std::string> >("PdfSetNames")) {
   if (!fixPOWHEG_.empty())
     pdfSetNames_.insert(pdfSetNames_.begin(), fixPOWHEG_);
  
   if (pdfSetNames_.size() > 3) {
     edm::LogWarning("") << pdfSetNames_.size()
                         << " PDF sets requested on input. Using only the first 3 sets and ignoring the rest!!";
     pdfSetNames_.erase(pdfSetNames_.begin() + 3, pdfSetNames_.end());
   }
  
   for (unsigned int k = 0; k < pdfSetNames_.size(); k++) {
     size_t dot = pdfSetNames_[k].find_first_of('.');
     size_t underscore = pdfSetNames_[k].find_first_of('_');
     if (underscore < dot) {
       pdfShortNames_.push_back(pdfSetNames_[k].substr(0, underscore));
     } else {
       pdfShortNames_.push_back(pdfSetNames_[k].substr(0, dot));
     }
     produces<std::vector<double> >(pdfShortNames_[k]);
   }
 }

References dot(), fixPOWHEG_, dqmdumpme::k, pdfSetNames_, and pdfShortNames_.

◆ ~PdfWeightProducer()

PdfWeightProducer::~PdfWeightProducer ( )

override

Definition at line 84 of file PdfWeightProducer.cc.

84 {}

Member Function Documentation

◆ beginJob()

void PdfWeightProducer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 87 of file PdfWeightProducer.cc.

                                  {
   for (unsigned int k = 1; k <= pdfSetNames_.size(); k++) {
     LHAPDF::initPDFSet(k, pdfSetNames_[k - 1]);
   }
 }

References LHAPDF::initPDFSet(), dqmdumpme::k, and pdfSetNames_.

◆ endJob()

void PdfWeightProducer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 94 of file PdfWeightProducer.cc.

94 {}

◆ produce()

void PdfWeightProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &
	)

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 97 of file PdfWeightProducer.cc.

                                                                       {
   if (iEvent.isRealData())
     return;
  
   edm::Handle<GenEventInfoProduct> pdfstuff;
   if (!iEvent.getByToken(pdfInfoToken_, pdfstuff)) {
     edm::LogError("PDFWeightProducer") << ">>> PdfInfo not found: " << pdfInfoTag_.encode() << " !!!";
     return;
   }
  
   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;
   if (useFirstAsDefault_ && pdf1 == -1. && pdf2 == -1.) {
     LHAPDF::usePDFMember(1, 0);
     pdf1 = LHAPDF::xfx(1, x1, Q, id1) / x1;
     pdf2 = LHAPDF::xfx(1, x2, Q, id2) / x2;
   }
  
   // Ad-hoc fix for POWHEG
   if (!fixPOWHEG_.empty()) {
     edm::Handle<reco::GenParticleCollection> genParticles;
     if (!iEvent.getByToken(genToken_, genParticles)) {
       edm::LogError("PDFWeightProducer") << ">>> genParticles  not found: " << genTag_.encode() << " !!!";
       return;
     }
     unsigned int gensize = genParticles->size();
     double mboson = 0.;
     for (unsigned int i = 0; i < gensize; ++i) {
       const reco::GenParticle& part = (*genParticles)[i];
       int status = part.status();
       if (status != 3)
         continue;
       int id = part.pdgId();
       if (id != 23 && abs(id) != 24)
         continue;
       mboson = part.mass();
       break;
     }
     Q = sqrt(mboson * mboson + Q * Q);
     LHAPDF::usePDFMember(1, 0);
     pdf1 = LHAPDF::xfx(1, x1, Q, id1) / x1;
     pdf2 = LHAPDF::xfx(1, x2, Q, id2) / x2;
   }
  
   // Put PDF weights in the event
   for (unsigned int k = 1; k <= pdfSetNames_.size(); ++k) {
     std::unique_ptr<std::vector<double> > weights(new std::vector<double>);
     unsigned int nweights = 1;
     if (LHAPDF::numberPDF(k) > 1)
       nweights += LHAPDF::numberPDF(k);
     weights->reserve(nweights);
  
     for (unsigned int i = 0; i < nweights; ++i) {
       LHAPDF::usePDFMember(k, i);
       double newpdf1 = LHAPDF::xfx(k, x1, Q, id1) / x1;
       double newpdf2 = LHAPDF::xfx(k, x2, Q, id2) / x2;
       weights->push_back(newpdf1 / pdf1 * newpdf2 / pdf2);
     }
     iEvent.put(std::move(weights), pdfShortNames_[k - 1]);
   }
 }

References funct::abs(), edm::InputTag::encode(), fixPOWHEG_, genParticles2HepMC_cfi::genParticles, genTag_, genToken_, mps_fire::i, gen::PdfInfo::id, globals_cff::id1, globals_cff::id2, iEvent, dqmdumpme::k, eostools::move(), LHAPDF::numberPDF(), GenEventInfoProduct::pdf(), pdfInfoTag_, pdfInfoToken_, pdfSetNames_, pdfShortNames_, class-composition::Q, gen::PdfInfo::scalePDF, mathSSE::sqrt(), mps_update::status, useFirstAsDefault_, LHAPDF::usePDFMember(), HLT_2018_cff::weights, gen::PdfInfo::x, testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, LHAPDF::xfx(), and gen::PdfInfo::xPDF.

Member Data Documentation

◆ fixPOWHEG_

std::string PdfWeightProducer::fixPOWHEG_

private

Definition at line 31 of file PdfWeightProducer.cc.

Referenced by PdfWeightProducer(), and produce().

◆ genTag_

edm::InputTag PdfWeightProducer::genTag_

private

Definition at line 33 of file PdfWeightProducer.cc.

Referenced by produce().

◆ genToken_

edm::EDGetTokenT<reco::GenParticleCollection> PdfWeightProducer::genToken_

private

Definition at line 34 of file PdfWeightProducer.cc.

Referenced by produce().

◆ pdfInfoTag_

edm::InputTag PdfWeightProducer::pdfInfoTag_

private

Definition at line 35 of file PdfWeightProducer.cc.

Referenced by produce().

◆ pdfInfoToken_

edm::EDGetTokenT<GenEventInfoProduct> PdfWeightProducer::pdfInfoToken_

private

Definition at line 36 of file PdfWeightProducer.cc.

Referenced by produce().

◆ pdfSetNames_

std::vector<std::string> PdfWeightProducer::pdfSetNames_

private

Definition at line 37 of file PdfWeightProducer.cc.

Referenced by beginJob(), PdfWeightProducer(), and produce().

◆ pdfShortNames_

std::vector<std::string> PdfWeightProducer::pdfShortNames_

private

Definition at line 38 of file PdfWeightProducer.cc.

Referenced by PdfWeightProducer(), and produce().

◆ useFirstAsDefault_

bool PdfWeightProducer::useFirstAsDefault_

private

Definition at line 32 of file PdfWeightProducer.cc.

Referenced by produce().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ PdfWeightProducer()

◆ ~PdfWeightProducer()

Member Function Documentation

◆ beginJob()

◆ endJob()

◆ produce()

Member Data Documentation

◆ fixPOWHEG_

◆ genTag_

◆ genToken_

◆ pdfInfoTag_

◆ pdfInfoToken_

◆ pdfSetNames_

◆ pdfShortNames_

◆ useFirstAsDefault_