/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/RecoTracker/MeasurementDet/plugins/TkMeasurementDetSet.cc

Go to the documentation of this file.

#include "TkMeasurementDetSet.h"
#include "TkStripMeasurementDet.h"
#include "DataFormats/SiStripDetId/interface/SiStripDetId.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"


void StMeasurementDetSet::init(std::vector<TkStripMeasurementDet> & stripDets) {
  // assume vector is full and ordered!
  int size = stripDets.size();
  
  empty_.resize(size,true);
  activeThisEvent_.resize(size,true);
  activeThisPeriod_.resize(size,true);
  id_.resize(size);
  subId_.resize(size);
  totalStrips_.resize(size);
  
  bad128Strip_.resize(size*6);
  hasAny128StripBad_.resize(size);
  badStripBlocks_.resize(size);

  if (isRegional()) {
    clusterI_.resize(2*size);
  }  else {
    detSet_.resize(size);
  }
  
  for (int i=0; i!=size; ++i) {
    auto & mdet =  stripDets[i]; 
    mdet.setIndex(i);
    //intialize the detId !
    id_[i] = mdet.specificGeomDet().geographicalId().rawId();
    subId_[i]=SiStripDetId(id_[i]).subdetId()-3;
    //initalize the total number of strips
    totalStrips_[i] =  mdet.specificGeomDet().specificTopology().nstrips();
  }
}


void StMeasurementDetSet::set128StripStatus(int i, bool good, int idx) { 
    int offset =  nbad128*i;
    if (idx == -1) {
      std::fill(bad128Strip_.begin()+offset, bad128Strip_.begin()+offset+6, !good);
      hasAny128StripBad_[i] = !good;
    } else {
      bad128Strip_[offset+idx] = !good;
      if (good == false) {
        hasAny128StripBad_[i] = false;
      } else { // this should not happen, as usually you turn on all fibers
        // and then turn off the bad ones, and not vice-versa,
        // so I don't care if it's not optimized
        hasAny128StripBad_[i] = true;
        for (int j = 0; i < (totalStrips_[j] >> 7); j++) {
          if (bad128Strip_[j+offset] == false) hasAny128StripBad_[i] = false; break;
        }
      }    
    } 
  }
  

void StMeasurementDetSet::initializeStripStatus(const SiStripQuality *quality, int qualityFlags, int qualityDebugFlags, const edm::ParameterSet& cutPset) {
  if (qualityFlags & BadStrips) {
     badStripCuts_[SiStripDetId::TIB-3] = BadStripCuts(cutPset.getParameter<edm::ParameterSet>("TIB"));
    badStripCuts_[SiStripDetId::TOB-3] = BadStripCuts(cutPset.getParameter<edm::ParameterSet>("TOB"));
    badStripCuts_[SiStripDetId::TID-3] = BadStripCuts(cutPset.getParameter<edm::ParameterSet>("TID"));
    badStripCuts_[SiStripDetId::TEC-3] = BadStripCuts(cutPset.getParameter<edm::ParameterSet>("TEC"));
  }
  setMaskBad128StripBlocks((qualityFlags & MaskBad128StripBlocks) != 0);
  
  
  if ((quality != 0) && (qualityFlags != 0))  {
    edm::LogInfo("MeasurementTracker") << "qualityFlags = " << qualityFlags;
    unsigned int on = 0, tot = 0; 
    unsigned int foff = 0, ftot = 0, aoff = 0, atot = 0; 
    for (int i=0; i!=nDet(); i++) {
      uint32_t detid = id(i);
      if (qualityFlags & BadModules) {
        bool isOn = quality->IsModuleUsable(detid);
        setActive(i,isOn);
        tot++; on += (unsigned int) isOn;
        if (qualityDebugFlags & BadModules) {
          edm::LogInfo("MeasurementTracker")<< "MeasurementTrackerImpl::initializeStripStatus : detid " << detid << " is " << (isOn ?  "on" : "off");
        }
      } else {
        setActive(i,true);
      }
      // first turn all APVs and fibers ON
      set128StripStatus(i,true); 
      if (qualityFlags & BadAPVFibers) {
        short badApvs   = quality->getBadApvs(detid);
        short badFibers = quality->getBadFibers(detid);
        for (int j = 0; j < 6; j++) {
          atot++;
          if (badApvs & (1 << j)) {
            set128StripStatus(i,false, j);
            aoff++;
          }
        }
        for (int j = 0; j < 3; j++) {
          ftot++;
             if (badFibers & (1 << j)) {
               set128StripStatus(i,false, 2*j);
               set128StripStatus(i,false, 2*j+1);
               foff++;
             }
        }
      } 
      auto & badStrips = getBadStripBlocks(i);
       badStrips.clear();
       if (qualityFlags & BadStrips) {
         SiStripBadStrip::Range range = quality->getRange(detid);
         for (SiStripBadStrip::ContainerIterator bit = range.first; bit != range.second; ++bit) {
           badStrips.push_back(quality->decode(*bit));
         }
       }
    }
    if (qualityDebugFlags & BadModules) {
      edm::LogInfo("MeasurementTracker StripModuleStatus") << 
        " Total modules: " << tot << ", active " << on <<", inactive " << (tot - on);
    }
    if (qualityDebugFlags & BadAPVFibers) {
      edm::LogInfo("MeasurementTracker StripAPVStatus") << 
        " Total APVs: " << atot << ", active " << (atot-aoff) <<", inactive " << (aoff);
        edm::LogInfo("MeasurementTracker StripFiberStatus") << 
          " Total Fibers: " << ftot << ", active " << (ftot-foff) <<", inactive " << (foff);
    }
  } else {
    for (int i=0; i!=nDet(); i++) {
      setActive(i,true);          // module ON
      set128StripStatus(i,true);  // all APVs and fibers ON
    }
  }
}