jptJetAnalysis::StripSignalOverNoiseCalculator Class Reference

Public Member Functions
double	operator() (const SiStripCluster &cluster) const
double	signalOverNoise (const SiStripCluster &cluster) const
	StripSignalOverNoiseCalculator (const std::string &theQualityLabel=std::string(""))
void	update (const edm::EventSetup &eventSetup)

Private Attributes
const SiStripGain *	gain_
uint32_t	gainCacheId_
const SiStripNoises *	noise_
uint32_t	noiseCacheId_
const SiStripQuality *	quality_
uint32_t	qualityCacheId_
const std::string	qualityLabel_

Detailed Description

Definition at line 55 of file JPTJetAnalyzer.cc.

Constructor & Destructor Documentation

jptJetAnalysis::StripSignalOverNoiseCalculator::StripSignalOverNoiseCalculator

(

const std::string &

theQualityLabel = std::string("")

)

Definition at line 737 of file JPTJetAnalyzer.cc.

    : qualityLabel_(theQualityLabel),
      quality_(NULL),
      noise_(NULL),
      gain_(NULL),
      qualityCacheId_(0),
      noiseCacheId_(0),
      gainCacheId_(0)
    {}

Member Function Documentation

double jptJetAnalysis::StripSignalOverNoiseCalculator::operator() ( const SiStripCluster &  cluster ) const

inline

Definition at line 61 of file JPTJetAnalyzer.cc.

References signalOverNoise().

      { return signalOverNoise(cluster); }

double jptJetAnalysis::StripSignalOverNoiseCalculator::signalOverNoise

(

const SiStripCluster &

cluster

)

const

Definition at line 778 of file JPTJetAnalyzer.cc.

References ecalMGPA::adc(), SiStripCluster::amplitudes(), SiStripCluster::firstStrip(), gain_, SiStripCluster::geographicalId(), SiStripNoises::getNoise(), SiStripNoises::getRange(), SiStripBadStrip::getRange(), SiStripGain::getRange(), SiStripGain::getStripGain(), SiStripQuality::IsStripBad(), noise_, quality_, mathSSE::sqrt(), and strip().

Referenced by operator()().

  {
    const uint32_t detId = cluster.geographicalId();
    const uint16_t firstStrip = cluster.firstStrip();
    const SiStripQuality::Range& qualityRange = quality_->getRange(detId);
    const SiStripNoises::Range& noiseRange = noise_->getRange(detId);
    const SiStripApvGain::Range& gainRange = gain_->getRange(detId);
    double signal = 0;
    double noise2 = 0;
    unsigned int nNonZeroStrips = 0;
    const std::vector<uint8_t>& clusterAmplitudes = cluster.amplitudes();
    const std::vector<uint8_t>::const_iterator clusterAmplitudesEnd = clusterAmplitudes.end();
    const std::vector<uint8_t>::const_iterator clusterAmplitudesBegin = clusterAmplitudes.begin();
    for (std::vector<uint8_t>::const_iterator iAmp = clusterAmplitudesBegin; iAmp != clusterAmplitudesEnd; ++iAmp) {
      const uint8_t adc = *iAmp;
      const uint16_t strip = iAmp-clusterAmplitudesBegin+firstStrip;
      const bool stripBad = quality_->IsStripBad(qualityRange,strip);
      const double noise = noise_->getNoise(strip,noiseRange);
      const double gain = gain_->getStripGain(strip,gainRange);
      signal += adc;
      if (adc) ++nNonZeroStrips;
      const double noiseContrib = (stripBad ? 0 : noise/gain);
      noise2 += noiseContrib*noiseContrib;
    }
    const double noise = sqrt(noise2/nNonZeroStrips);
    if (noise) return signal/noise;
    else return 0;
  }

void jptJetAnalysis::StripSignalOverNoiseCalculator::update

(

const edm::EventSetup &

eventSetup

)

Definition at line 747 of file JPTJetAnalyzer.cc.

References edm::eventsetup::EventSetupRecord::cacheIdentifier(), gain_, gainCacheId_, edm::EventSetup::get(), edm::eventsetup::EventSetupRecord::get(), noise_, noiseCacheId_, edm::ESHandle< class >::product(), quality_, qualityCacheId_, and qualityLabel_.

Referenced by python.Vispa.Gui.VispaWidget.VispaWidget::autosize(), python.Vispa.Views.LineDecayView.LineDecayContainer::createObject(), python.Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), python.Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), python.Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), python.Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), python.Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), python.Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), python.Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), python.Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), python.Vispa.Gui.FindDialog.FindDialog::reset(), python.Vispa.Gui.PortConnection.PointToPointConnection::select(), python.Vispa.Gui.VispaWidget.VispaWidget::select(), python.Vispa.Views.LineDecayView.LineDecayContainer::select(), python.Vispa.Gui.VispaWidget.VispaWidget::setText(), python.Vispa.Gui.VispaWidget.VispaWidget::setTitle(), python.Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), python.Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and python.Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

  {
    //update the quality if necessary
    const SiStripQualityRcd& qualityRecord = eventSetup.get<SiStripQualityRcd>();
    const uint32_t newQualityCacheId = qualityRecord.cacheIdentifier();
    if ((newQualityCacheId != qualityCacheId_) || !quality_) {
      edm::ESHandle<SiStripQuality> qualityHandle;
      qualityRecord.get(qualityLabel_,qualityHandle);
      quality_ = qualityHandle.product();
      qualityCacheId_ = newQualityCacheId;
    }
    //update the noise if necessary
    const SiStripNoisesRcd& noiseRecord = eventSetup.get<SiStripNoisesRcd>();
    const uint32_t newNoiseCacheId = noiseRecord.cacheIdentifier();
    if ((newNoiseCacheId != noiseCacheId_) || !noise_) {
      edm::ESHandle<SiStripNoises> noiseHandle;
      noiseRecord.get(noiseHandle);
      noise_ = noiseHandle.product();
      noiseCacheId_ = newNoiseCacheId;
    }
    //update the gain if necessary
    const SiStripGainRcd& gainRecord = eventSetup.get<SiStripGainRcd>();
    const uint32_t newGainCacheId = gainRecord.cacheIdentifier();
    if ((newGainCacheId != gainCacheId_) || !gain_) {
      edm::ESHandle<SiStripGain> gainHandle;
      gainRecord.get(gainHandle);
      gain_ = gainHandle.product();
      gainCacheId_ = newGainCacheId;
    }
  }

Member Data Documentation

const SiStripGain* jptJetAnalysis::StripSignalOverNoiseCalculator::gain_

private

Definition at line 67 of file JPTJetAnalyzer.cc.

Referenced by signalOverNoise(), and update().

uint32_t jptJetAnalysis::StripSignalOverNoiseCalculator::gainCacheId_

private

Definition at line 70 of file JPTJetAnalyzer.cc.

Referenced by update().

const SiStripNoises* jptJetAnalysis::StripSignalOverNoiseCalculator::noise_

private

Definition at line 66 of file JPTJetAnalyzer.cc.

Referenced by signalOverNoise(), and update().

uint32_t jptJetAnalysis::StripSignalOverNoiseCalculator::noiseCacheId_

private

Definition at line 69 of file JPTJetAnalyzer.cc.

Referenced by update().

const SiStripQuality* jptJetAnalysis::StripSignalOverNoiseCalculator::quality_

private

Definition at line 65 of file JPTJetAnalyzer.cc.

Referenced by signalOverNoise(), and update().

uint32_t jptJetAnalysis::StripSignalOverNoiseCalculator::qualityCacheId_

private

Definition at line 68 of file JPTJetAnalyzer.cc.

Referenced by update().

const std::string jptJetAnalysis::StripSignalOverNoiseCalculator::qualityLabel_

private

Definition at line 64 of file JPTJetAnalyzer.cc.

Referenced by update().