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 40 of file SiPixelStatusManager.h.

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

Definition at line 39 of file SiPixelStatusManager.h.

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

Definition at line 41 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 25 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 93 of file SiPixelStatusManager.cc.

References aveDigiOcc_, 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();

  if(outputBase_ == "nLumibased" && nLumi_>1){ // doesn't work for nLumi_=1 cos any integer can be completely divided by 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 241 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(), siPixelStatusVtr_, and jetUpdater_cfi::sort.

Referenced by SiPixelStatusHarvester::endRunProduce().

                                         {
  if(!siPixelStatusVtr_.empty()){ //only create std::map payloads when the number of non-zero DIGI lumi sections is greater than ZERO otherwise segmentation fault

   // 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 36 of file SiPixelStatusManager.h.

References siPixelStatusMap_.

Referenced by SiPixelStatusHarvester::endRunProduce().

{return siPixelStatusMap_; } 

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

inline

Definition at line 37 of file SiPixelStatusManager.h.

References FEDerror25Map_.

Referenced by SiPixelStatusHarvester::endRunProduce().

{return FEDerror25Map_;}

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

staticprivate

Definition at line 53 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 74 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 47 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().