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

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.

References indexes, and v_threshold.

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

~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 173 of file SiStripThreshold.cc.

References WZElectronSkims53X_cff::max, SiStripPI::min, isotrackApplyRegressor::range, and relativeConstraints::value.

Referenced by printSummary().

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

allThresholds()

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(), EcalPhiSymFlatTableProducers_cfi::fill, LaserClient_cfi::high, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, LaserClient_cfi::low, isotrackApplyRegressor::range, 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);
  }
}

compact()

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

private

Definition at line 33 of file SiStripThreshold.cc.

References input, and tier0::unique().

Referenced by put().

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

getData()

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, isotrackApplyRegressor::range, nano_mu_digi_cff::strip, and pfDeepBoostedJetPreprocessParams_cfi::upper_bound.

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

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

getDataVectorBegin()

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

getDataVectorEnd()

ContainerIterator SiStripThreshold::getDataVectorEnd ( ) const

inline

Definition at line 123 of file SiStripThreshold.h.

References v_threshold.

{ return v_threshold.end(); }

getDetIds()

void SiStripThreshold::getDetIds

(

std::vector< uint32_t > &

DetIds_

)

const

Definition at line 48 of file SiStripThreshold.cc.

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

getRange()

const SiStripThreshold::Range SiStripThreshold::getRange

(

const uint32_t &

detID

)

const

Definition at line 38 of file SiStripThreshold.cc.

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

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

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

getRegistryVectorBegin()

RegistryIterator SiStripThreshold::getRegistryVectorBegin ( ) const

inline

Definition at line 124 of file SiStripThreshold.h.

References indexes.

Referenced by printDebug(), and printSummary().

{ return indexes.begin(); }

getRegistryVectorEnd()

RegistryIterator SiStripThreshold::getRegistryVectorEnd ( ) const

inline

Definition at line 125 of file SiStripThreshold.h.

References indexes.

Referenced by printDebug(), and printSummary().

{ return indexes.end(); }

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.

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

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

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.

References addToStat(), getDataVectorBegin(), getRegistryVectorBegin(), getRegistryVectorEnd(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, create_idmaps::n, mathSSE::sqrt(), and contentValuesCheck::ss.

                                                                                       {
  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;
    for (; it != eit; ++it) {
      itp = it + 1;
      firstStrip = it->getFirstStrip();
      if (itp != eit)
        stripRange = (itp->getFirstStrip() - firstStrip);
      else
        //*FIXME In this class the number of strips in the detector is not known, therefore it is (wrongly) assumed that
        //       if the last firstStrip<511 then the detector has just 512 strips: this is clearly wrong, and it needs to be fixed
        stripRange = firstStrip > 511 ? 768 - firstStrip : 512 - firstStrip;
      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;
  }
}

put()

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, pfDeepBoostedJetPreprocessParams_cfi::lower_bound, AlCaHLTBitMon_ParallelJobs::p, and v_threshold.

Referenced by 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;
}

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,
		Container &	vthr
	)

Definition at line 57 of file SiStripThreshold.cc.

References a, and nano_mu_digi_cff::strip.

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

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

setData() [2/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.

References a, and nano_mu_digi_cff::strip.

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

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().