Inheritance diagram for LHECOMWeightProducer:

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

	~LHECOMWeightProducer ()

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

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

	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	beginJob () override

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

virtual 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
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)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, 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 ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

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

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

Detailed Description

Definition at line 24 of file LHECOMWeightProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 56 of file LHECOMWeightProducer.cc.

References _label, and _newECMS.

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

LHECOMWeightProducer::~LHECOMWeightProducer ( )

Definition at line 67 of file LHECOMWeightProducer.cc.

67 {}

Member Function Documentation

void LHECOMWeightProducer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 86 of file LHECOMWeightProducer.cc.

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

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

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 70 of file LHECOMWeightProducer.cc.

References _newECMS, _origECMS, _pdfmember, _pdfset, Exception, edm::Run::getByLabel(), LHAPDF::initPDFSet(), and lheTag_.

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

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

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 91 of file LHECOMWeightProducer.cc.

References _label, _newECMS, _origECMS, _pdfmember, _pdfset, gather_cfg::cout, edm::Event::getByLabel(), info(), edm::EventBase::isRealData(), lheTag_, LogTrace, edm::Event::put(), class-composition::Q, LHAPDF::usePDFMember(), histoStyle::weight, create_public_pileup_plots::weights, and LHAPDF::xfx().

                                                                          {
 
       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::auto_ptr<GenEventInfoProduct> info(new GenEventInfoProduct());
       info->setWeights(weights);
       iEvent.put(info, _label);
 }

Member Data Documentation

std::string LHECOMWeightProducer::_label

private

Definition at line 39 of file LHECOMWeightProducer.cc.

Referenced by cmsswVersionTools.PickRelValInputFiles::apply(), ElectronCalibrator.EmbeddedElectronCalibrator::correct(), LHECOMWeightProducer(), cmsswVersionTools.PickRelValInputFiles::messageEmptyList(), produce(), trigTools.SwitchOnTrigger::toolCode(), trigTools.SwitchOnTriggerStandAlone::toolCode(), trigTools.SwitchOnTriggerMatching::toolCode(), trigTools.SwitchOnTriggerMatchingStandAlone::toolCode(), and trigTools.SwitchOnTriggerMatchEmbedding::toolCode().

double LHECOMWeightProducer::_newECMS

private

Definition at line 38 of file LHECOMWeightProducer.cc.

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

double LHECOMWeightProducer::_origECMS

private

Definition at line 37 of file LHECOMWeightProducer.cc.

Referenced by beginRun(), and produce().

int LHECOMWeightProducer::_pdfmember

private

Definition at line 36 of file LHECOMWeightProducer.cc.

Referenced by beginRun(), and produce().

int LHECOMWeightProducer::_pdfset

private

Definition at line 35 of file LHECOMWeightProducer.cc.

Referenced by beginRun(), and produce().

edm::InputTag LHECOMWeightProducer::lheTag_

private

Definition at line 34 of file LHECOMWeightProducer.cc.

Referenced by beginRun(), and produce().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation