/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/CondFormats/SiPixelObjects/src/SiPixelQuality.cc

Go to the documentation of this file.

// SiPixelQuality.cc
//
// class implementation to hold a list of disabled pixel modules
//
// M. Eads
// Apr 2008

#include "CondFormats/SiPixelObjects/interface/SiPixelQuality.h"
#include "CondFormats/SiPixelObjects/interface/SiPixelFrameReverter.h"
#include "CondFormats/SiPixelObjects/interface/SiPixelFrameConverter.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/ESWatcher.h"
#include "CondFormats/DataRecord/interface/SiPixelFedCablingMapRcd.h"
#include "CondFormats/SiPixelObjects/interface/SiPixelFedCablingMap.h"
#include "CondFormats/SiPixelObjects/interface/SiPixelFedCablingTree.h"
#include "CondFormats/SiPixelObjects/interface/PixelROC.h"
#include "CondFormats/SiPixelObjects/interface/LocalPixel.h"
#include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetUnit.h"
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"


#include <algorithm>
#include <iostream>


using namespace sipixelobjects;


      //  errortype "whole" = int 0 in DB //
      //  errortype "tbmA" = int 1 in DB  //
      //  errortype "tbmB" = int 2 in DB  //
      //  errortype "none" = int 3 in DB  //
    
      //each bad roc correspond to a bit to 1: num=  //
      // 0 <-> all good rocs                         //
      // 1 <-> only roc 0 bad                        //
      // 2<-> only roc 1 bad                         //
      // 3<->  roc 0 and 1 bad                       //
      // 4 <-> only roc 2 bad                        //
      //  ...                                        //


// add a list of modules to the vector of disabled modules
void SiPixelQuality::addDisabledModule(std::vector<SiPixelQuality::disabledModuleType> & idVector) {
  theDisabledModules.insert(theDisabledModules.end(),
                            idVector.begin(),
                            idVector.end());

} 

void SiPixelQuality::add(const SiStripDetVOff* Voff) {
  
  //Get vector of Voff dets
  std::vector<uint32_t> vdets;
  Voff->getDetIds(vdets);

  std::vector<uint32_t>::const_iterator iter=vdets.begin();
  std::vector<uint32_t>::const_iterator iterEnd=vdets.end();

  for(;iter!=iterEnd;++iter){
    SiPixelQuality::disabledModuleType BadModule;  
    BadModule.DetID = *iter;
    BadModule.errorType = 0;
    BadModule.BadRocs = 65535; //NOTE: here all module is bad - no difference for half modules
    if(IsModuleUsable(BadModule.DetID))
      addDisabledModule(BadModule);
    
  }
}

// bool SiPixelQuality::removeDisabledModule(const uint32_t & detid){
//   std::vector<disabledModuleType>::const_iterator iter = std::lower_bound(theDisabledModules.begin(),theDisabledModules.end(),detid,SiPixelQuality::BadComponentStrictWeakOrdering());
//   if (iter != theDisabledModules.end() && iter->DetID==detid){
//     theDisabledModules.erase(iter);
//     return true;
//   }
//   return false;
// }

// ask if module is OK
bool SiPixelQuality::IsModuleUsable(const uint32_t& detid) const {
   if(IsFedBad(detid))
     return true;
   std::vector<SiPixelQuality::disabledModuleType>disabledModules = theDisabledModules;
   std::sort(disabledModules.begin(),disabledModules.end(),SiPixelQuality::BadComponentStrictWeakOrdering());
   std::vector<disabledModuleType>::const_iterator iter = std::lower_bound(disabledModules.begin(),disabledModules.end(),detid,SiPixelQuality::BadComponentStrictWeakOrdering());
   if (iter != disabledModules.end() && iter->DetID==detid && iter->errorType ==0)
    return false;
  return true;
}

int SiPixelQuality::BadModuleNumber() {
  std::vector<SiPixelQuality::disabledModuleType>disabledModules = theDisabledModules;
  //does not distinguish between partially dead or not 
  return disabledModules.size();
}

//ask if module is bad
 bool SiPixelQuality::IsModuleBad(const uint32_t & detid) const {
   if(IsFedBad(detid))
     return true;
   std::vector<SiPixelQuality::disabledModuleType>disabledModules = theDisabledModules;
   std::sort(disabledModules.begin(),disabledModules.end(),SiPixelQuality::BadComponentStrictWeakOrdering());
   std::vector<disabledModuleType>::const_iterator iter = std::lower_bound(disabledModules.begin(),disabledModules.end(),detid,SiPixelQuality::BadComponentStrictWeakOrdering());
   if (iter != disabledModules.end() && iter->DetID==detid && iter->errorType == 0) //errorType 0 corresponds to "whole" dead module
      return true;
   return false;
}

//ask if roc is bad
bool SiPixelQuality::IsRocBad(const uint32_t& detid, const short& rocNb) const {
  if(IsModuleBad(detid))
    return true;
   std::vector<SiPixelQuality::disabledModuleType>disabledModules = theDisabledModules;
   std::sort(disabledModules.begin(),disabledModules.end(),SiPixelQuality::BadComponentStrictWeakOrdering());
   std::vector<disabledModuleType>::const_iterator iter = std::lower_bound(disabledModules.begin(),disabledModules.end(),detid,SiPixelQuality::BadComponentStrictWeakOrdering());
   if (iter != disabledModules.end() && iter->DetID == detid){
     return ((iter->BadRocs >> rocNb)&0x1);}
  return false;
}

bool SiPixelQuality::IsAreaBad(uint32_t detid, sipixelobjects::GlobalPixel global, const edm::EventSetup& es, const SiPixelFedCabling* map ) const {
 
  SiPixelFrameReverter reverter(es,map);
  int rocfromarea = -1;  
  rocfromarea = reverter.findRocInDet(detid, global);

  return SiPixelQuality::IsRocBad(detid, rocfromarea);
  return false;
}

//ask for bad ROCS


short SiPixelQuality::getBadRocs(const uint32_t& detid) const{
  std::vector<SiPixelQuality::disabledModuleType>disabledModules = theDisabledModules;
  std::sort(disabledModules.begin(),disabledModules.end(),SiPixelQuality::BadComponentStrictWeakOrdering());
  std::vector<disabledModuleType>::const_iterator iter = std::lower_bound(disabledModules.begin(),disabledModules.end(),detid,SiPixelQuality::BadComponentStrictWeakOrdering());
  if (iter != disabledModules.end() && iter->DetID==detid)
    return iter->BadRocs;
  return 0;
}

const std::vector<LocalPoint> SiPixelQuality::getBadRocPositions(const uint32_t & detid, const TrackerGeometry& theTracker, const SiPixelFedCabling* map) const{
  std::vector<LocalPoint> badrocpositions (0);
  std::pair<uint8_t, uint8_t> coord(1,1);
   for(unsigned int i = 0; i < 16; i++){
     if (IsRocBad(detid, i) == true){
    std::vector<CablingPathToDetUnit> path = map->pathToDetUnit(detid);
    typedef  std::vector<CablingPathToDetUnit>::const_iterator IT;
       for  (IT it = path.begin(); it != path.end(); ++it) {
          const PixelROC * myroc = map->findItem(*it);
          if( myroc->idInDetUnit() == i) {
              LocalPixel::RocRowCol  local = { 39, 25};   //corresponding to center of ROC row, col
              GlobalPixel global = myroc->toGlobal( LocalPixel(local) );
              //       edm::ESHandle<TrackerGeometry> geom;
              //     es.get<TrackerDigiGeometryRecord>().get( geom );
              //    const TrackerGeometry& theTracker(*geom);
              const PixelGeomDetUnit * theGeomDet = dynamic_cast<const PixelGeomDetUnit*> (theTracker.idToDet(detid) );

              PixelTopology const * topology = &(theGeomDet->specificTopology());

              MeasurementPoint thepoint(global.row, global.col);
              LocalPoint localpoint = topology->localPosition(thepoint);
              badrocpositions.push_back(localpoint);
         break;
          }
       }
     }
   }
  return badrocpositions;
}


bool SiPixelQuality::IsFedBad(const uint32_t & detid) const{
  return false;
}