#include <RPCRecHitProbability.h>

Inheritance diagram for RPCRecHitProbability:

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

virtual void	beginJob ()

void	beginRun (const edm::Run &r, const edm::EventSetup &c)

virtual void	endJob (void)

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

	RPCRecHitProbability (const edm::ParameterSet &)
	Data Format. More...

	~RPCRecHitProbability ()

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

std::string	workerType () const

virtual	~EDAnalyzer ()

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

Private Attributes
int	counter

DQMStore *	dbe
	DQM store. More...

bool	dcs_

std::string	globalFolder_

float	muEtaCut_

std::string	muonFolder_

edm::InputTag	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

std::string	subsystemFolder_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

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::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 19 of file RPCRecHitProbability.h.

Constructor & Destructor Documentation

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

explicit

Data Format.

Definition at line 17 of file RPCRecHitProbability.cc.

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

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

RPCRecHitProbability::~RPCRecHitProbability ( )

Definition at line 31 of file RPCRecHitProbability.cc.

31 {}

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 137 of file RPCRecHitProbability.cc.

References counter, dcs_, reco::LeafCandidate::eta(), eta(), MonitorElement::Fill(), HcalObjRepresent::Fill(), i, reco::Muon::isGlobalMuon(), edm::HandleBase::isValid(), 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(), reco::LeafCandidate::pt(), recHitEta_, recHitPhi_B_, recHitPhi_EM_, recHitPhi_EP_, recHitPt_B_, recHitPt_EM_, recHitPt_EP_, reco::Track::recHitsBegin(), reco::Track::recHitsEnd(), MuonSubdetId::RPC, and RPCRecHitMuonEta_.

                                                                                    {
   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.getByLabel(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::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 32 of file RPCRecHitProbability.cc.

32 {}

void RPCRecHitProbability::beginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

virtual

get hold of back-end interface

Reimplemented from edm::EDAnalyzer.

Definition at line 34 of file RPCRecHitProbability.cc.

References DQMStore::book1D(), DQMStore::book2D(), dbe, dcs_, DQMStore::get(), MonitorElement::getName(), globalFolder_, i, muEtaCut_, muonFolder_, muPtCut_, mergeVDriftHistosByStation::name, NumberOfMuonEta_, NumberOfMuonPhi_B_, NumberOfMuonPhi_EM_, NumberOfMuonPhi_EP_, NumberOfMuonPt_B_, NumberOfMuonPt_EM_, NumberOfMuonPt_EP_, cppFunctionSkipper::operator, Pi, jptDQMConfig_cff::ptBins, recHitEta_, recHitPhi_B_, recHitPhi_EM_, recHitPhi_EP_, recHitPt_B_, recHitPt_EM_, recHitPt_EP_, DQMStore::removeElement(), RPCRecHitMuonEta_, DQMStore::setCurrentFolder(), subsystemFolder_, and indexGen::title.

                                                                                {
 
   edm::LogInfo ("rpcrechitprobability") <<"[RPCRecHitProbability]: Begin Run " ;
   
   dbe = edm::Service<DQMStore>().operator->();
 
   std::string  currentFolder = subsystemFolder_ +"/"+muonFolder_+"/"+ globalFolder_;
   dbe->setCurrentFolder(currentFolder); 
 
   int ptBins = 100 - (int)muPtCut_;
 
   //General part
   NumberOfMuonEta_ = dbe->get(currentFolder+"/NumberOfMuonEta");
   if(NumberOfMuonEta_) dbe->removeElement(NumberOfMuonEta_->getName());
   NumberOfMuonEta_ = dbe->book1D("NumberOfMuonEta", "Muons vs Eta", 20*muEtaCut_,  -muEtaCut_,  muEtaCut_ );
 
   NumberOfMuonPt_B_ = dbe->get(currentFolder+"/NumberOfMuonPt_Barrel");
   if(NumberOfMuonPt_B_) dbe->removeElement(NumberOfMuonPt_B_->getName());
   NumberOfMuonPt_B_ = dbe->book1D("NumberOfMuonPt_Barrel", "Muons vs Pt - Barrel", ptBins, muPtCut_, 100);
 
   NumberOfMuonPt_EP_ = dbe->get(currentFolder+"/NumberOfMuonPt_EndcapP");
   if(NumberOfMuonPt_EP_) dbe->removeElement(NumberOfMuonPt_EP_->getName());
   NumberOfMuonPt_EP_ = dbe->book1D("NumberOfMuonPt_EndcapP", "Muons vs Pt - Endcap+", ptBins , muPtCut_ , 100);
 
   NumberOfMuonPt_EM_ = dbe->get(currentFolder+"/NumberOfMuonPt_EndcapM");
   if(NumberOfMuonPt_EM_) dbe->removeElement(NumberOfMuonPt_EM_->getName());
   NumberOfMuonPt_EM_ = dbe->book1D("NumberOfMuonPt_EndcapM", "Muons vs Pt - Endcap-", ptBins , muPtCut_ , 100);
 
   NumberOfMuonPhi_B_ = dbe->get(currentFolder+"/NumberOfMuonPhi_Barrel");
   if(NumberOfMuonPhi_B_) dbe->removeElement(NumberOfMuonPhi_B_->getName());
   NumberOfMuonPhi_B_ = dbe->book1D("NumberOfMuonPhi_Barrel", "Muons vs Phi - Barrel", 144, -TMath::Pi(), TMath::Pi());
 
   NumberOfMuonPhi_EP_ = dbe->get(currentFolder+"/NumberOfMuonPhi_EndcapP");
   if(NumberOfMuonPhi_EP_) dbe->removeElement(NumberOfMuonPhi_EP_->getName());
   NumberOfMuonPhi_EP_ = dbe->book1D("NumberOfMuonPhi_EndcapP", "Muons vs Phi - Endcap+", 144,  -TMath::Pi(), TMath::Pi() );
 
   NumberOfMuonPhi_EM_ = dbe->get(currentFolder+"/NumberOfMuonPhi_EndcapM");
   if(NumberOfMuonPhi_EM_) dbe->removeElement(NumberOfMuonPhi_EM_->getName());
   NumberOfMuonPhi_EM_ = dbe->book1D("NumberOfMuonPhi_EndcapM", "Muons vs Phi - Endcap-", 144, -TMath::Pi(), TMath::Pi());
 
   //RPC part
   RPCRecHitMuonEta_ = dbe->get(currentFolder+"/RPCRecHitMuonEta");
   if(RPCRecHitMuonEta_) dbe->removeElement(RPCRecHitMuonEta_->getName());
   RPCRecHitMuonEta_ = dbe->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] = dbe->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] = dbe->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] = dbe->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] = dbe->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] = dbe->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] = dbe->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] = dbe->book1D(name.str(), title.str(), ptBins , muPtCut_ , 100);
 
   }
 
   dcs_ = true;
 }

void RPCRecHitProbability::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 129 of file RPCRecHitProbability.cc.

References dbe, RootFileName, DQMStore::save(), and saveRootFile.

                                      {
   if(saveRootFile) dbe->save(RootFileName); 
   dbe = 0;
 }

void RPCRecHitProbability::endLuminosityBlock	(	edm::LuminosityBlock const &	L,
		edm::EventSetup const &	E
	)

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 134 of file RPCRecHitProbability.cc.

134 {}

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

private

Definition at line 225 of file RPCRecHitProbability.cc.

References dcs_, edm::Event::getByLabel(), edm::HandleBase::isValid(), hitfit::return, and DcsStatus::RPC.

Referenced by analyze().

                                                          {
 
   edm::Handle<DcsStatusCollection> dcsStatus;
 
   if ( ! e.getByLabel("scalersRawToDigi", 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 ;
 }

Member Data Documentation

int RPCRecHitProbability::counter

private

Definition at line 40 of file RPCRecHitProbability.h.

Referenced by analyze().

DQMStore* RPCRecHitProbability::dbe

private

DQM store.

Definition at line 43 of file RPCRecHitProbability.h.

Referenced by beginRun(), and endJob().

bool RPCRecHitProbability::dcs_

private

Definition at line 44 of file RPCRecHitProbability.h.

Referenced by analyze(), beginRun(), and makeDcsInfo().

std::string RPCRecHitProbability::globalFolder_

private

Definition at line 47 of file RPCRecHitProbability.h.

Referenced by beginRun(), and RPCRecHitProbability().

float RPCRecHitProbability::muEtaCut_

private

Definition at line 45 of file RPCRecHitProbability.h.

Referenced by analyze(), beginRun(), and RPCRecHitProbability().

std::string RPCRecHitProbability::muonFolder_

private

Definition at line 39 of file RPCRecHitProbability.h.

Referenced by beginRun(), and RPCRecHitProbability().

edm::InputTag RPCRecHitProbability::muonLabel_

private

Definition at line 77 of file RPCRecHitProbability.h.

Referenced by analyze(), and RPCRecHitProbability().

float RPCRecHitProbability::muPtCut_

private

Definition at line 45 of file RPCRecHitProbability.h.

Referenced by analyze(), beginRun(), and RPCRecHitProbability().

MonitorElement* RPCRecHitProbability::NumberOfMuonEta_

private

Definition at line 62 of file RPCRecHitProbability.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCRecHitProbability::NumberOfMuonPhi_B_

private

Definition at line 54 of file RPCRecHitProbability.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCRecHitProbability::NumberOfMuonPhi_EM_

private

Definition at line 60 of file RPCRecHitProbability.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCRecHitProbability::NumberOfMuonPhi_EP_

private

Definition at line 57 of file RPCRecHitProbability.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCRecHitProbability::NumberOfMuonPt_B_

private

Definition at line 53 of file RPCRecHitProbability.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCRecHitProbability::NumberOfMuonPt_EM_

private

Definition at line 59 of file RPCRecHitProbability.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCRecHitProbability::NumberOfMuonPt_EP_

private

Definition at line 56 of file RPCRecHitProbability.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCRecHitProbability::recHitEta_[6]

private

Definition at line 66 of file RPCRecHitProbability.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCRecHitProbability::recHitPhi_B_[6]

private

Definition at line 69 of file RPCRecHitProbability.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCRecHitProbability::recHitPhi_EM_[6]

private

Definition at line 75 of file RPCRecHitProbability.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCRecHitProbability::recHitPhi_EP_[6]

private

Definition at line 72 of file RPCRecHitProbability.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCRecHitProbability::recHitPt_B_[6]

private

Definition at line 68 of file RPCRecHitProbability.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCRecHitProbability::recHitPt_EM_[6]

private

Definition at line 74 of file RPCRecHitProbability.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCRecHitProbability::recHitPt_EP_[6]

private

Definition at line 71 of file RPCRecHitProbability.h.

Referenced by analyze(), and beginRun().

std::string RPCRecHitProbability::RootFileName

private

Definition at line 51 of file RPCRecHitProbability.h.

Referenced by endJob(), and RPCRecHitProbability().

MonitorElement* RPCRecHitProbability::RPCRecHitMuonEta_

private

Definition at line 64 of file RPCRecHitProbability.h.

Referenced by analyze(), and beginRun().

bool RPCRecHitProbability::saveRootFile

private

Definition at line 50 of file RPCRecHitProbability.h.

Referenced by endJob(), and RPCRecHitProbability().

std::string RPCRecHitProbability::subsystemFolder_

private

Definition at line 48 of file RPCRecHitProbability.h.

Referenced by beginRun(), and RPCRecHitProbability().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation