SiPixelStatusManager Class Reference

#include <SiPixelStatusManager.h>

Public Types
typedef std::map < edm::LuminosityBlockNumber_t, std::map< int, std::vector < int > > >::iterator	FEDerror25Map_iterator
typedef std::map < edm::LuminosityBlockNumber_t, SiPixelDetectorStatus > ::iterator	siPixelStatusMap_iterator
typedef std::vector < SiPixelDetectorStatus > ::iterator	siPixelStatusVtr_iterator

Public Member Functions
void	createPayloads ()
const std::map < edm::LuminosityBlockNumber_t, SiPixelDetectorStatus > &	getBadComponents ()
const std::map < edm::LuminosityBlockNumber_t, std::map< int, std::vector < int > > > &	getFEDerror25Rocs ()
void	readLumi (const edm::LuminosityBlock &)
void	reset ()
	SiPixelStatusManager ()
	SiPixelStatusManager (const edm::ParameterSet &, edm::ConsumesCollector &&)
virtual	~SiPixelStatusManager ()

Private Member Functions
void	createBadComponents ()
void	createFEDerror25 ()

Static Private Member Functions
static bool	rankByLumi (SiPixelDetectorStatus status1, SiPixelDetectorStatus status2)

Private Attributes
int	aveDigiOcc_
std::map < edm::LuminosityBlockNumber_t, std::map< int, std::vector < int > > >	FEDerror25Map_
std::string	label_
std::string	moduleName_
int	nLumi_
std::string	outputBase_
std::map < edm::LuminosityBlockNumber_t, SiPixelDetectorStatus >	siPixelStatusMap_
edm::EDGetTokenT < SiPixelDetectorStatus >	siPixelStatusToken_
std::vector < SiPixelDetectorStatus >	siPixelStatusVtr_

Detailed Description

Author

helper class that set up IOV strcutre of SiPixelDetectorStatus

Author

Definition at line 24 of file SiPixelStatusManager.h.

Member Typedef Documentation

typedef std::map<edm::LuminosityBlockNumber_t, std::map<int, std::vector<int> > >::iterator SiPixelStatusManager::FEDerror25Map_iterator

Definition at line 41 of file SiPixelStatusManager.h.

typedef std::map<edm::LuminosityBlockNumber_t, SiPixelDetectorStatus>::iterator SiPixelStatusManager::siPixelStatusMap_iterator

Definition at line 40 of file SiPixelStatusManager.h.

typedef std::vector<SiPixelDetectorStatus>::iterator SiPixelStatusManager::siPixelStatusVtr_iterator

Definition at line 42 of file SiPixelStatusManager.h.

Constructor & Destructor Documentation

SiPixelStatusManager::SiPixelStatusManager

(

)

Definition at line 21 of file SiPixelStatusManager.cc.

{}

SiPixelStatusManager::SiPixelStatusManager	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &&	iC
	)

Definition at line 24 of file SiPixelStatusManager.cc.

References edm::InLumi, label_, moduleName_, outputBase_, reset(), and siPixelStatusToken_.

    : outputBase_(iConfig.getParameter<ParameterSet>("SiPixelStatusManagerParameters")
                      .getUntrackedParameter<std::string>("outputBase")),
      aveDigiOcc_(iConfig.getParameter<edm::ParameterSet>("SiPixelStatusManagerParameters")
                      .getUntrackedParameter<int>("aveDigiOcc")),
      nLumi_(iConfig.getParameter<edm::ParameterSet>("SiPixelStatusManagerParameters")
                 .getUntrackedParameter<int>("resetEveryNLumi")),
      moduleName_(iConfig.getParameter<ParameterSet>("SiPixelStatusManagerParameters")
                      .getUntrackedParameter<std::string>("moduleName")),
      label_(iConfig.getParameter<ParameterSet>("SiPixelStatusManagerParameters")
                 .getUntrackedParameter<std::string>("label")) {
  edm::InputTag siPixelStatusTag_(moduleName_, label_);
  siPixelStatusToken_ = iC.consumes<SiPixelDetectorStatus, edm::InLumi>(siPixelStatusTag_);

  LogInfo("SiPixelStatusManager") << "Output base: " << outputBase_ << std::endl;
  reset();
}

SiPixelStatusManager::~SiPixelStatusManager ( )

virtual

Definition at line 43 of file SiPixelStatusManager.cc.

{}

Member Function Documentation

void SiPixelStatusManager::createBadComponents ( )

private

Definition at line 88 of file SiPixelStatusManager.cc.

References aveDigiOcc_, counter, GetRecoTauVFromDQM_MC_cff::next, nLumi_, outputBase_, SiPixelDetectorStatus::perRocDigiOcc(), SiPixelDetectorStatus::setLSRange(), siPixelStatusMap_, siPixelStatusVtr_, and SiPixelDetectorStatus::updateDetectorStatus().

Referenced by createPayloads().

                                               {
  siPixelStatusVtr_iterator firstStatus = siPixelStatusVtr_.begin();
  siPixelStatusVtr_iterator lastStatus = siPixelStatusVtr_.end();

  siPixelStatusMap_.clear();

  // doesn't work for nLumi_=1 cos any integer can be completely divided by 1
  if (outputBase_ == "nLumibased" && nLumi_ > 1) {
    // if the total number of Lumi Blocks can't be completely divided by nLumi_,
    // the residual Lumi Blocks will be as the last IOV
    int iterationLumi = 0;

    LuminosityBlockNumber_t tmpLumi;
    SiPixelDetectorStatus tmpSiPixelStatus;
    for (siPixelStatusVtr_iterator it = firstStatus; it != lastStatus; it++) {
      // this is the begining of an IOV
      if (iterationLumi % nLumi_ == 0) {
        tmpLumi = edm::LuminosityBlockNumber_t(it->getLSRange().first);
        tmpSiPixelStatus = (*it);
      }

      // keep update detector status up to nLumi_ lumi sections
      if (iterationLumi % nLumi_ > 0) {
        tmpSiPixelStatus.updateDetectorStatus((*it));
        tmpSiPixelStatus.setLSRange(int(tmpLumi), (*it).getLSRange().second);
      }

      siPixelStatusVtr_iterator currentIt = it;
      siPixelStatusVtr_iterator nextIt = std::next(currentIt);
      // wirte out if current lumi is the last lumi-section in the IOV
      if (iterationLumi % nLumi_ == nLumi_ - 1 || nextIt == lastStatus) {
        // fill it into a new map (with IOV structured)
        siPixelStatusMap_[tmpLumi] = tmpSiPixelStatus;
      }

      iterationLumi = iterationLumi + 1;
    }

    // check whether there is not enough number of Lumi in the last IOV
    // (only when siPixelStatusVtr_.size() > nLumi_ or equivalently current siPixelStatusMap_.size()>1
    //            (otherwise there will be only one IOV, and not previous IOV before the last IOV)
    //            and the number of lumi can not be completely divided by the nLumi_.
    //                (then the number of lumis in the last IOV is equal to the residual, which is less than nLumi_)
    // if it is, combine last IOV with the IOV before it
    if (siPixelStatusVtr_.size() % nLumi_ != 0 && siPixelStatusMap_.size() > 1) {
      // start from the iterator of the end of std::map
      siPixelStatusMap_iterator iterEnd = siPixelStatusMap_.end();
      // the last IOV
      siPixelStatusMap_iterator iterLastIOV = std::prev(iterEnd);
      // the IOV before the last IOV
      siPixelStatusMap_iterator iterBeforeLastIOV = std::prev(iterLastIOV);

      // combine the last IOV data to the IOV before the last IOV
      (iterBeforeLastIOV->second).updateDetectorStatus(iterLastIOV->second);
      (iterBeforeLastIOV->second)
          .setLSRange((iterBeforeLastIOV->second).getLSRange().first, (iterLastIOV->second).getLSRange().second);

      // delete the last IOV, so the IOV before the last IOV becomes the new last IOV
      siPixelStatusMap_.erase(iterLastIOV);
    }

  } else if (outputBase_ == "dynamicLumibased") {
    double aveDigiOcc = 1.0 * aveDigiOcc_;

    edm::LuminosityBlockNumber_t tmpLumi;
    SiPixelDetectorStatus tmpSiPixelStatus;
    bool isNewIOV = true;

    int counter = 0;
    for (siPixelStatusVtr_iterator it = firstStatus; it != lastStatus; it++) {
      if (isNewIOV) {  // if it is new IOV, init with the current data
        tmpLumi = edm::LuminosityBlockNumber_t(it->getLSRange().first);
        tmpSiPixelStatus = (*it);
      } else {  // if it is not new IOV, append current data
        tmpSiPixelStatus.updateDetectorStatus((*it));
        tmpSiPixelStatus.setLSRange(int(tmpLumi), (*it).getLSRange().second);
      }

      // if reaching the end of data, write the last IOV to the map whatsoevec
      siPixelStatusVtr_iterator currentIt = it;
      siPixelStatusVtr_iterator nextIt = std::next(currentIt);
      if (tmpSiPixelStatus.perRocDigiOcc() < aveDigiOcc && nextIt != lastStatus) {
        isNewIOV = false;  // if digi occ is not enough, next data will not belong to new IOV
      } else {             // if (accunumated) digi occ is enough, write the data to the map
        isNewIOV = true;
        siPixelStatusMap_[tmpLumi] = tmpSiPixelStatus;
        // so next loop is the begining of a new IOV
      }
      counter++;

    }  // end of siPixelStatusMap

    // check whether last IOV has enough statistics
    // (ONLY when there are more than oneIOV(otherwise there is NO previous IOV before the last IOV) )
    // if not, combine with previous IOV
    if (siPixelStatusMap_.size() > 1) {
      // start from the end iterator of the std::map
      siPixelStatusMap_iterator iterEnd = siPixelStatusMap_.end();
      // the last IOV
      siPixelStatusMap_iterator iterLastIOV = std::prev(iterEnd);
      // if the statistics of the last IOV is not enough
      if ((iterLastIOV->second).perRocDigiOcc() < aveDigiOcc) {
        // the IOV before the last IOV of the map
        siPixelStatusMap_iterator iterBeforeLastIOV = std::prev(iterLastIOV);
        // combine the last IOV data to the IOV before the last IOV
        (iterBeforeLastIOV->second).updateDetectorStatus(iterLastIOV->second);
        (iterBeforeLastIOV->second)
            .setLSRange((iterBeforeLastIOV->second).getLSRange().first, (iterLastIOV->second).getLSRange().second);
        // erase the last IOV, so the IOV before the last IOV becomes the new last IOV
        siPixelStatusMap_.erase(iterLastIOV);
      }
    }

  } else if (outputBase_ == "runbased" || ((int(siPixelStatusVtr_.size()) <= nLumi_ && outputBase_ == "nLumibased"))) {
    edm::LuminosityBlockNumber_t tmpLumi = edm::LuminosityBlockNumber_t(firstStatus->getLSRange().first);
    SiPixelDetectorStatus tmpSiPixelStatus = (*firstStatus);

    siPixelStatusVtr_iterator nextStatus = ++siPixelStatusVtr_.begin();
    for (siPixelStatusVtr_iterator it = nextStatus; it != lastStatus; it++) {
      tmpSiPixelStatus.updateDetectorStatus((*it));
      tmpSiPixelStatus.setLSRange(int(tmpLumi), (*it).getLSRange().second);
    }

    siPixelStatusMap_[tmpLumi] = tmpSiPixelStatus;

  } else {
    LogInfo("SiPixelStatusManager") << "Unrecognized payload outputBase parameter: " << outputBase_ << endl;
  }
}

void SiPixelStatusManager::createFEDerror25 ( )

private

Definition at line 218 of file SiPixelStatusManager.cc.

References FEDerror25Map_, SiPixelDetectorStatus::getFEDerror25Rocs(), GetRecoTauVFromDQM_MC_cff::next, and siPixelStatusVtr_.

Referenced by createPayloads().

                                            {
  // initialize the first IOV and SiPixelDetector status (in the first IOV)
  siPixelStatusVtr_iterator firstStatus = siPixelStatusVtr_.begin();
  edm::LuminosityBlockNumber_t firstLumi = edm::LuminosityBlockNumber_t(firstStatus->getLSRange().first);
  SiPixelDetectorStatus firstFEDerror25 = (*firstStatus);
  FEDerror25Map_[firstLumi] = firstFEDerror25.getFEDerror25Rocs();

  siPixelStatusVtr_iterator lastStatus = siPixelStatusVtr_.end();

  bool sameAsLastIOV = true;
  edm::LuminosityBlockNumber_t previousLumi = firstLumi;

  siPixelStatusVtr_iterator secondStatus = std::next(siPixelStatusVtr_.begin());
  for (siPixelStatusVtr_iterator it = secondStatus; it != lastStatus; it++) {
    // init for each lumi section (iterator)
    edm::LuminosityBlockNumber_t tmpLumi = edm::LuminosityBlockNumber_t(it->getLSRange().first);
    SiPixelDetectorStatus tmpFEDerror25 = (*it);

    std::map<int, std::vector<int> > tmpBadRocLists = tmpFEDerror25.getFEDerror25Rocs();

    std::map<int, SiPixelModuleStatus>::iterator itModEnd = tmpFEDerror25.end();
    for (std::map<int, SiPixelModuleStatus>::iterator itMod = tmpFEDerror25.begin(); itMod != itModEnd; ++itMod) {
      int detid = itMod->first;
      // if the badroc list differs for any detid, update the payload
      if (tmpBadRocLists[detid] != (FEDerror25Map_[previousLumi])[detid]) {
        sameAsLastIOV = false;
        break;  // jump out of the loop once a new payload is found
      }
    }

    if (sameAsLastIOV == false) {
      //only write new IOV when this Lumi's FEDerror25 ROC list is not equal to the previous one
      FEDerror25Map_[tmpLumi] = tmpBadRocLists;
      // and reset
      previousLumi = tmpLumi;
      sameAsLastIOV = true;
    }
  }
}

void SiPixelStatusManager::createPayloads

(

)

Definition at line 56 of file SiPixelStatusManager.cc.

References createBadComponents(), createFEDerror25(), rankByLumi(), and siPixelStatusVtr_.

Referenced by SiPixelStatusHarvester::dqmEndRun().

                                          {
  //only create std::map payloads when the number of non-zero DIGI lumi sections is greater than ZERO otherwise segmentation fault
  if (!siPixelStatusVtr_.empty()) {
    // sort the vector according to lumi
    std::sort(siPixelStatusVtr_.begin(), siPixelStatusVtr_.end(), SiPixelStatusManager::rankByLumi);

    // create FEDerror25 ROCs and bad ROCs from PCL
    SiPixelStatusManager::createFEDerror25();
    SiPixelStatusManager::createBadComponents();

    // realse the cost of siPixelStatusVtr_ since it is not needed anymore
    siPixelStatusVtr_.clear();
  }
}

const std::map<edm::LuminosityBlockNumber_t, SiPixelDetectorStatus>& SiPixelStatusManager::getBadComponents ( )

inline

Definition at line 35 of file SiPixelStatusManager.h.

References siPixelStatusMap_.

Referenced by SiPixelStatusHarvester::dqmEndRun().

{ return siPixelStatusMap_; }

const std::map<edm::LuminosityBlockNumber_t, std::map<int, std::vector<int> > >& SiPixelStatusManager::getFEDerror25Rocs ( )

inline

Definition at line 36 of file SiPixelStatusManager.h.

References FEDerror25Map_.

Referenced by SiPixelStatusHarvester::dqmEndRun().

                                                                                             {
    return FEDerror25Map_;
  }

bool SiPixelStatusManager::rankByLumi ( SiPixelDetectorStatus  status1, SiPixelDetectorStatus  status2  )

staticprivate

Definition at line 52 of file SiPixelStatusManager.cc.

References SiPixelDetectorStatus::getLSRange().

Referenced by createPayloads().

                                                                                                  {
  return (status1.getLSRange().first < status2.getLSRange().first);
}

void SiPixelStatusManager::readLumi

(

const edm::LuminosityBlock &

iLumi

)

Definition at line 72 of file SiPixelStatusManager.cc.

References SiPixelDetectorStatus::digiOccDET(), edm::LuminosityBlock::getByToken(), edm::HandleBase::isValid(), edm::LuminosityBlockBase::luminosityBlock(), edm::LuminosityBlockBase::run(), siPixelStatusToken_, and siPixelStatusVtr_.

Referenced by SiPixelStatusHarvester::endLuminosityBlock().

                                                                {
  edm::Handle<SiPixelDetectorStatus> siPixelStatusHandle;
  iLumi.getByToken(siPixelStatusToken_, siPixelStatusHandle);

  if (siPixelStatusHandle.isValid()) {  // check the product
    SiPixelDetectorStatus tmpStatus = (*siPixelStatusHandle);
    if (tmpStatus.digiOccDET() > 0) {  // only put in SiPixelDetectorStatus with non zero digi (pixel hit)
      siPixelStatusVtr_.push_back(tmpStatus);
    }
  } else {
    edm::LogWarning("SiPixelStatusManager") << " SiPixelDetectorStatus is not valid for run " << iLumi.run() << " lumi "
                                            << iLumi.luminosityBlock() << std::endl;
  }
}

void SiPixelStatusManager::reset

(

void

)

Definition at line 46 of file SiPixelStatusManager.cc.

References siPixelStatusMap_, and siPixelStatusVtr_.

Referenced by MatrixReader.MatrixReader::__init__(), data_sources.json_list::__next__(), MatrixReader.MatrixReader::showRaw(), SiPixelStatusHarvester::SiPixelStatusHarvester(), and SiPixelStatusManager().

                                 {
  siPixelStatusMap_.clear();
  siPixelStatusVtr_.clear();
}

Member Data Documentation

int SiPixelStatusManager::aveDigiOcc_

private

Definition at line 54 of file SiPixelStatusManager.h.

Referenced by createBadComponents().

std::map<edm::LuminosityBlockNumber_t, std::map<int, std::vector<int> > > SiPixelStatusManager::FEDerror25Map_

private

Definition at line 51 of file SiPixelStatusManager.h.

Referenced by createFEDerror25(), and getFEDerror25Rocs().

std::string SiPixelStatusManager::label_

private

Definition at line 57 of file SiPixelStatusManager.h.

Referenced by SequenceTypes.DummyModule::__repr__(), Mixins._Labelable::_findDependencies(), Mixins._Labelable::setLabel(), and SiPixelStatusManager().

std::string SiPixelStatusManager::moduleName_

private

Definition at line 56 of file SiPixelStatusManager.h.

Referenced by SiPixelStatusManager().

int SiPixelStatusManager::nLumi_

private

Definition at line 55 of file SiPixelStatusManager.h.

Referenced by createBadComponents().

std::string SiPixelStatusManager::outputBase_

private

Definition at line 53 of file SiPixelStatusManager.h.

Referenced by createBadComponents(), and SiPixelStatusManager().

std::map<edm::LuminosityBlockNumber_t, SiPixelDetectorStatus> SiPixelStatusManager::siPixelStatusMap_

private

Definition at line 50 of file SiPixelStatusManager.h.

Referenced by createBadComponents(), getBadComponents(), and reset().

edm::EDGetTokenT<SiPixelDetectorStatus> SiPixelStatusManager::siPixelStatusToken_

private

Definition at line 59 of file SiPixelStatusManager.h.

Referenced by readLumi(), and SiPixelStatusManager().

std::vector<SiPixelDetectorStatus> SiPixelStatusManager::siPixelStatusVtr_

private

Definition at line 49 of file SiPixelStatusManager.h.

Referenced by createBadComponents(), createFEDerror25(), createPayloads(), readLumi(), and reset().