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, const float &cTh, Container &vthr)
void	setData (const uint16_t &strip, const float &lTh, const float &hTh, 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

Container

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

Definition at line 104 of file SiStripThreshold.h.

ContainerIterator

typedef Container::const_iterator SiStripThreshold::ContainerIterator

Definition at line 105 of file SiStripThreshold.h.

InputVector

typedef Container SiStripThreshold::InputVector

Definition at line 109 of file SiStripThreshold.h.

Range

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

Definition at line 106 of file SiStripThreshold.h.

Registry

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

Definition at line 107 of file SiStripThreshold.h.

RegistryIterator

typedef Registry::const_iterator SiStripThreshold::RegistryIterator

Definition at line 108 of file SiStripThreshold.h.

Constructor & Destructor Documentation

SiStripThreshold() [1/2]

SiStripThreshold::SiStripThreshold ( )

inline

Definition at line 111 of file SiStripThreshold.h.

{};

SiStripThreshold() [2/2]

SiStripThreshold::SiStripThreshold ( const SiStripThreshold &  orig )

inline

Definition at line 112 of file SiStripThreshold.h.

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

References indexes, and v_threshold.

~SiStripThreshold()

virtual SiStripThreshold::~SiStripThreshold ( )

inlinevirtual

Definition at line 116 of file SiStripThreshold.h.

{};

Member Function Documentation

addToStat()

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.

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

References SiStripPI::max, min(), FastTimerService_cff::range, combinedConstraintHelpers::sum2(), and relativeConstraints::value.

Referenced by printSummary().

allThresholds()

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

Definition at line 82 of file SiStripThreshold.cc.

                                                               {
  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);
  }
}

References cms::cuda::assert(), ntuplemaker::fill, LaserClient_cfi::high, LaserClient_cfi::low, FastTimerService_cff::range, and DigiDM_cff::strips.

compact()

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

private

Definition at line 33 of file SiStripThreshold.cc.

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

References input, and tier0::unique().

Referenced by put().

getData()

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

Definition at line 70 of file SiStripThreshold.cc.

                                                                                              {
  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();
  }
}

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

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

getDataVectorBegin()

ContainerIterator SiStripThreshold::getDataVectorBegin ( ) const

inline

Definition at line 122 of file SiStripThreshold.h.

{ return v_threshold.begin(); }

References v_threshold.

Referenced by printDebug(), and printSummary().

getDataVectorEnd()

ContainerIterator SiStripThreshold::getDataVectorEnd ( ) const

inline

Definition at line 123 of file SiStripThreshold.h.

{ return v_threshold.end(); }

References v_threshold.

getDetIds()

void SiStripThreshold::getDetIds

(

std::vector< uint32_t > &

DetIds_

)

const

Definition at line 48 of file SiStripThreshold.cc.

                                                                   {
  // 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);
  }
}

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

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

getRange()

const SiStripThreshold::Range SiStripThreshold::getRange

(

const uint32_t &

detID

)

const

Definition at line 38 of file SiStripThreshold.cc.

                                                                                  {
  // 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);
}

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

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

getRegistryVectorBegin()

RegistryIterator SiStripThreshold::getRegistryVectorBegin ( ) const

inline

Definition at line 124 of file SiStripThreshold.h.

{ return indexes.begin(); }

References indexes.

Referenced by printDebug(), and printSummary().

getRegistryVectorEnd()

RegistryIterator SiStripThreshold::getRegistryVectorEnd ( ) const

inline

Definition at line 125 of file SiStripThreshold.h.

{ return indexes.end(); }

References indexes.

Referenced by printDebug(), and printSummary().

printDebug()

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.

                                                                                     {
  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);
    }
  }
}

References getDataVectorBegin(), getRegistryVectorBegin(), getRegistryVectorEnd(), and contentValuesCheck::ss.

printSummary()

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.

                                                                                       {
  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;
  }
}

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

put()

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

Definition at line 8 of file SiStripThreshold.cc.

                                                                          {
  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. Skippig 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;
}

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 SiStripCondObjBuilderFromDb::storeThreshold().

serialize()

template<class Archive >

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

private

setData() [1/2]

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.

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

References a, and digitizers_cfi::strip.

setData() [2/2]

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

Definition at line 57 of file SiStripThreshold.cc.

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

References a, and digitizers_cfi::strip.

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

Friends And Related Function Documentation

boost::serialization::access

friend class boost::serialization::access

friend

Definition at line 146 of file SiStripThreshold.h.

cond::serialization::access

template<typename CondSerializationT , typename Enabled >

friend struct cond::serialization::access

friend

Definition at line 146 of file SiStripThreshold.h.

Member Data Documentation

indexes

Registry SiStripThreshold::indexes

private

Definition at line 144 of file SiStripThreshold.h.

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

v_threshold

Container SiStripThreshold::v_threshold

private

Definition at line 143 of file SiStripThreshold.h.

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