CMS 3D CMS Logo

List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
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 (const uint32_t &detID, const uint32_t index) const
 
const SiStripApvGain::Range getRange (uint32_t detID) 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 SiStripGainoperator= (const SiStripGain &)=delete
 
void printDebug (std::stringstream &ss, const TrackerTopology *trackerTopo) const
 
void printSummary (std::stringstream &ss, const TrackerTopology *trackerTopo) const
 
 SiStripGain ()
 
 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)
 
 SiStripGain (const SiStripGain &)=delete
 

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=nullptr, const double &factor2=1.)
 

Private Attributes

const SiStripApvGainapvgain_
 
std::unique_ptr< SiStripApvGainapvgainAutoPtr_
 
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 44 of file SiStripGain.h.

Constructor & Destructor Documentation

◆ SiStripGain() [1/4]

SiStripGain::SiStripGain ( )
inline

Definition at line 46 of file SiStripGain.h.

46 {}

◆ SiStripGain() [2/4]

SiStripGain::SiStripGain ( const SiStripGain )
delete

◆ SiStripGain() [3/4]

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

Kept for compatibility.

Definition at line 51 of file SiStripGain.h.

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

References DQMScaleToClient_cfi::factor, and multiply().

◆ SiStripGain() [4/4]

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

Definition at line 55 of file SiStripGain.h.

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

References DQMScaleToClient_cfi::factor, and multiply().

Member Function Documentation

◆ fillNewGain()

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

Definition at line 39 of file SiStripGain.cc.

42  {
43  SiStripApvGain *newApvGain = new SiStripApvGain;
44  edm::FileInPath fp("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat");
46  const std::map<uint32_t, SiStripDetInfoFileReader::DetInfo> &DetInfos = reader.getAllData();
47 
48  // Loop on the apvgain in input and fill the newApvGain with the
49  // values/factor.
50  std::vector<uint32_t> detIds;
51  apvgain->getDetIds(detIds);
52  std::vector<uint32_t>::const_iterator it = detIds.begin();
53  for (; it != detIds.end(); ++it) {
54  std::map<uint32_t, SiStripDetInfoFileReader::DetInfo>::const_iterator detInfoIt = DetInfos.find(*it);
55  if (detInfoIt != DetInfos.end()) {
56  std::vector<float> theSiStripVector;
57 
58  // Loop on all the apvs and then on the strips
59  SiStripApvGain::Range range = apvgain->getRange(*it);
60 
61  SiStripApvGain::Range range2;
62  if (apvgain2 != nullptr) {
63  range2 = apvgain2->getRange(*it);
64  }
65 
66  for (int apv = 0; apv < detInfoIt->second.nApvs; ++apv) {
67  float apvGainValue = apvgain->getApvGain(apv, range) / factor;
68 
69  if ((apvgain2 != nullptr) && (factor2 != 0.)) {
70  apvGainValue *= apvgain2->getApvGain(apv, range2) / factor2;
71  }
72 
73  theSiStripVector.push_back(apvGainValue);
74  }
75  SiStripApvGain::Range inputRange(theSiStripVector.begin(), theSiStripVector.end());
76  if (!newApvGain->put(*it, inputRange)) {
77  edm::LogError("SiStripGain") << "detid already exists" << std::endl;
78  }
79  }
80  }
81  apvgain_ = newApvGain;
82  // Deletes the managed object and replaces it with the new one
83  apvgainAutoPtr_.reset(newApvGain);
84 }

References apvgain_, apvgainAutoPtr_, DQMScaleToClient_cfi::factor, personalPlayback::fp, SiStripApvGain::getApvGain(), SiStripApvGain::getDetIds(), SiStripApvGain::getRange(), pileupCalc::inputRange, SiStripApvGain::put(), FastTimerService_cff::range, and DQM::reader.

Referenced by multiply().

◆ getApvGain() [1/2]

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

◆ 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.

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

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

◆ 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.

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

References apvgain_, and SiStripApvGain::getDetIds().

Referenced by printDebug(), and printSummary().

◆ getLabelName()

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

Definition at line 97 of file SiStripGain.h.

97 { return recordLabelPair_[index].second; }

References recordLabelPair_.

◆ getNumberOfTags()

size_t SiStripGain::getNumberOfTags ( ) const
inline

Definition at line 95 of file SiStripGain.h.

95 { return apvgainVector_.size(); }

References apvgainVector_.

◆ getRange() [1/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.

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

References apvgainVector_.

◆ getRange() [2/2]

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

◆ getRangeByPos()

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

Definition at line 72 of file SiStripGain.h.

72 { return apvgain_->getRangeByPos(pos); }

References apvgain_, and SiStripApvGain::getRangeByPos().

◆ getRcdName()

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

Definition at line 96 of file SiStripGain.h.

96 { return recordLabelPair_[index].first; }

References recordLabelPair_.

◆ getStripGain() [1/2]

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

◆ 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.

86  {
87  if (!(apvgainVector_.empty())) {
89  }
90  edm::LogError("SiStripGain::getStripGain") << "ERROR: no gain available. Returning gain = 1." << std::endl;
91  return 1.;
92 }

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

◆ getTagNorm()

double SiStripGain::getTagNorm ( const uint32_t  index) const
inline

Definition at line 98 of file SiStripGain.h.

98 { return normVector_[index]; }

References normVector_.

◆ multiply()

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.

19  {
20  // When inserting the first ApvGain
21  if (apvgain_ == nullptr) {
22  if ((factor != 1) && (factor != 0)) {
23  fillNewGain(&apvgain, factor);
24  } else {
25  // If the normalization factor is one, no need to create a new
26  // SiStripApvGain
27  apvgain_ = &apvgain;
28  }
29  } else {
30  // There is already an ApvGain inside the SiStripGain. Multiply it by the
31  // new one and save the new pointer.
32  fillNewGain(apvgain_, 1., &apvgain, factor);
33  }
34  recordLabelPair_.push_back(recordLabelPair);
35  apvgainVector_.push_back(&apvgain);
36  normVector_.push_back(factor);
37 }

References apvgain_, apvgainVector_, DQMScaleToClient_cfi::factor, fillNewGain(), normVector_, and recordLabelPair_.

Referenced by SiStripGain().

◆ 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.

111  {
112  std::vector<unsigned int> detIds;
113  getDetIds(detIds);
114  std::vector<unsigned int>::const_iterator detid = detIds.begin();
115  ss << "Number of detids " << detIds.size() << std::endl;
116 
117  for (; detid != detIds.end(); ++detid) {
119  int apv = 0;
120  for (int it = 0; it < range.second - range.first; ++it) {
121  ss << "detid " << *detid << " \t"
122  << " apv " << apv++ << " \t" << getApvGain(it, range) << " \t" << std::endl;
123  }
124  }
125 }

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

◆ printSummary()

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

Definition at line 127 of file SiStripGain.cc.

127  {
128  SiStripDetSummary summaryGain{trackerTopo};
129 
130  std::vector<unsigned int> detIds;
131  getDetIds(detIds);
132  std::vector<uint32_t>::const_iterator detid = detIds.begin();
133  for (; detid != detIds.end(); ++detid) {
135  for (int it = 0; it < range.second - range.first; ++it) {
136  summaryGain.add(*detid, getApvGain(it, range));
137  }
138  }
139  ss << "Summary of gain values:" << std::endl;
140  summaryGain.print(ss, true);
141 }

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

Member Data Documentation

◆ apvgain_

const SiStripApvGain* SiStripGain::apvgain_
private

Definition at line 113 of file SiStripGain.h.

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

◆ apvgainAutoPtr_

std::unique_ptr<SiStripApvGain> SiStripGain::apvgainAutoPtr_
private

Definition at line 114 of file SiStripGain.h.

Referenced by fillNewGain().

◆ apvgainVector_

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

Definition at line 111 of file SiStripGain.h.

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

◆ normVector_

std::vector<double> SiStripGain::normVector_
private

Definition at line 112 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 115 of file SiStripGain.h.

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

FastTimerService_cff.range
range
Definition: FastTimerService_cff.py:34
SiStripApvGain::getStripGain
static float getStripGain(uint16_t strip, const Range &range)
Definition: SiStripApvGain.h:78
SiStripGain::recordLabelPair_
std::vector< std::pair< std::string, std::string > > recordLabelPair_
Definition: SiStripGain.h:115
digitizers_cfi.strip
strip
Definition: digitizers_cfi.py:19
pos
Definition: PixelAliasList.h:18
SiStripGain::multiply
void multiply(const SiStripApvGain &apvgain, const double &factor, const std::pair< std::string, std::string > &recordLabelPair)
Definition: SiStripGain.cc:17
SiStripGain::normVector_
std::vector< double > normVector_
Definition: SiStripGain.h:112
SiStripGain::apvgainAutoPtr_
std::unique_ptr< SiStripApvGain > apvgainAutoPtr_
Definition: SiStripGain.h:114
personalPlayback.fp
fp
Definition: personalPlayback.py:523
SiStripApvGain::getRange
const Range getRange(const uint32_t detID) const
Definition: SiStripApvGain.cc:30
SiStripDetInfoFileReader
Definition: SiStripDetInfoFileReader.h:30
SiStripGain::apvgainVector_
std::vector< const SiStripApvGain * > apvgainVector_
Definition: SiStripGain.h:111
SiStripGain::fillNewGain
void fillNewGain(const SiStripApvGain *apvgain, const double &factor, const SiStripApvGain *apvgain2=nullptr, const double &factor2=1.)
Definition: SiStripGain.cc:39
SiStripGain::getStripGain
static float getStripGain(const uint16_t &strip, const SiStripApvGain::Range &range)
Definition: SiStripGain.h:73
contentValuesCheck.ss
ss
Definition: contentValuesCheck.py:33
pileupCalc.inputRange
inputRange
Definition: pileupCalc.py:218
DetId
Definition: DetId.h:17
edm::FileInPath
Definition: FileInPath.h:64
SiStripGain::apvgain_
const SiStripApvGain * apvgain_
Definition: SiStripGain.h:113
DQM.reader
reader
Definition: DQM.py:105
SiStripGain::getRange
const SiStripApvGain::Range getRange(uint32_t detID) const
Definition: SiStripGain.h:71
SiStripApvGain::put
bool put(const uint32_t &detID, Range input)
Definition: SiStripApvGain.cc:8
DQMScaleToClient_cfi.factor
factor
Definition: DQMScaleToClient_cfi.py:8
SiStripGain::getDetIds
void getDetIds(std::vector< uint32_t > &DetIds_) const
ATTENTION: we assume the detIds are the same as those from the first gain.
Definition: SiStripGain.cc:102
SiStripGain::getApvGain
static float getApvGain(const uint16_t &apv, const SiStripApvGain::Range &range)
Definition: SiStripGain.h:76
SiStripApvGain::getRangeByPos
Range getRangeByPos(unsigned short pos) const
Definition: SiStripApvGain.cc:44
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
SiStripApvGain::Range
std::pair< ContainerIterator, ContainerIterator > Range
Definition: SiStripApvGain.h:28
SiStripDetSummary
Definition: SiStripDetSummary.h:28
AlignmentPI::index
index
Definition: AlignmentPayloadInspectorHelper.h:46
SiStripApvGain::getDetIds
void getDetIds(std::vector< uint32_t > &DetIds_) const
Definition: SiStripApvGain.cc:53
SiStripApvGain
Definition: SiStripApvGain.h:25
SiStripApvGain::getApvGain
static float getApvGain(uint16_t apv, const Range &range)
Definition: SiStripApvGain.h:82