#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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

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

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
	)

overrideprivate

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