#include <SiStripPedestals.h>

Classes
struct	DetRegistry

class	StrictWeakOrdering

Public Types
typedef std::vector< char >	Container

typedef std::vector< char >::const_iterator	ContainerIterator

typedef std::vector< uint16_t >	InputVector

typedef std::pair< ContainerIterator, ContainerIterator >	Range

typedef std::vector< DetRegistry >	Registry

typedef Registry::const_iterator	RegistryIterator

Public Member Functions
void	allPeds (std::vector< int > &pefs, const Range &range) const

ContainerIterator	getDataVectorBegin () const

ContainerIterator	getDataVectorEnd () const

void	getDetIds (std::vector< uint32_t > &DetIds_) const

float	getPed (const uint16_t &strip, const Range &range) const

const Range	getRange (const uint32_t &detID) const

RegistryIterator	getRegistryVectorBegin () const

RegistryIterator	getRegistryVectorEnd () const

void	printDebug (std::stringstream &ss, const TrackerTopology *trackerTopo) const
	Prints all pedestals. More...

void	printSummary (std::stringstream &ss, const TrackerTopology *trackerTopo) const
	Prints mean pedestal value divided for subdet, layer and mono/stereo. More...

bool	put (const uint32_t &detID, InputVector &input)

void	setData (float ped, InputVector &vped)

	SiStripPedestals ()

	~SiStripPedestals ()

Private Member Functions
uint16_t	decode (const uint16_t &strip, const Range &range) const

void	encode (InputVector &Vi, std::vector< unsigned char > &Vo_CHAR)

uint16_t	get10bits (const uint8_t *&ptr, int8_t skip) const

template<class Archive >
void	serialize (Archive &ar, const unsigned int version)

Private Attributes
Registry	indexes

Container	v_pedestals

Friends
class	boost::serialization::access

template<typename CondSerializationT , typename Enabled >
struct	cond::serialization::access

Detailed Description

Stores the pedestal of every strip.
Encodes the information in a vector<char> and uses a vector<DetRegistry> to connect each range of values to the corresponding detId.
The DetRegistry struct contains the detId and two uint32_t giving the index of begin and end of the corresponding range in the vector<char>.
Has methods to return the pedestal of a given strip of of all the strips.

The printSummary method uses SiStripDetSummary. See description therein.
The printDebug method prints the pedestal value for every strip of every detId.

Definition at line 27 of file SiStripPedestals.h.

Member Typedef Documentation

◆ Container

typedef std::vector<char> SiStripPedestals::Container

Definition at line 49 of file SiStripPedestals.h.

◆ ContainerIterator

typedef std::vector<char>::const_iterator SiStripPedestals::ContainerIterator

Definition at line 50 of file SiStripPedestals.h.

◆ InputVector

typedef std::vector<uint16_t> SiStripPedestals::InputVector

Definition at line 54 of file SiStripPedestals.h.

◆ Range

typedef std::pair<ContainerIterator, ContainerIterator> SiStripPedestals::Range

Definition at line 51 of file SiStripPedestals.h.

◆ Registry

typedef std::vector<DetRegistry> SiStripPedestals::Registry

Definition at line 52 of file SiStripPedestals.h.

◆ RegistryIterator

typedef Registry::const_iterator SiStripPedestals::RegistryIterator

Definition at line 53 of file SiStripPedestals.h.

Constructor & Destructor Documentation

◆ SiStripPedestals()

SiStripPedestals::SiStripPedestals ( )

inline

Definition at line 56 of file SiStripPedestals.h.

56 {};

◆ ~SiStripPedestals()

SiStripPedestals::~SiStripPedestals ( )

inline

Definition at line 57 of file SiStripPedestals.h.

57 {};

Member Function Documentation

◆ allPeds()

void SiStripPedestals::allPeds	(	std::vector< int > &	pefs,
		const Range &	range
	)		const

Definition at line 143 of file SiStripPedestals.cc.

                                                                              {
   size_t mysize = ((range.second - range.first) << 3) / 10;
   size_t size = peds.size();
   if (mysize < size)
     throw cms::Exception("CorruptedData") << "[SiStripPedestals::allPeds] Requested pedestals for " << peds.size()
                                           << " strips, I have it only for " << mysize << " strips\n";
   size_t size4 = size & (~0x3), carry = size & 0x3;  // we have an optimized way of unpacking 4 strips
   const uint8_t* ptr = reinterpret_cast<const uint8_t*>(&*range.second) - 1;
   std::vector<int>::iterator out = peds.begin(), end4 = peds.begin() + size4;
   // we do it this baroque way instead of just loopin on all the strips because it's faster
   // as the value of 'skip' is a constant, so the compiler can compute the masks directly
   while (out < end4) {
     *out = static_cast<int>(get10bits(ptr, 0));
     ++out;
     *out = static_cast<int>(get10bits(ptr, 2));
     ++out;
     *out = static_cast<int>(get10bits(ptr, 4));
     ++out;
     *out = static_cast<int>(get10bits(ptr, 6));
     ++out;
     --ptr;  // every 4 strips we have to skip one more bit
   }
   for (size_t rem = 0; rem < carry; ++rem) {
     *out = static_cast<int>(get10bits(ptr, 2 * rem));
     ++out;
   }
 }

References Exception, get10bits(), MillePedeFileConverter_cfg::out, FastTimerService_cff::range, and findQualityFiles::size.

Referenced by sistrip::FEDEmulator::retrievePedestals(), and SiStripPedestalsSubtractor::subtract_().

◆ decode()

uint16_t SiStripPedestals::decode	(	const uint16_t &	strip,
		const Range &	range
	)		const

private

Definition at line 101 of file SiStripPedestals.cc.

                                                                                  {
   const char* data = &*(range.second - 1);  // pointer to the last byte of data
   static const uint16_t BITS_PER_STRIP = 10;
  
   uint32_t lowBit = strip * BITS_PER_STRIP;
   uint8_t firstByteBit = (lowBit & 6);  //module
   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);
  
   /*
     if(strip  < 25){
     std::cout       << "***************DECODE*********************"<<"\n"
     << "strip "<<strip << " " 
     << value 
     <<"\t   :"<<print_as_binary(value) 
     <<"\t  :"<<print_as_binary(    ((uint16_t(*(data-lowBit/8  )) & firstByteMask) >>   firstByteBit)       )
     << "-"<<print_as_binary(  ((uint16_t(*(data-lowBit/8-1)) & secondByteMask) <<firstByteNBits)    )
     << "\t *(data-lowBit/8) " << print_as_binary(    *(data-lowBit/8 ))
     << "\t *(data-lowBit/8-1) " << print_as_binary(    *(data-lowBit/8 -1 ))
     << "\tlowBit:"<< lowBit
     << "\tfirstByteMask :"<<print_as_binary(firstByteMask)
     << "\tsecondByteMask:"<<print_as_binary(secondByteMask)
     << "\tfirstByteBit:"<<print_as_binary(firstByteBit)
     << std::endl;
     }
   */
   return value;
 }

References data, FastTimerService_cff::range, digitizers_cfi::strip, and relativeConstraints::value.

Referenced by ztail.Decoder::follow(), getPed(), and ztail.Decoder::initial_synchronize().

◆ encode()

void SiStripPedestals::encode	(	InputVector &	Vi,
		std::vector< unsigned char > &	Vo_CHAR
	)

private

Definition at line 59 of file SiStripPedestals.cc.

                                                                            {
   static const uint16_t BITS_PER_STRIP = 10;
   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 & 0x6);
     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);
  
     /*
       if(stripIndex   < 25 ){
       std::cout       << "***************ENCODE*********************"<<std::endl
       << "\tdata-lowBit/8     :"<<print_as_binary((*(data-lowBit/8)   & ~(firstByteMask)))
       << "-"<<print_as_binary(((Vi[stripIndex] & 0xffu) <<firstByteBit))
       << "\tdata-lowBit/8-1   :"<<print_as_binary((*(data-lowBit/8-1)   & ~(secondByteMask)))
       << "-"<<print_as_binary((((Vi[stripIndex]>> firstByteNBits) & secondByteMask)))
       << std::endl;
       std::cout       << "strip "<<stripIndex<<"\tvi: " << Vi[stripIndex] <<"\t"
       << print_short_as_binary(Vi[stripIndex])
       << "\tvo1:"<< print_char_as_binary(*(data-lowBit/8))
       << "\tvo2:"<< print_char_as_binary(*(data-lowBit/8-1))
       << "\tlowBit:"<< lowBit
       << "\tfirstByteMask :"<<print_as_binary(firstByteMask)
       << "\tsecondByteMask:"<<print_as_binary(secondByteMask)
       << "\tfirstByteBit:"<<print_as_binary(firstByteBit)
       << std::endl;
       }
     */
   }
 }

References data, and mps_fire::i.

Referenced by put().

◆ get10bits()

uint16_t SiStripPedestals::get10bits	(	const uint8_t *&	ptr,
		int8_t	skip
	)		const

inlineprivate

Get 9 bit words from a bit stream, starting from the right, skipping the first 'skip' bits (0 < skip < 7). Ptr must point to the rightmost byte that has some bits of this word, and is updated by this function

Definition at line 135 of file SiStripPedestals.cc.

                                                                                   {
   uint8_t maskThis = (0xFF << skip);
   uint8_t maskThat = ((4 << skip) - 1);
   uint16_t ret = (((*ptr) & maskThis) >> skip);
   --ptr;
   return ret | (((*ptr) & maskThat) << (8 - skip));
 }

References runTheMatrix::ret, and optionsL1T::skip.

Referenced by allPeds().

◆ getDataVectorBegin()

ContainerIterator SiStripPedestals::getDataVectorBegin ( ) const

inline

Definition at line 64 of file SiStripPedestals.h.

64 { return v_pedestals.begin(); }

References v_pedestals.

◆ getDataVectorEnd()

ContainerIterator SiStripPedestals::getDataVectorEnd ( ) const

inline

Definition at line 65 of file SiStripPedestals.h.

65 { return v_pedestals.end(); }

References v_pedestals.

◆ getDetIds()

void SiStripPedestals::getDetIds ( std::vector< uint32_t > & DetIds_ ) const

Definition at line 38 of file SiStripPedestals.cc.

                                                                    {
   // returns vector of DetIds in map
   SiStripPedestals::RegistryIterator begin = indexes.begin();
   SiStripPedestals::RegistryIterator end = indexes.end();
   for (SiStripPedestals::RegistryIterator p = begin; p != end; ++p) {
     DetIds_.push_back(p->detid);
   }
 }

References mps_fire::end, indexes, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by SiStripPedestalsDQM::getActiveDetIds(), printDebug(), and printSummary().

◆ getPed()

float SiStripPedestals::getPed	(	const uint16_t &	strip,
		const Range &	range
	)		const

Definition at line 51 of file SiStripPedestals.cc.

                                                                               {
   if (10 * strip >= (range.second - range.first) * 8) {
     throw cms::Exception("CorruptedData")
         << "[SiStripPedestals::getPed] looking for SiStripPedestals for a strip out of range: strip " << strip;
   }
   return static_cast<float>(decode(strip, range));
 }

References decode(), Exception, FastTimerService_cff::range, and digitizers_cfi::strip.

Referenced by CalibrationScanTask::CalibrationScanTask(), CalibrationTask::CalibrationTask(), SiStripAPVRestorer::createCMMapRealPed(), SiStripDigitizerAlgorithm::digitize(), SiStripPedestalsDQM::fillMEsForDet(), SiStripPedestalsDQM::fillMEsForLayer(), LaserAlignment::fillPedestalProfiles(), printDebug(), printSummary(), and DigiSimLinkAlgorithm::run().

◆ getRange()

const SiStripPedestals::Range SiStripPedestals::getRange ( const uint32_t & detID ) const

Definition at line 28 of file SiStripPedestals.cc.

                                                                                   {
   // get SiStripPedestals Range of DetId
  
   RegistryIterator p = std::lower_bound(indexes.begin(), indexes.end(), DetId, SiStripPedestals::StrictWeakOrdering());
   if (p == indexes.end() || p->detid != DetId)
     return SiStripPedestals::Range(v_pedestals.end(), v_pedestals.end());
   else
     return SiStripPedestals::Range(v_pedestals.begin() + p->ibegin, v_pedestals.begin() + p->iend);
 }

References indexes, pfDeepBoostedJetPreprocessParams_cfi::lower_bound, AlCaHLTBitMon_ParallelJobs::p, and v_pedestals.

Referenced by CalibrationScanTask::CalibrationScanTask(), CalibrationTask::CalibrationTask(), SiStripAPVRestorer::createCMMapRealPed(), SiStripDigitizerAlgorithm::digitize(), SiStripPedestalsDQM::fillMEsForDet(), SiStripPedestalsDQM::fillMEsForLayer(), LaserAlignment::fillPedestalProfiles(), printDebug(), printSummary(), sistrip::FEDEmulator::retrievePedestals(), DigiSimLinkAlgorithm::run(), and SiStripPedestalsSubtractor::subtract_().

◆ getRegistryVectorBegin()

RegistryIterator SiStripPedestals::getRegistryVectorBegin ( ) const

inline

Definition at line 66 of file SiStripPedestals.h.

66 { return indexes.begin(); }

References indexes.

◆ getRegistryVectorEnd()

RegistryIterator SiStripPedestals::getRegistryVectorEnd ( ) const

inline

Definition at line 67 of file SiStripPedestals.h.

67 { return indexes.end(); }

References indexes.

◆ printDebug()

void SiStripPedestals::printDebug	(	std::stringstream &	ss,
		const TrackerTopology *	trackerTopo
	)		const

Prints all pedestals.

Definition at line 185 of file SiStripPedestals.cc.

                                                                                      {
   std::vector<uint32_t> detid;
   getDetIds(detid);
  
   ss << "Number of detids = " << detid.size() << std::endl;
  
   for (size_t id = 0; id < detid.size(); ++id) {
     SiStripPedestals::Range range = getRange(detid[id]);
  
     int strip = 0;
     ss << "detid" << std::setw(15) << "strip" << std::setw(10) << "pedestal" << std::endl;
     int detId = 0;
     int oldDetId = 0;
     for (int it = 0; it < (range.second - range.first) * 8 / 10; ++it) {
       detId = detid[id];
       if (detId != oldDetId) {
         oldDetId = detId;
         ss << detid[id];
       } else
         ss << "   ";
       ss << std::setw(15) << strip++ << std::setw(10) << getPed(it, range) << std::endl;
     }
   }
 }

References getDetIds(), getPed(), getRange(), triggerObjects_cff::id, FastTimerService_cff::range, contentValuesCheck::ss, and digitizers_cfi::strip.

◆ printSummary()

void SiStripPedestals::printSummary	(	std::stringstream &	ss,
		const TrackerTopology *	trackerTopo
	)		const

Prints mean pedestal value divided for subdet, layer and mono/stereo.

Definition at line 171 of file SiStripPedestals.cc.

                                                                                                  {
   std::vector<uint32_t> detid;
   getDetIds(detid);
   SiStripDetSummary summary{trackerTopo};
   for (size_t id = 0; id < detid.size(); ++id) {
     SiStripPedestals::Range range = getRange(detid[id]);
     for (int it = 0; it < (range.second - range.first) * 8 / 10; ++it) {
       summary.add(detid[id], getPed(it, range));
     }
   }
   ss << "Summary of pedestals:" << std::endl;
   summary.print(ss);
 }

References getDetIds(), getPed(), getRange(), triggerObjects_cff::id, FastTimerService_cff::range, contentValuesCheck::ss, and edmLumisInFiles::summary.

◆ put()

bool SiStripPedestals::put	(	const uint32_t &	detID,
		InputVector &	input
	)

Definition at line 5 of file SiStripPedestals.cc.

                                                                     {
   // put in SiStripPedestals of DetId
   std::vector<unsigned char> Vo_CHAR;
   encode(input, Vo_CHAR);
  
   Registry::iterator p =
       std::lower_bound(indexes.begin(), indexes.end(), DetId, SiStripPedestals::StrictWeakOrdering());
   if (p != indexes.end() && p->detid == DetId)
     return false;
  
   //size_t sd= input.second-input.first;
   size_t sd = Vo_CHAR.end() - Vo_CHAR.begin();
   DetRegistry detregistry;
   detregistry.detid = DetId;
   detregistry.ibegin = v_pedestals.size();
   detregistry.iend = v_pedestals.size() + sd;
   indexes.insert(p, detregistry);
  
   //v_pedestals.insert(v_pedestals.end(),input.first,input.second);
   v_pedestals.insert(v_pedestals.end(), Vo_CHAR.begin(), Vo_CHAR.end());
   return true;
 }

References SiStripPedestals::DetRegistry::detid, encode(), SiStripPedestals::DetRegistry::ibegin, SiStripPedestals::DetRegistry::iend, indexes, input, pfDeepBoostedJetPreprocessParams_cfi::lower_bound, AlCaHLTBitMon_ParallelJobs::p, sd, and v_pedestals.

Referenced by SiStripPopConPedestalsHandlerFromDQM::dqmEndJob(), and SiStripCondObjBuilderFromDb::storePedestals().

◆ serialize()

template<class Archive >

void SiStripPedestals::serialize	(	Archive &	ar,
		const unsigned int	version
	)

private

◆ setData()

void SiStripPedestals::setData	(	float	ped,
		SiStripPedestals::InputVector &	vped
	)

Definition at line 47 of file SiStripPedestals.cc.

                                                                            {
   vped.push_back((static_cast<uint16_t>(ped) & 0x3FF));
 }

Referenced by SiStripPopConPedestalsHandlerFromDQM::dqmEndJob(), SiStripCondObjBuilderFromDb::setDefaultValuesCabling(), and SiStripCondObjBuilderFromDb::setValuesCabling().

Friends And Related Function Documentation

◆ boost::serialization::access

friend class boost::serialization::access

friend

Definition at line 86 of file SiStripPedestals.h.

◆ cond::serialization::access

template<typename CondSerializationT , typename Enabled >

friend struct cond::serialization::access

friend

Definition at line 86 of file SiStripPedestals.h.

Member Data Documentation

◆ indexes

Registry SiStripPedestals::indexes

private

Definition at line 84 of file SiStripPedestals.h.

Referenced by getDetIds(), getRange(), getRegistryVectorBegin(), getRegistryVectorEnd(), and put().

◆ v_pedestals

Container SiStripPedestals::v_pedestals

private

Definition at line 83 of file SiStripPedestals.h.

Referenced by getDataVectorBegin(), getDataVectorEnd(), getRange(), and put().

Classes

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Friends

Detailed Description

Member Typedef Documentation

◆ Container

◆ ContainerIterator

◆ InputVector

◆ Range

◆ Registry

◆ RegistryIterator

Constructor & Destructor Documentation

◆ SiStripPedestals()

◆ ~SiStripPedestals()

Member Function Documentation

◆ allPeds()

◆ decode()

◆ encode()

◆ get10bits()

◆ getDataVectorBegin()

◆ getDataVectorEnd()

◆ getDetIds()

◆ getPed()

◆ getRange()

◆ getRegistryVectorBegin()

◆ getRegistryVectorEnd()

◆ printDebug()

◆ printSummary()

◆ put()

◆ serialize()

◆ setData()

Friends And Related Function Documentation

◆ boost::serialization::access

◆ cond::serialization::access

Member Data Documentation

◆ indexes

◆ v_pedestals