StripClusterizerAlgorithm Class Reference

#include <StripClusterizerAlgorithm.h>

Inheritance diagram for StripClusterizerAlgorithm:

Classes
struct	Index
struct	InvalidChargeException

Public Types
typedef edmNew::DetSetVector < SiStripCluster >	output_t

Public Member Functions
virtual void	addFed (sistrip::FEDZSChannelUnpacker &unpacker, uint16_t ipair, std::vector< SiStripCluster > &out)
virtual void	addFed (sistrip::FEDZSChannelUnpacker &unpacker, uint16_t ipair, output_t::FastFiller &out)
std::vector< uint32_t > const &	allDetIds () const
SiStripDetCabling const *	cabling () const
void	clusterize (const edm::DetSetVector< SiStripDigi > &, output_t &)
void	clusterize (const edmNew::DetSetVector< SiStripDigi > &, output_t &)
virtual void	clusterizeDetUnit (const edm::DetSet< SiStripDigi > &, output_t::FastFiller &)=0
virtual void	clusterizeDetUnit (const edmNew::DetSet< SiStripDigi > &, output_t::FastFiller &)=0
std::vector< const FedChannelConnection * > const &	currentConnection () const
virtual void	initialize (const edm::EventSetup &)
virtual void	stripByStripAdd (uint16_t strip, uint8_t adc, std::vector< SiStripCluster > &out)
virtual void	stripByStripAdd (uint16_t strip, uint8_t adc, output_t::FastFiller &out)
virtual bool	stripByStripBegin (uint32_t id)=0
virtual void	stripByStripEnd (std::vector< SiStripCluster > &out)
virtual void	stripByStripEnd (output_t::FastFiller &out)
virtual	~StripClusterizerAlgorithm ()

Protected Member Functions
bool	allBadBetween (uint16_t L, const uint16_t &R) const
bool	bad (const uint16_t &strip) const
uint32_t	currentId ()
float	gain (const uint16_t &strip) const
bool	isModuleBad (const uint32_t &id) const
bool	isModuleUsable (const uint32_t &id) const
float	noise (const uint16_t &strip) const
bool	setDetId (const uint32_t)
	StripClusterizerAlgorithm ()

Protected Attributes
bool	_setDetId
std::string	qualityLabel

Private Member Functions
template<class T >
void	clusterize_ (const T &input, output_t &output)

Private Attributes
std::vector< std::vector < const FedChannelConnection * > >	connections
uint32_t	detId =0
std::vector< uint32_t >	detIds
uint32_t	gain_cache_id
edm::ESHandle< SiStripGain >	gainHandle
SiStripApvGain::Range	gainRange
unsigned short	ind =invalidI
std::vector< Index >	indices
uint32_t	noise_cache_id
edm::ESHandle< SiStripNoises >	noiseHandle
SiStripNoises::Range	noiseRange
uint32_t	quality_cache_id
edm::ESHandle< SiStripQuality >	qualityHandle
SiStripQuality::Range	qualityRange
SiStripDetCabling const *	theCabling = nullptr

Static Private Attributes
static constexpr unsigned short	invalidI = std::numeric_limits<unsigned short>::max()

Detailed Description

Definition at line 19 of file StripClusterizerAlgorithm.h.

Member Typedef Documentation

typedef edmNew::DetSetVector<SiStripCluster> StripClusterizerAlgorithm::output_t

Definition at line 27 of file StripClusterizerAlgorithm.h.

Constructor & Destructor Documentation

virtual StripClusterizerAlgorithm::~StripClusterizerAlgorithm ( )

inlinevirtual

Definition at line 23 of file StripClusterizerAlgorithm.h.

{}

StripClusterizerAlgorithm::StripClusterizerAlgorithm ( )

inlineprotected

Definition at line 54 of file StripClusterizerAlgorithm.h.

: qualityLabel(""), noise_cache_id(0), gain_cache_id(0), quality_cache_id(0) {}

Member Function Documentation

virtual void StripClusterizerAlgorithm::addFed ( sistrip::FEDZSChannelUnpacker &  unpacker, uint16_t  ipair, std::vector< SiStripCluster > &  out  )

inlinevirtual

Reimplemented in ThreeThresholdAlgorithm.

Definition at line 36 of file StripClusterizerAlgorithm.h.

Referenced by sistrip::RawToClustersLazyUnpacker::fill().

{}

virtual void StripClusterizerAlgorithm::addFed ( sistrip::FEDZSChannelUnpacker &  unpacker, uint16_t  ipair, output_t::FastFiller &  out  )

inlinevirtual

Reimplemented in ThreeThresholdAlgorithm.

Definition at line 40 of file StripClusterizerAlgorithm.h.

{}

bool StripClusterizerAlgorithm::allBadBetween ( uint16_t  L, const uint16_t &  R  ) const

inlineprotected

Definition at line 63 of file StripClusterizerAlgorithm.h.

References bad(), and dttmaxenums::R.

Referenced by ThreeThresholdAlgorithm::candidateEnded().

{ while( ++L < R  &&  bad(L) ); return L == R; }

std::vector<uint32_t> const& StripClusterizerAlgorithm::allDetIds ( ) const

inline

Definition at line 49 of file StripClusterizerAlgorithm.h.

References detIds.

{ return detIds;}

bool StripClusterizerAlgorithm::bad ( const uint16_t &  strip ) const

inlineprotected

Definition at line 60 of file StripClusterizerAlgorithm.h.

References qualityHandle, and qualityRange.

Referenced by ThreeThresholdAlgorithm::addToCandidate(), allBadBetween(), and ThreeThresholdAlgorithm::appendBadNeighbors().

{ return qualityHandle->IsStripBad( qualityRange, strip ); }

SiStripDetCabling const* StripClusterizerAlgorithm::cabling ( ) const

inline

Definition at line 48 of file StripClusterizerAlgorithm.h.

References theCabling.

Referenced by initialize().

{ return theCabling;} 

void StripClusterizerAlgorithm::clusterize	(	const edm::DetSetVector< SiStripDigi > &	input,
		output_t &	output
	)

Definition at line 191 of file StripClusterizerAlgorithm.cc.

References clusterize_().

{clusterize_(input, output);}

void StripClusterizerAlgorithm::clusterize	(	const edmNew::DetSetVector< SiStripDigi > &	input,
		output_t &	output
	)

Definition at line 192 of file StripClusterizerAlgorithm.cc.

References clusterize_().

{clusterize_(input, output);}

template<class T >

void StripClusterizerAlgorithm::clusterize_ ( const T &  input, output_t &  output  )

inlineprivate

Definition at line 70 of file StripClusterizerAlgorithm.h.

References clusterizeDetUnit(), and alignCSCRings::ff.

Referenced by clusterize().

                                                                       {
    for(typename T::const_iterator it = input.begin(); it!=input.end(); it++) {
      output_t::FastFiller ff(output, it->detId()); 
      clusterizeDetUnit(*it, ff);   
      if(ff.empty()) ff.abort();    
    }   
  }

virtual void StripClusterizerAlgorithm::clusterizeDetUnit ( const edm::DetSet< SiStripDigi > &  , output_t::FastFiller &    )

pure virtual

Implemented in ThreeThresholdAlgorithm.

Referenced by clusterize_().

virtual void StripClusterizerAlgorithm::clusterizeDetUnit ( const edmNew::DetSet< SiStripDigi > &  , output_t::FastFiller &    )

pure virtual

Implemented in ThreeThresholdAlgorithm.

std::vector<const FedChannelConnection *> const& StripClusterizerAlgorithm::currentConnection ( ) const

inline

Definition at line 50 of file StripClusterizerAlgorithm.h.

References connections, and ind.

{ return connections[ind]; }

uint32_t StripClusterizerAlgorithm::currentId ( )

inlineprotected

Definition at line 56 of file StripClusterizerAlgorithm.h.

References detId.

Referenced by ThreeThresholdAlgorithm::endCandidate().

{return detId;}

float StripClusterizerAlgorithm::gain ( const uint16_t &  strip ) const

inlineprotected

Definition at line 59 of file StripClusterizerAlgorithm.h.

References gainRange, and SiStripGain::getStripGain().

Referenced by ThreeThresholdAlgorithm::applyGains().

{ return SiStripGain::getStripGain( strip, gainRange ); }

void StripClusterizerAlgorithm::initialize ( const edm::EventSetup &  es )

virtual

Reimplemented in OldThreeThresholdAlgorithm.

Definition at line 25 of file StripClusterizerAlgorithm.cc.

References trackerHits::c, cabling(), results_mgr::conn, SiStripDetCabling::connected(), connections, COUT, detId, detIds, gain_cache_id, gainHandle, gainRange, edm::EventSetup::get(), SiStripDetCabling::getDetCabling(), i, ind, indices, invalidI, isModuleBad(), j, gen::k, mod(), noise_cache_id, noiseHandle, noiseRange, quality_cache_id, qualityHandle, qualityLabel, qualityRange, and theCabling.

                                    {
  uint32_t n_cache_id = es.get<SiStripNoisesRcd>().cacheIdentifier();
  uint32_t g_cache_id = es.get<SiStripGainRcd>().cacheIdentifier();
  uint32_t q_cache_id = es.get<SiStripQualityRcd>().cacheIdentifier();

  bool mod=false;
  if(g_cache_id != gain_cache_id) {
    es.get<SiStripGainRcd>().get( gainHandle );
    gain_cache_id = g_cache_id;
    mod=true;
  }
  if(n_cache_id != noise_cache_id) {
    es.get<SiStripNoisesRcd>().get( noiseHandle );
    noise_cache_id = n_cache_id;
    mod=true;
  }
  if(q_cache_id != quality_cache_id) {
    es.get<SiStripQualityRcd>().get( qualityLabel, qualityHandle );
    quality_cache_id = q_cache_id;
    mod=true;
  }

  if (mod) { 
    // redo indexing!
    theCabling = qualityHandle->cabling();
    assert(theCabling); 
    auto const & conn = cabling()->connected();
    COUT << "cabling " << conn.size() << std::endl;
    detIds.clear();
    detIds.reserve(conn.size());
    for (auto const & c : conn) { if (!isModuleBad(c.first)) detIds.push_back(c.first);}
    indices.clear();
    indices.resize(detIds.size());
    COUT << "good detIds " << detIds.size() << std::endl;

    if (0==detIds.size()) {
       ind=0; detId=0; return;
    }

    {
      connections.clear();
      connections.resize(detIds.size());
      //connections (slow, not a big deal)
      auto const & conns = cabling()->getDetCabling();
      for (auto i=0U; i<detIds.size(); ++i) {
    auto c = conns.find(detIds[i]);
    if (c!=conns.end()) connections[i]=(*c).second;
      }
      

    }


    { // quality
      std::vector<uint32_t> dum; qualityHandle->getDetIds(dum); 
      assert(dum.size()<invalidI);
      unsigned short j=0, i=0;
      while (i<dum.size() && j<detIds.size()) {
    if (dum[i]<detIds[j]) ++i;
    else if (detIds[j]<dum[i]) {indices[j].qi=invalidI; ++j;}
    else {
      indices[j].qi=i; ++i; ++j;
    }
      }
      unsigned int  nn=0;
      for(auto k=0U; k<detIds.size();++k) { if (indices[k].qi<invalidI) {++nn; assert(dum[indices[k].qi]==detIds[k]);}}
      assert(nn<=dum.size());
      COUT << "quality " << dum.size() << " " <<nn<< std::endl;
    }
    { //noise
      std::vector<uint32_t> dum; noiseHandle->getDetIds(dum); 
      assert(dum.size()<invalidI);
      unsigned short j=0, i=0;
      while (i<dum.size() && j<detIds.size()) {
    if (dum[i]<detIds[j]) ++i;
    else if (detIds[j]<dum[i]) {indices[j].ni=invalidI; ++j;}
    else {
      indices[j].ni=i; ++i; ++j;
    }
      }
      unsigned int  nn=0;
      for(auto k=0U; k<detIds.size();++k) { if (indices[k].ni<invalidI) {++nn; assert(dum[indices[k].ni]==detIds[k]);}}
      assert(nn<=dum.size());
      COUT << "noise " << dum.size() << " " <<nn<< std::endl;
    }
    { //gain
      std::vector<uint32_t> dum; gainHandle->getDetIds(dum); 
      assert(dum.size()<invalidI);
      unsigned short j=0, i=0;
      while (i<dum.size() && j<detIds.size()) {
    if (dum[i]<detIds[j]) ++i;
    else if (detIds[j]<dum[i]) {indices[j].gi=invalidI; ++j;}
    else {
      indices[j].gi=i; ++i; ++j;
    }
      }
      unsigned int  nn=0;
      for(auto k=0U; k<detIds.size();++k) { if (indices[k].gi<invalidI) {++nn; assert(dum[indices[k].gi]==detIds[k]);}}
      assert(nn<=dum.size());
      COUT << "gain " << dum.size() << " " <<nn<< std::endl;
    }
  }

  if (0==detIds.size()) {
       ind=0; detId=0; return;
  }

  // initalize first det
  ind = 0;
  detId = detIds[ind];
  
  noiseRange = noiseHandle->getRangeByPos(indices[ind].ni);
  gainRange = gainHandle->getRangeByPos(indices[ind].gi);
  qualityRange = qualityHandle->getRangeByPos(indices[ind].qi);
  
}

bool StripClusterizerAlgorithm::isModuleBad ( const uint32_t &  id ) const

inlineprotected

Definition at line 61 of file StripClusterizerAlgorithm.h.

References qualityHandle.

Referenced by ThreeThresholdAlgorithm::clusterizeDetUnit_(), and initialize().

{ return qualityHandle->IsModuleBad( id ); }

bool StripClusterizerAlgorithm::isModuleUsable ( const uint32_t &  id ) const

inlineprotected

Definition at line 62 of file StripClusterizerAlgorithm.h.

References qualityHandle.

Referenced by ThreeThresholdAlgorithm::clusterizeDetUnit_(), and ThreeThresholdAlgorithm::stripByStripBegin().

{ return qualityHandle->IsModuleUsable( id ); }

float StripClusterizerAlgorithm::noise ( const uint16_t &  strip ) const

inlineprotected

Definition at line 58 of file StripClusterizerAlgorithm.h.

References SiStripNoises::getNoise(), and noiseRange.

Referenced by ThreeThresholdAlgorithm::addToCandidate().

{ return SiStripNoises::getNoise( strip, noiseRange ); }

bool StripClusterizerAlgorithm::setDetId ( const uint32_t  id )

protected

Definition at line 144 of file StripClusterizerAlgorithm.cc.

References b, detId, detIds, alignCSCRings::e, gainHandle, gainRange, ind, indices, noiseHandle, noiseRange, AlCaHLTBitMon_ParallelJobs::p, qualityHandle, qualityRange, and unlikely.

Referenced by ThreeThresholdAlgorithm::clusterizeDetUnit_(), OldThreeThresholdAlgorithm::OldThreeThresholdAlgorithm(), ThreeThresholdAlgorithm::stripByStripBegin(), and ThreeThresholdAlgorithm::ThreeThresholdAlgorithm().

                            {
  if (id==detId) return true; // rare....
  // priority to sequential scan
  if (id==detIds[ind+1]) {
    ++ind;
  } else if unlikely( (id>detId) &  (id<detIds[ind+1]) ) {
      /*
      edm::LogWarning("StripClusterizerAlgorithm") 
    <<"id " << id << " not connected. this is impossible on data "
    << "old id " << detId << std::endl;
      */
      return false; 
    }   else{
    auto b = detIds.begin();
    auto e = detIds.end();
    if (id>detId) b+=ind;
    else e=b+ind;
    auto p = std::lower_bound(b,e,id);
    if (p==e || id!=(*p)) {
#ifdef NOT_ON_MONTECARLO
      edm::LogWarning("StripClusterizerAlgorithm") 
    <<"id " << id << " not connected. this is impossible on data "
    << "old id " << detId << std::endl;
#endif
      return false;
    }
    ind = p-detIds.begin();
  }

  detId = id;
  noiseRange = noiseHandle->getRangeByPos(indices[ind].ni);
  gainRange = gainHandle->getRangeByPos(indices[ind].gi);
  qualityRange = qualityHandle->getRangeByPos(indices[ind].qi);

#ifdef EDM_ML_DEBUG
  assert(detIds[ind]==detId); 
  auto oldg =  gainHandle->getRange(id);
  assert(oldg==gainRange);
  auto oldn = noiseHandle->getRange(id);
  assert(oldn==noiseRange);
  auto oldq = qualityHandle->getRange(id);
  assert(oldq==qualityRange);
#endif

  return true;
}

virtual void StripClusterizerAlgorithm::stripByStripAdd ( uint16_t  strip, uint8_t  adc, std::vector< SiStripCluster > &  out  )

inlinevirtual

Reimplemented in ThreeThresholdAlgorithm.

Definition at line 37 of file StripClusterizerAlgorithm.h.

Referenced by sistrip::RawToClustersLazyUnpacker::fill().

{}

virtual void StripClusterizerAlgorithm::stripByStripAdd ( uint16_t  strip, uint8_t  adc, output_t::FastFiller &  out  )

inlinevirtual

Reimplemented in ThreeThresholdAlgorithm.

Definition at line 41 of file StripClusterizerAlgorithm.h.

{}

virtual bool StripClusterizerAlgorithm::stripByStripBegin ( uint32_t  id )

pure virtual

Implemented in OldThreeThresholdAlgorithm, and ThreeThresholdAlgorithm.

Referenced by sistrip::RawToClustersLazyUnpacker::fill().

virtual void StripClusterizerAlgorithm::stripByStripEnd ( std::vector< SiStripCluster > &  out )

inlinevirtual

Reimplemented in ThreeThresholdAlgorithm.

Definition at line 38 of file StripClusterizerAlgorithm.h.

Referenced by sistrip::RawToClustersLazyUnpacker::fill().

{}

virtual void StripClusterizerAlgorithm::stripByStripEnd ( output_t::FastFiller &  out )

inlinevirtual

Reimplemented in ThreeThresholdAlgorithm.

Definition at line 42 of file StripClusterizerAlgorithm.h.

{}

Member Data Documentation

bool StripClusterizerAlgorithm::_setDetId

protected

Definition at line 66 of file StripClusterizerAlgorithm.h.

Referenced by OldThreeThresholdAlgorithm::clusterizeDetUnit_(), ThreeThresholdAlgorithm::endCandidate(), OldThreeThresholdAlgorithm::OldThreeThresholdAlgorithm(), and ThreeThresholdAlgorithm::ThreeThresholdAlgorithm().

std::vector<std::vector<const FedChannelConnection *> > StripClusterizerAlgorithm::connections

private

Definition at line 86 of file StripClusterizerAlgorithm.h.

Referenced by currentConnection(), and initialize().

uint32_t StripClusterizerAlgorithm::detId =0

private

Definition at line 95 of file StripClusterizerAlgorithm.h.

Referenced by currentId(), initialize(), and setDetId().

std::vector<uint32_t> StripClusterizerAlgorithm::detIds

private

Definition at line 85 of file StripClusterizerAlgorithm.h.

Referenced by allDetIds(), initialize(), and setDetId().

uint32_t StripClusterizerAlgorithm::gain_cache_id

private

Definition at line 95 of file StripClusterizerAlgorithm.h.

Referenced by initialize().

edm::ESHandle<SiStripGain> StripClusterizerAlgorithm::gainHandle

private

Definition at line 91 of file StripClusterizerAlgorithm.h.

Referenced by initialize(), and setDetId().

SiStripApvGain::Range StripClusterizerAlgorithm::gainRange

private

Definition at line 88 of file StripClusterizerAlgorithm.h.

Referenced by gain(), initialize(), and setDetId().

unsigned short StripClusterizerAlgorithm::ind =invalidI

private

Definition at line 96 of file StripClusterizerAlgorithm.h.

Referenced by currentConnection(), initialize(), and setDetId().

std::vector<Index> StripClusterizerAlgorithm::indices

private

Definition at line 87 of file StripClusterizerAlgorithm.h.

Referenced by initialize(), and setDetId().

constexpr unsigned short StripClusterizerAlgorithm::invalidI = std::numeric_limits<unsigned short>::max()

staticprivate

Definition at line 78 of file StripClusterizerAlgorithm.h.

Referenced by initialize().

uint32_t StripClusterizerAlgorithm::noise_cache_id

private

Definition at line 95 of file StripClusterizerAlgorithm.h.

Referenced by initialize().

edm::ESHandle<SiStripNoises> StripClusterizerAlgorithm::noiseHandle

private

Definition at line 92 of file StripClusterizerAlgorithm.h.

Referenced by initialize(), and setDetId().

SiStripNoises::Range StripClusterizerAlgorithm::noiseRange

private

Definition at line 89 of file StripClusterizerAlgorithm.h.

Referenced by initialize(), noise(), and setDetId().

uint32_t StripClusterizerAlgorithm::quality_cache_id

private

Definition at line 95 of file StripClusterizerAlgorithm.h.

Referenced by initialize().

edm::ESHandle<SiStripQuality> StripClusterizerAlgorithm::qualityHandle

private

Definition at line 93 of file StripClusterizerAlgorithm.h.

Referenced by bad(), initialize(), isModuleBad(), isModuleUsable(), and setDetId().

std::string StripClusterizerAlgorithm::qualityLabel

protected

Definition at line 65 of file StripClusterizerAlgorithm.h.

Referenced by initialize(), and ThreeThresholdAlgorithm::ThreeThresholdAlgorithm().

SiStripQuality::Range StripClusterizerAlgorithm::qualityRange

private

Definition at line 90 of file StripClusterizerAlgorithm.h.

Referenced by bad(), initialize(), and setDetId().

SiStripDetCabling const* StripClusterizerAlgorithm::theCabling = nullptr

private

Definition at line 94 of file StripClusterizerAlgorithm.h.

Referenced by cabling(), and initialize().