Inheritance diagram for PDFWeightsTest:

Public Member Functions
	PDFWeightsTest (const ParameterSet &cfg)

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () 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::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	analyze (const Event &evt, const EventSetup &es) override

Private Attributes
unsigned int	nPdfEigWeights_

unsigned int	nPdfWeights_

std::vector< float >	pdfeigweights_

unsigned int	pdfWeightOffset_

std::vector< float >	pdfweights_

PDFWeightsHelper	pdfweightshelper_

EDGetTokenT< LHEEventProduct >	srcToken_

EDGetTokenT< LHEEventProduct >	srcTokenAlt_

EDGetTokenT< GenEventInfoProduct >	srcTokenGen_

TTree *	tree_

float	weight_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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 21 of file PDFWeightsTest.cc.

Constructor & Destructor Documentation

◆ PDFWeightsTest()

PDFWeightsTest::PDFWeightsTest ( const ParameterSet & cfg )

inlineexplicit

Definition at line 38 of file PDFWeightsTest.cc.

       : srcToken_(consumes<LHEEventProduct>(edm::InputTag("externalLHEProducer"))),
         srcTokenAlt_(consumes<LHEEventProduct>(edm::InputTag("source"))),
         srcTokenGen_(consumes<GenEventInfoProduct>(edm::InputTag("generator"))),
         pdfWeightOffset_(cfg.getParameter<unsigned int>("pdfWeightOffset")),
         nPdfWeights_(cfg.getParameter<unsigned int>("nPdfWeights")),
         nPdfEigWeights_(cfg.getParameter<unsigned int>("nPdfEigWeights")),
         pdfweights_(nPdfWeights_),
         pdfeigweights_(nPdfEigWeights_) {
     edm::Service<TFileService> fs;
     tree_ = fs->make<TTree>("tree", "");
  
     tree_->Branch("pdfrep", &pdfweights_);
     tree_->Branch("pdfeig", &pdfeigweights_);
     tree_->Branch("weight", &weight_);
  
     edm::FileInPath mc2hessianCSV = cfg.getParameter<edm::FileInPath>("mc2hessianCSV");
     pdfweightshelper_.Init(nPdfWeights_, nPdfEigWeights_, mc2hessianCSV);
   }

References looper::cfg, PDFWeightsHelper::Init(), and TFileService::make().

Member Function Documentation

◆ analyze()

void PDFWeightsTest::analyze	(	const Event &	evt,
		const EventSetup &	es
	)

inlineoverrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 59 of file PDFWeightsTest.cc.

                                                                 {
     Handle<LHEEventProduct> lheInfo;
     evt.getByToken(srcToken_, lheInfo);
  
     if (!lheInfo.isValid()) {
       evt.getByToken(srcTokenAlt_, lheInfo);
     }
  
     double nomlheweight = lheInfo->weights()[0].wgt;
  
     Handle<GenEventInfoProduct> genInfo;
     evt.getByToken(srcTokenGen_, genInfo);
  
     weight_ = genInfo->weight();
  
     //get the original mc replica weights
     std::vector<double> inpdfweights(nPdfWeights_);
     for (unsigned int ipdf = 0; ipdf < nPdfWeights_; ++ipdf) {
       unsigned int iwgt = ipdf + pdfWeightOffset_;
  
       //this is the weight to be used for evaluating uncertainties with mc replica weights
       pdfweights_[ipdf] = lheInfo->weights()[iwgt].wgt * weight_ / nomlheweight;
  
       //this is the raw weight to be fed to the mc2hessian convertor
       inpdfweights[ipdf] = lheInfo->weights()[iwgt].wgt;
     }
  
     std::vector<double> outpdfweights(nPdfEigWeights_);
     //do the actual conversion, where the nominal lhe weight is needed as the reference point for the linearization
     pdfweightshelper_.DoMC2Hessian(nomlheweight, inpdfweights.data(), outpdfweights.data());
  
     for (unsigned int iwgt = 0; iwgt < nPdfEigWeights_; ++iwgt) {
       double wgtval = outpdfweights[iwgt];
  
       //the is the weight to be used for evaluating uncertainties with hessian weights
       pdfeigweights_[iwgt] = wgtval * weight_ / nomlheweight;
     }
  
     tree_->Fill();
   }

References PDFWeightsHelper::DoMC2Hessian(), edm::Event::getByToken(), genWeightsTable_cfi::lheInfo, and GenEventInfoProduct::weight().

Member Data Documentation

◆ nPdfEigWeights_

unsigned int PDFWeightsTest::nPdfEigWeights_

private

Definition at line 30 of file PDFWeightsTest.cc.

◆ nPdfWeights_

unsigned int PDFWeightsTest::nPdfWeights_

private

Definition at line 29 of file PDFWeightsTest.cc.

◆ pdfeigweights_

std::vector<float> PDFWeightsTest::pdfeigweights_

private

Definition at line 34 of file PDFWeightsTest.cc.

◆ pdfWeightOffset_

unsigned int PDFWeightsTest::pdfWeightOffset_

private

Definition at line 28 of file PDFWeightsTest.cc.

◆ pdfweights_

std::vector<float> PDFWeightsTest::pdfweights_

private

Definition at line 33 of file PDFWeightsTest.cc.

◆ pdfweightshelper_

PDFWeightsHelper PDFWeightsTest::pdfweightshelper_

private

Definition at line 23 of file PDFWeightsTest.cc.

◆ srcToken_

EDGetTokenT<LHEEventProduct> PDFWeightsTest::srcToken_

private

Definition at line 24 of file PDFWeightsTest.cc.

◆ srcTokenAlt_

EDGetTokenT<LHEEventProduct> PDFWeightsTest::srcTokenAlt_

private

Definition at line 25 of file PDFWeightsTest.cc.

◆ srcTokenGen_

EDGetTokenT<GenEventInfoProduct> PDFWeightsTest::srcTokenGen_

private

Definition at line 26 of file PDFWeightsTest.cc.

◆ tree_

TTree* PDFWeightsTest::tree_

private

Definition at line 32 of file PDFWeightsTest.cc.

◆ weight_

float PDFWeightsTest::weight_

private

Definition at line 35 of file PDFWeightsTest.cc.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ PDFWeightsTest()

Member Function Documentation

◆ analyze()

Member Data Documentation

◆ nPdfEigWeights_

◆ nPdfWeights_

◆ pdfeigweights_

◆ pdfWeightOffset_

◆ pdfweights_

◆ pdfweightshelper_

◆ srcToken_

◆ srcTokenAlt_

◆ srcTokenGen_

◆ tree_

◆ weight_