/data/doxygen/doxygen-1.7.3/gen/CMSSW_4_2_8/src/DQM/RPCMonitorClient/src/RPCMon_SS_Dbx_Global.cc

Go to the documentation of this file.

#include <sstream>
#include <iomanip>
#include <DQM/RPCMonitorClient/interface/RPCMon_SS_Dbx_Global.h>
#include "DQM/RPCMonitorDigi/interface/utils.h"
// Framework
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
//DQM Services
#include "DQMServices/Core/interface/MonitorElement.h"
//DataFormats
#include <DataFormats/MuonDetId/interface/RPCDetId.h>
#include "DataFormats/RPCDigi/interface/RPCDigiCollection.h"
//Geometry
#include "Geometry/RPCGeometry/interface/RPCGeomServ.h"
#include "Geometry/RPCGeometry/interface/RPCGeometry.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"


RPCMon_SS_Dbx_Global::RPCMon_SS_Dbx_Global(const edm::ParameterSet& iConfig ){


  edm::LogVerbatim ("rpcmonitorerror") << "[RPCMon_SS_Dbx_Global]: Constructor";


  globalFolder_ = iConfig.getUntrackedParameter<std::string>("GlobalHistogramsFolder","RPC/RecHits/SummaryHistograms");

  saveRootFile_ =  iConfig.getUntrackedParameter<bool>("SaveRootFile", false);
  rootFileName_ = iConfig.getUntrackedParameter<std::string>("RootFileName","out.root");
  digiLabel_ = iConfig.getUntrackedParameter<std::string>("DigiLabel","muonRPCDigis");
  numberOfRings_ = iConfig.getUntrackedParameter<int>("NumberOfEndcapRings", 2);

}

RPCMon_SS_Dbx_Global::~RPCMon_SS_Dbx_Global(){
  edm::LogVerbatim ("rpcmonitorerror") << "[RPCMon_SS_Dbx_Global]: Destructor ";
  dbe_=0;
}

void RPCMon_SS_Dbx_Global::beginJob(){
 edm::LogVerbatim ("rpcmonitorerror") << "[RPCMon_SS_Dbx_Global]: Begin job ";
 dbe_ = edm::Service<DQMStore>().operator->();
}

void RPCMon_SS_Dbx_Global::endJob(){}

void RPCMon_SS_Dbx_Global::beginRun(const edm::Run& r, const edm::EventSetup& c){
  edm::LogVerbatim ("rpcmonitorerror") << "[RPCMon_SS_Dbx_Global]: Begin run";

  dbe_->setCurrentFolder(globalFolder_);

  MonitorElement* me;
  me =0;
  me = dbe_->get(globalFolder_ + "/AfterPulseBxDiff");
  if ( 0!=me) 
    dbe_->removeElement(me->getName());
  me = dbe_->book1D("AfterPulseBxDiff","After Pulse Bx Difference",13,-6.5,6.5);
}

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

 edm::Handle<RPCDigiCollection> rpcDigis;
  iEvent.getByType(rpcDigis);

 edm::ESHandle<RPCGeometry> pDD;
 iSetup.get<MuonGeometryRecord>().get( pDD );
 rpcdqm::utils rpcUtils;
 
 RPCDigiCollection::DigiRangeIterator detUnitIt;
 
 dbe_->setCurrentFolder(globalFolder_);

 // Loop over DetUnit's
 for (detUnitIt=rpcDigis->begin();detUnitIt!=rpcDigis->end();++detUnitIt){      
   const RPCDetId& id = (*detUnitIt).first;
   const RPCRoll* roll = dynamic_cast<const RPCRoll* >( pDD->roll(id));
   const RPCDigiCollection::Range& range = (*detUnitIt).second;
   //     std::cout <<" detector "<< id.region()<< std::endl;
   std::ostringstream tag;
   std::ostringstream name;     
   MonitorElement* me;

   //get roll name
   RPCGeomServ RPCname(id);
   std::string nameRoll = RPCname.name();
   //get roll number      
   RPCGeomServ RPCnumber(id);
   int nr = RPCnumber.chambernr();
   //get segment (valid only if endcap)
   RPCGeomServ RPCsegment(id);
   int seg = RPCsegment.segment();
   
   // Loop over the digis of this DetUnit 
   for (RPCDigiCollection::const_iterator digiIt = range.first;digiIt!=range.second;++digiIt) {
     if (digiIt->strip() < 1 || digiIt->strip() > roll->nstrips() )
       edm::LogVerbatim ("rpcmonitorerror") <<" XXXX  Problem with detector "<< id << ", strip # " << digiIt->strip();
     else {  //Loop on digi2
       for (RPCDigiCollection::const_iterator digiIt2 = digiIt; digiIt2!=range.second; ++digiIt2) {
         int dstrip = digiIt->strip() - digiIt2->strip();
         int dbx = digiIt->bx() - digiIt2->bx();
         if ( dstrip == 0 && abs(dbx) != 0 ) {
           //      std::cout <<" detector 3333 "<< id.region()<<" " << id.ring()<<std::endl;
           me = dbe_->get(globalFolder_ + "/AfterPulseBxDiff");
           if(!me) continue;
           me->Fill(dbx);

           if (id.region()!=0){//Endcap
             name.str("");
             tag.str("");
             tag << "AfterPulse_Ring_vs_Segment_Disk_" << id.region()*id.station();
             name << "Endcap, After Pulse, Diff bx, Disk # " << id.region()*id.station();
             me = dbe_->get(globalFolder_ + "/"+tag.str());
             if (!me){
               me = dbe_->book2D (tag.str(),name.str(),36,0.5,36.5, 3*numberOfRings_, 0.0, 3*numberOfRings_);
               
               rpcUtils.labelXAxisSegment(me);
               rpcUtils.labelYAxisRing(me, numberOfRings_);
               
             }
             int yBin;
             (numberOfRings_ == 3 ? yBin= id.ring()*3-id.roll()+1 : yBin= (id.ring()-1)*3-id.roll()+1);      
             me->Fill(seg, yBin); 
           } else { // Barrel ( region == 0 ).
             
             name.str("");
             tag.str("");
             tag << "AfterPulse_Roll_vs_Sector_Wheel_" << std::setw(2) << std::setfill('+') << id.ring();
             name << "Barrel, After Pulse, Diff. bx, Wheel # " << std::setw(2) << std::setfill('+') << id.ring();
             me = dbe_->get(globalFolder_ + "/"+tag.str());
             if (!me){
               me = dbe_->book2D (tag.str(),name.str(),12,0.5,12.5,21,0.5,21.5);
               rpcUtils.labelXAxisSector( me);
               rpcUtils.labelYAxisRoll(me, 0, id.ring());
             }
             me->Fill(id.sector(),nr);
           }//Barrel
         }
       }
     }
   }
 }
}