Inheritance diagram for PdfSystematicsAnalyzer:

Public Member Functions
void	beginJob () override

void	endJob () override

bool	filter (edm::Event &, const edm::EventSetup &) override

	PdfSystematicsAnalyzer (const edm::ParameterSet &pset)

	~PdfSystematicsAnalyzer () override

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

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDFilter () 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

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)

	~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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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 Attributes
unsigned int	originalEvents_

std::vector< int >	pdfStart_

std::vector< edm::InputTag >	pdfWeightTags_

std::vector< edm::EDGetTokenT< std::vector< double > > >	pdfWeightTokens_

unsigned int	selectedEvents_

std::string	selectorPath_

edm::EDGetTokenT< edm::TriggerResults >	triggerResultsToken_

std::vector< double >	weighted2SelectedEvents_

std::vector< double >	weightedEvents_

std::vector< double >	weightedSelectedEvents_

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	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::EDFilter
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

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

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 6 of file PdfSystematicsAnalyzer.cc.

Constructor & Destructor Documentation

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

Definition at line 40 of file PdfSystematicsAnalyzer.cc.

References GlobalPosition_Frontier_DevDB_cff::tag, and triggerResultsToken_.

     : selectorPath_(pset.getUntrackedParameter<std::string>("SelectorPath", "")),
       pdfWeightTags_(pset.getUntrackedParameter<std::vector<edm::InputTag> >("PdfWeightTags")),
       pdfWeightTokens_(edm::vector_transform(
           pdfWeightTags_, [this](edm::InputTag const& tag) { return consumes<std::vector<double> >(tag); })),
       triggerResultsToken_(consumes<edm::TriggerResults>(edm::InputTag("TriggerResults"))) {}

PdfSystematicsAnalyzer::~PdfSystematicsAnalyzer ( )

override

Definition at line 48 of file PdfSystematicsAnalyzer.cc.

48 {}

Member Function Documentation

void PdfSystematicsAnalyzer::beginJob ( void )

overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 51 of file PdfSystematicsAnalyzer.cc.

References mps_fire::i, originalEvents_, pdfStart_, pdfWeightTags_, and selectedEvents_.

                                       {
   originalEvents_ = 0;
   selectedEvents_ = 0;
   edm::LogVerbatim("PDFAnalysis") << "PDF uncertainties will be determined for the following sets: ";
   for (unsigned int i = 0; i < pdfWeightTags_.size(); ++i) {
     edm::LogVerbatim("PDFAnalysis") << "\t" << pdfWeightTags_[i].instance();
     pdfStart_.push_back(-1);
   }
 }

void PdfSystematicsAnalyzer::endJob ( void )

overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 62 of file PdfSystematicsAnalyzer.cc.

References mps_fire::i, createfilelist::int, dqmiolumiharvest::j, originalEvents_, pdfStart_, pdfWeightTags_, selectedEvents_, mathSSE::sqrt(), weighted2SelectedEvents_, weightedEvents_, weightedSelectedEvents_, and hybridSuperClusters_cfi::xi.

                                     {
   if (originalEvents_ == 0) {
     edm::LogVerbatim("PDFAnalysis") << "NO EVENTS => NO RESULTS";
     return;
   }
   if (selectedEvents_ == 0) {
     edm::LogVerbatim("PDFAnalysis") << "NO SELECTED EVENTS => NO RESULTS";
     return;
   }
 
   edm::LogVerbatim("PDFAnalysis") << "\n>>>> Begin of PDF weight systematics summary >>>>";
   edm::LogVerbatim("PDFAnalysis") << "Total number of analyzed data: " << originalEvents_ << " [events]";
   double originalAcceptance = double(selectedEvents_) / originalEvents_;
   edm::LogVerbatim("PDFAnalysis") << "Total number of selected data: " << selectedEvents_
                                   << " [events], corresponding to acceptance: [" << originalAcceptance * 100 << " +- "
                                   << 100 * sqrt(originalAcceptance * (1. - originalAcceptance) / originalEvents_)
                                   << "] %";
 
   edm::LogVerbatim("PDFAnalysis") << "\n>>>>> PDF UNCERTAINTIES ON RATE >>>>>>";
   for (unsigned int i = 0; i < pdfWeightTags_.size(); ++i) {
     bool nnpdfFlag = (pdfWeightTags_[i].instance().substr(0, 5) == "NNPDF");
     unsigned int nmembers = weightedSelectedEvents_.size() - pdfStart_[i];
     if (i < pdfWeightTags_.size() - 1)
       nmembers = pdfStart_[i + 1] - pdfStart_[i];
     unsigned int npairs = (nmembers - 1) / 2;
     edm::LogVerbatim("PDFAnalysis") << "RATE Results for PDF set " << pdfWeightTags_[i].instance() << " ---->";
 
     double events_central = weightedSelectedEvents_[pdfStart_[i]];
     edm::LogVerbatim("PDFAnalysis") << "\tEstimate for central PDF member: " << int(events_central) << " [events]";
     double events2_central = weighted2SelectedEvents_[pdfStart_[i]];
     edm::LogVerbatim("PDFAnalysis")
         << "\ti.e. [" << std::setprecision(4) << 100 * (events_central - selectedEvents_) / selectedEvents_ << " +- "
         << 100 * sqrt(events2_central - events_central + selectedEvents_ * (1 - originalAcceptance)) / selectedEvents_
         << "] % relative variation with respect to original PDF";
 
     if (npairs > 0) {
       edm::LogVerbatim("PDFAnalysis") << "\tNumber of eigenvectors for uncertainty estimation: " << npairs;
       double wplus = 0.;
       double wminus = 0.;
       unsigned int nplus = 0;
       unsigned int nminus = 0;
       for (unsigned int j = 0; j < npairs; ++j) {
         double wa = weightedSelectedEvents_[pdfStart_[i] + 2 * j + 1] / events_central - 1.;
         double wb = weightedSelectedEvents_[pdfStart_[i] + 2 * j + 2] / events_central - 1.;
         if (nnpdfFlag) {
           if (wa > 0.) {
             wplus += wa * wa;
             nplus++;
           } else {
             wminus += wa * wa;
             nminus++;
           }
           if (wb > 0.) {
             wplus += wb * wb;
             nplus++;
           } else {
             wminus += wb * wb;
             nminus++;
           }
         } else {
           if (wa > wb) {
             if (wa < 0.)
               wa = 0.;
             if (wb > 0.)
               wb = 0.;
             wplus += wa * wa;
             wminus += wb * wb;
           } else {
             if (wb < 0.)
               wb = 0.;
             if (wa > 0.)
               wa = 0.;
             wplus += wb * wb;
             wminus += wa * wa;
           }
         }
       }
       if (wplus > 0)
         wplus = sqrt(wplus);
       if (wminus > 0)
         wminus = sqrt(wminus);
       if (nnpdfFlag) {
         if (nplus > 0)
           wplus /= sqrt(nplus);
         if (nminus > 0)
           wminus /= sqrt(nminus);
       }
       edm::LogVerbatim("PDFAnalysis") << "\tRelative uncertainty with respect to central member: +"
                                       << std::setprecision(4) << 100. * wplus << " / -" << std::setprecision(4)
                                       << 100. * wminus << " [%]";
     } else {
       edm::LogVerbatim("PDFAnalysis") << "\tNO eigenvectors for uncertainty estimation";
     }
   }
 
   edm::LogVerbatim("PDFAnalysis") << "\n>>>>> PDF UNCERTAINTIES ON ACCEPTANCE >>>>>>";
   for (unsigned int i = 0; i < pdfWeightTags_.size(); ++i) {
     bool nnpdfFlag = (pdfWeightTags_[i].instance().substr(0, 5) == "NNPDF");
     unsigned int nmembers = weightedEvents_.size() - pdfStart_[i];
     if (i < pdfWeightTags_.size() - 1)
       nmembers = pdfStart_[i + 1] - pdfStart_[i];
     unsigned int npairs = (nmembers - 1) / 2;
     edm::LogVerbatim("PDFAnalysis") << "ACCEPTANCE Results for PDF set " << pdfWeightTags_[i].instance() << " ---->";
 
     double acc_central = 0.;
     double acc2_central = 0.;
     if (weightedEvents_[pdfStart_[i]] > 0) {
       acc_central = weightedSelectedEvents_[pdfStart_[i]] / weightedEvents_[pdfStart_[i]];
       acc2_central = weighted2SelectedEvents_[pdfStart_[i]] / weightedEvents_[pdfStart_[i]];
     }
     double waverage = weightedEvents_[pdfStart_[i]] / originalEvents_;
     edm::LogVerbatim("PDFAnalysis") << "\tEstimate for central PDF member acceptance: [" << acc_central * 100 << " +- "
                                     << 100 *
                                            sqrt((acc2_central / waverage - acc_central * acc_central) / originalEvents_)
                                     << "] %";
     double xi = acc_central - originalAcceptance;
     double deltaxi = (acc2_central - (originalAcceptance + 2 * xi + xi * xi)) / originalEvents_;
     if (deltaxi > 0)
       deltaxi = sqrt(deltaxi);  //else deltaxi = 0.;
     edm::LogVerbatim("PDFAnalysis") << "\ti.e. [" << std::setprecision(4) << 100 * xi / originalAcceptance << " +- "
                                     << std::setprecision(4) << 100 * deltaxi / originalAcceptance
                                     << "] % relative variation with respect to the original PDF";
 
     if (npairs > 0) {
       edm::LogVerbatim("PDFAnalysis") << "\tNumber of eigenvectors for uncertainty estimation: " << npairs;
       double wplus = 0.;
       double wminus = 0.;
       unsigned int nplus = 0;
       unsigned int nminus = 0;
       for (unsigned int j = 0; j < npairs; ++j) {
         double wa = 0.;
         if (weightedEvents_[pdfStart_[i] + 2 * j + 1] > 0)
           wa = (weightedSelectedEvents_[pdfStart_[i] + 2 * j + 1] / weightedEvents_[pdfStart_[i] + 2 * j + 1]) /
                    acc_central -
                1.;
         double wb = 0.;
         if (weightedEvents_[pdfStart_[i] + 2 * j + 2] > 0)
           wb = (weightedSelectedEvents_[pdfStart_[i] + 2 * j + 2] / weightedEvents_[pdfStart_[i] + 2 * j + 2]) /
                    acc_central -
                1.;
         if (nnpdfFlag) {
           if (wa > 0.) {
             wplus += wa * wa;
             nplus++;
           } else {
             wminus += wa * wa;
             nminus++;
           }
           if (wb > 0.) {
             wplus += wb * wb;
             nplus++;
           } else {
             wminus += wb * wb;
             nminus++;
           }
         } else {
           if (wa > wb) {
             if (wa < 0.)
               wa = 0.;
             if (wb > 0.)
               wb = 0.;
             wplus += wa * wa;
             wminus += wb * wb;
           } else {
             if (wb < 0.)
               wb = 0.;
             if (wa > 0.)
               wa = 0.;
             wplus += wb * wb;
             wminus += wa * wa;
           }
         }
       }
       if (wplus > 0)
         wplus = sqrt(wplus);
       if (wminus > 0)
         wminus = sqrt(wminus);
       if (nnpdfFlag) {
         if (nplus > 0)
           wplus /= sqrt(nplus);
         if (nminus > 0)
           wminus /= sqrt(nminus);
       }
       edm::LogVerbatim("PDFAnalysis") << "\tRelative uncertainty with respect to central member: +"
                                       << std::setprecision(4) << 100. * wplus << " / -" << std::setprecision(4)
                                       << 100. * wminus << " [%]";
     } else {
       edm::LogVerbatim("PDFAnalysis") << "\tNO eigenvectors for uncertainty estimation";
     }
   }
   edm::LogVerbatim("PDFAnalysis") << ">>>> End of PDF weight systematics summary >>>>";
 }

bool PdfSystematicsAnalyzer::filter	(	edm::Event &	ev,
		const edm::EventSetup &
	)

override

Definition at line 256 of file PdfSystematicsAnalyzer.cc.

References edm::HLTGlobalStatus::accept(), DEFINE_FWK_MODULE, edm::Event::getByToken(), mps_fire::i, dqmiolumiharvest::j, originalEvents_, pdfStart_, pdfWeightTags_, pdfWeightTokens_, selectedEvents_, selectorPath_, edm::TriggerNames::size(), edm::TriggerNames::triggerIndex(), edm::Event::triggerNames(), triggerResults, triggerResultsToken_, trigNames, weighted2SelectedEvents_, weightedEvents_, weightedSelectedEvents_, and HGCalRecHit_cfi::weights.

Referenced by Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::filter(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::setDataAccessor(), and Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::switchCenterView().

                                                                       {
   edm::Handle<std::vector<double> > weightHandle;
   for (unsigned int i = 0; i < pdfWeightTags_.size(); ++i) {
     if (!ev.getByToken(pdfWeightTokens_[i], weightHandle)) {
       if (originalEvents_ == 0) {
         edm::LogError("PDFAnalysis") << ">>> WARNING: some weights not found!";
         edm::LogError("PDFAnalysis") << ">>> But maybe OK, if you are prefiltering!";
         edm::LogError("PDFAnalysis") << ">>> If things are OK, this warning should disappear after a while!";
       }
       return false;
     }
   }
 
   originalEvents_++;
 
   bool selectedEvent = false;
   edm::Handle<edm::TriggerResults> triggerResults;
   if (!ev.getByToken(triggerResultsToken_, triggerResults)) {
     edm::LogError("PDFAnalysis") << ">>> TRIGGER collection does not exist !!!";
     return false;
   }
 
   const edm::TriggerNames& trigNames = ev.triggerNames(*triggerResults);
   unsigned int pathIndex = trigNames.triggerIndex(selectorPath_);
   bool pathFound = (pathIndex < trigNames.size());  // pathIndex >= 0, since pathIndex is unsigned
   if (pathFound) {
     if (triggerResults->accept(pathIndex))
       selectedEvent = true;
   }
   //edm::LogVerbatim("PDFAnalysis") << ">>>> Path Name: " << selectorPath_ << ", selected? " << selectedEvent;
 
   if (selectedEvent)
     selectedEvents_++;
 
   for (unsigned int i = 0; i < pdfWeightTags_.size(); ++i) {
     if (!ev.getByToken(pdfWeightTokens_[i], weightHandle))
       return false;
     std::vector<double> weights = (*weightHandle);
     unsigned int nmembers = weights.size();
     // Set up arrays the first time wieghts are read
     if (pdfStart_[i] < 0) {
       pdfStart_[i] = weightedEvents_.size();
       for (unsigned int j = 0; j < nmembers; ++j) {
         weightedEvents_.push_back(0.);
         weightedSelectedEvents_.push_back(0.);
         weighted2SelectedEvents_.push_back(0.);
       }
     }
 
     for (unsigned int j = 0; j < nmembers; ++j) {
       weightedEvents_[pdfStart_[i] + j] += weights[j];
       if (selectedEvent) {
         weightedSelectedEvents_[pdfStart_[i] + j] += weights[j];
         weighted2SelectedEvents_[pdfStart_[i] + j] += weights[j] * weights[j];
       }
     }
 
     /*
             printf("\n>>>>>>>>> Run %8d Event %d, members %3d PDF set %s : Weights >>>> \n", ev.id().run(), ev.id().event(), nmembers, pdfWeightTags_[i].instance().data());
             for (unsigned int i=0; i<nmembers; i+=5) {
                   for (unsigned int j=0; ((j<5)&&(i+j<nmembers)); ++j) {
                         printf(" %2d: %7.4f", i+j, weights[i+j]);
                   }
                   safe_printf("\n");
             }
             */
   }
 
   return true;
 }

Member Data Documentation

unsigned int PdfSystematicsAnalyzer::originalEvents_

private

Definition at line 19 of file PdfSystematicsAnalyzer.cc.

Referenced by beginJob(), endJob(), and filter().

std::vector<int> PdfSystematicsAnalyzer::pdfStart_

private

Definition at line 21 of file PdfSystematicsAnalyzer.cc.

Referenced by beginJob(), endJob(), and filter().

std::vector<edm::InputTag> PdfSystematicsAnalyzer::pdfWeightTags_

private

Definition at line 16 of file PdfSystematicsAnalyzer.cc.

Referenced by beginJob(), endJob(), and filter().

std::vector<edm::EDGetTokenT<std::vector<double> > > PdfSystematicsAnalyzer::pdfWeightTokens_

private

Definition at line 17 of file PdfSystematicsAnalyzer.cc.

Referenced by filter().

unsigned int PdfSystematicsAnalyzer::selectedEvents_

private

Definition at line 20 of file PdfSystematicsAnalyzer.cc.

Referenced by beginJob(), endJob(), and filter().

std::string PdfSystematicsAnalyzer::selectorPath_

private

Definition at line 15 of file PdfSystematicsAnalyzer.cc.

Referenced by filter().

edm::EDGetTokenT<edm::TriggerResults> PdfSystematicsAnalyzer::triggerResultsToken_

private

Definition at line 18 of file PdfSystematicsAnalyzer.cc.

Referenced by filter(), and PdfSystematicsAnalyzer().

std::vector<double> PdfSystematicsAnalyzer::weighted2SelectedEvents_

private

Definition at line 23 of file PdfSystematicsAnalyzer.cc.

Referenced by endJob(), and filter().

std::vector<double> PdfSystematicsAnalyzer::weightedEvents_

private

Definition at line 24 of file PdfSystematicsAnalyzer.cc.

Referenced by endJob(), and filter().

std::vector<double> PdfSystematicsAnalyzer::weightedSelectedEvents_

private

Definition at line 22 of file PdfSystematicsAnalyzer.cc.

Referenced by endJob(), and filter().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation