StripClusterizerAlgorithm Class Reference

#include <StripClusterizerAlgorithm.h>

Inheritance diagram for StripClusterizerAlgorithm:

Classes
struct	Det
struct	Index
struct	InvalidChargeException
struct	State

Public Types
typedef edmNew::DetSetVector< SiStripCluster >	output_t

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

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

Protected Member Functions
Det	findDetId (const uint32_t) const
bool	isModuleBad (const uint32_t &id) const
bool	isModuleUsable (const uint32_t &id) const
	StripClusterizerAlgorithm ()

Protected Attributes
std::string	qualityLabel

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

Private Attributes
std::vector< std::vector< const FedChannelConnection * > >	connections
std::vector< uint32_t >	detIds
uint32_t	gain_cache_id
edm::ESHandle< SiStripGain >	gainHandle
std::vector< Index >	indices
uint32_t	noise_cache_id
edm::ESHandle< SiStripNoises >	noiseHandle
uint32_t	quality_cache_id
edm::ESHandle< SiStripQuality >	qualityHandle
SiStripDetCabling const *	theCabling = 0

Detailed Description

Definition at line 20 of file StripClusterizerAlgorithm.h.

Member Typedef Documentation

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

Definition at line 60 of file StripClusterizerAlgorithm.h.

Constructor & Destructor Documentation

virtual StripClusterizerAlgorithm::~StripClusterizerAlgorithm ( )

inlinevirtual

Definition at line 56 of file StripClusterizerAlgorithm.h.

References initialize().

{}

StripClusterizerAlgorithm::StripClusterizerAlgorithm ( )

inlineprotected

Definition at line 96 of file StripClusterizerAlgorithm.h.

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

Member Function Documentation

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

inlinevirtual

Definition at line 69 of file StripClusterizerAlgorithm.h.

Referenced by ThreeThresholdAlgorithm::addFed().

                                                            {}

virtual void StripClusterizerAlgorithm::addFed ( State &  state, sistrip::FEDZSChannelUnpacker &  unpacker, uint16_t  ipair, output_t::TSFastFiller &  out  ) const

inlinevirtual

Reimplemented in ThreeThresholdAlgorithm.

Definition at line 76 of file StripClusterizerAlgorithm.h.

                                                       {}

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

inline

Definition at line 89 of file StripClusterizerAlgorithm.h.

{ return detIds; }

SiStripDetCabling const* StripClusterizerAlgorithm::cabling ( ) const

inline

Definition at line 88 of file StripClusterizerAlgorithm.h.

Referenced by initialize().

{ return theCabling; }

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

Definition at line 196 of file StripClusterizerAlgorithm.cc.

References clusterize_().

                                                                                                            {
  clusterize_(input, output);
}

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

Definition at line 199 of file StripClusterizerAlgorithm.cc.

References clusterize_().

                                                                                                               {
  clusterize_(input, output);
}

template<class T >

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

inlineprivate

Definition at line 106 of file StripClusterizerAlgorithm.h.

References alignCSCRings::ff.

Referenced by clusterize().

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

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

pure virtual

Implemented in ThreeThresholdAlgorithm.

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

pure virtual

Implemented in ThreeThresholdAlgorithm.

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

inline

Definition at line 91 of file StripClusterizerAlgorithm.h.

References StripClusterizerAlgorithm::Det::ind.

                                                                                        {
    return connections[det.ind];
  }

StripClusterizerAlgorithm::Det StripClusterizerAlgorithm::findDetId ( const uint32_t  id ) const

protected

Definition at line 161 of file StripClusterizerAlgorithm.cc.

References b, StripClusterizerAlgorithm::Det::detId, detIds, MillePedeFileConverter_cfg::e, gainHandle, StripClusterizerAlgorithm::Det::gainRange, SiStripNoises::getRange(), SiStripBadStrip::getRange(), SiStripGain::getRange(), SiStripNoises::getRangeByPos(), SiStripBadStrip::getRangeByPos(), SiStripGain::getRangeByPos(), triggerObjects_cff::id, StripClusterizerAlgorithm::Det::ind, indices, isModuleUsable(), noiseHandle, StripClusterizerAlgorithm::Det::noiseRange, AlCaHLTBitMon_ParallelJobs::p, edm::ESHandle< T >::product(), StripClusterizerAlgorithm::Det::quality, qualityHandle, and StripClusterizerAlgorithm::Det::qualityRange.

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

                                                                                         {
  auto b = detIds.begin();
  auto e = detIds.end();
  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 Det();
  }
  Det det;
  det.ind = p - detIds.begin();

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

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

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

virtual

Definition at line 24 of file StripClusterizerAlgorithm.cc.

References HltBtagPostValidation_cff::c, cabling(), SiStripQuality::cabling(), SiStripCommissioningClient_cfg::conn, SiStripDetCabling::connected(), connections, COUT, detIds, gain_cache_id, gainHandle, edm::EventSetup::get(), SiStripDetCabling::getDetCabling(), SiStripNoises::getDetIds(), SiStripBadStrip::getDetIds(), SiStripGain::getDetIds(), mps_fire::i, indices, invalidI, isModuleBad(), dqmiolumiharvest::j, dqmdumpme::k, mod(), groupFilesInBlocks::nn, noise_cache_id, noiseHandle, quality_cache_id, qualityHandle, qualityLabel, theCabling, and mitigatedMETSequence_cff::U.

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

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

inlineprotected

Definition at line 99 of file StripClusterizerAlgorithm.h.

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

{ return qualityHandle->IsModuleBad(id); }

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

inlineprotected

Definition at line 100 of file StripClusterizerAlgorithm.h.

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

{ return qualityHandle->IsModuleUsable(id); }

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

inlinevirtual

Reimplemented in ThreeThresholdAlgorithm.

Definition at line 73 of file StripClusterizerAlgorithm.h.

{}

virtual void StripClusterizerAlgorithm::stripByStripAdd ( State &  state, uint16_t  strip, uint8_t  adc, output_t::TSFastFiller &  out  ) const

inlinevirtual

Reimplemented in ThreeThresholdAlgorithm.

Definition at line 80 of file StripClusterizerAlgorithm.h.

{}

virtual Det StripClusterizerAlgorithm::stripByStripBegin ( uint32_t  id ) const

pure virtual

Implemented in ThreeThresholdAlgorithm.

virtual void StripClusterizerAlgorithm::stripByStripEnd ( State &  state, std::vector< SiStripCluster > &  out  ) const

inlinevirtual

Reimplemented in ThreeThresholdAlgorithm.

Definition at line 74 of file StripClusterizerAlgorithm.h.

{}

virtual void StripClusterizerAlgorithm::stripByStripEnd ( State &  state, output_t::TSFastFiller &  out  ) const

inlinevirtual

Reimplemented in ThreeThresholdAlgorithm.

Definition at line 81 of file StripClusterizerAlgorithm.h.

{}

Member Data Documentation

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

private

Definition at line 119 of file StripClusterizerAlgorithm.h.

Referenced by initialize().

std::vector<uint32_t> StripClusterizerAlgorithm::detIds

private

Definition at line 118 of file StripClusterizerAlgorithm.h.

Referenced by findDetId(), and initialize().

uint32_t StripClusterizerAlgorithm::gain_cache_id

private

Definition at line 125 of file StripClusterizerAlgorithm.h.

Referenced by initialize().

edm::ESHandle<SiStripGain> StripClusterizerAlgorithm::gainHandle

private

Definition at line 121 of file StripClusterizerAlgorithm.h.

Referenced by findDetId(), and initialize().

std::vector<Index> StripClusterizerAlgorithm::indices

private

Definition at line 120 of file StripClusterizerAlgorithm.h.

Referenced by findDetId(), and initialize().

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

static

Definition at line 22 of file StripClusterizerAlgorithm.h.

Referenced by initialize().

uint32_t StripClusterizerAlgorithm::noise_cache_id

private

Definition at line 125 of file StripClusterizerAlgorithm.h.

Referenced by initialize().

edm::ESHandle<SiStripNoises> StripClusterizerAlgorithm::noiseHandle

private

Definition at line 122 of file StripClusterizerAlgorithm.h.

Referenced by findDetId(), and initialize().

uint32_t StripClusterizerAlgorithm::quality_cache_id

private

Definition at line 125 of file StripClusterizerAlgorithm.h.

Referenced by initialize().

edm::ESHandle<SiStripQuality> StripClusterizerAlgorithm::qualityHandle

private

Definition at line 123 of file StripClusterizerAlgorithm.h.

Referenced by findDetId(), and initialize().

std::string StripClusterizerAlgorithm::qualityLabel

protected

Definition at line 102 of file StripClusterizerAlgorithm.h.

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

SiStripDetCabling const* StripClusterizerAlgorithm::theCabling = 0

private

Definition at line 124 of file StripClusterizerAlgorithm.h.

Referenced by initialize().