/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/CondFormats/Calibration/src/mySiStripNoises.cc

Go to the documentation of this file.

#include "CondFormats/Calibration/interface/mySiStripNoises.h"
#include <algorithm>
bool mySiStripNoises::put(const uint32_t DetId, InputVector& input) {
  // put in SiStripNoises of DetId
  std::vector<unsigned char>    Vo_CHAR;
  encode(input, Vo_CHAR);
  Registry::iterator p=std::lower_bound(indexes.begin(),indexes.end(),DetId,mySiStripNoises::StrictWeakOrdering());
  if (p!=indexes.end()  &&      p->detid==DetId)
    return false;
  size_t sd = Vo_CHAR.end() - Vo_CHAR.begin();
  DetRegistry detregistry;
  detregistry.detid=DetId;
  detregistry.ibegin=v_noises.size();
  detregistry.iend=v_noises.size()+sd;
  indexes.insert(p,detregistry);
  v_noises.insert(v_noises.end(),Vo_CHAR.begin(),Vo_CHAR.end());
  return true;
}

const mySiStripNoises::Range mySiStripNoises::getRange(const uint32_t& DetId) const {
  // get SiStripNoises Range of DetId
  RegistryIterator p = std::lower_bound(indexes.begin(),indexes.end(),DetId,StrictWeakOrdering());
  if (p==indexes.end()|| p->detid!=DetId) 
    return mySiStripNoises::Range(v_noises.end(),v_noises.end()); 
  else 
    return mySiStripNoises::Range(v_noises.begin()+p->ibegin,v_noises.begin()+p->iend);
}

void mySiStripNoises::getDetIds(std::vector<uint32_t>& DetIds) const {
  // returns vector of DetIds in map
  mySiStripNoises::RegistryIterator begin = indexes.begin();
  mySiStripNoises::RegistryIterator end   = indexes.end();
  for (mySiStripNoises::RegistryIterator p=begin; p != end; ++p) {
    DetIds.push_back(p->detid);
  }
}

float mySiStripNoises::getNoise( const uint16_t& strip, const Range& range ) const{
  return   static_cast<float> (decode(strip,range)/10.0);
} 

void mySiStripNoises::setData(float noise_, std::vector<short>& v){
  v.push_back((static_cast<int16_t>  (noise_*10.0 + 0.5) & 0x01FF)) ;
}

void mySiStripNoises::encode(InputVector& Vi, std::vector<unsigned char>& Vo){
  static const uint16_t  BITS_PER_STRIP  = 9;
  const size_t           VoSize          = (size_t)((Vi.size() *       BITS_PER_STRIP)/8+.999);
  Vo.resize(VoSize);
  for(size_t i = 0; i<Vo.size(); ++i)
    Vo[i]   &=      0x00u;
  
  for(unsigned int stripIndex =0; stripIndex<Vi.size(); ++stripIndex){
    unsigned char*  data    =       &Vo[Vo.size()-1];
    uint32_t lowBit         =       stripIndex * BITS_PER_STRIP;
    uint8_t firstByteBit    =       (lowBit & 0x7);
    uint8_t firstByteNBits  =       8 - firstByteBit;
    uint8_t firstByteMask   =       0xffu << firstByteBit;
    uint8_t secondByteNbits =       (BITS_PER_STRIP - firstByteNBits);
    uint8_t secondByteMask  =       ~(0xffu << secondByteNbits);

    *(data-lowBit/8)        =       (*(data-lowBit/8)   & ~(firstByteMask))         | ((Vi[stripIndex] & 0xffu) <<firstByteBit);
    *(data-lowBit/8-1)      =       (*(data-lowBit/8-1) & ~(secondByteMask))        | ((Vi[stripIndex] >> firstByteNBits) & secondByteMask);
  }
}

uint16_t mySiStripNoises::decode (const uint16_t& strip, const Range& range) const{
  const unsigned char *data = &*(range.second -1);  // pointer to the last byte of data
  static const uint16_t BITS_PER_STRIP = 9;
  
  uint32_t lowBit        = strip * BITS_PER_STRIP;
  uint8_t firstByteBit   = (lowBit & 7);//module 8
  uint8_t firstByteNBits = 8 - firstByteBit;
  uint8_t firstByteMask  = 0xffu << firstByteBit;
  uint8_t secondByteMask = ~(0xffu << (BITS_PER_STRIP - firstByteNBits));
  uint16_t value         =   ((uint16_t(*(data-lowBit/8  )) & firstByteMask) >> firstByteBit) | ((uint16_t(*(data-lowBit/8-1)) & secondByteMask) << firstByteNBits);
  return value;
}