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

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_

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

Constructor & Destructor Documentation

◆ SimpleSystematicsAnalyzer()

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

Definition at line 39 of file SimpleSystematicsAnalyzer.cc.

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

References GlobalPosition_Frontier_DevDB_cff::tag.

◆ ~SimpleSystematicsAnalyzer()

SimpleSystematicsAnalyzer::~SimpleSystematicsAnalyzer ( )

override

Definition at line 47 of file SimpleSystematicsAnalyzer.cc.

47 {}

Member Function Documentation

◆ beginJob()

void SimpleSystematicsAnalyzer::beginJob ( void )

overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 50 of file SimpleSystematicsAnalyzer.cc.

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

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

◆ endJob()

void SimpleSystematicsAnalyzer::endJob ( void )

overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 63 of file SimpleSystematicsAnalyzer.cc.

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

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

◆ filter()

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

overridevirtual

Implements edm::EDFilter.

Definition at line 112 of file SimpleSystematicsAnalyzer.cc.

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