#include <RPCMonitorSync.h>

Inheritance diagram for RPCMonitorSync:

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

std::map< std::string, MonitorElement * >	bookDetUnitME (RPCDetId &detId)
	Booking of MonitoringElemnt for one RPCDetId (= roll) More...

virtual void	endJob (void)

	RPCMonitorSync (const edm::ParameterSet &)

	~RPCMonitorSync ()

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

std::string	workerType () const

virtual	~EDAnalyzer ()

Private Member Functions
MonitorElement *	barrelOffsetHist (std::string name, std::string title)
	Log messages. More...

MonitorElement *	barrelWidthHist (std::string name, std::string title)

MonitorElement *	endcapOffsetHist (std::string name, std::string title)

MonitorElement *	endcapWidthHist (std::string name, std::string title)

void	readRPCDAQStrips (const edm::Event &iEvent)

Private Attributes
int	counter

DQMStore *	dbe
	back-end interface More...

MonitorElement *	h1

std::map< uint32_t, std::map < std::string, MonitorElement * > >	meCollection

std::string	nameInLog

std::string	RootFileName

bool	saveRootFile

int	saveRootFileEventsInterval

std::map< uint32_t, timing >	synchroMap

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

RPC Synchronization Monitoring Class

Date:: 2009/08/27 09:42:37

Revision:: 1.6

Author: Piotr Traczyk (SINS)

Anna Cimmino (INFN)

Definition at line 91 of file RPCMonitorSync.h.

Constructor & Destructor Documentation

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

explicit

get hold of back-end interface

Definition at line 64 of file RPCMonitorSync.cc.

References counter, dbe, edm::ParameterSet::getUntrackedParameter(), nameInLog, cppFunctionSkipper::operator, RootFileName, saveRootFile, and saveRootFileEventsInterval.

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

Definition at line 78 of file RPCMonitorSync.cc.

                                {
 
 }

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 270 of file RPCMonitorSync.cc.

References counter, nameInLog, and readRPCDAQStrips().

                                                                                   {
   
   edm::LogInfo (nameInLog) << "Beginning analyzing event " << counter++;
 
    //   time_t aTime = iEvent.time().value();
 
    readRPCDAQStrips(iEvent);
 
    return;
 }

MonitorElement * RPCMonitorSync::barrelOffsetHist	(	std::string	name,
		std::string	title
	)

private

Log messages.

Definition at line 29 of file RPCMonitorSync.cc.

References DQMStore::book2D(), dbe, estimatePileup::hist, DQMStore::setCurrentFolder(), and indexGen::title.

Referenced by endJob().

                                                                                  {
   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::barrelWidthHist	(	std::string	name,
		std::string	title
	)

private

Definition at line 45 of file RPCMonitorSync.cc.

References DQMStore::book2D(), dbe, estimatePileup::hist, DQMStore::setCurrentFolder(), and indexGen::title.

Referenced by endJob().

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

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

Booking of MonitoringElemnt for one RPCDetId (= roll)

Definition at line 85 of file RPCMonitorSync.cc.

References DQMStore::book1D(), dbe, RPCDetId::layer(), RPCDetId::region(), RPCDetId::ring(), RPCDetId::roll(), RPCDetId::sector(), DQMStore::setCurrentFolder(), RPCDetId::station(), and RPCDetId::subsector().

Referenced by readRPCDAQStrips().

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

MonitorElement * RPCMonitorSync::endcapOffsetHist	(	std::string	name,
		std::string	title
	)

private

Definition at line 37 of file RPCMonitorSync.cc.

References DQMStore::book2D(), dbe, estimatePileup::hist, DQMStore::setCurrentFolder(), and indexGen::title.

Referenced by endJob().

                                                                                   {
   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::endcapWidthHist	(	std::string	name,
		std::string	title
	)

private

Definition at line 53 of file RPCMonitorSync.cc.

References DQMStore::book2D(), dbe, estimatePileup::hist, DQMStore::setCurrentFolder(), and indexGen::title.

Referenced by endJob().

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

void RPCMonitorSync::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 125 of file RPCMonitorSync.cc.

References barrelOffsetHist(), barrelWidthHist(), dbe, endcapOffsetHist(), endcapWidthHist(), MonitorElement::Fill(), RPCDetId::layer(), evf::evtn::offset(), RPCDetId::region(), RPCDetId::ring(), RPCDetId::roll(), RootFileName, DQMStore::save(), saveRootFile, RPCDetId::sector(), RPCDetId::station(), RPCDetId::subsector(), synchroMap, and tablePrinter::width.

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

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

private

Definition at line 217 of file RPCMonitorSync.cc.

References abs, bookDetUnitME(), counter, dbe, timing::early_all, edm::Event::getByType(), i, errorMatrix2Lands_multiChannel::id, timing::inTime, timing::late_all, RPCDetId::layer(), meCollection, RPCDetId::roll(), RootFileName, DQMStore::save(), saveRootFile, saveRootFileEventsInterval, RPCDetId::subsector(), and synchroMap.

Referenced by analyze().

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