CMSSW_4_4_3_patch1/src/DQM/RPCMonitorDigi/src/RPCMonitorSync.cc

Go to the documentation of this file.

#include <DQM/RPCMonitorDigi/interface/RPCMonitorSync.h>



// RPC Digi
#include "DataFormats/MuonDetId/interface/RPCDetId.h"
#include "DataFormats/RPCDigi/interface/RPCDigiCollection.h"

#include <FWCore/MessageLogger/interface/MessageLogger.h>

#include <vector>
#include <ctime>
#include <sstream>


// Detector overview historam booking routines

MonitorElement *RPCMonitorSync::barrelOffsetHist(std::string name, std::string title ) {
  std::stringstream htit;
  htit<<"Synchronization offset - "<<title;
  dbe->setCurrentFolder("Global");
  MonitorElement *hist = dbe->book2D(name,htit.str(),18,1,7,24,1,13);
  return hist;
} 

MonitorElement *RPCMonitorSync::endcapOffsetHist( std::string name, std::string title ) {
  std::stringstream htit;
  htit<<"Synchronization offset - "<<title;
  dbe->setCurrentFolder("Global");
  MonitorElement *hist = dbe->book2D(name,htit.str(),12,1,4,36,1,7);
  return hist;
} 

MonitorElement *RPCMonitorSync::barrelWidthHist( std::string name, std::string title ) {
  std::stringstream htit;
  htit<<"Synchronization width - "<<title;
  dbe->setCurrentFolder("Global");
  MonitorElement *hist = dbe->book2D(name,htit.str(),18,1,7,24,1,13);
  return hist;
} 

MonitorElement *RPCMonitorSync::endcapWidthHist( std::string name, std::string title ) {
  std::stringstream htit;
  htit<<"Synchronization width - "<<title;
  dbe->setCurrentFolder("Global");
  MonitorElement *hist = dbe->book2D(name,htit.str(),12,1,4,36,1,7);
  return hist;
} 


// Class constructor and destructor

RPCMonitorSync::RPCMonitorSync( const edm::ParameterSet& pset )
{
  counter=0;

  nameInLog = pset.getUntrackedParameter<std::string>("moduleLogName", "RPC_DQM");

  saveRootFile  = pset.getUntrackedParameter<bool>("SyncDQMSaveRootFile", false); 
  saveRootFileEventsInterval  = pset.getUntrackedParameter<int>("SyncEventsInterval", 10000); 
  RootFileName  = pset.getUntrackedParameter<std::string>("RootFileNameSync", "RPCMonitorSync.root"); 
  dbe = edm::Service<DQMStore>().operator->();
  dbe->showDirStructure();
}

RPCMonitorSync::~RPCMonitorSync(){

}


// Histogram for individual DetUnit - booking routine

std::map<std::string, MonitorElement*> RPCMonitorSync::bookDetUnitME(RPCDetId & detId) {

  std::map<std::string, MonitorElement*> meMap;
  std::string regionName;
  std::string ringType;

  if(detId.region() ==  0) {
    regionName="Barrel";
    ringType="Wheel";
  } else {
    ringType="Disk";
    if(detId.region() == -1) regionName="Encap-";
    if(detId.region() ==  1) regionName="Encap+";
  }
 
  char folder[120];
  sprintf(folder,"%s/%s_%d/station_%d/sector_%d",regionName.c_str(),ringType.c_str(),
                                detId.ring(),detId.station(),detId.sector());
 
  dbe->setCurrentFolder(folder);

  // Name components common to current RPDDetId  
  char detUnitLabel[128];
  char layerLabel[128];
  sprintf(detUnitLabel ,"%d",detId());
  sprintf(layerLabel ,"layer%d_subsector%d_roll%d",detId.layer(),detId.subsector(),detId.roll());

  char meId [128];
  char meTitle [128];
  
  sprintf(meId,"Sync_%s",detUnitLabel);
  sprintf(meTitle,"Sychronization_for_%s",layerLabel);
  meMap[meId] = dbe->book1D(meId, meTitle, 9, -4.5, 4.5);
  
  return meMap;
}


// Fill histograms and save data

void RPCMonitorSync::endJob(void)
{
  int station_map[8]={1,2,3,4,5,0,6,0}; // map RPC layer numbers to histogram bins
    
  MonitorElement *hOffsetMBp1 = barrelOffsetHist("hOffsetMBp1","Barrell Wheel +1");
  MonitorElement *hOffsetMBp2 = barrelOffsetHist("hOffsetMBp2","Barrell Wheel +2");
  MonitorElement *hOffsetMB0  = barrelOffsetHist("hOffsetMB0" ,"Barrell Wheel 0");
  MonitorElement *hOffsetMBm1 = barrelOffsetHist("hOffsetMBm1","Barrell Wheel -1");
  MonitorElement *hOffsetMBm2 = barrelOffsetHist("hOffsetMBm2","Barrell Wheel -2");
  MonitorElement *hOffsetMEp1 = endcapOffsetHist("hOffsetMEp1","Endcap Disk +1");
  MonitorElement *hOffsetMEp2 = endcapOffsetHist("hOffsetMEp2","Endcap Disk +2");
  MonitorElement *hOffsetMEp3 = endcapOffsetHist("hOffsetMEp3","Endcap Disk +3");
  MonitorElement *hOffsetMEm1 = endcapOffsetHist("hOffsetMEm1","Endcap Disk -1");
  MonitorElement *hOffsetMEm2 = endcapOffsetHist("hOffsetMEm2","Endcap Disk -2");
  MonitorElement *hOffsetMEm3 = endcapOffsetHist("hOffsetMEm3","Endcap Disk -3");

  MonitorElement *hWidthMBp1 = barrelWidthHist("hWidthMBp1","Barrell Wheel +1");
  MonitorElement *hWidthMBp2 = barrelWidthHist("hWidthMBp2","Barrell Wheel +2");
  MonitorElement *hWidthMB0  = barrelWidthHist("hWidthMB0" ,"Barrell Wheel 0");
  MonitorElement *hWidthMBm1 = barrelWidthHist("hWidthMBm1","Barrell Wheel -1");
  MonitorElement *hWidthMBm2 = barrelWidthHist("hWidthMBm2","Barrell Wheel -2");
  MonitorElement *hWidthMEp1 = endcapWidthHist("hWidthMEp1","Endcap Disk +1");
  MonitorElement *hWidthMEp2 = endcapWidthHist("hWidthMEp2","Endcap Disk +2");
  MonitorElement *hWidthMEp3 = endcapWidthHist("hWidthMEp3","Endcap Disk +3");
  MonitorElement *hWidthMEm1 = endcapWidthHist("hWidthMEm1","Endcap Disk -1");
  MonitorElement *hWidthMEm2 = endcapWidthHist("hWidthMEm2","Endcap Disk -2");
  MonitorElement *hWidthMEm3 = endcapWidthHist("hWidthMEm3","Endcap Disk -3");
 
  std::map <uint32_t,timing>::const_iterator ci;
  float offset,width;
  float xf=0,yf=0;

  for(ci=synchroMap.begin();ci!=synchroMap.end();ci++){
    uint32_t id = ci->first;
    RPCDetId detId(id);
    
    offset = ci->second.offset();
    width = ci->second.width();

    //    RPCDetId *tempDetId=new RPCDetId(ci->first); 
        
      
    if( detId.region()==0 ) {
      // Fill barrel histogram
      xf=station_map[detId.station()*2+detId.layer()-3]+((float)(detId.roll()-0.5)/3.);
      yf=detId.sector()+((float)(detId.subsector()-0.5)/2.);
      if ((detId.sector()==4) && (detId.station()==4)) 
        yf=detId.sector()+((float)(detId.subsector()-0.5)/4.);

      if (detId.ring()==1) hOffsetMBp1->Fill(xf,yf,offset);
      if (detId.ring()==2) hOffsetMBp2->Fill(xf,yf,offset);
      if (detId.ring()==-1) hOffsetMBm1->Fill(xf,yf,offset);
      if (detId.ring()==-2) hOffsetMBm2->Fill(xf,yf,offset);
      if (detId.ring()==0) hOffsetMB0->Fill(xf,yf,offset);
      if (detId.ring()==1) hWidthMBp1->Fill(xf,yf,width);
      if (detId.ring()==2) hWidthMBp2->Fill(xf,yf,width);
      if (detId.ring()==-1) hWidthMBm1->Fill(xf,yf,width);
      if (detId.ring()==-2) hWidthMBm2->Fill(xf,yf,width);
      if (detId.ring()==0) hWidthMB0->Fill(xf,yf,width);
    } else {
      // Fill endcap histogram
      xf=detId.ring()  +((float)(detId.roll()-0.5)/4.);
      yf=detId.sector()+((float)(detId.subsector()-0.5)/6.);

      if (detId.region()==1) {
        if (detId.station()==1) hOffsetMEp1->Fill(xf,yf,offset);
        if (detId.station()==2) hOffsetMEp2->Fill(xf,yf,offset);
        if (detId.station()==3) hOffsetMEp3->Fill(xf,yf,offset);
        if (detId.station()==1) hWidthMEp1->Fill(xf,yf,width);
        if (detId.station()==2) hWidthMEp2->Fill(xf,yf,width);
        if (detId.station()==3) hWidthMEp3->Fill(xf,yf,width);
      }
      if (detId.region()==-1) {
        if (detId.station()==1) hOffsetMEm1->Fill(xf,yf,offset);
        if (detId.station()==2) hOffsetMEm2->Fill(xf,yf,offset);
        if (detId.station()==3) hOffsetMEm3->Fill(xf,yf,offset);
        if (detId.station()==1) hWidthMEm1->Fill(xf,yf,width);
        if (detId.station()==2) hWidthMEm2->Fill(xf,yf,width);
        if (detId.station()==3) hWidthMEm3->Fill(xf,yf,width);
      }
    } // end histogram filling

  } // end loop over synchroMap

  if(saveRootFile) 
    dbe->save(RootFileName);

}


// Loop over Digis and read synchronization information

void RPCMonitorSync::readRPCDAQStrips(const edm::Event& iEvent) {

  timing aTiming;
  aTiming.inTime = 0;
  for (int i=0;i<4;i++) {
    aTiming.early_all[i] = 0;
    aTiming.late_all[i] = 0;
  }

  char detUnitLabel[128];
  char layerLabel[128];
  char meId [128];

  edm::Handle<RPCDigiCollection> rpcDigis;
//  LogInfo (nameInLog) << "   RPC Digis: " << rpcDigis->size() << endl;
  iEvent.getByType(rpcDigis);
  RPCDigiCollection::DigiRangeIterator rpcDigiCI;
  for(rpcDigiCI = rpcDigis->begin();rpcDigiCI!=rpcDigis->end();rpcDigiCI++){

    RPCDetId detId=(*rpcDigiCI).first; 
    uint32_t id=detId(); 

    sprintf(detUnitLabel ,"%d",detId());
    sprintf(layerLabel ,"layer%d_subsector%d_roll%d",detId.layer(),detId.subsector(),detId.roll());
 
    std::map<uint32_t, std::map<std::string,MonitorElement*> >::iterator meItr = meCollection.find(id);
    if (meItr == meCollection.end() || (meCollection.size()==0)) {
      meCollection[id]=bookDetUnitME(detId);
    }
    std::map<std::string, MonitorElement*> meMap=meCollection[id];

    if(synchroMap.find(id)==synchroMap.end()) { 
      synchroMap[id] = aTiming;
    }
    const RPCDigiCollection::Range& range = (*rpcDigiCI).second;
    int aBX=4;    
    for (RPCDigiCollection::const_iterator digiIt = range.first; digiIt!=range.second; ++digiIt){
      if( digiIt->bx()<aBX) aBX=digiIt->bx();
    }
    if (abs(aBX)<4) {
      if(aBX<0) synchroMap[id].early_all[-aBX]++;
      if(aBX==0) synchroMap[id].inTime++;
      if(aBX>0) synchroMap[id].late_all[aBX]++;
    }
    sprintf(meId,"Sync_%s",detUnitLabel);
    meMap[meId]->Fill(aBX);
  }
  if((!(counter%saveRootFileEventsInterval))&&(saveRootFile) ) dbe->save(RootFileName);
}


// Process event

void RPCMonitorSync::analyze( const edm::Event& iEvent, const edm::EventSetup& iSetup ) {
  
  edm::LogInfo (nameInLog) << "Beginning analyzing event " << counter++;

   //   time_t aTime = iEvent.time().value();

   
   readRPCDAQStrips(iEvent);

   return;
}