SiStripGain Class Reference

#include <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, const SiStripDetInfo &detInfo)
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, const SiStripDetInfo &detInfo)
	Kept for compatibility. More...
	SiStripGain (const SiStripApvGain &apvgain, const double &factor, const std::pair< std::string, std::string > &recordLabelPair, const SiStripDetInfo &detInfo)

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, SiStripDetInfo const &detInfo, const SiStripApvGain *apvgain2=nullptr, 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

CalibFormats/SiStripObjects/interface/SiStripGain.h

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 45 of file SiStripGain.h.

Constructor & Destructor Documentation

SiStripGain() [1/4]

SiStripGain::SiStripGain ( )

inline

Definition at line 47 of file SiStripGain.h.

{}

SiStripGain() [2/4]

SiStripGain::SiStripGain ( const SiStripGain &  )

delete

SiStripGain() [3/4]

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

inline

Kept for compatibility.

Definition at line 52 of file SiStripGain.h.

References DQMScaleToClient_cfi::factor, and multiply().

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

SiStripGain() [4/4]

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

inline

Definition at line 57 of file SiStripGain.h.

References DQMScaleToClient_cfi::factor, and multiply().

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

Member Function Documentation

fillNewGain()

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

private

Definition at line 40 of file SiStripGain.cc.

References apvgain_, apvgainAutoPtr_, DQMScaleToClient_cfi::factor, SiStripDetInfo::getAllData(), SiStripApvGain::getApvGain(), SiStripApvGain::getDetIds(), SiStripApvGain::getRange(), pileupCalc::inputRange, SiStripApvGain::put(), and FastTimerService_cff::range.

Referenced by multiply().

                                                     {
  SiStripApvGain *newApvGain = new SiStripApvGain;
  const auto &DetInfos = detInfo.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) {
    auto 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);
}

getApvGain() [1/2]

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

inlinestatic

Definition at line 80 of file SiStripGain.h.

References SiStripApvGain::getApvGain(), and FastTimerService_cff::range.

Referenced by getApvGain(), printDebug(), printSummary(), and ShallowGainCalibration::produce().

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

getApvGain() [2/2]

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

Definition at line 94 of file SiStripGain.cc.

References apvgainVector_, getApvGain(), normVector_, and FastTimerService_cff::range.

                                                                                                               {
  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.;
}

getDetIds()

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 102 of file SiStripGain.cc.

References apvgain_, and SiStripApvGain::getDetIds().

Referenced by printDebug(), and printSummary().

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

getLabelName()

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

inline

Definition at line 101 of file SiStripGain.h.

References recordLabelPair_.

{ return recordLabelPair_[index].second; }

getNumberOfTags()

size_t SiStripGain::getNumberOfTags ( ) const

inline

Definition at line 99 of file SiStripGain.h.

References apvgainVector_.

{ return apvgainVector_.size(); }

getRange() [1/2]

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

inline

Definition at line 75 of file SiStripGain.h.

References apvgain_, and SiStripApvGain::getRange().

Referenced by SiStripMonitorCluster::analyze(), SiStripMonitorTrack::clusterInfos(), printDebug(), printSummary(), ShallowGainCalibration::produce(), DigiSimLinkAlgorithm::run(), SiStripClusterToDigiProducer::setDetId(), SiStripClusterInfo::stripGains(), and SiStripClusterInfo::stripNoisesRescaledByGain().

{ return apvgain_->getRange(detID); }

getRange() [2/2]

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 107 of file SiStripGain.cc.

References apvgainVector_.

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

getRangeByPos()

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

inline

Definition at line 76 of file SiStripGain.h.

References apvgain_, and SiStripApvGain::getRangeByPos().

{ return apvgain_->getRangeByPos(pos); }

getRcdName()

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

inline

Definition at line 100 of file SiStripGain.h.

References recordLabelPair_.

{ return recordLabelPair_[index].first; }

getStripGain() [1/2]

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

inlinestatic

Definition at line 77 of file SiStripGain.h.

References SiStripApvGain::getStripGain(), FastTimerService_cff::range, and digitizers_cfi::strip.

Referenced by SiStripMonitorCluster::analyze(), SiStripMonitorTrack::clusterInfos(), SiStripClusterToDigiProducer::gain(), getStripGain(), DigiSimLinkAlgorithm::run(), SiStripClusterInfo::stripGains(), and SiStripClusterInfo::stripNoisesRescaledByGain().

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

getStripGain() [2/2]

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

Definition at line 86 of file SiStripGain.cc.

References apvgainVector_, getStripGain(), FastTimerService_cff::range, and digitizers_cfi::strip.

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

getTagNorm()

double SiStripGain::getTagNorm ( const uint32_t  index ) const

inline

Definition at line 102 of file SiStripGain.h.

References normVector_.

{ return normVector_[index]; }

multiply()

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

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_, DQMScaleToClient_cfi::factor, fillNewGain(), normVector_, and recordLabelPair_.

Referenced by SiStripGain().

                                                          {
  // When inserting the first ApvGain
  if (apvgain_ == nullptr) {
    if ((factor != 1) && (factor != 0)) {
      fillNewGain(&apvgain, factor, detInfo);
    } 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., detInfo, &apvgain, factor);
  }
  recordLabelPair_.push_back(recordLabelPair);
  apvgainVector_.push_back(&apvgain);
  normVector_.push_back(factor);
}

operator=()

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

delete

printDebug()

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

Definition at line 111 of file SiStripGain.cc.

References getApvGain(), getDetIds(), getRange(), FastTimerService_cff::range, and contentValuesCheck::ss.

                                                                                 {
  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;
    }
  }
}

printSummary()

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

Definition at line 127 of file SiStripGain.cc.

References getApvGain(), getDetIds(), getRange(), FastTimerService_cff::range, and contentValuesCheck::ss.

                                                                                            {
  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

apvgain_

const SiStripApvGain* SiStripGain::apvgain_

private

Definition at line 118 of file SiStripGain.h.

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

apvgainAutoPtr_

std::unique_ptr<SiStripApvGain> SiStripGain::apvgainAutoPtr_

private

Definition at line 119 of file SiStripGain.h.

Referenced by fillNewGain().

apvgainVector_

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

private

Definition at line 116 of file SiStripGain.h.

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

normVector_

std::vector<double> SiStripGain::normVector_

private

Definition at line 117 of file SiStripGain.h.

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

recordLabelPair_

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

private

Definition at line 120 of file SiStripGain.h.

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