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 19 of file StripClusterizerAlgorithm.h.

Member Typedef Documentation

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

Definition at line 57 of file StripClusterizerAlgorithm.h.

Constructor & Destructor Documentation

virtual StripClusterizerAlgorithm::~StripClusterizerAlgorithm ( )

inlinevirtual

Definition at line 52 of file StripClusterizerAlgorithm.h.

References initialize().

{}

StripClusterizerAlgorithm::StripClusterizerAlgorithm ( )

inlineprotected

Definition at line 85 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 66 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 70 of file StripClusterizerAlgorithm.h.

{}

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

inline

Definition at line 79 of file StripClusterizerAlgorithm.h.

{ return detIds;}

SiStripDetCabling const* StripClusterizerAlgorithm::cabling ( ) const

inline

Definition at line 78 of file StripClusterizerAlgorithm.h.

Referenced by initialize().

{ return theCabling;} 

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

Definition at line 166 of file StripClusterizerAlgorithm.cc.

References clusterize_().

{clusterize_(input, output);}

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

Definition at line 167 of file StripClusterizerAlgorithm.cc.

References clusterize_(), and StripClusterizerAlgorithm::InvalidChargeException::InvalidChargeException().

{clusterize_(input, output);}

template<class T >

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

inlineprivate

Definition at line 95 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 81 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 130 of file StripClusterizerAlgorithm.cc.

References b, StripClusterizerAlgorithm::Det::detId, detIds, MillePedeFileConverter_cfg::e, gainHandle, StripClusterizerAlgorithm::Det::gainRange, SiStripNoises::getRange(), SiStripGain::getRange(), SiStripBadStrip::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_(), initialize(), 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 25 of file StripClusterizerAlgorithm.cc.

References EnergyCorrector::c, cabling(), SiStripQuality::cabling(), SiStripDetCabling::connected(), connections, COUT, detIds, findDetId(), gain_cache_id, gainHandle, edm::EventSetup::get(), SiStripDetCabling::getDetCabling(), SiStripNoises::getDetIds(), SiStripBadStrip::getDetIds(), SiStripGain::getDetIds(), mps_fire::i, indices, invalidI, isModuleBad(), gen::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 88 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 89 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 67 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 71 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 68 of file StripClusterizerAlgorithm.h.

{}

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

inlinevirtual

Reimplemented in ThreeThresholdAlgorithm.

Definition at line 72 of file StripClusterizerAlgorithm.h.

{}

Member Data Documentation

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

private

Definition at line 110 of file StripClusterizerAlgorithm.h.

Referenced by initialize().

std::vector<uint32_t> StripClusterizerAlgorithm::detIds

private

Definition at line 109 of file StripClusterizerAlgorithm.h.

Referenced by findDetId(), and initialize().

uint32_t StripClusterizerAlgorithm::gain_cache_id

private

Definition at line 116 of file StripClusterizerAlgorithm.h.

Referenced by initialize().

edm::ESHandle<SiStripGain> StripClusterizerAlgorithm::gainHandle

private

Definition at line 112 of file StripClusterizerAlgorithm.h.

Referenced by findDetId(), and initialize().

std::vector<Index> StripClusterizerAlgorithm::indices

private

Definition at line 111 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 116 of file StripClusterizerAlgorithm.h.

Referenced by initialize().

edm::ESHandle<SiStripNoises> StripClusterizerAlgorithm::noiseHandle

private

Definition at line 113 of file StripClusterizerAlgorithm.h.

Referenced by findDetId(), and initialize().

uint32_t StripClusterizerAlgorithm::quality_cache_id

private

Definition at line 116 of file StripClusterizerAlgorithm.h.

Referenced by initialize().

edm::ESHandle<SiStripQuality> StripClusterizerAlgorithm::qualityHandle

private

Definition at line 114 of file StripClusterizerAlgorithm.h.

Referenced by findDetId(), and initialize().

std::string StripClusterizerAlgorithm::qualityLabel

protected

Definition at line 91 of file StripClusterizerAlgorithm.h.

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

SiStripDetCabling const* StripClusterizerAlgorithm::theCabling = 0

private

Definition at line 115 of file StripClusterizerAlgorithm.h.

Referenced by initialize().