#include <HLTVertexPerformanceAnalyzer.h>

Inheritance diagram for HLTVertexPerformanceAnalyzer:

Public Member Functions
void	dqmBeginRun (const edm::Run &iRun, const edm::EventSetup &iSetup) override

	HLTVertexPerformanceAnalyzer (const edm::ParameterSet &)

	~HLTVertexPerformanceAnalyzer () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () 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 Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

Private Attributes
std::vector< bool >	_isfoundHLTs

std::vector< std::string >	folders

std::vector< std::map< std::string, MonitorElement * > >	H1_

HLTConfigProvider	hltConfigProvider_

std::vector< int >	hltPathIndexs_

std::vector< std::string >	hltPathNames_

edm::EDGetTokenT< edm::TriggerResults >	hlTriggerResults_

std::string	hlTriggerResults_Label

edm::EDConsumerBase::Labels	label

std::string	mainFolder_

edm::EDGetTokenT< std::vector< SimVertex > >	simVertexCollection_

bool	triggerConfChanged_

std::vector< edm::EDGetTokenT< reco::VertexCollection > >	VertexCollection_

std::vector< std::string >	VertexCollection_Label

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	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::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Code used to produce DQM validation plots for b-tag at HLT. It plots the distribution of recoVertex.z() - simVertex.z() of the primary vertex.

Definition at line 33 of file HLTVertexPerformanceAnalyzer.h.

Constructor & Destructor Documentation

HLTVertexPerformanceAnalyzer::HLTVertexPerformanceAnalyzer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 6 of file HLTVertexPerformanceAnalyzer.cc.

References edm::ParameterSet::getParameter(), mps_fire::i, diffTwoXMLs::label, edm::EDConsumerBase::labelsForToken(), AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, and edm::vector_transform().

 {
     mainFolder_                     = iConfig.getParameter<std::string>("mainFolder");
     hlTriggerResults_           = consumes<TriggerResults>(iConfig.getParameter<InputTag> ("TriggerResults"));
     VertexCollection_           =   edm::vector_transform(iConfig.getParameter<std::vector<edm::InputTag> >( "Vertex" ), [this](edm::InputTag const & tag){return mayConsume< reco::VertexCollection>(tag);});
     hltPathNames_               = iConfig.getParameter< std::vector<std::string> > ("HLTPathNames");
     simVertexCollection_ = consumes<std::vector<SimVertex> >(iConfig.getParameter<edm::InputTag> ("SimVertexCollection"));
 
     EDConsumerBase::labelsForToken(hlTriggerResults_,label);
     hlTriggerResults_Label = label.module;
     
     for(unsigned int i=0; i<VertexCollection_.size() ; i++){
         EDConsumerBase::labelsForToken(VertexCollection_[i],label);
         VertexCollection_Label.push_back(label.module);
     }
 }

HLTVertexPerformanceAnalyzer::~HLTVertexPerformanceAnalyzer ( )

override

Definition at line 24 of file HLTVertexPerformanceAnalyzer.cc.

25 {

26 }

Member Function Documentation

void HLTVertexPerformanceAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 63 of file HLTVertexPerformanceAnalyzer.cc.

References edm::HLTGlobalStatus::accept(), coll, HcalObjRepresent::Fill(), edm::Event::getByToken(), edm::HandleBase::isValid(), RecoTauDiscriminantConfiguration::maxValue, CoreSimVertex::position(), edm::Handle< T >::product(), caloTruthProducer_cfi::simVertexCollection, and edm::triggerResults().

 {
     bool trigRes=false;
     using namespace edm;
 
     //get triggerResults
     Handle<TriggerResults> TriggerResulsHandler;
     if ( hlTriggerResults_Label == "" || hlTriggerResults_Label == "NULL" ) {
         edm::LogInfo("NoTriggerResults") << "TriggerResults ==> Empty";
         return;
     }
     iEvent.getByToken(hlTriggerResults_, TriggerResulsHandler);
     if (TriggerResulsHandler.isValid())   trigRes=true;
     if ( !trigRes ) { edm::LogInfo("NoTriggerResults") << "TriggerResults ==> not readable"; return;}
     const TriggerResults & triggerResults = *(TriggerResulsHandler.product());
 
     //get simVertex
     float simPV=0;
 
     Handle<std::vector<SimVertex> > simVertexCollection;
     iEvent.getByToken(simVertexCollection_, simVertexCollection);
     if (!simVertexCollection.isValid()) { edm::LogInfo("NoSimVertex") << "SimVertex collection is invalid"; return;}
     if (simVertexCollection->empty()) { edm::LogInfo("NoSimVertex") << "SimVertex collection is empty"; return;}
     const SimVertex simPVh = *(simVertexCollection->begin());
     simPV=simPVh.position().z();
     
     //fill the DQM plot
     Handle<VertexCollection> VertexHandler;
     for (unsigned int ind=0; ind<hltPathNames_.size();ind++) {
         for (unsigned int coll=0; coll<VertexCollection_.size();coll++) {
             bool VertexOK=false;
             if ( !_isfoundHLTs[ind]) continue;  //if the hltPath is not in the event, skip the event
             if ( !triggerResults.accept(hltPathIndexs_[ind]) ) continue;    //if the hltPath was not accepted skip the event
             
             //get the recoVertex
             if (VertexCollection_Label.at(coll) != "" && VertexCollection_Label.at(coll) != "NULL" )
             {
                 iEvent.getByToken(VertexCollection_.at(coll), VertexHandler);
                 if (VertexHandler.isValid()>0)   VertexOK=true;
             }
             
             if (VertexOK){
                 //calculate the variable (RecoVertex - SimVertex)
                 float value=VertexHandler->begin()->z()-simPV;
                 
                 //if value is over/under flow, assign the extreme value
                 float maxValue=H1_.at(ind)["Vertex_"+VertexCollection_Label.at(coll)]->getTH1F()->GetXaxis()->GetXmax();
                 if(value>maxValue)  value=maxValue-0.0001; 
                 if(value<-maxValue) value=-maxValue+0.0001; 
                 //fill the histo
                 H1_.at(ind)["Vertex_"+VertexCollection_Label.at(coll)] -> Fill(value);
             }
         }// for on VertexCollection_
     }//for on hltPathNames_
 }

void HLTVertexPerformanceAnalyzer::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 120 of file HLTVertexPerformanceAnalyzer.cc.

References DQMStore::IBooker::book1D(), coll, DEFINE_FWK_MODULE, RecoTauDiscriminantConfiguration::maxValue, DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

 {
     //book the DQM plots
     using namespace std;
     std::string dqmFolder;
     for (unsigned int ind=0; ind<hltPathNames_.size();ind++) {
         dqmFolder = Form("%s/Vertex/%s",mainFolder_.c_str(),hltPathNames_[ind].c_str());
         H1_.push_back(std::map<std::string, MonitorElement *>());
         ibooker.setCurrentFolder(dqmFolder);
         for (unsigned int coll=0; coll<VertexCollection_.size();coll++) {
             float maxValue = 0.02;
             if(VertexCollection_Label.at(coll)==("hltFastPrimaryVertex")) maxValue = 2.; //for the hltFastPrimaryVertex use a larger scale (res ~ 1 cm)
             float vertexU = maxValue;
             float vertexL = -maxValue;
             int   vertexBins = 100;
             if ( VertexCollection_Label.at(coll) != "" && VertexCollection_Label.at(coll) != "NULL" ) { 
                 H1_.back()["Vertex_"+VertexCollection_Label.at(coll)]       = ibooker.book1D("Vertex_"+VertexCollection_Label.at(coll),      VertexCollection_Label.at(coll).c_str(),  vertexBins, vertexL, vertexU );
                 H1_.back()["Vertex_"+VertexCollection_Label.at(coll)]      -> setAxisTitle("vertex error (cm)",1);
             }
         }
     }
 }

void HLTVertexPerformanceAnalyzer::dqmBeginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 29 of file HLTVertexPerformanceAnalyzer.cc.

References runEdmFileComparison::found, createfilelist::int, diffTwoXMLs::label, and edm::EDConsumerBase::labelsForToken().

                                                                                               {
     triggerConfChanged_ = true;
     EDConsumerBase::labelsForToken(hlTriggerResults_,label);
     
     hltConfigProvider_.init(iRun, iSetup, label.process, triggerConfChanged_);
     const std::vector< std::string > & allHltPathNames = hltConfigProvider_.triggerNames();
 
     //fill hltPathIndexs_ with the trigger number of each hltPathNames_
     for ( size_t trgs=0; trgs<hltPathNames_.size(); trgs++) {
         unsigned int found = 1;
         int it_mem = -1;
         for (size_t it=0 ; it < allHltPathNames.size() ; ++it )
         {
             found = allHltPathNames.at(it).find(hltPathNames_[trgs]);
             if ( found == 0 )
             {
                 it_mem= (int) it;
             }
         }//for allallHltPathNames
         hltPathIndexs_.push_back(it_mem);
     }//for hltPathNames_
     
     //fill _isfoundHLTs for each hltPathNames_
     for ( size_t trgs=0; trgs<hltPathNames_.size(); trgs++) {
         if ( hltPathIndexs_[trgs] < 0 ) {
             _isfoundHLTs.push_back(false);
         } 
         else {
             _isfoundHLTs.push_back(true);
         }
     }
 }