PDFWeightsTest Class Reference

Inheritance diagram for PDFWeightsTest:

Public Member Functions
	PDFWeightsTest (const ParameterSet &cfg)
Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
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 const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
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 const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
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
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 22 of file PDFWeightsTest.cc.

Constructor & Destructor Documentation

PDFWeightsTest::PDFWeightsTest ( const ParameterSet &  cfg )

inlineexplicit

Definition at line 39 of file PDFWeightsTest.cc.

References edm::ParameterSet::getParameter(), PDFWeightsHelper::Init(), and TFileService::make().

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

Member Function Documentation

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

inlineoverrideprivate

Definition at line 63 of file PDFWeightsTest.cc.

References DEFINE_FWK_MODULE, PDFWeightsHelper::DoMC2Hessian(), edm::Event::getByToken(), edm::HandleBase::isValid(), GenEventInfoProduct::weight(), and LHEEventProduct::weights().

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

Member Data Documentation

unsigned int PDFWeightsTest::nPdfEigWeights_

private

Definition at line 31 of file PDFWeightsTest.cc.

unsigned int PDFWeightsTest::nPdfWeights_

private

Definition at line 30 of file PDFWeightsTest.cc.

std::vector<float> PDFWeightsTest::pdfeigweights_

private

Definition at line 35 of file PDFWeightsTest.cc.

unsigned int PDFWeightsTest::pdfWeightOffset_

private

Definition at line 29 of file PDFWeightsTest.cc.

std::vector<float> PDFWeightsTest::pdfweights_

private

Definition at line 34 of file PDFWeightsTest.cc.

PDFWeightsHelper PDFWeightsTest::pdfweightshelper_

private

Definition at line 24 of file PDFWeightsTest.cc.

EDGetTokenT<LHEEventProduct> PDFWeightsTest::srcToken_

private

Definition at line 25 of file PDFWeightsTest.cc.

EDGetTokenT<LHEEventProduct> PDFWeightsTest::srcTokenAlt_

private

Definition at line 26 of file PDFWeightsTest.cc.

EDGetTokenT<GenEventInfoProduct> PDFWeightsTest::srcTokenGen_

private

Definition at line 27 of file PDFWeightsTest.cc.

TTree* PDFWeightsTest::tree_

private

Definition at line 33 of file PDFWeightsTest.cc.

float PDFWeightsTest::weight_

private

Definition at line 36 of file PDFWeightsTest.cc.