|
|
|
#include <RPCMonitorSync.h>
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) | |
| virtual void | endJob (void) |
| RPCMonitorSync (const edm::ParameterSet &) | |
| ~RPCMonitorSync () | |
Private Member Functions | |
| MonitorElement * | barrelOffsetHist (std::string name, std::string title) |
| Log messages. | |
| 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 | |
| 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 |
RPC Synchronization Monitoring Class
Anna Cimmino (INFN)
Definition at line 91 of file RPCMonitorSync.h.
| 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.
{
}
| 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);
}
int RPCMonitorSync::counter [private] |
Definition at line 113 of file RPCMonitorSync.h.
Referenced by analyze(), readRPCDAQStrips(), and RPCMonitorSync().
DQMStore* RPCMonitorSync::dbe [private] |
back-end interface
Definition at line 115 of file RPCMonitorSync.h.
Referenced by barrelOffsetHist(), barrelWidthHist(), bookDetUnitME(), endcapOffsetHist(), endcapWidthHist(), endJob(), readRPCDAQStrips(), and RPCMonitorSync().
MonitorElement* RPCMonitorSync::h1 [private] |
Definition at line 116 of file RPCMonitorSync.h.
std::map<uint32_t, std::map<std::string, MonitorElement*> > RPCMonitorSync::meCollection [private] |
Definition at line 117 of file RPCMonitorSync.h.
Referenced by readRPCDAQStrips().
std::string RPCMonitorSync::nameInLog [private] |
Definition at line 119 of file RPCMonitorSync.h.
Referenced by analyze(), and RPCMonitorSync().
std::string RPCMonitorSync::RootFileName [private] |
Definition at line 122 of file RPCMonitorSync.h.
Referenced by endJob(), readRPCDAQStrips(), and RPCMonitorSync().
bool RPCMonitorSync::saveRootFile [private] |
Definition at line 120 of file RPCMonitorSync.h.
Referenced by endJob(), readRPCDAQStrips(), and RPCMonitorSync().
int RPCMonitorSync::saveRootFileEventsInterval [private] |
Definition at line 121 of file RPCMonitorSync.h.
Referenced by readRPCDAQStrips(), and RPCMonitorSync().
std::map<uint32_t,timing> RPCMonitorSync::synchroMap [private] |
Definition at line 112 of file RPCMonitorSync.h.
Referenced by endJob(), and readRPCDAQStrips().
1.7.3