SimpleSystematicsAnalyzer Class Reference

Inheritance diagram for SimpleSystematicsAnalyzer:

Public Member Functions
void	beginJob () override
void	endJob () override
bool	filter (edm::Event &, const edm::EventSetup &) override
	SimpleSystematicsAnalyzer (const edm::ParameterSet &pset)
	~SimpleSystematicsAnalyzer () 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_
unsigned int	selectedEvents_
std::string	selectorPath_
edm::EDGetTokenT< edm::TriggerResults >	triggerResultsToken_
std::vector< double >	weighted2SelectedEvents_
std::vector< double >	weightedEvents_
std::vector< double >	weightedSelectedEvents_
std::vector< edm::InputTag >	weightTags_
std::vector< edm::EDGetTokenT< double > >	weightTokens_

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 SimpleSystematicsAnalyzer.cc.

Constructor & Destructor Documentation

SimpleSystematicsAnalyzer::SimpleSystematicsAnalyzer

(

const edm::ParameterSet &

pset

)

Definition at line 39 of file SimpleSystematicsAnalyzer.cc.

References GlobalPosition_Frontier_DevDB_cff::tag, and triggerResultsToken_.

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

SimpleSystematicsAnalyzer::~SimpleSystematicsAnalyzer ( )

override

Definition at line 47 of file SimpleSystematicsAnalyzer.cc.

{}

Member Function Documentation

void SimpleSystematicsAnalyzer::beginJob ( void  )

overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 50 of file SimpleSystematicsAnalyzer.cc.

References mps_fire::i, originalEvents_, selectedEvents_, weighted2SelectedEvents_, weightedEvents_, weightedSelectedEvents_, and weightTags_.

                                         {
  originalEvents_ = 0;
  selectedEvents_ = 0;
  edm::LogVerbatim("SimpleSystematicsAnalysis") << "Uncertainties will be determined for the following tags: ";
  for (unsigned int i = 0; i < weightTags_.size(); ++i) {
    edm::LogVerbatim("SimpleSystematicsAnalysis") << "\t" << weightTags_[i].encode();
    weightedEvents_.push_back(0.);
    weightedSelectedEvents_.push_back(0.);
    weighted2SelectedEvents_.push_back(0.);
  }
}

void SimpleSystematicsAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 63 of file SimpleSystematicsAnalyzer.cc.

References mps_fire::i, originalEvents_, selectedEvents_, mathSSE::sqrt(), weighted2SelectedEvents_, weightedEvents_, weightedSelectedEvents_, weightTags_, and hybridSuperClusters_cfi::xi.

                                       {
  if (originalEvents_ == 0) {
    edm::LogVerbatim("SimpleSystematicsAnalysis") << "NO EVENTS => NO RESULTS";
    return;
  }
  if (selectedEvents_ == 0) {
    edm::LogVerbatim("SimpleSystematicsAnalysis") << "NO SELECTED EVENTS => NO RESULTS";
    return;
  }

  edm::LogVerbatim("SimpleSystematicsAnalysis") << "\n>>>> Begin of Weight systematics summary >>>>";
  edm::LogVerbatim("SimpleSystematicsAnalysis") << "Total number of analyzed data: " << originalEvents_ << " [events]";
  double originalAcceptance = double(selectedEvents_) / originalEvents_;
  edm::LogVerbatim("SimpleSystematicsAnalysis")
      << "Total number of selected data: " << selectedEvents_ << " [events], corresponding to acceptance: ["
      << originalAcceptance * 100 << " +- "
      << 100 * sqrt(originalAcceptance * (1. - originalAcceptance) / originalEvents_) << "] %";

  for (unsigned int i = 0; i < weightTags_.size(); ++i) {
    edm::LogVerbatim("SimpleSystematicsAnalysis") << "Results for Weight Tag: " << weightTags_[i].encode() << " ---->";

    double acc_central = 0.;
    double acc2_central = 0.;
    if (weightedEvents_[i] > 0) {
      acc_central = weightedSelectedEvents_[i] / weightedEvents_[i];
      acc2_central = weighted2SelectedEvents_[i] / weightedEvents_[i];
    }
    double waverage = weightedEvents_[i] / originalEvents_;
    edm::LogVerbatim("SimpleSystematicsAnalysis")
        << "\tTotal Events after reweighting: " << weightedEvents_[i] << " [events]";
    edm::LogVerbatim("SimpleSystematicsAnalysis")
        << "\tEvents selected after reweighting: " << weightedSelectedEvents_[i] << " [events]";
    edm::LogVerbatim("SimpleSystematicsAnalysis")
        << "\tAcceptance after reweighting: [" << 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("SimpleSystematicsAnalysis")
        << "\ti.e. [" << std::setprecision(4) << 100 * xi / originalAcceptance << " +- " << std::setprecision(4)
        << 100 * deltaxi / originalAcceptance << "] % relative variation with respect to the original acceptance";
  }
  edm::LogVerbatim("SimpleSystematicsAnalysis") << ">>>> End of Weight systematics summary >>>>";
}

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

override

Definition at line 112 of file SimpleSystematicsAnalyzer.cc.

References edm::HLTGlobalStatus::accept(), DEFINE_FWK_MODULE, edm::Event::getByToken(), mps_fire::i, originalEvents_, funct::pow(), selectedEvents_, selectorPath_, edm::TriggerNames::size(), edm::TriggerNames::triggerIndex(), edm::Event::triggerNames(), triggerResults, triggerResultsToken_, trigNames, weighted2SelectedEvents_, weightedEvents_, weightedSelectedEvents_, weightTags_, and weightTokens_.

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().

                                                                         {
  originalEvents_++;

  bool selectedEvent = false;
  edm::Handle<edm::TriggerResults> triggerResults;
  if (!ev.getByToken(triggerResultsToken_, triggerResults)) {
    edm::LogError("SimpleSystematicsAnalysis") << ">>> 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("SimpleSystematicsAnalysis") << ">>>> Path Name: " << selectorPath_ << ", selected? " << selectedEvent;

  if (selectedEvent)
    selectedEvents_++;

  for (unsigned int i = 0; i < weightTags_.size(); ++i) {
    edm::Handle<double> weightHandle;
    ev.getByToken(weightTokens_[i], weightHandle);
    weightedEvents_[i] += (*weightHandle);
    if (selectedEvent) {
      weightedSelectedEvents_[i] += (*weightHandle);
      weighted2SelectedEvents_[i] += pow((*weightHandle), 2);
    }
  }

  return true;
}

Member Data Documentation

unsigned int SimpleSystematicsAnalyzer::originalEvents_

private

Definition at line 19 of file SimpleSystematicsAnalyzer.cc.

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

unsigned int SimpleSystematicsAnalyzer::selectedEvents_

private

Definition at line 21 of file SimpleSystematicsAnalyzer.cc.

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

std::string SimpleSystematicsAnalyzer::selectorPath_

private

Definition at line 15 of file SimpleSystematicsAnalyzer.cc.

Referenced by filter().

edm::EDGetTokenT<edm::TriggerResults> SimpleSystematicsAnalyzer::triggerResultsToken_

private

Definition at line 18 of file SimpleSystematicsAnalyzer.cc.

Referenced by filter(), and SimpleSystematicsAnalyzer().

std::vector<double> SimpleSystematicsAnalyzer::weighted2SelectedEvents_

private

Definition at line 23 of file SimpleSystematicsAnalyzer.cc.

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

std::vector<double> SimpleSystematicsAnalyzer::weightedEvents_

private

Definition at line 20 of file SimpleSystematicsAnalyzer.cc.

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

std::vector<double> SimpleSystematicsAnalyzer::weightedSelectedEvents_

private

Definition at line 22 of file SimpleSystematicsAnalyzer.cc.

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

std::vector<edm::InputTag> SimpleSystematicsAnalyzer::weightTags_

private

Definition at line 16 of file SimpleSystematicsAnalyzer.cc.

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

std::vector<edm::EDGetTokenT<double> > SimpleSystematicsAnalyzer::weightTokens_

private

Definition at line 17 of file SimpleSystematicsAnalyzer.cc.

Referenced by filter().