PdfSystematicsAnalyzer Class Reference

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
std::vector< bool > const &	recordProvenanceList () 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)
TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...
	~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 &&)=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 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	wantsInputProcessBlocks ()
static bool	wantsProcessBlocks ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)
template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)
ProducesCollector	producesCollector ()
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 6 of file PdfSystematicsAnalyzer.cc.

Constructor & Destructor Documentation

PdfSystematicsAnalyzer()

PdfSystematicsAnalyzer::PdfSystematicsAnalyzer

(

const edm::ParameterSet &

pset

)

Definition at line 40 of file PdfSystematicsAnalyzer.cc.

    : 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"))) {}

References GlobalPosition_Frontier_DevDB_cff::tag.

~PdfSystematicsAnalyzer()

PdfSystematicsAnalyzer::~PdfSystematicsAnalyzer ( )

override

Definition at line 48 of file PdfSystematicsAnalyzer.cc.

{}

Member Function Documentation

beginJob()

void PdfSystematicsAnalyzer::beginJob ( void  )

overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 51 of file PdfSystematicsAnalyzer.cc.

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

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

endJob()

void PdfSystematicsAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 62 of file PdfSystematicsAnalyzer.cc.

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

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

filter()

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

overridevirtual

Implements edm::EDFilter.

Definition at line 256 of file PdfSystematicsAnalyzer.cc.

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

References ev, mps_fire::i, dqmiolumiharvest::j, originalEvents_, pdfStart_, pdfWeightTags_, pdfWeightTokens_, selectedEvents_, selectorPath_, triggerResults, triggerResultsToken_, trigNames, weighted2SelectedEvents_, weightedEvents_, weightedSelectedEvents_, and HLT_FULL_cff::weights.

Member Data Documentation

originalEvents_

unsigned int PdfSystematicsAnalyzer::originalEvents_

private

Definition at line 19 of file PdfSystematicsAnalyzer.cc.

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

pdfStart_

std::vector<int> PdfSystematicsAnalyzer::pdfStart_

private

Definition at line 21 of file PdfSystematicsAnalyzer.cc.

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

pdfWeightTags_

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

private

Definition at line 16 of file PdfSystematicsAnalyzer.cc.

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

pdfWeightTokens_

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

private

Definition at line 17 of file PdfSystematicsAnalyzer.cc.

Referenced by filter().

selectedEvents_

unsigned int PdfSystematicsAnalyzer::selectedEvents_

private

Definition at line 20 of file PdfSystematicsAnalyzer.cc.

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

selectorPath_

std::string PdfSystematicsAnalyzer::selectorPath_

private

Definition at line 15 of file PdfSystematicsAnalyzer.cc.

Referenced by filter().

triggerResultsToken_

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

private

Definition at line 18 of file PdfSystematicsAnalyzer.cc.

Referenced by filter().

weighted2SelectedEvents_

std::vector<double> PdfSystematicsAnalyzer::weighted2SelectedEvents_

private

Definition at line 23 of file PdfSystematicsAnalyzer.cc.

Referenced by endJob(), and filter().

weightedEvents_

std::vector<double> PdfSystematicsAnalyzer::weightedEvents_

private

Definition at line 24 of file PdfSystematicsAnalyzer.cc.

Referenced by endJob(), and filter().

weightedSelectedEvents_

std::vector<double> PdfSystematicsAnalyzer::weightedSelectedEvents_

private

Definition at line 22 of file PdfSystematicsAnalyzer.cc.

Referenced by endJob(), and filter().