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 &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)

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

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Data Format.

Definition at line 21 of file RPCRecHitProbability.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 14 of file RPCRecHitProbability.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), globalFolder_, muEtaCut_, muonFolder_, muonLabel_, muPtCut_, RootFileName, saveRootFile, scalersRawToDigiLabel_, AlCaHLTBitMon_QueryRunRegistry::string, and subsystemFolder_.

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

RPCRecHitProbability::~RPCRecHitProbability ( )

override

Definition at line 29 of file RPCRecHitProbability.cc.

29 {}

Member Function Documentation

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

overrideprotected

Definition at line 103 of file RPCRecHitProbability.cc.

References counter, dcs_, PVValHelper::eta, reco::LeafCandidate::eta(), MonitorElement::Fill(), HcalObjRepresent::Fill(), mps_fire::i, reco::Muon::isGlobalMuon(), edm::HandleBase::isValid(), makeDcsInfo(), muEtaCut_, muonLabel_, muPtCut_, NumberOfMuonEta_, NumberOfMuonPhi_B_, NumberOfMuonPhi_EM_, NumberOfMuonPhi_EP_, NumberOfMuonPt_B_, NumberOfMuonPt_EM_, NumberOfMuonPt_EP_, reco::Muon::outerTrack(), phi, reco::LeafCandidate::phi(), EnergyCorrector::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().

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

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 31 of file RPCRecHitProbability.cc.

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

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

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

private

Definition at line 191 of file RPCRecHitProbability.cc.

References dcs_, edm::Event::getByToken(), edm::HandleBase::isValid(), mathSSE::return(), DcsStatus::RPC, and scalersRawToDigiLabel_.

Referenced by analyze().

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