LHECOMWeightProducer Class Reference

Inheritance diagram for LHECOMWeightProducer:

Public Member Functions
	LHECOMWeightProducer (const edm::ParameterSet &)
	~LHECOMWeightProducer () 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
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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	beginRun (edm::Run const &run, const edm::EventSetup &es) override
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
std::string	_label
double	_newECMS
double	_origECMS
int	_pdfmember
int	_pdfset
edm::InputTag	lheTag_

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	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
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 = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
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 24 of file LHECOMWeightProducer.cc.

Constructor & Destructor Documentation

LHECOMWeightProducer()

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

explicit

Definition at line 56 of file LHECOMWeightProducer.cc.

    : lheTag_(pset.getParameter<edm::InputTag>("lheSrc")), _newECMS(pset.getParameter<double>("NewECMS")) {
  std::stringstream labelStr;
  labelStr << "com"
           << "To" << _newECMS;
  _label = labelStr.str();
  produces<GenEventInfoProduct>(_label);
}

References _label, and _newECMS.

~LHECOMWeightProducer()

LHECOMWeightProducer::~LHECOMWeightProducer ( )

override

Definition at line 66 of file LHECOMWeightProducer.cc.

{}

Member Function Documentation

beginJob()

void LHECOMWeightProducer::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 85 of file LHECOMWeightProducer.cc.

                                    {
  //LHAPDF::initPDFSet(1,pdfSetName_);
}

beginRun()

void LHECOMWeightProducer::beginRun ( edm::Run const &  run, const edm::EventSetup &  es  )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 69 of file LHECOMWeightProducer.cc.

                                                                              {
  using namespace edm;
  Handle<LHERunInfoProduct> lheRun;
  run.getByLabel(lheTag_, lheRun);
  //assumes the same pdf is used for both beams
  _pdfset = lheRun->heprup().PDFSUP.first;
  _pdfmember = lheRun->heprup().PDFGUP.first;
  _origECMS = lheRun->heprup().EBMUP.first + lheRun->heprup().EBMUP.second;
  edm::LogInfo("LHECOMWeightProducer") << "PDFSET: " << _pdfset << "; member: " << _pdfmember
                                       << "; COM energy: " << _origECMS;
  if (_newECMS > _origECMS)
    throw cms::Exception("LHECOMWeightProducer") << "You cannot reweight COM energy to a higher than original energy ";
  LHAPDF::initPDFSet(1, _pdfset, _pdfmember);
}

References _newECMS, _origECMS, _pdfmember, _pdfset, lhef::HEPRUP::EBMUP, Exception, LHERunInfoProduct::heprup(), LHAPDF::initPDFSet(), lheTag_, lhef::HEPRUP::PDFGUP, lhef::HEPRUP::PDFSUP, and writedatasetfile::run.

produce()

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

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 90 of file LHECOMWeightProducer.cc.

                                                                         {
  using namespace std;
  bool verbose = false;
 
  if (iEvent.isRealData())
    return;
 
  edm::Handle<LHEEventProduct> lheevent;
  iEvent.getByLabel(lheTag_, lheevent);
 
  float Q = lheevent->hepeup().SCALUP;
 
  int id1 = lheevent->hepeup().IDUP[0];
  double x1 = fabs(lheevent->hepeup().PUP[0][2] / (_origECMS / 2));
  double x1prime = fabs(lheevent->hepeup().PUP[0][2] / (_newECMS / 2));
 
  int id2 = lheevent->hepeup().IDUP[1];
  double x2 = fabs(lheevent->hepeup().PUP[1][2] / (_origECMS / 2));
  double x2prime = fabs(lheevent->hepeup().PUP[1][2] / (_newECMS / 2));
 
  LogTrace("LHECOMWeightProducer") << "*******LHECOMWeightProducer*******\n"
                                   << " Q  : " << Q << "\n"
                                   << " id1: " << id1 << "\n"
                                   << " x1 : " << x1 << "\n"
                                   << " x1': " << x1prime << "\n"
                                   << " id2: " << id2 << "\n"
                                   << " x2 : " << x2 << "\n"
                                   << " x2': " << x2prime;
  //gluon is 0 in the LHAPDF numbering
  if (id1 == 21)
    id1 = 0;
  if (id2 == 21)
    id2 = 0;
 
  // Put PDF weights in the event
  if (verbose)
    cout << " Set : " << _pdfset << "  member : " << _pdfmember << endl;
 
  LHAPDF::usePDFMember(1, _pdfmember);
  double oldpdf1 = LHAPDF::xfx(1, x1, Q, id1) / x1;
  double oldpdf2 = LHAPDF::xfx(1, x2, Q, id2) / x2;
  double newpdf1 = LHAPDF::xfx(1, x1prime, Q, id1) / x1prime;
  double newpdf2 = LHAPDF::xfx(1, x2prime, Q, id2) / x2prime;
  LogTrace("LHECOMWeightProducer") << "     xfx1 : " << oldpdf1 << "\n"
                                   << "     xfx2 : " << oldpdf2 << "\n"
                                   << "     xfx1': " << newpdf1 << "\n"
                                   << "     xfx2': " << newpdf2 << "\n"
                                   << "     weight:" << (newpdf1 / oldpdf1) * (newpdf2 / oldpdf2);
  double weight = (newpdf1 / oldpdf1) * (newpdf2 / oldpdf2);
  std::vector<double> weights;
  weights.push_back(weight);
  std::unique_ptr<GenEventInfoProduct> info(new GenEventInfoProduct());
  info->setWeights(weights);
  iEvent.put(std::move(info), _label);
}

References _label, _newECMS, _origECMS, _pdfmember, _pdfset, gather_cfg::cout, LHEEventProduct::hepeup(), globals_cff::id1, globals_cff::id2, lhef::HEPEUP::IDUP, iEvent, info(), lheTag_, LogTrace, eostools::move(), lhef::HEPEUP::PUP, class-composition::Q, lhef::HEPEUP::SCALUP, LHAPDF::usePDFMember(), verbose, HLT_FULL_cff::weights, testProducerWithPsetDescEmpty_cfi::x1, testProducerWithPsetDescEmpty_cfi::x2, and LHAPDF::xfx().

Member Data Documentation

_label

std::string LHECOMWeightProducer::_label

private

Definition at line 39 of file LHECOMWeightProducer.cc.

Referenced by cmsswVersionTools.PickRelValInputFiles::apply(), validation.SimpleSample::digest(), ConfigToolBase.ConfigToolBase::dumpPython(), validation.SimpleSample::label(), LHECOMWeightProducer(), cmsswVersionTools.PickRelValInputFiles::messageEmptyList(), produce(), trigTools.SwitchOnTrigger::toolCode(), trigTools.SwitchOnTriggerStandAlone::toolCode(), trigTools.SwitchOnTriggerMatching::toolCode(), trigTools.SwitchOnTriggerMatchingStandAlone::toolCode(), and trigTools.SwitchOnTriggerMatchEmbedding::toolCode().

_newECMS

double LHECOMWeightProducer::_newECMS

private

Definition at line 38 of file LHECOMWeightProducer.cc.

Referenced by beginRun(), LHECOMWeightProducer(), and produce().

_origECMS

double LHECOMWeightProducer::_origECMS

private

Definition at line 37 of file LHECOMWeightProducer.cc.

Referenced by beginRun(), and produce().

_pdfmember

int LHECOMWeightProducer::_pdfmember

private

Definition at line 36 of file LHECOMWeightProducer.cc.

Referenced by beginRun(), and produce().

_pdfset

int LHECOMWeightProducer::_pdfset

private

Definition at line 35 of file LHECOMWeightProducer.cc.

Referenced by beginRun(), and produce().

lheTag_

edm::InputTag LHECOMWeightProducer::lheTag_

private

Definition at line 34 of file LHECOMWeightProducer.cc.

Referenced by beginRun(), and produce().