Data Format. More...

#include <RPCRecHitProbability.h>

Inheritance diagram for RPCRecHitProbability:

Public Member Functions
	RPCRecHitProbability (const edm::ParameterSet &)

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

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

	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

Protected Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

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

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
void	makeDcsInfo (const edm::Event &)

Private Attributes
int	counter

bool	dcs_

std::string	globalFolder_

float	muEtaCut_

std::string	muonFolder_

edm::EDGetTokenT< reco::CandidateView >	muonLabel_

float	muPtCut_

MonitorElement *	NumberOfMuonEta_

MonitorElement *	NumberOfMuonPhi_B_

MonitorElement *	NumberOfMuonPhi_EM_

MonitorElement *	NumberOfMuonPhi_EP_

MonitorElement *	NumberOfMuonPt_B_

MonitorElement *	NumberOfMuonPt_EM_

MonitorElement *	NumberOfMuonPt_EP_

MonitorElement *	recHitEta_ [6]

MonitorElement *	recHitPhi_B_ [6]

MonitorElement *	recHitPhi_EM_ [6]

MonitorElement *	recHitPhi_EP_ [6]

MonitorElement *	recHitPt_B_ [6]

MonitorElement *	recHitPt_EM_ [6]

MonitorElement *	recHitPt_EP_ [6]

std::string	RootFileName

MonitorElement *	RPCRecHitMuonEta_

bool	saveRootFile

edm::EDGetTokenT< DcsStatusCollection >	scalersRawToDigiLabel_

std::string	subsystemFolder_

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 Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Data Format.

Definition at line 19 of file RPCRecHitProbability.h.

Constructor & Destructor Documentation

◆ RPCRecHitProbability()

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

explicit

Definition at line 14 of file RPCRecHitProbability.cc.

                                                                       : counter(0) {
   saveRootFile = pset.getUntrackedParameter<bool>("SaveRootFile", false);
   RootFileName = pset.getUntrackedParameter<std::string>("RootFileName", "RPCRecHitProbabilityDQM.root");
  
   muonLabel_ = consumes<reco::CandidateView>(pset.getParameter<edm::InputTag>("MuonLabel"));
   muPtCut_ = pset.getUntrackedParameter<double>("MuonPtCut", 3.0);
   muEtaCut_ = pset.getUntrackedParameter<double>("MuonEtaCut", 1.9);
   scalersRawToDigiLabel_ = consumes<DcsStatusCollection>(pset.getParameter<edm::InputTag>("ScalersRawToDigiLabel"));
  
   subsystemFolder_ = pset.getUntrackedParameter<std::string>("RPCFolder", "RPC");
   globalFolder_ = pset.getUntrackedParameter<std::string>("GlobalFolder", "SummaryHistograms");
   muonFolder_ = pset.getUntrackedParameter<std::string>("MuonFolder", "Muon");
 }

References globalFolder_, muEtaCut_, muonFolder_, muonLabel_, muPtCut_, muonDTDigis_cfi::pset, RootFileName, saveRootFile, scalersRawToDigiLabel_, AlCaHLTBitMon_QueryRunRegistry::string, and subsystemFolder_.

◆ ~RPCRecHitProbability()

RPCRecHitProbability::~RPCRecHitProbability ( )

override

Definition at line 28 of file RPCRecHitProbability.cc.

28 {}

Member Function Documentation

◆ analyze()

void RPCRecHitProbability::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	setup
	)

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 104 of file RPCRecHitProbability.cc.

                                                                                     {
   dcs_ = true;
   //Check HV status
   this->makeDcsInfo(event);
   if (!dcs_) {
     edm::LogWarning("rpcrechitprobability") << "[RPCRecHitProbability]: DCS bit OFF";
     return;  //if RPC not ON there's no need to continue
   }
  
   counter++;
   edm::LogInfo("rpcrechitprobability") << "[RPCRecHitProbability]: Beginning analyzing event " << counter;
  
   //Muons
   edm::Handle<reco::CandidateView> muonCands;
   event.getByToken(muonLabel_, muonCands);
   std::map<RPCDetId, std::vector<RPCRecHit> > rechitMuon;
  
   if (muonCands.isValid()) {
     int nStaMuons = muonCands->size();
  
     for (int i = 0; i < nStaMuons; i++) {
       const reco::Candidate& goodMuon = (*muonCands)[i];
       const reco::Muon* muCand = dynamic_cast<const reco::Muon*>(&goodMuon);
  
       if (!muCand->isGlobalMuon())
         continue;
       float eta = muCand->eta();
       float pt = muCand->pt();
       if (pt < muPtCut_ || fabs(eta) > muEtaCut_)
         continue;
  
       float phi = muCand->phi();
  
       NumberOfMuonEta_->Fill(eta);
  
       if (eta > 0.8) {
         NumberOfMuonPt_EP_->Fill(pt);
         NumberOfMuonPhi_EP_->Fill(phi);
       } else if (eta < -0.8) {
         NumberOfMuonPt_EM_->Fill(pt);
         NumberOfMuonPhi_EM_->Fill(phi);
       } else {
         NumberOfMuonPt_B_->Fill(pt);
         NumberOfMuonPhi_B_->Fill(phi);
       }
  
       reco::Track muTrack = (*(muCand->outerTrack()));
       std::vector<TrackingRecHitRef> rpcTrackRecHits;
  
       //loop on mu rechits
  
       int recHitCounter = 0;
       for (trackingRecHit_iterator it = muTrack.recHitsBegin(); it != muTrack.recHitsEnd(); it++) {
         if (!(*it)->isValid())
           continue;
         int muSubDetId = (*it)->geographicalId().subdetId();
         if (muSubDetId == MuonSubdetId::RPC) {
           recHitCounter++;
         }
       }  // end loop on mu rechits
  
       RPCRecHitMuonEta_->Fill(eta, recHitCounter);
  
       int j = 0;
       while (recHitCounter >= j + 1 && j < 6) {
         if (recHitEta_[j])
           recHitEta_[j]->Fill(eta);
         if (eta > 0.8) {
           recHitPt_EP_[j]->Fill(pt);
           recHitPhi_EP_[j]->Fill(phi);
         } else if (eta < -0.8) {
           recHitPt_EM_[j]->Fill(pt);
           recHitPhi_EM_[j]->Fill(phi);
         } else {
           recHitPt_B_[j]->Fill(pt);
           recHitPhi_B_[j]->Fill(phi);
         }
  
         j++;
       }
     }
   } else {
     edm::LogError("rpcrechitprobability") << "[RPCRecHitProbability]: Muons - Product not valid for event" << counter;
   }
 }

References counter, dcs_, PVValHelper::eta, reco::LeafCandidate::eta(), dqm::impl::MonitorElement::Fill(), mps_fire::i, reco::Muon::isGlobalMuon(), edm::HandleBase::isValid(), dqmiolumiharvest::j, makeDcsInfo(), muEtaCut_, muonLabel_, muPtCut_, NumberOfMuonEta_, NumberOfMuonPhi_B_, NumberOfMuonPhi_EM_, NumberOfMuonPhi_EP_, NumberOfMuonPt_B_, NumberOfMuonPt_EM_, NumberOfMuonPt_EP_, reco::Muon::outerTrack(), phi, reco::LeafCandidate::phi(), DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), recHitEta_, recHitPhi_B_, recHitPhi_EM_, recHitPhi_EP_, recHitPt_B_, recHitPt_EM_, recHitPt_EP_, reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), MuonSubdetId::RPC, RPCRecHitMuonEta_, and edm::View< T >::size().

◆ bookHistograms()

void RPCRecHitProbability::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	r,
		edm::EventSetup const &	iSetup
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 30 of file RPCRecHitProbability.cc.

                                                                        {
   edm::LogInfo("rpcrechitprobability") << "[RPCRecHitProbability]: Begin Run ";
  
   std::string currentFolder = subsystemFolder_ + "/" + muonFolder_ + "/" + globalFolder_;
   ibooker.setCurrentFolder(currentFolder);
  
   int ptBins = 100 - (int)muPtCut_;
  
   //General part
   NumberOfMuonEta_ = ibooker.book1D("NumberOfMuonEta", "Muons vs Eta", 20 * muEtaCut_, -muEtaCut_, muEtaCut_);
   NumberOfMuonPt_B_ = ibooker.book1D("NumberOfMuonPt_Barrel", "Muons vs Pt - Barrel", ptBins, muPtCut_, 100);
   NumberOfMuonPt_EP_ = ibooker.book1D("NumberOfMuonPt_EndcapP", "Muons vs Pt - Endcap+", ptBins, muPtCut_, 100);
   NumberOfMuonPt_EM_ = ibooker.book1D("NumberOfMuonPt_EndcapM", "Muons vs Pt - Endcap-", ptBins, muPtCut_, 100);
   NumberOfMuonPhi_B_ =
       ibooker.book1D("NumberOfMuonPhi_Barrel", "Muons vs Phi - Barrel", 144, -TMath::Pi(), TMath::Pi());
   NumberOfMuonPhi_EP_ =
       ibooker.book1D("NumberOfMuonPhi_EndcapP", "Muons vs Phi - Endcap+", 144, -TMath::Pi(), TMath::Pi());
   NumberOfMuonPhi_EM_ =
       ibooker.book1D("NumberOfMuonPhi_EndcapM", "Muons vs Phi - Endcap-", 144, -TMath::Pi(), TMath::Pi());
  
   //RPC part
   RPCRecHitMuonEta_ = ibooker.book2D(
       "RPCRecHitMuonEta", "Number Of RecHits per Muons vs Eta", 20 * muEtaCut_, -muEtaCut_, muEtaCut_, 7, 0.5, 7.5);
  
   std::stringstream name, title;
   for (int i = 0; i < 6; i++) {
     name.str("");
     title.str("");
     name << (i + 1) << "RecHitMuonEta";
     title << "At least " << (i + 1) << " Cluster vs Eta";
     recHitEta_[i] = ibooker.book1D(name.str(), title.str(), 20 * muEtaCut_, -muEtaCut_, muEtaCut_);
  
     name.str("");
     title.str("");
     name << (i + 1) << "RecHitMuonPhiB";
     title << "At least " << (i + 1) << " Cluster vs Phi-Barrel";
     recHitPhi_B_[i] = ibooker.book1D(name.str(), title.str(), 144, -TMath::Pi(), TMath::Pi());
  
     name.str("");
     title.str("");
     name << (i + 1) << "RecHitMuonPtB";
     title << "At least " << (i + 1) << " Cluster vs Pt-Barrel";
     recHitPt_B_[i] = ibooker.book1D(name.str(), title.str(), ptBins, muPtCut_, 100);
  
     name.str("");
     title.str("");
     name << (i + 1) << "RecHitMuonPhiEP";
     title << "At least " << (i + 1) << " Cluster vs Phi-Endcap+";
     recHitPhi_EP_[i] = ibooker.book1D(name.str(), title.str(), 144, -TMath::Pi(), TMath::Pi());
  
     name.str("");
     title.str("");
     name << (i + 1) << "RecHitMuonPtEP";
     title << "At least " << (i + 1) << " Cluster vs Pt-Endcap+";
     recHitPt_EP_[i] = ibooker.book1D(name.str(), title.str(), ptBins, muPtCut_, 100);
  
     name.str("");
     title.str("");
     name << (i + 1) << "RecHitMuonPhiEM";
     title << "At least " << (i + 1) << " Cluster vs Phi-Endcap-";
     recHitPhi_EM_[i] = ibooker.book1D(name.str(), title.str(), 144, -TMath::Pi(), TMath::Pi());
  
     name.str("");
     title.str("");
     name << (i + 1) << "RecHitMuonPtEM";
     title << "At least " << (i + 1) << " Cluster vs Pt-Endcap-";
     recHitPt_EM_[i] = ibooker.book1D(name.str(), title.str(), ptBins, muPtCut_, 100);
   }
  
   dcs_ = true;
 }

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dcs_, globalFolder_, mps_fire::i, createfilelist::int, muEtaCut_, muonFolder_, muPtCut_, Skims_PA_cff::name, NumberOfMuonEta_, NumberOfMuonPhi_B_, NumberOfMuonPhi_EM_, NumberOfMuonPhi_EP_, NumberOfMuonPt_B_, NumberOfMuonPt_EM_, NumberOfMuonPt_EP_, Pi, L1EGammaCrystalsEmulatorProducer_cfi::ptBins, recHitEta_, recHitPhi_B_, recHitPhi_EM_, recHitPhi_EP_, recHitPt_B_, recHitPt_EM_, recHitPt_EP_, RPCRecHitMuonEta_, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, subsystemFolder_, and overlapproblemtsosanalyzer_cfi::title.

◆ makeDcsInfo()

void RPCRecHitProbability::makeDcsInfo ( const edm::Event & e )

private

Definition at line 190 of file RPCRecHitProbability.cc.

                                                         {
   edm::Handle<DcsStatusCollection> dcsStatus;
  
   if (!e.getByToken(scalersRawToDigiLabel_, dcsStatus)) {
     dcs_ = true;
     return;
   }
  
   if (!dcsStatus.isValid()) {
     edm::LogWarning("RPCDcsInfo") << "scalersRawToDigi not found";
     dcs_ = true;  // info not available: set to true
     return;
   }
  
   for (DcsStatusCollection::const_iterator dcsStatusItr = dcsStatus->begin(); dcsStatusItr != dcsStatus->end();
        ++dcsStatusItr) {
     if (!dcsStatusItr->ready(DcsStatus::RPC))
       dcs_ = false;
   }
  
   return;
 }