SiStripGain Class Reference

#include <CalibFormats/SiStripObjects/interface/SiStripGain.h>

Public Member Functions
float	getApvGain (const uint16_t &apv, const SiStripApvGain::Range &range, const uint32_t index) const
void	getDetIds (std::vector< uint32_t > &DetIds_) const
	ATTENTION: we assume the detIds are the same as those from the first gain. More...
std::string	getLabelName (const uint32_t index) const
size_t	getNumberOfTags () const
const SiStripApvGain::Range	getRange (uint32_t detID) const
const SiStripApvGain::Range	getRange (const uint32_t &detID, const uint32_t index) const
SiStripApvGain::Range	getRangeByPos (unsigned short pos) const
std::string	getRcdName (const uint32_t index) const
float	getStripGain (const uint16_t &strip, const SiStripApvGain::Range &range, const uint32_t index) const
double	getTagNorm (const uint32_t index) const
void	multiply (const SiStripApvGain &apvgain, const double &factor, const std::pair< std::string, std::string > &recordLabelPair)
	Used to input additional gain values that will be multiplied to the first one. More...
const SiStripGain &	operator= (const SiStripGain &)=delete
void	printDebug (std::stringstream &ss, const TrackerTopology *trackerTopo) const
void	printSummary (std::stringstream &ss, const TrackerTopology *trackerTopo) const
	SiStripGain ()
	SiStripGain (const SiStripGain &)=delete
	SiStripGain (const SiStripApvGain &apvgain, const double &factor)
	Kept for compatibility. More...
	SiStripGain (const SiStripApvGain &apvgain, const double &factor, const std::pair< std::string, std::string > &recordLabelPair)

Static Public Member Functions
static float	getApvGain (const uint16_t &apv, const SiStripApvGain::Range &range)
static float	getStripGain (const uint16_t &strip, const SiStripApvGain::Range &range)

Private Member Functions
void	fillNewGain (const SiStripApvGain *apvgain, const double &factor, const SiStripApvGain *apvgain2=0, const double &factor2=1.)

Private Attributes
const SiStripApvGain *	apvgain_
std::unique_ptr< SiStripApvGain >	apvgainAutoPtr_
std::vector< const SiStripApvGain * >	apvgainVector_
std::vector< double >	normVector_
std::vector< std::pair< std::string, std::string > >	recordLabelPair_

Detailed Description

Description: give detector view for the cabling classes

Usage: <usage>

Original Author: gbruno Created: Wed Mar 22 12:24:20 CET 2006

Modifications by M. De Mattia (demat.nosp@m.tia@.nosp@m.pd.in.nosp@m.fn.i.nosp@m.t) on 11/11/2009: It now holds a std::vector of pointers to ApvGain and a std::vector of corresponding normalization factors.
It returns the product of all the Gain/norm ratios.
The multiply method allows to input additional gain records.
ATTENTION: the code assumes that the second tag has at least the same DetIds that the first tag and only the DetIds present in the first tag will be used.

There are two set of methods to access the gain value. The first one returns the products of all ApvGain/norm. The second set of methods take an additional integer paramter and return the corresponding ApvGain (without normalization). Note that no check is done inside these methods to see if the ApvGain really exists. It is responsibility of the user to not pass an index value that exceeds the number of ApvGains.
The normalization factors for each of the stored ApvGains are also accessible passing the corresponding index.
Additional method are provided to get the number of ApvGains used to build this object, the names of the records that stored those ApvGains and the labels (they can be used to go back to the tags looking in the cfg).

Definition at line 42 of file SiStripGain.h.

Constructor & Destructor Documentation

SiStripGain::SiStripGain ( )

inline

Definition at line 45 of file SiStripGain.h.

References operator=().

{}

SiStripGain::SiStripGain ( const SiStripGain &  )

delete

SiStripGain::SiStripGain ( const SiStripApvGain &  apvgain, const double &  factor  )

inline

Kept for compatibility.

Definition at line 50 of file SiStripGain.h.

References multiply().

                                                                           :
    apvgain_(nullptr)
  {
    multiply(apvgain, factor, std::make_pair("", ""));
  }

SiStripGain::SiStripGain ( const SiStripApvGain &  apvgain, const double &  factor, const std::pair< std::string, std::string > &  recordLabelPair  )

inline

Definition at line 56 of file SiStripGain.h.

References multiply().

                                                                    :
    apvgain_(nullptr)
  {
    multiply(apvgain, factor, recordLabelPair);
  }

Member Function Documentation

void SiStripGain::fillNewGain ( const SiStripApvGain *  apvgain, const double &  factor, const SiStripApvGain *  apvgain2 = 0, const double &  factor2 = 1.  )

private

Definition at line 39 of file SiStripGain.cc.

References apvgain_, apvgainAutoPtr_, edm::FileInPath::fullPath(), SiStripApvGain::getApvGain(), SiStripApvGain::getDetIds(), SiStripApvGain::getRange(), estimatePileup::inputRange, SiStripApvGain::put(), and matplotRender::reader.

Referenced by getTagNorm(), and multiply().

{
  SiStripApvGain * newApvGain = new SiStripApvGain;
  edm::FileInPath fp("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat");
  SiStripDetInfoFileReader reader(fp.fullPath());
  const std::map<uint32_t, SiStripDetInfoFileReader::DetInfo>& DetInfos = reader.getAllData();

  // Loop on the apvgain in input and fill the newApvGain with the values/factor.
  std::vector<uint32_t> detIds;
  apvgain->getDetIds(detIds);
  std::vector<uint32_t>::const_iterator it = detIds.begin();
  for( ; it != detIds.end(); ++it ) {

    std::map<uint32_t, SiStripDetInfoFileReader::DetInfo>::const_iterator detInfoIt = DetInfos.find(*it);
    if( detInfoIt != DetInfos.end() ) {

      std::vector<float> theSiStripVector;

      // Loop on all the apvs and then on the strips
      SiStripApvGain::Range range = apvgain->getRange(*it);

      SiStripApvGain::Range range2;
      if( apvgain2 != nullptr ) {
    range2 = apvgain2->getRange(*it);
      }

      for( int apv = 0; apv < detInfoIt->second.nApvs; ++apv ) {
    float apvGainValue = apvgain->getApvGain( apv, range )/factor;

    if( (apvgain2 != nullptr) && (factor2 != 0.) ) {
      apvGainValue *= apvgain2->getApvGain( apv, range2 )/factor2;
    }

    theSiStripVector.push_back(apvGainValue);
      }
      SiStripApvGain::Range inputRange(theSiStripVector.begin(), theSiStripVector.end());
      if( ! newApvGain->put(*it, inputRange) ) {
    edm::LogError("SiStripGain") << "detid already exists" << std::endl;
      }
    }
  }
  apvgain_ = newApvGain;
  // Deletes the managed object and replaces it with the new one
  apvgainAutoPtr_.reset(newApvGain);
}

static float SiStripGain::getApvGain ( const uint16_t &  apv, const SiStripApvGain::Range &  range  )

inlinestatic

Definition at line 73 of file SiStripGain.h.

References SiStripApvGain::getApvGain(), getDetIds(), getRange(), getStripGain(), and digitizers_cfi::strip.

Referenced by SiStripGainFromCalibTree::algoBeginRun(), SiStripApvGainReader::analyze(), SiStripApvGainRescaler::analyze(), SiStripChannelGainFromDBMiscalibrator::analyze(), SiStripGainsPCLHarvester::beginRun(), SiStripGainsPCLWorker::dqmBeginRun(), getApvGain(), printDebug(), printSummary(), and ShallowGainCalibration::produce().

{  return SiStripApvGain::getApvGain(apv, range); }

float SiStripGain::getApvGain	(	const uint16_t &	apv,
		const SiStripApvGain::Range &	range,
		const uint32_t	index
	)		const

Definition at line 99 of file SiStripGain.cc.

References apvgainVector_, getApvGain(), and normVector_.

{
  if( !(apvgainVector_.empty()) ) {
    return (apvgainVector_[index]->getApvGain(apv, range))/(normVector_[index]);
  }
  edm::LogError("SiStripGain::getApvGain") << "ERROR: no gain available. Returning gain = 1." << std::endl;
  return 1.;
}

void SiStripGain::getDetIds

(

std::vector< uint32_t > &

DetIds_

)

const

ATTENTION: we assume the detIds are the same as those from the first gain.

Definition at line 108 of file SiStripGain.cc.

References apvgain_, and SiStripApvGain::getDetIds().

Referenced by SiStripApvGainReader::analyze(), SiStripApvGainRescaler::analyze(), SiStripChannelGainFromDBMiscalibrator::analyze(), getApvGain(), StripClusterizerAlgorithm::initialize(), printDebug(), and printSummary().

{
  // ATTENTION: we assume the detIds are the same as those from the first gain
  return apvgain_->getDetIds(DetIds_);
}

std::string SiStripGain::getLabelName ( const uint32_t  index ) const

inline

Definition at line 99 of file SiStripGain.h.

References recordLabelPair_.

  {
    return recordLabelPair_[index].second;
  }

size_t SiStripGain::getNumberOfTags ( ) const

inline

Definition at line 91 of file SiStripGain.h.

References apvgainVector_.

Referenced by SiStripGainFromCalibTree::algoBeginRun(), SiStripGainsPCLHarvester::beginRun(), and SiStripGainsPCLWorker::dqmBeginRun().

  {
    return apvgainVector_.size();
  }

const SiStripApvGain::Range SiStripGain::getRange ( uint32_t  detID ) const

inline

Definition at line 70 of file SiStripGain.h.

References apvgain_, and SiStripApvGain::getRange().

Referenced by SiStripGainFromData::algoBeginRun(), SiStripGainFromCalibTree::algoBeginRun(), SiStripApvGainReader::analyze(), SiStripMonitorCluster::analyze(), SiStripApvGainRescaler::analyze(), SiStripChannelGainFromDBMiscalibrator::analyze(), SiStripGainsPCLHarvester::beginRun(), SiStripMonitorTrack::clusterInfos(), SiStripProcessedRawDigiProducer::common_process(), SiTrivialDigitalConverter::convert(), SiTrivialDigitalConverter::convertRaw(), SiStripDigitizerAlgorithm::digitize(), SiStripGainsPCLWorker::dqmBeginRun(), StripClusterizerAlgorithm::findDetId(), getApvGain(), printDebug(), printSummary(), ShallowGainCalibration::produce(), edm::DataMixingSiStripMCDigiWorker::putSiStrip(), DigiSimLinkAlgorithm::run(), SiStripClusterToDigiProducer::setDetId(), SiStripClusterInfo::stripGains(), and SiStripClusterInfo::stripNoisesRescaledByGain().

{ return apvgain_->getRange(detID);}

const SiStripApvGain::Range SiStripGain::getRange	(	const uint32_t &	detID,
		const uint32_t	index
	)		const

The second parameter allows to specify which gain to retrieve, considering that they are in input order. NOTE that no protection is inside the method (because we want to keep it very light) therefore it is the caller duty to check that the index is in the correct range.

Definition at line 115 of file SiStripGain.cc.

References apvgainVector_.

{
  return apvgainVector_[index]->getRange(DetId);
}

SiStripApvGain::Range SiStripGain::getRangeByPos ( unsigned short  pos ) const

inline

Definition at line 71 of file SiStripGain.h.

References apvgain_, and SiStripApvGain::getRangeByPos().

Referenced by StripClusterizerAlgorithm::findDetId().

{ return apvgain_->getRangeByPos(pos);}

std::string SiStripGain::getRcdName ( const uint32_t  index ) const

inline

Definition at line 95 of file SiStripGain.h.

References recordLabelPair_.

  {
    return recordLabelPair_[index].first;
  }

static float SiStripGain::getStripGain ( const uint16_t &  strip, const SiStripApvGain::Range &  range  )

inlinestatic

Definition at line 72 of file SiStripGain.h.

References SiStripApvGain::getStripGain().

Referenced by SiStripMonitorCluster::analyze(), SiStripMonitorTrack::clusterInfos(), SiStripProcessedRawDigiProducer::common_process(), SiTrivialDigitalConverter::convert(), SiTrivialDigitalConverter::convertRaw(), SiStripDigitizerAlgorithm::digitize(), StripClusterizerAlgorithm::Det::gain(), SiStripClusterToDigiProducer::gain(), getApvGain(), getStripGain(), edm::DataMixingSiStripMCDigiWorker::putSiStrip(), DigiSimLinkAlgorithm::run(), SiStripClusterInfo::stripGains(), and SiStripClusterInfo::stripNoisesRescaledByGain().

{ return  SiStripApvGain::getStripGain(strip, range);}

float SiStripGain::getStripGain	(	const uint16_t &	strip,
		const SiStripApvGain::Range &	range,
		const uint32_t	index
	)		const

Definition at line 88 of file SiStripGain.cc.

References apvgainVector_, and getStripGain().

{
  if( !(apvgainVector_.empty()) ) {
    return( apvgainVector_[index]->getStripGain(strip, range) );
  }
  edm::LogError("SiStripGain::getStripGain") << "ERROR: no gain available. Returning gain = 1." << std::endl;
  return 1.;
}

double SiStripGain::getTagNorm ( const uint32_t  index ) const

inline

Definition at line 103 of file SiStripGain.h.

References fillNewGain(), normVector_, printDebug(), and printSummary().

  {
    return normVector_[index];
  }

void SiStripGain::multiply	(	const SiStripApvGain &	apvgain,
		const double &	factor,
		const std::pair< std::string, std::string > &	recordLabelPair
	)

Used to input additional gain values that will be multiplied to the first one.

Definition at line 17 of file SiStripGain.cc.

References apvgain_, apvgainVector_, fillNewGain(), normVector_, and recordLabelPair_.

Referenced by SiStripGain(), and SiStripGainESProducerTemplate< TDependentRecord, TInputRecord >::SiStripGainNormalizationFunction().

{
  // When inserting the first ApvGain
  if( apvgain_ == nullptr ) {
    if( (factor != 1) && (factor != 0) ) {
      fillNewGain( &apvgain, factor );
    }
    else {
      // If the normalization factor is one, no need to create a new SiStripApvGain
      apvgain_ = &apvgain;
    }
  }
  else {
    // There is already an ApvGain inside the SiStripGain. Multiply it by the new one and save the new pointer.
    fillNewGain( apvgain_, 1., &apvgain, factor ); 
  }
  recordLabelPair_.push_back(recordLabelPair);
  apvgainVector_.push_back(&apvgain);
  normVector_.push_back(factor);
}

const SiStripGain& SiStripGain::operator= ( const SiStripGain &  )

delete

Referenced by SiStripGain().

void SiStripGain::printDebug	(	std::stringstream &	ss,
		const TrackerTopology *	trackerTopo
	)		const

Definition at line 120 of file SiStripGain.cc.

References getApvGain(), getDetIds(), and getRange().

Referenced by getTagNorm().

{
  std::vector<unsigned int> detIds;
  getDetIds(detIds);
  std::vector<unsigned int>::const_iterator detid = detIds.begin();
  ss << "Number of detids " << detIds.size() << std::endl;

  for( ; detid != detIds.end(); ++detid ) {
    SiStripApvGain::Range range = getRange(*detid);
    int apv=0;
    for( int it=0; it < range.second - range.first; ++it ) {
      ss << "detid " << *detid << " \t"
     << " apv " << apv++ << " \t"
     << getApvGain(it,range) << " \t"
     << std::endl;
    }
  }
}

void SiStripGain::printSummary	(	std::stringstream &	ss,
		const TrackerTopology *	trackerTopo
	)		const

Definition at line 139 of file SiStripGain.cc.

References getApvGain(), getDetIds(), and getRange().

Referenced by getTagNorm().

{
  SiStripDetSummary summaryGain{trackerTopo};

  std::vector<unsigned int> detIds;
  getDetIds(detIds);
  std::vector<uint32_t>::const_iterator detid = detIds.begin();
  for( ; detid != detIds.end(); ++detid ) {
    SiStripApvGain::Range range = getRange(*detid);
    for( int it=0; it < range.second - range.first; ++it ) {
      summaryGain.add(*detid, getApvGain(it, range));
    }
  }
  ss << "Summary of gain values:" << std::endl;
  summaryGain.print(ss, true);
}

Member Data Documentation

const SiStripApvGain* SiStripGain::apvgain_

private

Definition at line 120 of file SiStripGain.h.

Referenced by fillNewGain(), getDetIds(), getRange(), getRangeByPos(), and multiply().

std::unique_ptr<SiStripApvGain> SiStripGain::apvgainAutoPtr_

private

Definition at line 121 of file SiStripGain.h.

Referenced by fillNewGain().

std::vector<const SiStripApvGain *> SiStripGain::apvgainVector_

private

Definition at line 118 of file SiStripGain.h.

Referenced by getApvGain(), getNumberOfTags(), getRange(), getStripGain(), and multiply().

std::vector<double> SiStripGain::normVector_

private

Definition at line 119 of file SiStripGain.h.

Referenced by getApvGain(), getTagNorm(), and multiply().

std::vector<std::pair<std::string, std::string> > SiStripGain::recordLabelPair_

private

Definition at line 122 of file SiStripGain.h.

Referenced by getLabelName(), getRcdName(), and multiply().