SiStripThreshold Class Reference

#include <SiStripThreshold.h>

Classes
struct	Data
class	dataStrictWeakOrdering
struct	DetRegistry
class	StrictWeakOrdering

Public Types
typedef std::vector< Data >	Container
typedef Container::const_iterator	ContainerIterator
typedef Container	InputVector
typedef std::pair < ContainerIterator, ContainerIterator >	Range
typedef std::vector< DetRegistry >	Registry
typedef Registry::const_iterator	RegistryIterator

Public Member Functions
void	allThresholds (std::vector< float > &lowThs, std::vector< float > &highThs, const Range &range) const
SiStripThreshold::Data	getData (const uint16_t &strip, const Range &range) const
ContainerIterator	getDataVectorBegin () const
ContainerIterator	getDataVectorEnd () const
void	getDetIds (std::vector< uint32_t > &DetIds_) 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 the thresholds for all DetIds. More...
void	printSummary (std::stringstream &ss, const TrackerTopology *trackerTopo) const
	Prints mean, rms, min and max threshold values for each DetId. More...
bool	put (const uint32_t &detID, const InputVector &vect)
void	setData (const uint16_t &strip, const float &lTh, const float &hTh, Container &vthr)
void	setData (const uint16_t &strip, const float &lTh, const float &hTh, const float &cTh, Container &vthr)
	SiStripThreshold ()
	SiStripThreshold (const SiStripThreshold &orig)
virtual	~SiStripThreshold ()

Private Member Functions
void	addToStat (float value, uint16_t &range, float &sum, float &sum2, float &min, float &max) const
Container::iterator	compact (Container &input)
template<class Archive >
void	serialize (Archive &ar, const unsigned int version)

Private Attributes
Registry	indexes
Container	v_threshold

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

Detailed Description

Holds the thresholds:

• High threshold
  
• Low threshold
  
• Cluster threshold
  The values are stored as bitsets, in particular:
• One uint16_t stores the first strip and the high threshold.
• One uint8_t stores the low threshold.
• One uint8_t stores the cluster threshold. The information are stored in a Data struct, which is also responsible for the encoding and decoding.
  To fill the SiStripThreshold object:
• create and empty vector<Data>: SiStripThreshold::Container theSiStripVector;
• use the setData method to fill it with the thresholds
• use the put method to save the values in the object. The put method is used to actually fill the object. It receives the DetId and the vector<Data> with the threshold values.
  Before being saved the vector<Data> is sorted in the FirstStrip and consecutive entries having all the same thresholds are removed. This way it still stores the same information using less space.
  To retrieve the information:
• getDetIds: fills a vector with all detIds for which the thresholds have been stored.
• getDataVectorBegin and getDataVectorEnd: return the begin and end iterators to the vector<Data>.
• getData can be used to get the thresholds for a single strip. The printSummary method prints mean, rms, min and max threshold values for each DetId. The printDebug method prints all the thresholds for all DetIds.

Definition at line 42 of file SiStripThreshold.h.

Member Typedef Documentation

typedef std::vector<Data> SiStripThreshold::Container

Definition at line 104 of file SiStripThreshold.h.

typedef Container::const_iterator SiStripThreshold::ContainerIterator

Definition at line 105 of file SiStripThreshold.h.

typedef Container SiStripThreshold::InputVector

Definition at line 109 of file SiStripThreshold.h.

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

Definition at line 106 of file SiStripThreshold.h.

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

Definition at line 107 of file SiStripThreshold.h.

typedef Registry::const_iterator SiStripThreshold::RegistryIterator

Definition at line 108 of file SiStripThreshold.h.

Constructor & Destructor Documentation

SiStripThreshold::SiStripThreshold ( )

inline

Definition at line 111 of file SiStripThreshold.h.

{};

SiStripThreshold::SiStripThreshold ( const SiStripThreshold &  orig )

inline

Definition at line 112 of file SiStripThreshold.h.

References indexes, and v_threshold.

                                                 {
    v_threshold = orig.v_threshold;
    indexes = orig.indexes;
  }

virtual SiStripThreshold::~SiStripThreshold ( )

inlinevirtual

Definition at line 116 of file SiStripThreshold.h.

{};

Member Function Documentation

void SiStripThreshold::addToStat ( float  value, uint16_t &  range, float &  sum, float &  sum2, float &  min, float &  max  ) const

private

Definition at line 177 of file SiStripThreshold.cc.

References sistrip::SpyUtilities::range(), and relativeConstraints::value.

Referenced by printSummary().

                                                                                                                    {
  sum += value * range;
  sum2 += value * value * range;
  if (value < min)
    min = value;
  if (value > max)
    max = value;
}

void SiStripThreshold::allThresholds	(	std::vector< float > &	lowThs,
		std::vector< float > &	highThs,
		const Range &	range
	)		const

Definition at line 82 of file SiStripThreshold.cc.

References cms::cuda::assert(), PVValHelper::fill(), and DigiDM_cff::strips.

                                                               {
  ContainerIterator it = range.first;
  size_t strips = lowThs.size();
  assert(strips == highThs.size());
  while (it != range.second) {
    size_t firstStrip = it->getFirstStrip();
    //std::cout << "First strip is " << firstStrip << std::endl;
    float high = it->getHth(), low = it->getLth();
    //std::cout << "High is " << high << ", low is " << low << std::endl;
    ++it;  // increment the pointer
    size_t lastStrip = (it == range.second ? strips : it->getFirstStrip());
    //std::cout << "Last strip is " << lastStrip << std::endl;
    if (lastStrip > strips) {
      it = range.second;   // I should stop here,
      lastStrip = strips;  // and fill only 'strips' strips
    }
    std::fill(&lowThs[firstStrip], &lowThs[lastStrip], low);
    std::fill(&highThs[firstStrip], &highThs[lastStrip], high);
  }
}

SiStripThreshold::Container::iterator SiStripThreshold::compact ( Container &  input )

private

Definition at line 33 of file SiStripThreshold.cc.

References tier0::unique().

Referenced by put().

                                                                            {
  std::stable_sort(input.begin(), input.end());
  return std::unique(input.begin(), input.end());
}

SiStripThreshold::Data SiStripThreshold::getData	(	const uint16_t &	strip,
		const Range &	range
	)		const

Definition at line 70 of file SiStripThreshold.cc.

References Exception, sistrip::FirstThStripMask_, sistrip::FirstThStripShift_, sistrip::HighThStripMask_, AlCaHLTBitMon_ParallelJobs::p, and cuda_std::upper_bound().

Referenced by SiStripThresholdDQM::fillMEsForDet(), SiStripThresholdDQM::fillMEsForLayer(), and SiStripFedZeroSuppression::suppress().

                                                                                              {
  uint16_t estrip =
      (strip & sistrip::FirstThStripMask_) << sistrip::FirstThStripShift_ | (63 & sistrip::HighThStripMask_);
  ContainerIterator p = std::upper_bound(range.first, range.second, estrip, SiStripThreshold::dataStrictWeakOrdering());
  if (p != range.first) {
    return *(--p);
  } else {
    throw cms::Exception("CorruptedData") << "[SiStripThreshold::getData] asking for data for a strip " << strip
                                          << " lower then the first stored strip " << p->getFirstStrip();
  }
}

ContainerIterator SiStripThreshold::getDataVectorBegin ( ) const

inline

Definition at line 122 of file SiStripThreshold.h.

References v_threshold.

Referenced by printDebug(), and printSummary().

{ return v_threshold.begin(); }

ContainerIterator SiStripThreshold::getDataVectorEnd ( ) const

inline

Definition at line 123 of file SiStripThreshold.h.

References v_threshold.

{ return v_threshold.end(); }

void SiStripThreshold::getDetIds

(

std::vector< uint32_t > &

DetIds_

)

const

Definition at line 48 of file SiStripThreshold.cc.

References SplitLinear::begin, dataset::end, indexes, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by SiStripThresholdDQM::fillMEsForDet(), and SiStripThresholdDQM::getActiveDetIds().

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

const SiStripThreshold::Range SiStripThreshold::getRange

(

const uint32_t &

detID

)

const

Definition at line 38 of file SiStripThreshold.cc.

References indexes, cuda_std::lower_bound(), AlCaHLTBitMon_ParallelJobs::p, and v_threshold.

Referenced by SiStripThresholdDQM::fillMEsForDet(), SiStripThresholdDQM::fillMEsForLayer(), and SiStripFedZeroSuppression::suppress().

                                                                                  {
  // get SiStripThreshold Range of DetId

  RegistryIterator p = std::lower_bound(indexes.begin(), indexes.end(), DetId, SiStripThreshold::StrictWeakOrdering());
  if (p == indexes.end() || p->detid != DetId)
    return SiStripThreshold::Range(v_threshold.end(), v_threshold.end());
  else
    return SiStripThreshold::Range(v_threshold.begin() + p->ibegin, v_threshold.begin() + p->iend);
}

RegistryIterator SiStripThreshold::getRegistryVectorBegin ( ) const

inline

Definition at line 124 of file SiStripThreshold.h.

References indexes.

Referenced by printDebug(), and printSummary().

{ return indexes.begin(); }

RegistryIterator SiStripThreshold::getRegistryVectorEnd ( ) const

inline

Definition at line 125 of file SiStripThreshold.h.

References indexes.

Referenced by printDebug(), and printSummary().

{ return indexes.end(); }

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

Prints all the thresholds for all DetIds.

Definition at line 105 of file SiStripThreshold.cc.

References getDataVectorBegin(), getRegistryVectorBegin(), and getRegistryVectorEnd().

                                                                                     {
  RegistryIterator rit = getRegistryVectorBegin(), erit = getRegistryVectorEnd();
  ContainerIterator it, eit;
  for (; rit != erit; ++rit) {
    it = getDataVectorBegin() + rit->ibegin;
    eit = getDataVectorBegin() + rit->iend;
    ss << "\ndetid: " << rit->detid << " \t ";
    for (; it != eit; ++it) {
      ss << "\n \t ";
      it->print(ss);
    }
  }
}

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

Prints mean, rms, min and max threshold values for each DetId.

Definition at line 119 of file SiStripThreshold.cc.

References addToStat(), getDataVectorBegin(), getRegistryVectorBegin(), getRegistryVectorEnd(), dqmiodumpmetadata::n, and mathSSE::sqrt().

                                                                                       {
  RegistryIterator rit = getRegistryVectorBegin(), erit = getRegistryVectorEnd();
  ContainerIterator it, eit, itp;
  float meanLth, meanHth, meanCth;  //mean value
  float rmsLth, rmsHth, rmsCth;     //rms value
  float maxLth, maxHth, maxCth;     //max value
  float minLth, minHth, minCth;     //min value
  uint16_t n;
  uint16_t firstStrip, stripRange;
  for (; rit != erit; ++rit) {
    it = getDataVectorBegin() + rit->ibegin;
    eit = getDataVectorBegin() + rit->iend;
    ss << "\ndetid: " << rit->detid << " \t ";

    meanLth = 0;
    meanHth = 0;
    meanCth = 0;  //mean value
    rmsLth = 0;
    rmsHth = 0;
    rmsCth = 0;  //rms value
    maxLth = 0;
    maxHth = 0;
    maxCth = 0;  //max value
    minLth = 10000;
    minHth = 10000;
    minCth = 10000;  //min value
    n = 0;
    firstStrip = 0;
    for (; it != eit; ++it) {
      itp = it + 1;
      firstStrip = it->getFirstStrip();
      if (itp != eit)
        stripRange = (itp->getFirstStrip() - firstStrip);
      else
        stripRange =
            firstStrip > 511
                ? 768 - firstStrip
                : 512 -
                      firstStrip;  //*FIXME, I dont' know ithis class the strip number of a detector, so I assume wrongly that if the last firstStrip<511 the detector has only 512 strips. Clearly wrong. to be fixed

      addToStat(it->getLth(), stripRange, meanLth, rmsLth, minLth, maxLth);
      addToStat(it->getHth(), stripRange, meanHth, rmsHth, minHth, maxHth);
      addToStat(it->getClusth(), stripRange, meanCth, rmsCth, minCth, maxCth);
      n += stripRange;
    }
    meanLth /= n;
    meanHth /= n;
    meanCth /= n;
    rmsLth = sqrt(rmsLth / n - meanLth * meanLth);
    rmsHth = sqrt(rmsHth / n - meanHth * meanHth);
    rmsCth = sqrt(rmsCth / n - meanCth * meanCth);
    ss << "\nn " << n << " \tmeanLth " << meanLth << " \t rmsLth " << rmsLth << " \t minLth " << minLth << " \t maxLth "
       << maxLth;
    ss << "\n\tmeanHth " << meanHth << " \t rmsHth " << rmsHth << " \t minHth " << minHth << " \t maxHth " << maxHth;
    ss << "\n\tmeanCth " << meanCth << " \t rmsCth " << rmsCth << " \t minCth " << minCth << " \t maxCth " << maxCth;
  }
}

bool SiStripThreshold::put	(	const uint32_t &	detID,
		const InputVector &	vect
	)

Definition at line 8 of file SiStripThreshold.cc.

References compact(), SiStripThreshold::DetRegistry::detid, SiStripThreshold::DetRegistry::ibegin, SiStripThreshold::DetRegistry::iend, indexes, cuda_std::lower_bound(), AlCaHLTBitMon_ParallelJobs::p, sd, and v_threshold.

Referenced by SiStripThresholdBuilder::analyze(), and SiStripCondObjBuilderFromDb::storeThreshold().

                                                                          {
  InputVector vect = _vect;
  // put in SiStripThreshold::v_threshold of DetId
  Registry::iterator p =
      std::lower_bound(indexes.begin(), indexes.end(), DetId, SiStripThreshold::StrictWeakOrdering());
  if (p != indexes.end() && p->detid == DetId) {
    edm::LogError("SiStripThreshold") << "[" << __PRETTY_FUNCTION__ << "] SiStripThreshold for DetID " << DetId
                                      << " is already stored. Skipping this put" << std::endl;
    return false;
  }

  SiStripThreshold::Container::iterator new_end = compact(vect);

  size_t sd = new_end - vect.begin();
  DetRegistry detregistry;
  detregistry.detid = DetId;
  detregistry.ibegin = v_threshold.size();
  detregistry.iend = v_threshold.size() + sd;
  indexes.insert(p, detregistry);

  v_threshold.insert(v_threshold.end(), vect.begin(), new_end);

  return true;
}

template<class Archive >

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

private

void SiStripThreshold::setData	(	const uint16_t &	strip,
		const float &	lTh,
		const float &	hTh,
		Container &	vthr
	)

Definition at line 57 of file SiStripThreshold.cc.

References a, and SiStripThreshold::Data::encode().

Referenced by SiStripThresholdBuilder::analyze(), SiStripCondObjBuilderFromDb::setDefaultValuesCabling(), and SiStripCondObjBuilderFromDb::setValuesCabling().

                                                                                                         {
  Data a;
  a.encode(strip, lTh, hTh);
  vthr.push_back(a);
}

void SiStripThreshold::setData	(	const uint16_t &	strip,
		const float &	lTh,
		const float &	hTh,
		const float &	cTh,
		Container &	vthr
	)

Definition at line 63 of file SiStripThreshold.cc.

References a, and SiStripThreshold::Data::encode().

                                                                                                  {
  Data a;
  a.encode(strip, lTh, hTh, cTh);
  vthr.push_back(a);
}

Friends And Related Function Documentation

friend class boost::serialization::access

friend

Definition at line 146 of file SiStripThreshold.h.

template<typename CondSerializationT , typename Enabled >

friend struct cond::serialization::access

friend

Definition at line 146 of file SiStripThreshold.h.

Member Data Documentation

Registry SiStripThreshold::indexes

private

Definition at line 144 of file SiStripThreshold.h.

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

Container SiStripThreshold::v_threshold

private

Definition at line 143 of file SiStripThreshold.h.

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