#include <HLTVertexPerformanceAnalyzer.h>

Inheritance diagram for HLTVertexPerformanceAnalyzer:

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

	HLTVertexPerformanceAnalyzer (const edm::ParameterSet &)

	~HLTVertexPerformanceAnalyzer ()

Public Member Functions inherited from thread_unsafe::DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

	DQMEDAnalyzer (void)

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

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< TriggerResults >	hlTriggerResults_

std::string	hlTriggerResults_Label

edm::EDConsumerBase::Labels	label

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::EDAnalyzer
typedef EDAnalyzer	ModuleType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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

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 35 of file HLTVertexPerformanceAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 3 of file HLTVertexPerformanceAnalyzer.cc.

References edm::ParameterSet::getParameter(), hlTriggerResults_, GlobalPosition_Frontier_DevDB_cff::tag, edm::vector_transform(), and VertexCollection_.

 {
     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> >(edm::InputTag("g4SimHits"));
 
     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 ( )

Definition at line 20 of file HLTVertexPerformanceAnalyzer.cc.

21 {

22 }

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 59 of file HLTVertexPerformanceAnalyzer.cc.

References _isfoundHLTs, coll, HcalObjRepresent::Fill(), edm::Event::getByToken(), H1_, hltPathIndexs_, hltPathNames_, hlTriggerResults_, hlTriggerResults_Label, edm::errors::LogicError, CoreSimVertex::position(), cms::Exception::raise(), simVertexCollection_, triggerResults, relativeConstraints::value, VertexCollection_, and VertexCollection_Label.

 {
     bool trigRes=false;
     using namespace edm;
 
     //get triggerResults
     Handle<TriggerResults> TriggerResulsHandler;
     Exception excp(errors::LogicError);
     if ( hlTriggerResults_Label == "" || hlTriggerResults_Label == "NULL" ) {
         excp << "TriggerResults ==> Empty";
         excp.raise();
     }
     try {
         iEvent.getByToken(hlTriggerResults_, TriggerResulsHandler);
         if (TriggerResulsHandler.isValid())   trigRes=true;
     }  catch (...) { }
     if ( !trigRes ) {    excp << "TriggerResults ==> not readable";            excp.raise(); }
     const TriggerResults & triggerResults = *(TriggerResulsHandler.product());
 
     //get simVertex
     float simPV=0;
 
     Handle<std::vector<SimVertex> > simVertexCollection;
     
     try {
         iEvent.getByToken(simVertexCollection_, simVertexCollection);
         const SimVertex simPVh = *(simVertexCollection->begin());
         simPV=simPVh.position().z();
     }
     catch (...) {}
     
     //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" )
             {
                 try {
                     iEvent.getByToken(VertexCollection_.at(coll), VertexHandler);
                     if (VertexHandler.isValid())   VertexOK=true;                       
                 }  catch (...) { }          
             }
             
             //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
             if (VertexOK) 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 thread_unsafe::DQMEDAnalyzer.

Definition at line 122 of file HLTVertexPerformanceAnalyzer.cc.

References DQMStore::IBooker::book1D(), coll, H1_, hltPathNames_, DQMStore::IBooker::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, VertexCollection_, and VertexCollection_Label.

 {
     //book the DQM plots
     using namespace std;
     std::string dqmFolder;
     for (unsigned int ind=0; ind<hltPathNames_.size();ind++) {
         dqmFolder = Form("HLT/BTag/Vertex/%s",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
	)

virtual

Reimplemented from thread_unsafe::DQMEDAnalyzer.

Definition at line 25 of file HLTVertexPerformanceAnalyzer.cc.

References _isfoundHLTs, newFWLiteAna::found, hltConfigProvider_, hltPathIndexs_, hltPathNames_, hlTriggerResults_, HLTConfigProvider::init(), label, edm::EDConsumerBase::Labels::process, triggerConfChanged_, and HLTConfigProvider::triggerNames().

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