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

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)

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

Constructor & Destructor Documentation

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

Definition at line 39 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 47 of file PdfSystematicsAnalyzer.cc.

47 {}

Member Function Documentation

void PdfSystematicsAnalyzer::beginJob ( void )

overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 50 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 61 of file PdfSystematicsAnalyzer.cc.

References mps_fire::i, createfilelist::int, 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 221 of file PdfSystematicsAnalyzer.cc.

References edm::HLTGlobalStatus::accept(), DEFINE_FWK_MODULE, edm::Event::getByToken(), mps_fire::i, originalEvents_, pdfStart_, pdfWeightTags_, pdfWeightTokens_, selectedEvents_, selectorPath_, edm::TriggerNames::size(), edm::TriggerNames::triggerIndex(), edm::Event::triggerNames(), triggerResults, triggerResultsToken_, trigNames, weighted2SelectedEvents_, weightedEvents_, weightedSelectedEvents_, and create_public_lumi_plots::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 18 of file PdfSystematicsAnalyzer.cc.

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

std::vector<int> PdfSystematicsAnalyzer::pdfStart_

private

Definition at line 20 of file PdfSystematicsAnalyzer.cc.

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

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

private

Definition at line 15 of file PdfSystematicsAnalyzer.cc.

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

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

private

Definition at line 16 of file PdfSystematicsAnalyzer.cc.

Referenced by filter().

unsigned int PdfSystematicsAnalyzer::selectedEvents_

private

Definition at line 19 of file PdfSystematicsAnalyzer.cc.

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

std::string PdfSystematicsAnalyzer::selectorPath_

private

Definition at line 14 of file PdfSystematicsAnalyzer.cc.

Referenced by filter().

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

private

Definition at line 17 of file PdfSystematicsAnalyzer.cc.

Referenced by filter(), and PdfSystematicsAnalyzer().

std::vector<double> PdfSystematicsAnalyzer::weighted2SelectedEvents_

private

Definition at line 22 of file PdfSystematicsAnalyzer.cc.

Referenced by endJob(), and filter().

std::vector<double> PdfSystematicsAnalyzer::weightedEvents_

private

Definition at line 23 of file PdfSystematicsAnalyzer.cc.

Referenced by endJob(), and filter().

std::vector<double> PdfSystematicsAnalyzer::weightedSelectedEvents_

private

Definition at line 21 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