SiStripClusterInfo Class Reference

#include <SiStripClusterInfo.h>

Public Member Functions
float	baryStrip () const
uint16_t	charge () const
std::pair< uint16_t, uint16_t >	chargeLR () const
const SiStripCluster *	cluster () const
uint32_t	detId () const
uint16_t	firstStrip () const
void	initEvent (const edm::EventSetup &iSetup)
bool	IsAnythingBad () const
bool	IsApvBad () const
bool	IsFiberBad () const
bool	IsModuleBad () const
bool	IsModuleUsable () const
uint8_t	maxCharge () const
uint16_t	maxIndex () const
uint16_t	maxStrip () const
float	noise () const
float	noiseRescaledByGain () const
void	setCluster (const SiStripCluster &cluster, int detId)
float	signalOverNoise () const
	SiStripClusterInfo (edm::ConsumesCollector &&, const std::string &qualityLabel="")
const SiStripGain *	siStripGain () const
const SiStripQuality *	siStripQuality () const
auto	stripCharges () const -> decltype(cluster() ->amplitudes())
std::vector< float >	stripGains () const
std::vector< float >	stripNoises () const
std::vector< float >	stripNoisesRescaledByGain () const
std::vector< bool >	stripQualitiesBad () const
float	variance () const
uint16_t	width () const

Private Member Functions
float	calculate_noise (const std::vector< float > &) const

Private Attributes
const SiStripCluster *	cluster_ptr = nullptr
uint32_t	detId_ = 0
const SiStripGain *	siStripGain_ = nullptr
edm::ESGetToken< SiStripGain, SiStripGainRcd >	siStripGainToken_
const SiStripNoises *	siStripNoises_ = nullptr
edm::ESGetToken< SiStripNoises, SiStripNoisesRcd >	siStripNoisesToken_
const SiStripQuality *	siStripQuality_ = nullptr
edm::ESGetToken< SiStripQuality, SiStripQualityRcd >	siStripQualityToken_

Detailed Description

Definition at line 21 of file SiStripClusterInfo.h.

Constructor & Destructor Documentation

SiStripClusterInfo()

SiStripClusterInfo::SiStripClusterInfo	(	edm::ConsumesCollector &&	iC,
		const std::string &	qualityLabel = ""
	)

Definition at line 5 of file SiStripClusterInfo.cc.

    : siStripNoisesToken_{iC.esConsumes<SiStripNoises, SiStripNoisesRcd>()},
      siStripGainToken_{iC.esConsumes<SiStripGain, SiStripGainRcd>()},
      siStripQualityToken_{iC.esConsumes<SiStripQuality, SiStripQualityRcd>(edm::ESInputTag("", qualityLabel))} {}

Member Function Documentation

baryStrip()

float SiStripClusterInfo::baryStrip ( ) const

inline

Definition at line 33 of file SiStripClusterInfo.h.

References SiStripCluster::barycenter(), and cluster().

Referenced by TrackerDpgAnalysis::analyze(), SiStripMonitorTrack::clusterInfos(), and ApeEstimator::fillHitVariables().

{ return cluster()->barycenter(); }

calculate_noise()

float SiStripClusterInfo::calculate_noise ( const std::vector< float > &  noise ) const

private

Definition at line 82 of file SiStripClusterInfo.cc.

References mps_fire::i, noise(), mathSSE::sqrt(), stripCharges(), and width().

Referenced by noise(), and noiseRescaledByGain().

                                                                             {
  float noiseSumInQuadrature = 0;
  int numberStripsOverThreshold = 0;
  for (int i = 0; i < width(); i++) {
    if (stripCharges()[i] != 0) {
      noiseSumInQuadrature += noise.at(i) * noise.at(i);
      numberStripsOverThreshold++;
    }
  }
  return std::sqrt(noiseSumInQuadrature / numberStripsOverThreshold);
}

charge()

uint16_t SiStripClusterInfo::charge ( void  ) const

inline

Definition at line 43 of file SiStripClusterInfo.h.

References mps_fire::end, and stripCharges().

Referenced by TrackerDpgAnalysis::analyze(), SiStripMonitorTrack::clusterInfos(), ApeEstimator::fillHitVariables(), StandaloneTrackMonitor::processClusters(), StandaloneTrackMonitor::processHit(), and signalOverNoise().

{ return std::accumulate(stripCharges().begin(), stripCharges().end(), uint16_t(0)); }

chargeLR()

std::pair< uint16_t, uint16_t > SiStripClusterInfo::chargeLR

(

)

const

Definition at line 21 of file SiStripClusterInfo.cc.

References mps_fire::end, SiStripPI::max, and stripCharges().

Referenced by ApeEstimator::fillHitVariables().

                                                               {
  std::vector<uint8_t>::const_iterator begin(stripCharges().begin()), end(stripCharges().end()), max;
  max = max_element(begin, end);
  return std::make_pair(accumulate(begin, max, uint16_t(0)), accumulate(max + 1, end, uint16_t(0)));
}

cluster()

const SiStripCluster* SiStripClusterInfo::cluster ( ) const

inline

Definition at line 28 of file SiStripClusterInfo.h.

References cluster_ptr.

Referenced by baryStrip(), firstStrip(), setCluster(), stripCharges(), and width().

{ return cluster_ptr; }

detId()

uint32_t SiStripClusterInfo::detId ( ) const

inline

Definition at line 30 of file SiStripClusterInfo.h.

References detId_.

Referenced by SiStripMonitorTrack::clusterInfos(), and setCluster().

{ return detId_; }

firstStrip()

uint16_t SiStripClusterInfo::firstStrip ( ) const

inline

Definition at line 32 of file SiStripClusterInfo.h.

References cluster(), and SiStripCluster::firstStrip().

Referenced by SiStripMonitorTrack::clusterInfos(), IsApvBad(), IsFiberBad(), maxStrip(), stripGains(), stripNoises(), stripNoisesRescaledByGain(), and stripQualitiesBad().

{ return cluster()->firstStrip(); }

initEvent()

void SiStripClusterInfo::initEvent

(

const edm::EventSetup &

iSetup

)

Definition at line 10 of file SiStripClusterInfo.cc.

References edm::EventSetup::getData(), siStripGain_, siStripGainToken_, siStripNoises_, siStripNoisesToken_, siStripQuality_, and siStripQualityToken_.

Referenced by StandaloneTrackMonitor::analyze(), HitEff::analyze(), SiStripMonitorTrack::analyze(), TrackerDpgAnalysis::analyze(), ApeEstimator::analyze(), ShallowClustersProducer::produce(), and reco::modules::TrackerTrackHitFilter::produce().

                                                              {
  siStripNoises_ = &iSetup.getData(siStripNoisesToken_);
  siStripGain_ = &iSetup.getData(siStripGainToken_);
  siStripQuality_ = &iSetup.getData(siStripQualityToken_);
}

IsAnythingBad()

bool SiStripClusterInfo::IsAnythingBad

(

)

const

Definition at line 94 of file SiStripClusterInfo.cc.

References IsApvBad(), IsFiberBad(), IsModuleBad(), and stripQualitiesBad().

                                             {
  std::vector<bool> stripBad = stripQualitiesBad();
  return IsApvBad() || IsFiberBad() || IsModuleBad() ||
         accumulate(stripBad.begin(), stripBad.end(), false, std::logical_or<bool>());
}

IsApvBad()

bool SiStripClusterInfo::IsApvBad

(

)

const

Definition at line 100 of file SiStripClusterInfo.cc.

References detId_, firstStrip(), SiStripQuality::IsApvBad(), siStripQuality_, and width().

Referenced by IsAnythingBad().

                                        {
  return siStripQuality_->IsApvBad(detId_, firstStrip() / 128) ||
         siStripQuality_->IsApvBad(detId_, (firstStrip() + width()) / 128);
}

IsFiberBad()

bool SiStripClusterInfo::IsFiberBad

(

)

const

Definition at line 105 of file SiStripClusterInfo.cc.

References detId_, firstStrip(), SiStripQuality::IsFiberBad(), siStripQuality_, and width().

Referenced by IsAnythingBad().

                                          {
  return siStripQuality_->IsFiberBad(detId_, firstStrip() / 256) ||
         siStripQuality_->IsFiberBad(detId_, (firstStrip() + width()) / 256);
}

IsModuleBad()

bool SiStripClusterInfo::IsModuleBad

(

)

const

Definition at line 110 of file SiStripClusterInfo.cc.

References detId_, SiStripQuality::IsModuleBad(), and siStripQuality_.

Referenced by IsAnythingBad().

{ return siStripQuality_->IsModuleBad(detId_); }

IsModuleUsable()

bool SiStripClusterInfo::IsModuleUsable

(

)

const

Definition at line 112 of file SiStripClusterInfo.cc.

References detId_, SiStripQuality::IsModuleUsable(), and siStripQuality_.

Referenced by ApeEstimator::fillHitVariables().

{ return siStripQuality_->IsModuleUsable(detId_); }

maxCharge()

uint8_t SiStripClusterInfo::maxCharge ( ) const

inline

Definition at line 44 of file SiStripClusterInfo.h.

References mps_fire::end, and stripCharges().

Referenced by TrackerDpgAnalysis::analyze(), SiStripMonitorTrack::clusterInfos(), and ApeEstimator::fillHitVariables().

{ return *std::max_element(stripCharges().begin(), stripCharges().end()); }

maxIndex()

uint16_t SiStripClusterInfo::maxIndex ( ) const

inline

Definition at line 45 of file SiStripClusterInfo.h.

References mps_fire::end, and stripCharges().

Referenced by SiStripMonitorTrack::clusterInfos(), ApeEstimator::fillHitVariables(), and maxStrip().

                            {
    return std::max_element(stripCharges().begin(), stripCharges().end()) - stripCharges().begin();
  }

maxStrip()

uint16_t SiStripClusterInfo::maxStrip ( ) const

inline

Definition at line 34 of file SiStripClusterInfo.h.

References firstStrip(), and maxIndex().

Referenced by ApeEstimator::fillHitVariables().

{ return firstStrip() + maxIndex(); }

noise()

float SiStripClusterInfo::noise ( ) const

inline

Definition at line 50 of file SiStripClusterInfo.h.

References calculate_noise(), and stripNoises().

Referenced by TrackerDpgAnalysis::analyze(), and calculate_noise().

{ return calculate_noise(stripNoises()); }

noiseRescaledByGain()

float SiStripClusterInfo::noiseRescaledByGain ( ) const

inline

Definition at line 51 of file SiStripClusterInfo.h.

References calculate_noise(), and stripNoisesRescaledByGain().

Referenced by TrackerDpgAnalysis::analyze(), SiStripMonitorTrack::clusterInfos(), and signalOverNoise().

{ return calculate_noise(stripNoisesRescaledByGain()); }

setCluster()

void SiStripClusterInfo::setCluster	(	const SiStripCluster &	cluster,
		int	detId
	)

Definition at line 16 of file SiStripClusterInfo.cc.

References cluster(), cluster_ptr, detId(), and detId_.

Referenced by SiStripMonitorTrack::AllClusters(), HitEff::analyze(), TrackerDpgAnalysis::analyze(), reco::modules::TrackerTrackHitFilter::checkStoN(), SiStripMonitorTrack::fillControlViewHistos(), ApeEstimator::fillHitVariables(), SiStripPositionCorrectionsTableProducer::fillTable(), StandaloneTrackMonitor::processClusters(), StandaloneTrackMonitor::processHit(), ShallowClustersProducer::produce(), and SiStripMonitorTrack::RecHitInfo().

                                                                            {
  cluster_ptr = &cluster;
  detId_ = detId;
}

signalOverNoise()

float SiStripClusterInfo::signalOverNoise ( ) const

inline

Definition at line 53 of file SiStripClusterInfo.h.

References charge(), and noiseRescaledByGain().

Referenced by HitEff::analyze(), TrackerDpgAnalysis::analyze(), reco::modules::TrackerTrackHitFilter::checkStoN(), SiStripMonitorTrack::clusterInfos(), SiStripMonitorTrack::fillControlViewHistos(), and ApeEstimator::fillHitVariables().

{ return charge() / noiseRescaledByGain(); }

siStripGain()

const SiStripGain* SiStripClusterInfo::siStripGain ( ) const

inline

Definition at line 61 of file SiStripClusterInfo.h.

References siStripGain_.

Referenced by SiStripMonitorTrack::AllClusters(), and SiStripMonitorTrack::RecHitInfo().

{ return siStripGain_; }

siStripQuality()

const SiStripQuality* SiStripClusterInfo::siStripQuality ( ) const

inline

Definition at line 62 of file SiStripClusterInfo.h.

References siStripQuality_.

Referenced by SiStripMonitorTrack::AllClusters(), and SiStripMonitorTrack::RecHitInfo().

{ return siStripQuality_; }

stripCharges()

auto SiStripClusterInfo::stripCharges ( ) const -> decltype(cluster()->amplitudes())

inline

Definition at line 37 of file SiStripClusterInfo.h.

References SiStripCluster::amplitudes(), and cluster().

Referenced by calculate_noise(), charge(), chargeLR(), SiStripMonitorTrack::clusterInfos(), ApeEstimator::fillHitVariables(), maxCharge(), maxIndex(), and variance().

{ return cluster()->amplitudes(); }

stripGains()

std::vector< float > SiStripClusterInfo::stripGains

(

)

const

Definition at line 62 of file SiStripClusterInfo.cc.

References detId_, firstStrip(), ecalph2::gains, SiStripGain::getRange(), SiStripGain::getStripGain(), mps_fire::i, siStripGain_, and width().

                                                      {
  SiStripApvGain::Range detGainRange = siStripGain_->getRange(detId_);

  std::vector<float> gains;
  gains.reserve(width());
  for (size_t i = 0; i < width(); i++) {
    gains.push_back(siStripGain_->getStripGain(firstStrip() + i, detGainRange));
  }
  return gains;
}

stripNoises()

std::vector< float > SiStripClusterInfo::stripNoises

(

)

const

Definition at line 51 of file SiStripClusterInfo.cc.

References detId_, firstStrip(), SiStripNoises::getNoise(), SiStripNoises::getRange(), mps_fire::i, hgcalLayerClusters_cfi::noises, siStripNoises_, and width().

Referenced by noise().

                                                       {
  SiStripNoises::Range detNoiseRange = siStripNoises_->getRange(detId_);

  std::vector<float> noises;
  noises.reserve(width());
  for (size_t i = 0; i < width(); i++) {
    noises.push_back(siStripNoises_->getNoise(firstStrip() + i, detNoiseRange));
  }
  return noises;
}

stripNoisesRescaledByGain()

std::vector< float > SiStripClusterInfo::stripNoisesRescaledByGain

(

)

const

Definition at line 38 of file SiStripClusterInfo.cc.

References detId_, MillePedeFileConverter_cfg::e, firstStrip(), SiStripNoises::getNoise(), SiStripNoises::getRange(), SiStripGain::getRange(), SiStripGain::getStripGain(), mps_fire::i, bookConverter::results, siStripGain_, siStripNoises_, and width().

Referenced by noiseRescaledByGain().

                                                                     {
  SiStripNoises::Range detNoiseRange = siStripNoises_->getRange(detId_);
  SiStripApvGain::Range detGainRange = siStripGain_->getRange(detId_);

  std::vector<float> results;
  results.reserve(width());
  for (size_t i = 0, e = width(); i < e; i++) {
    results.push_back(siStripNoises_->getNoise(firstStrip() + i, detNoiseRange) /
                      siStripGain_->getStripGain(firstStrip() + i, detGainRange));
  }
  return results;
}

stripQualitiesBad()

std::vector< bool > SiStripClusterInfo::stripQualitiesBad

(

)

const

Definition at line 73 of file SiStripClusterInfo.cc.

References detId_, firstStrip(), mps_fire::i, SiStripQuality::IsStripBad(), siStripQuality_, and width().

Referenced by IsAnythingBad().

                                                            {
  std::vector<bool> isBad;
  isBad.reserve(width());
  for (int i = 0; i < width(); i++) {
    isBad.push_back(siStripQuality_->IsStripBad(detId_, firstStrip() + i));
  }
  return isBad;
}

variance()

float SiStripClusterInfo::variance

(

)

const

Definition at line 27 of file SiStripClusterInfo.cc.

References mps_fire::end, mps_fire::i, submitPVResolutionJobs::q, stripCharges(), testProducerWithPsetDescEmpty_cfi::x1, and testProducerWithPsetDescEmpty_cfi::x2.

Referenced by SiStripPositionCorrectionsTableProducer::fillTable().

                                         {
  float q(0), x1(0), x2(0);
  for (auto begin(stripCharges().begin()), end(stripCharges().end()), it(begin); it != end; ++it) {
    unsigned i = it - begin;
    q += (*it);
    x1 += (*it) * (i + 0.5);
    x2 += (*it) * (i * i + i + 1. / 3);
  }
  return (x2 - x1 * x1 / q) / q;
}

width()

uint16_t SiStripClusterInfo::width ( ) const

inline

Definition at line 31 of file SiStripClusterInfo.h.

References SiStripCluster::amplitudes(), cluster(), and SiStripCluster::size().

Referenced by TrackerDpgAnalysis::analyze(), calculate_noise(), SiStripMonitorTrack::clusterInfos(), ApeEstimator::fillHitVariables(), IsApvBad(), IsFiberBad(), StandaloneTrackMonitor::processClusters(), StandaloneTrackMonitor::processHit(), stripGains(), stripNoises(), stripNoisesRescaledByGain(), and stripQualitiesBad().

{ return cluster()->amplitudes().size(); }

Member Data Documentation

cluster_ptr

const SiStripCluster* SiStripClusterInfo::cluster_ptr = nullptr

private

Definition at line 67 of file SiStripClusterInfo.h.

Referenced by cluster(), and setCluster().

detId_

uint32_t SiStripClusterInfo::detId_ = 0

private

Definition at line 77 of file SiStripClusterInfo.h.

Referenced by detId(), IsApvBad(), IsFiberBad(), IsModuleBad(), IsModuleUsable(), setCluster(), stripGains(), stripNoises(), stripNoisesRescaledByGain(), and stripQualitiesBad().

siStripGain_

const SiStripGain* SiStripClusterInfo::siStripGain_ = nullptr

private

Definition at line 74 of file SiStripClusterInfo.h.

Referenced by initEvent(), siStripGain(), stripGains(), and stripNoisesRescaledByGain().

siStripGainToken_

edm::ESGetToken<SiStripGain, SiStripGainRcd> SiStripClusterInfo::siStripGainToken_

private

Definition at line 70 of file SiStripClusterInfo.h.

Referenced by initEvent().

siStripNoises_

const SiStripNoises* SiStripClusterInfo::siStripNoises_ = nullptr

private

Definition at line 73 of file SiStripClusterInfo.h.

Referenced by initEvent(), stripNoises(), and stripNoisesRescaledByGain().

siStripNoisesToken_

edm::ESGetToken<SiStripNoises, SiStripNoisesRcd> SiStripClusterInfo::siStripNoisesToken_

private

Definition at line 69 of file SiStripClusterInfo.h.

Referenced by initEvent().

siStripQuality_

const SiStripQuality* SiStripClusterInfo::siStripQuality_ = nullptr

private

Definition at line 75 of file SiStripClusterInfo.h.

Referenced by initEvent(), IsApvBad(), IsFiberBad(), IsModuleBad(), IsModuleUsable(), siStripQuality(), and stripQualitiesBad().

siStripQualityToken_

edm::ESGetToken<SiStripQuality, SiStripQualityRcd> SiStripClusterInfo::siStripQualityToken_

private

Definition at line 71 of file SiStripClusterInfo.h.

Referenced by initEvent().