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

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)

virtual	~ProducerBase () noexcept(false)

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_

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

Constructor & Destructor Documentation

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

Definition at line 38 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 46 of file SimpleSystematicsAnalyzer.cc.

46 {}

Member Function Documentation

void SimpleSystematicsAnalyzer::beginJob ( void )

overridevirtual

Reimplemented from edm::EDFilter.

Definition at line 49 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 62 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 103 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;
 }