SiPixelModuleStatus.cc

Go to the documentation of this file.

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

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

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

// ----------------------------------------------------------------------
SiPixelModuleStatus::SiPixelModuleStatus(int detId, int nrocs) : fDetid(detId), fNrocs(nrocs) {
  for (int i = 0; i < fNrocs; ++i) {
    SiPixelRocStatus a;
    fRocs.push_back(a);
  }
};

// ----------------------------------------------------------------------
SiPixelModuleStatus::~SiPixelModuleStatus(){};

// ----------------------------------------------------------------------
void SiPixelModuleStatus::fillDIGI(int iroc) {
  if (iroc < fNrocs)
    fRocs[iroc].fillDIGI();
}

// ----------------------------------------------------------------------
void SiPixelModuleStatus::updateDIGI(int iroc, unsigned int nhit) {
  if (iroc < fNrocs)
    fRocs[iroc].updateDIGI(nhit);
}

// ----------------------------------------------------------------------
void SiPixelModuleStatus::fillFEDerror25(PixelFEDChannel ch) {
  int roc_first = int(ch.roc_first);
  int roc_last = int(ch.roc_last);
  for (int iroc = 0; iroc < fNrocs; iroc++) {
    if (iroc >= roc_first && iroc <= roc_last) {
      fRocs[iroc].fillFEDerror25();
    }
  }
}

// ----------------------------------------------------------------------
unsigned int SiPixelModuleStatus::digiOccROC(int iroc) { return (iroc < fNrocs ? fRocs[iroc].digiOccROC() : -1); }

// ----------------------------------------------------------------------
unsigned int SiPixelModuleStatus::digiOccMOD() {
  unsigned int count(0);
  for (int iroc = 0; iroc < fNrocs; ++iroc) {
    count += digiOccROC(iroc);
  }
  return count;
}

// ----------------------------------------------------------------------
int SiPixelModuleStatus::detid() { return fDetid; }

// ----------------------------------------------------------------------
int SiPixelModuleStatus::nrocs() { return fNrocs; }

// ----------------------------------------------------------------------
void SiPixelModuleStatus::setNrocs(int iroc) { fNrocs = iroc; }

// ----------------------------------------------------------------------
double SiPixelModuleStatus::perRocDigiOcc() {
  unsigned int ave(0);
  for (int iroc = 0; iroc < fNrocs; ++iroc) {
    unsigned int inc = digiOccROC(iroc);
    ave += inc;
  }
  return (1.0 * ave) / fNrocs;
}

double SiPixelModuleStatus::perRocDigiOccVar() {
  double fModAverage = SiPixelModuleStatus::perRocDigiOcc();

  double sig = 1.0;
  for (int iroc = 0; iroc < fNrocs; ++iroc) {
    unsigned int inc = digiOccROC(iroc);
    sig += (fModAverage - inc) * (fModAverage - inc);
  }

  double fModSigma = sig / (fNrocs - 1);
  return TMath::Sqrt(fModSigma);
}

// ----------------------------------------------------------------------
// Be careful : return the address not the value of ROC status
SiPixelRocStatus* SiPixelModuleStatus::getRoc(int i) { return &fRocs[i]; }

// ----------------------------------------------------------------------
void SiPixelModuleStatus::updateModuleDIGI(int iroc, unsigned int nhits) { fRocs[iroc].updateDIGI(nhits); }

void SiPixelModuleStatus::updateModuleStatus(SiPixelModuleStatus newData) {
  bool isSameModule = true;
  if (fDetid != newData.detid() || fNrocs != newData.nrocs()) {
    isSameModule = false;
  }

  if (isSameModule) {
    for (int iroc = 0; iroc < fNrocs; ++iroc) {
      //update occupancy
      fRocs[iroc].updateDIGI(newData.digiOccROC(iroc));
      //update FEDerror25
      SiPixelRocStatus* rocStatus = newData.getRoc(iroc);
      if (rocStatus->isFEDerror25()) {
        fRocs[iroc].fillFEDerror25();
      }
    }

  }  // if same module
}