SiPixelDetectorStatus.cc

Go to the documentation of this file.

#include <iostream>
#include <fstream>
#include <map>
#include <cmath>
#include <vector>

#include <TString.h>
#include <TMath.h>
#include <TH1.h>

#include "CalibTracker/SiPixelQuality/interface/SiPixelDetectorStatus.h"

// ----------------------------------------------------------------------
SiPixelDetectorStatus::SiPixelDetectorStatus(): fLS0(99999999), fLS1(0), fRun0(99999999), fRun1(0) {

  fDetHits = 0;
  fNevents = 0;

}

// ----------------------------------------------------------------------
SiPixelDetectorStatus::~SiPixelDetectorStatus() {

}


// ----------------------------------------------------------------------
void SiPixelDetectorStatus::readFromFile(std::string filename) {

  std::ifstream INS;
  std::string sline;
  INS.open(filename.c_str());

  int oldDetId(-1);
  int detid(0), roc(0), hits(0), nroc(0);
  SiPixelModuleStatus *pMod(nullptr);
  bool readOK(false);
  while (std::getline(INS, sline)) {

    if (std::string::npos != sline.find("# SiPixelDetectorStatus START")) {
      readOK = true;
      continue;
    }
    if (!readOK) continue;

    if (std::string::npos != sline.find("# SiPixelDetectorStatus END")) {
      pMod->setNrocs(nroc+1);
      break;
    }

    if (sline.find("# SiPixelDetectorStatus for LS") != std::string::npos) {
      std::sscanf(sline.c_str(), "# SiPixelDetectorStatus for LS  %d .. %d", &fLS0, &fLS1);
      continue;
    }
    if (sline.find("# SiPixelDetectorStatus for run") != std::string::npos) {
      std::sscanf(sline.c_str(), "# SiPixelDetectorStatus for run %d .. %d", &fRun0, &fRun1);
      continue;
    }
    if (sline.find("# SiPixelDetectorStatus total hits = ") != std::string::npos) {
      std::sscanf(sline.c_str(), "# SiPixelDetectorStatus total hits = %ld", &fDetHits);
      continue;
    }

    std::sscanf(sline.c_str(), "%d %d %d", &detid, &roc, &hits);
    if (roc > nroc) nroc = roc;
    if (detid != oldDetId) {
      if (pMod) {
    pMod->setNrocs(nroc+1);
      }

      oldDetId = detid;
      if (getModule(detid) == nullptr) {
    addModule(detid,nroc+1);
      } 

      pMod = getModule(detid);
      nroc = 0;
    }
    if (pMod) {
      fDetHits += hits;
      pMod->updateModuleDIGI(roc, hits);
    }

  }

  INS.close();

}


// ----------------------------------------------------------------------
void SiPixelDetectorStatus::dumpToFile(std::string filename) {

  std::ofstream OD(filename.c_str());
  OD << "# SiPixelDetectorStatus START" << std::endl;
  OD << "# SiPixelDetectorStatus for LS  " << fLS0 << " .. " << fLS1 << std::endl;
  OD << "# SiPixelDetectorStatus for run " << fRun0 << " .. " << fRun1 << std::endl;
  OD << "# SiPixelDetectorStatus total hits = " << fDetHits << std::endl;

  for (std::map<int, SiPixelModuleStatus>::iterator it = SiPixelDetectorStatus::begin(); it != SiPixelDetectorStatus::end(); ++it) {
    for (int iroc = 0; iroc < it->second.nrocs(); ++iroc) {
      for (int idc = 0; idc < 26; ++idc) {
    OD << Form("%10d %2d %3d", it->first, iroc, int(it->second.getRoc(iroc)->digiOccROC())) << std::endl;
      }
    }
  }
  OD << "# SiPixelDetectorStatus END" << std::endl;
  OD.close();

}


// ----------------------------------------------------------------------
void SiPixelDetectorStatus::addModule(int detid, int nrocs) {

     SiPixelModuleStatus a(detid, nrocs);
     fModules.insert(std::make_pair(detid, a));

}

// ----------------------------------------------------------------------
void SiPixelDetectorStatus::addModule(int detid, SiPixelModuleStatus a) {

     fModules.insert(std::make_pair(detid, a));

}


// ----------------------------------------------------------------------
void SiPixelDetectorStatus::fillDIGI(int detid, int roc) {

     ++fDetHits;
     fModules[detid].fillDIGI(roc);

}

// ----------------------------------------------------------------------
void SiPixelDetectorStatus::fillFEDerror25(int detid,PixelFEDChannel ch){

   if (fModules.find(detid) != fModules.end()){
        fModules[detid].fillFEDerror25(ch);
   }

}

// FEDerror25 effected ROCs in for each module
std::map<int, std::vector<int>> SiPixelDetectorStatus::getFEDerror25Rocs(){

    std::map<int, std::vector<int>> badRocLists_;

    for(std::map<int, SiPixelModuleStatus>::iterator itMod = SiPixelDetectorStatus::begin(); itMod != SiPixelDetectorStatus::end(); ++itMod)
    {
         int detid = itMod->first;
         // FEDerror25 effected ROCs in a given module
         std::vector<int> list;
         SiPixelModuleStatus modStatus = itMod->second;
         for (int iroc = 0; iroc < modStatus.nrocs(); ++iroc) {

              SiPixelRocStatus* roc = modStatus.getRoc(iroc);
              if(roc->isFEDerror25()){
                 list.push_back(iroc);
              }
              badRocLists_[detid]=list;
         }

    }

    return badRocLists_;
}

// ----------------------------------------------------------------------
std::map<int, SiPixelModuleStatus>::iterator SiPixelDetectorStatus::begin() {

  return fModules.begin();

}

// ----------------------------------------------------------------------
//map<int, SiPixelModuleStatus>::iterator SiPixelDetectorStatus::next() {
//  return fNext++;
//}

// ----------------------------------------------------------------------
std::map<int, SiPixelModuleStatus>::iterator SiPixelDetectorStatus::end() {

  return fModules.end();

}

// ----------------------------------------------------------------------
int SiPixelDetectorStatus::nmodules() {

  return static_cast<int>(fModules.size());

}

// ----------------------------------------------------------------------
SiPixelModuleStatus* SiPixelDetectorStatus::getModule(int detid) {

  if (fModules.find(detid) == fModules.end()) {
    return nullptr;
  }
  return &(fModules[detid]);

}

bool SiPixelDetectorStatus::findModule(int detid) {

  if (fModules.find(detid) == fModules.end()) 
    return false;
  else
    return true;

}

// ----------------------------------------------------------------------
double SiPixelDetectorStatus::perRocDigiOcc() {

  unsigned long int ave(0);
  int nrocs(0);
  for (std::map<int, SiPixelModuleStatus>::iterator it = SiPixelDetectorStatus::begin(); it != SiPixelDetectorStatus::end(); ++it) {
    unsigned long int inc = it->second.digiOccMOD();
    ave   += inc;
    nrocs += it->second.nrocs();
  }
  return (1.0*ave)/nrocs;

}

double SiPixelDetectorStatus::perRocDigiOccVar(){

  double fDetAverage = SiPixelDetectorStatus::perRocDigiOcc();

  double sig = 0.0;
  int nrocs(0);
  for(std::map<int, SiPixelModuleStatus>::iterator it = SiPixelDetectorStatus::begin(); it != SiPixelDetectorStatus::end(); ++it) {
    unsigned long int inc = it->second.digiOccMOD();
    sig += (fDetAverage - inc) * (fDetAverage - inc);
    nrocs += it->second.nrocs();
  }

  double fDetSigma   = sig/(nrocs - 1);
  return TMath::Sqrt(fDetSigma);
}

// combine status from different data (coming from different run/lumi)
void SiPixelDetectorStatus::updateDetectorStatus(SiPixelDetectorStatus newData){

  // loop over new data status
  for(std::map<int, SiPixelModuleStatus>::iterator it = newData.begin(); it != newData.end(); ++it) {
       
       int detid = it->first;
       if(fModules.find(detid) != fModules.end()){// if the detid is in the module lists
          fModules[detid].updateModuleStatus( *(newData.getModule(detid)) );
       }
       else{
          fModules.insert(std::make_pair(detid, *(newData.getModule(detid))));
       }

  }

  fDetHits = fDetHits + newData.digiOccDET();
  fNevents = fNevents + newData.getNevents();

}