#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.
	~RPCRecHitProbability ()
Private Member Functions
void	makeDcsInfo (const edm::Event &)
Private Attributes
int	counter
DQMStore *	dbe
	DQM store.
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_

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.

{}

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_, eta(), reco::LeafCandidate::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(), reco::LeafCandidate::phi(), 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.

{}

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.

{}

void RPCRecHitProbability::makeDcsInfo ( const edm::Event & e ) [private]

Definition at line 225 of file RPCRecHitProbability.cc.

References dcs_, edm::Event::getByLabel(), edm::HandleBase::isValid(), 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 ;
}