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 103 of file SiStripThreshold.h.

typedef Container::const_iterator SiStripThreshold::ContainerIterator

Definition at line 104 of file SiStripThreshold.h.

typedef Container SiStripThreshold::InputVector

Definition at line 108 of file SiStripThreshold.h.

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

Definition at line 105 of file SiStripThreshold.h.

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

Definition at line 106 of file SiStripThreshold.h.

typedef Registry::const_iterator SiStripThreshold::RegistryIterator

Definition at line 107 of file SiStripThreshold.h.

Constructor & Destructor Documentation

SiStripThreshold::SiStripThreshold ( )

inline

Definition at line 110 of file SiStripThreshold.h.

{};

SiStripThreshold::SiStripThreshold ( const SiStripThreshold &  orig )

inline

Definition at line 111 of file SiStripThreshold.h.

References indexes, and v_threshold.

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

virtual SiStripThreshold::~SiStripThreshold ( )

inlinevirtual

Definition at line 115 of file SiStripThreshold.h.

References getDetIds(), getRange(), and put().

{};

Member Function Documentation

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

private

Definition at line 160 of file SiStripThreshold.cc.

References relativeConstraints::value.

Referenced by getRegistryVectorEnd(), and 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 80 of file SiStripThreshold.cc.

References lumiContext::fill, and RecoTauPiZeroBuilderPlugins_cfi::strips.

Referenced by getRegistryVectorEnd().

                                                                                                                   {
    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 31 of file SiStripThreshold.cc.

References tier0::unique().

Referenced by getRegistryVectorEnd(), and 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 68 of file SiStripThreshold.cc.

References Exception, sistrip::FirstThStripMask_, sistrip::FirstThStripShift_, sistrip::HighThStripMask_, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by SiStripThresholdReader::analyze(), SiStripThresholdDQM::fillMEsForDet(), SiStripThresholdDQM::fillMEsForLayer(), getRegistryVectorEnd(), 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 47 of file SiStripThreshold.cc.

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

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

                                                                   {
  // 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 36 of file SiStripThreshold.cc.

References indexes, AlCaHLTBitMon_ParallelJobs::p, and v_threshold.

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

                                                                                  {
  // 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 addToStat(), allThresholds(), compact(), getData(), indexes, input, SiStripPI::max, min(), printDebug(), printSummary(), setData(), digi_MixPreMix_cfi::strip, and combinedConstraintHelpers::sum2().

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

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

Referenced by 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 115 of file SiStripThreshold.cc.

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

Referenced by getRegistryVectorEnd().

                                                                                      {
  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, AlCaHLTBitMon_ParallelJobs::p, sd, and v_threshold.

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

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

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

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

Referenced by SiStripThresholdBuilder::analyze(), getRegistryVectorEnd(), 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 62 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 147 of file SiStripThreshold.h.

template<typename CondSerializationT , typename Enabled >

friend struct cond::serialization::access

friend

Definition at line 147 of file SiStripThreshold.h.

Member Data Documentation

Registry SiStripThreshold::indexes

private

Definition at line 145 of file SiStripThreshold.h.

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

Container SiStripThreshold::v_threshold

private

Definition at line 144 of file SiStripThreshold.h.

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