#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 &event, edm::EventSetup const &setup) final

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

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

void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer ()

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

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

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

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 DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

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

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

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::HLTVertexPerformanceAnalyzer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 6 of file HLTVertexPerformanceAnalyzer.cc.

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

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

◆ ~HLTVertexPerformanceAnalyzer()

HLTVertexPerformanceAnalyzer::~HLTVertexPerformanceAnalyzer ( )

override

Definition at line 24 of file HLTVertexPerformanceAnalyzer.cc.

24 {}

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 56 of file HLTVertexPerformanceAnalyzer.cc.

                                                                                               {
   bool trigRes = false;
   using namespace edm;
  
   // get triggerResults
   Handle<TriggerResults> TriggerResulsHandler;
   if (hlTriggerResults_Label.empty() || 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).empty() && 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_
 }

References iEvent, edm::HandleBase::isValid(), seedmultiplicitymonitor_newtracking_cfi::maxValue, CoreSimVertex::position(), edm::Handle< T >::product(), caloTruthProducer_cfi::simVertexCollection, and edm::triggerResults().

◆ bookHistograms()

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

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 125 of file HLTVertexPerformanceAnalyzer.cc.

                                                                                {
   // 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).empty() && 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);
       }
     }
   }
 }

References dqm::implementation::IBooker::book1D(), L1TMonitor_cff::dqmFolder, seedmultiplicitymonitor_newtracking_cfi::maxValue, dqm::impl::MonitorElement::setAxisTitle(), dqm::implementation::NavigatorBase::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

◆ dqmBeginRun()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 26 of file HLTVertexPerformanceAnalyzer.cc.

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