SiPixelHistoricInfoEDAClient Class Reference

#include <SiPixelHistoricInfoEDAClient.h>

Inheritance diagram for SiPixelHistoricInfoEDAClient:

Public Member Functions
	SiPixelHistoricInfoEDAClient (const edm::ParameterSet &)
	~SiPixelHistoricInfoEDAClient ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
virtual void	beginLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &)
virtual void	beginRun (const edm::Run &, const edm::EventSetup &)
virtual void	endJob ()
virtual void	endLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &)
virtual void	endRun (const edm::Run &, const edm::EventSetup &)
void	fillPerformanceSummary () const
void	retrieveMEs ()
void	saveFile (std::string filename) const
void	writeDB () const

Private Attributes
DQMStore *	dbe_
bool	firstEventInRun
SiPixelHistogramId	histogramManager
std::map< uint32_t, std::vector< MonitorElement * > >	mapOfdetIDtoMEs
int	nEventsInRun
std::string	outputDir_
edm::ParameterSet	parameterSet_
SiPixelPerformanceSummary *	performanceSummary
bool	printDebug_
bool	writeHisto_

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

Definition at line 19 of file SiPixelHistoricInfoEDAClient.h.

Constructor & Destructor Documentation

SiPixelHistoricInfoEDAClient::SiPixelHistoricInfoEDAClient ( const edm::ParameterSet &  parameterSet )

explicit

Definition at line 18 of file SiPixelHistoricInfoEDAClient.cc.

References dbe_, edm::ParameterSet::getUntrackedParameter(), and cppFunctionSkipper::operator.

                                                                                           { 
  dbe_ = Service<DQMStore>().operator->(); 
  
  parameterSet_ = parameterSet;  
  
  printDebug_ = parameterSet_.getUntrackedParameter<bool>("printDebug", false); 
  writeHisto_ = parameterSet_.getUntrackedParameter<bool>("writeHisto", false); 
  
  outputDir_ = parameterSet_.getUntrackedParameter<string>("outputDir", ".");
}

SiPixelHistoricInfoEDAClient::~SiPixelHistoricInfoEDAClient

(

)

Definition at line 30 of file SiPixelHistoricInfoEDAClient.cc.

{}

Member Function Documentation

void SiPixelHistoricInfoEDAClient::analyze ( const edm::Event &  event, const edm::EventSetup &  eventSetup  )

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 51 of file SiPixelHistoricInfoEDAClient.cc.

References edm::EventBase::time(), and edm::Timestamp::value().

                                                                                           {
  nEventsInRun++; 
  if (firstEventInRun) { 
    firstEventInRun = false; 

    performanceSummary->setTimeStamp(event.time().value()); 
  }
}

void SiPixelHistoricInfoEDAClient::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 33 of file SiPixelHistoricInfoEDAClient.cc.

{}

void SiPixelHistoricInfoEDAClient::beginLuminosityBlock ( const edm::LuminosityBlock &  luminosityBlock, const edm::EventSetup &  eventSetup  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 45 of file SiPixelHistoricInfoEDAClient.cc.

References edm::LuminosityBlockBase::luminosityBlock().

                                                                                         {
  performanceSummary->setLuminosityBlock(luminosityBlock.luminosityBlock());
}

void SiPixelHistoricInfoEDAClient::beginRun ( const edm::Run &  run, const edm::EventSetup &  eventSetup  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 36 of file SiPixelHistoricInfoEDAClient.cc.

References edm::RunBase::run().

                                                                                        {
  performanceSummary = new SiPixelPerformanceSummary(); 
  performanceSummary->setRunNumber(run.run());

  nEventsInRun = 0; 
  firstEventInRun = true; 
}

void SiPixelHistoricInfoEDAClient::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 80 of file SiPixelHistoricInfoEDAClient.cc.

References dbe_, and DQMStore::save().

                                          {
  if (writeHisto_) {
    ostringstream endJobOutputFile; 
    endJobOutputFile << outputDir_ << "/SiPixelHistoricInfoEDAClient_endJob.root"; 
    dbe_->save(endJobOutputFile.str());
  }
}

void SiPixelHistoricInfoEDAClient::endLuminosityBlock ( const edm::LuminosityBlock &  luminosityBlock, const edm::EventSetup &  eventSetup  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 61 of file SiPixelHistoricInfoEDAClient.cc.

                                                                                       {}

void SiPixelHistoricInfoEDAClient::endRun ( const edm::Run &  run, const edm::EventSetup &  eventSetup  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 65 of file SiPixelHistoricInfoEDAClient.cc.

References dbe_, edm::RunBase::run(), and DQMStore::save().

                                                                                      {
  performanceSummary->setNumberOfEvents(nEventsInRun);

  retrieveMEs();
  fillPerformanceSummary();
  writeDB(); 
  
  if (writeHisto_) {
    ostringstream endRunOutputFile; 
    endRunOutputFile << outputDir_ << "/SiPixelHistoricInfoEDAClient_" << run.run() <<".root"; 
    dbe_->save(endRunOutputFile.str()); 
  }
}

void SiPixelHistoricInfoEDAClient::fillPerformanceSummary ( ) const

private

Definition at line 163 of file SiPixelHistoricInfoEDAClient.cc.

References b, gather_cfg::cout, python.multivaluedict::map(), mathSSE::sqrt(), and v.

                                                                {
  for (map< uint32_t, vector<MonitorElement*> >::const_iterator iMEvec=mapOfdetIDtoMEs.begin(); 
       iMEvec!=mapOfdetIDtoMEs.end(); iMEvec++) {
    uint32_t theMEdetID = iMEvec->first;
    vector<MonitorElement*> theMEvector = iMEvec->second;
    
    if (printDebug_) { 
      cout << theMEdetID << ":"; for (vector<MonitorElement*>::const_iterator iME = theMEvector.begin(); 
                          iME!=theMEvector.end(); iME++) cout << (*iME)->getName() << ","; cout << endl; 
    }     
    for (vector<MonitorElement*>::const_iterator iME = theMEvector.begin(); iME!=theMEvector.end(); iME++) {
      string theMEname = (*iME)->getName();

      
      if (theMEdetID<50) {
        // from SiPixelMonitorRawData
        if (theMEname.find("errorType_siPixelDigis")!=string::npos) { 
      for (int v=25; v<40; v++) {
        int b = (*iME)->getTH1()->GetXaxis()->FindBin(v); 
        performanceSummary->setRawDataErrorType(theMEdetID, v-25, (*iME)->getBinContent(b)); 
      }
        } 
      } 
      // from SiPixelMonitorTrack
      else if (theMEdetID==77) {
        if (theMEname.find("ntracks")!=string::npos && theMEname.find("InPixVol")==string::npos) { 
          float trkMean=0.0, trkRMS=0.0; 
      if ((*iME)->getBinContent(1)>0.0) {
            trkMean = float((*iME)->getBinContent(3))/(*iME)->getBinContent(1); // Barrel regionID: 80
        trkRMS = sqrt(trkMean*(trkMean+1.0)/(*iME)->getBinContent(1)); 
            performanceSummary->setFractionOfTracks(80, trkMean, trkRMS); 

            trkMean = float((*iME)->getBinContent(4))/(*iME)->getBinContent(1); // Endcap regionID: 81
        trkRMS = sqrt(trkMean*(trkMean+1.0)/(*iME)->getBinContent(1)); 
            performanceSummary->setFractionOfTracks(81, trkMean, trkRMS); 
          }
    }
      } 
      else if (theMEdetID==78) { // OnTrack
        if (theMEname.find("nclusters_siPixelClusters")!=string::npos) { 
          performanceSummary->setNumberOfOnTrackClusters(80, (*iME)->getBinContent(2)); 
          performanceSummary->setNumberOfOnTrackClusters(81, (*iME)->getBinContent(3)); 
        }
        if (theMEname.find("charge_siPixelClusters_Barrel")!=string::npos) {
          performanceSummary->setClusterChargeOnTrack(80, (*iME)->getMean(), (*iME)->getRMS()); 
        }
        if (theMEname.find("charge_siPixelClusters_Endcap")!=string::npos) {
          performanceSummary->setClusterChargeOnTrack(81, (*iME)->getMean(), (*iME)->getRMS()); 
        }
        if (theMEname.find("size_siPixelClusters_Barrel")!=string::npos) {
          performanceSummary->setClusterSizeOnTrack(80, (*iME)->getMean(), (*iME)->getRMS()); 
        }
        if (theMEname.find("size_siPixelClusters_Endcap")!=string::npos) {
          performanceSummary->setClusterSizeOnTrack(81, (*iME)->getMean(), (*iME)->getRMS()); 
        }
      }
      else if (theMEdetID==79) { // OffTrack
        if (theMEname.find("nclusters_siPixelClusters")!=string::npos) { 
          performanceSummary->setNumberOfOffTrackClusters(80, (*iME)->getBinContent(2)); 
          performanceSummary->setNumberOfOffTrackClusters(81, (*iME)->getBinContent(3)); 
        }
        if (theMEname.find("charge_siPixelClusters_Barrel")!=string::npos) {
          performanceSummary->setClusterChargeOffTrack(80, (*iME)->getMean(), (*iME)->getRMS()); 
        }
        if (theMEname.find("charge_siPixelClusters_Endcap")!=string::npos) {
          performanceSummary->setClusterChargeOffTrack(81, (*iME)->getMean(), (*iME)->getRMS()); 
        }
        if (theMEname.find("size_siPixelClusters_Barrel")!=string::npos) {
          performanceSummary->setClusterSizeOffTrack(80, (*iME)->getMean(), (*iME)->getRMS()); 
        }
        if (theMEname.find("size_siPixelClusters_Endcap")!=string::npos) {
          performanceSummary->setClusterSizeOffTrack(81, (*iME)->getMean(), (*iME)->getRMS()); 
        }
      }
      else {
        // from SiPixelMonitorDigi 
        if (theMEname.find("ndigis_siPixelDigis")!=string::npos) { 
          performanceSummary->setNumberOfDigis(theMEdetID, (*iME)->getMean(), (*iME)->getRMS(), (*iME)->getEntries()==0.0);
        } 
        if (theMEname.find("adc_siPixelDigis")!=string::npos) { 
          performanceSummary->setADC(theMEdetID, (*iME)->getMean(), (*iME)->getRMS(), (*iME)->getEntries()==0.0);
        } 
        // from SiPixelMonitorCluster
        if (theMEname.find("nclusters_siPixelClusters")!=string::npos) {
          performanceSummary->setNumberOfClusters(theMEdetID, (*iME)->getMean(), (*iME)->getRMS(), (*iME)->getEntries()==0.0);
        }
        if (theMEname.find("charge_siPixelClusters")!=string::npos) {
          performanceSummary->setClusterCharge(theMEdetID, (*iME)->getMean(), (*iME)->getRMS(), (*iME)->getEntries()==0.0);
        } 
        if (theMEname.find("size_siPixelClusters")!=string::npos) {
          performanceSummary->setClusterSize(theMEdetID, (*iME)->getMean(), (*iME)->getRMS(), (*iME)->getEntries()==0.0);
        } 
        if (theMEname.find("sizeX_siPixelClusters")!=string::npos) {
          performanceSummary->setClusterSizeX(theMEdetID, (*iME)->getMean(), (*iME)->getRMS(), (*iME)->getEntries()==0.0);
        }
        if (theMEname.find("sizeY_siPixelClusters")!=string::npos) {
          performanceSummary->setClusterSizeY(theMEdetID, (*iME)->getMean(), (*iME)->getRMS(), (*iME)->getEntries()==0.0);
        }
        // from SiPixelMonitorRecHit
        if (theMEname.find("nRecHits_siPixelRecHits")!=string::npos) {
          performanceSummary->setNumberOfRecHits(theMEdetID, (*iME)->getMean(), (*iME)->getRMS(), (*iME)->getEntries()==0.0);
        }
        // from SiPixelMonitorTrack
        if (theMEname.find("residualX_siPixelTrackResiduals")!=string::npos) {
          performanceSummary->setResidualX(theMEdetID, (*iME)->getMean(), (*iME)->getRMS(), (*iME)->getEntries()==0.0);
        }
        if (theMEname.find("residualY_siPixelTrackResiduals")!=string::npos) {
          performanceSummary->setResidualY(theMEdetID, (*iME)->getMean(), (*iME)->getRMS(), (*iME)->getEntries()==0.0);
        } 
        // temporary solutions
        if (theMEname.find("hitmap_siPixelClusters")!=std::string::npos || 
            theMEname.find("hitmap_siPixelDigis")!=std::string::npos) { // if digi map exists, use it; else use cluster map
      int nNoisyCells=0, nEmptyCells=0;             //             not use xypos_siPixelRecHits
      for (int xBin=0; xBin<(*iME)->getNbinsX(); xBin++) {      // 4-1 pixel per bin
        for (int yBin=0; yBin<(*iME)->getNbinsY(); yBin++) { 
          if ((*iME)->getBinContent(xBin+1, yBin+1)>0.01*(*iME)->getEntries()) nNoisyCells++; 
          if ((*iME)->getBinContent(xBin+1, yBin+1)==.0 && (*iME)->getBinError(xBin+1, yBin+1)==0.0) nEmptyCells++; 
        } 
      } 
          performanceSummary->setNumberOfNoisCells(theMEdetID, float(nNoisyCells)); 
          performanceSummary->setNumberOfDeadCells(theMEdetID, float(nEmptyCells)); 
        } 
        // performanceSummary->setNumberOfPixelHitsInTrackFit(theMEdetId, float(nPixelHits)); 
      }
    }
  }
}

void SiPixelHistoricInfoEDAClient::retrieveMEs ( )

private

Definition at line 89 of file SiPixelHistoricInfoEDAClient.cc.

References gather_cfg::cout, dbe_, DQMStore::get(), DQMStore::getContents(), MonitorElement::getName(), and edm::second().

                                               {
  mapOfdetIDtoMEs.clear(); 
  
  vector<string> listOfMEswithFullPath;
  dbe_->getContents(listOfMEswithFullPath);
    
  for (vector<string>::const_iterator iMEstr = listOfMEswithFullPath.begin(); 
       iMEstr!=listOfMEswithFullPath.end(); iMEstr++) {
    if (printDebug_) cout << iMEstr->data() << endl; 
      
    size_t pathLength = iMEstr->find(":",0);     
    string thePath = iMEstr->substr(0, pathLength); 
    string allHists = iMEstr->substr(pathLength+1); 
        
    if (thePath.find("Pixel",0)!=string::npos) { 
      if (thePath.find("Track",0)!=string::npos) { // for Pixel/Tracks, Pixel/Clusters/On,OffTrack
                            uint32_t newMEdetID = 77; 
    if (thePath.find("On", 0)!=string::npos) newMEdetID = 78; 
    if (thePath.find("Off",0)!=string::npos) newMEdetID = 79; 
    
    size_t histnameLength;
    do {
      histnameLength = allHists.find(",",0);
      string theHist;
      if (histnameLength!=string::npos) {
        theHist = allHists.substr(0, histnameLength);
        allHists.erase(0, histnameLength+1);
      }
      else theHist = allHists; 

      string fullPathHist = thePath + "/" + theHist; 
      MonitorElement* newME = dbe_->get(fullPathHist); 
      if (newME) {
        if (mapOfdetIDtoMEs.find(newMEdetID)==mapOfdetIDtoMEs.end()) {
          vector<MonitorElement*> newMEvector;
          newMEvector.push_back(newME);
          mapOfdetIDtoMEs.insert(make_pair(newMEdetID, newMEvector));
        }
        else (mapOfdetIDtoMEs.find(newMEdetID)->second).push_back(newME);
      } 
    } 
    while (histnameLength!=string::npos); 
      }
      else {
        size_t histnameLength; 
        do {
          histnameLength = allHists.find(",",0);
          string theHist;
          if (histnameLength!=string::npos) {
            theHist = allHists.substr(0, histnameLength);
            allHists.erase(0, histnameLength+1);
          }
          else theHist = allHists; 

          string fullPathHist = thePath + "/" + theHist;
          MonitorElement* newME = dbe_->get(fullPathHist);
          if (newME) {
            uint32_t newMEdetID = histogramManager.getRawId(newME->getName());

            if (mapOfdetIDtoMEs.find(newMEdetID)==mapOfdetIDtoMEs.end()) {
              vector<MonitorElement*> newMEvector;
              newMEvector.push_back(newME);
              mapOfdetIDtoMEs.insert(make_pair(newMEdetID, newMEvector));
            }
            else ((mapOfdetIDtoMEs.find(newMEdetID))->second).push_back(newME);
          }
        } 
        while (histnameLength!=string::npos); 
      }
    } 
  } 
}

void SiPixelHistoricInfoEDAClient::saveFile ( std::string  filename ) const

inlineprivate

Definition at line 36 of file SiPixelHistoricInfoEDAClient.h.

References dbe_, and DQMStore::save().

{ dbe_->save(filename); }

void SiPixelHistoricInfoEDAClient::writeDB ( ) const

private

Definition at line 292 of file SiPixelHistoricInfoEDAClient.cc.

References cond::service::PoolDBOutputService::appendSinceTime(), cond::service::PoolDBOutputService::beginOfTime(), gather_cfg::cout, cond::service::PoolDBOutputService::createNewIOV(), cond::service::PoolDBOutputService::currentTime(), cond::service::PoolDBOutputService::endOfTime(), edm::Service< T >::isAvailable(), and cond::service::PoolDBOutputService::isNewTagRequest().

                                                 {  
  if (printDebug_) performanceSummary->printAll(); 
  else performanceSummary->print(); 
  cout << "SiPixelHistoricInfoEDAClient::writeDB()" << endl; 

  Service<cond::service::PoolDBOutputService> mydbservice; 
  if (mydbservice.isAvailable()) {
    if (mydbservice->isNewTagRequest("SiPixelPerformanceSummaryRcd")) {
      mydbservice->createNewIOV<SiPixelPerformanceSummary>(performanceSummary, 
                                                           mydbservice->beginOfTime(), 
                               mydbservice->endOfTime(), 
                              "SiPixelPerformanceSummaryRcd"); 
    } 
    else {
      mydbservice->appendSinceTime<SiPixelPerformanceSummary>(performanceSummary, 
                                                              mydbservice->currentTime(), 
                                 "SiPixelPerformanceSummaryRcd"); 
    }
  }
  else LogError("writeDB") << "service unavailable" << endl; 
}

Member Data Documentation

DQMStore* SiPixelHistoricInfoEDAClient::dbe_

private

Definition at line 44 of file SiPixelHistoricInfoEDAClient.h.

Referenced by saveFile().

bool SiPixelHistoricInfoEDAClient::firstEventInRun

private

Definition at line 46 of file SiPixelHistoricInfoEDAClient.h.

SiPixelHistogramId SiPixelHistoricInfoEDAClient::histogramManager

private

Definition at line 49 of file SiPixelHistoricInfoEDAClient.h.

std::map< uint32_t, std::vector<MonitorElement*> > SiPixelHistoricInfoEDAClient::mapOfdetIDtoMEs

private

Definition at line 50 of file SiPixelHistoricInfoEDAClient.h.

int SiPixelHistoricInfoEDAClient::nEventsInRun

private

Definition at line 47 of file SiPixelHistoricInfoEDAClient.h.

std::string SiPixelHistoricInfoEDAClient::outputDir_

private

Definition at line 41 of file SiPixelHistoricInfoEDAClient.h.

Referenced by valtools.webpage::readCaptions().

edm::ParameterSet SiPixelHistoricInfoEDAClient::parameterSet_

private

Definition at line 43 of file SiPixelHistoricInfoEDAClient.h.

SiPixelPerformanceSummary* SiPixelHistoricInfoEDAClient::performanceSummary

private

Definition at line 51 of file SiPixelHistoricInfoEDAClient.h.

bool SiPixelHistoricInfoEDAClient::printDebug_

private

Definition at line 39 of file SiPixelHistoricInfoEDAClient.h.

bool SiPixelHistoricInfoEDAClient::writeHisto_

private

Definition at line 40 of file SiPixelHistoricInfoEDAClient.h.