SiStripApvSimulationParameters Class Reference

#include <SiStripApvSimulationParameters.h>

Public Types
using	layerid = unsigned int
using	LayerParameters = PhysicsTools::Calibration::HistogramF3D

Public Member Functions
void	calculateIntegrals ()
const LayerParameters &	getTEC (layerid wheel) const
const LayerParameters &	getTIB (layerid layer) const
const LayerParameters &	getTID (layerid wheel) const
const LayerParameters &	getTOB (layerid layer) const
bool	putTEC (layerid wheel, const LayerParameters &params)
bool	putTEC (layerid wheel, LayerParameters &&params)
bool	putTIB (layerid layer, const LayerParameters &params)
bool	putTIB (layerid layer, LayerParameters &&params)
bool	putTID (layerid wheel, const LayerParameters &params)
bool	putTID (layerid wheel, LayerParameters &&params)
bool	putTOB (layerid layer, const LayerParameters &params)
bool	putTOB (layerid layer, LayerParameters &&params)
float	sampleTEC (layerid wheel, float r, float pu, CLHEP::HepRandomEngine *engine) const
float	sampleTIB (layerid layer, float z, float pu, CLHEP::HepRandomEngine *engine) const
float	sampleTID (layerid wheel, float r, float pu, CLHEP::HepRandomEngine *engine) const
float	sampleTOB (layerid layer, float z, float pu, CLHEP::HepRandomEngine *engine) const
	SiStripApvSimulationParameters ()
	SiStripApvSimulationParameters (layerid nTIB, layerid nTOB, layerid nTID, layerid nTEC)
	~SiStripApvSimulationParameters ()

Private Member Functions
float	sampleBarrel (layerid layerIdx, float z, float pu, CLHEP::HepRandomEngine *engine) const
float	sampleEndcap (layerid wheelIdx, float r, float pu, CLHEP::HepRandomEngine *engine) const
template<class Archive >
void	serialize (Archive &ar, const unsigned int version)

Private Attributes
std::vector< PhysicsTools::Calibration::HistogramF3D >	m_barrelParam
std::vector< PhysicsTools::Calibration::HistogramF2D >	m_barrelParam_xInt
std::vector< PhysicsTools::Calibration::HistogramF3D >	m_endcapParam
std::vector< PhysicsTools::Calibration::HistogramF2D >	m_endcapParam_xInt
layerid	m_nTEC
layerid	m_nTIB
layerid	m_nTID
layerid	m_nTOB

Friends
class	boost::serialization::access
template<typename CondSerializationT , typename Enabled >
struct	cond::serialization::access

Detailed Description

Stores a histogram binned in PU, z and baseline voltage, for every barrel layer of the strip tracker

Definition at line 17 of file SiStripApvSimulationParameters.h.

Member Typedef Documentation

layerid

using SiStripApvSimulationParameters::layerid = unsigned int

Definition at line 19 of file SiStripApvSimulationParameters.h.

LayerParameters

using SiStripApvSimulationParameters::LayerParameters = PhysicsTools::Calibration::HistogramF3D

Definition at line 20 of file SiStripApvSimulationParameters.h.

Constructor & Destructor Documentation

SiStripApvSimulationParameters() [1/2]

SiStripApvSimulationParameters::SiStripApvSimulationParameters ( layerid  nTIB, layerid  nTOB, layerid  nTID, layerid  nTEC  )

inline

Definition at line 22 of file SiStripApvSimulationParameters.h.

      : m_nTIB(nTIB), m_nTOB(nTOB), m_nTID(nTID), m_nTEC(nTEC) {
    m_barrelParam.resize(m_nTIB + m_nTOB);
    m_barrelParam_xInt.resize(m_nTIB + m_nTOB);
    m_endcapParam.resize(m_nTID + m_nTEC);
    m_endcapParam_xInt.resize(m_nTID + m_nTEC);
  }

References m_barrelParam, m_barrelParam_xInt, m_endcapParam, m_endcapParam_xInt, m_nTEC, m_nTIB, m_nTID, and m_nTOB.

SiStripApvSimulationParameters() [2/2]

SiStripApvSimulationParameters::SiStripApvSimulationParameters ( )

inline

Definition at line 29 of file SiStripApvSimulationParameters.h.

{}

~SiStripApvSimulationParameters()

SiStripApvSimulationParameters::~SiStripApvSimulationParameters ( )

inline

Definition at line 30 of file SiStripApvSimulationParameters.h.

{}

Member Function Documentation

calculateIntegrals()

void SiStripApvSimulationParameters::calculateIntegrals

(

)

Definition at line 41 of file SiStripApvSimulationParameters.cc.

                                                        {
  if (m_barrelParam.size() != m_barrelParam_xInt.size()) {
    m_barrelParam_xInt.resize(m_barrelParam.size());
    for (unsigned int i{0}; i != m_barrelParam.size(); ++i) {
      m_barrelParam_xInt[i] = calculateXInt(m_barrelParam[i]);
    }
  }
  if (m_endcapParam.size() != m_endcapParam_xInt.size()) {
    m_endcapParam_xInt.resize(m_endcapParam.size());
    for (unsigned int i{0}; i != m_endcapParam.size(); ++i) {
      m_endcapParam_xInt[i] = calculateXInt(m_endcapParam[i]);
    }
  }
}

References mps_fire::i, m_barrelParam, m_barrelParam_xInt, m_endcapParam, and m_endcapParam_xInt.

getTEC()

const LayerParameters& SiStripApvSimulationParameters::getTEC ( layerid  wheel ) const

inline

Definition at line 49 of file SiStripApvSimulationParameters.h.

{ return m_endcapParam[m_nTID + wheel - 1]; }

References m_endcapParam, m_nTID, and makeMuonMisalignmentScenario::wheel.

getTIB()

const LayerParameters& SiStripApvSimulationParameters::getTIB ( layerid  layer ) const

inline

Definition at line 46 of file SiStripApvSimulationParameters.h.

{ return m_barrelParam[layer - 1]; }

References m_barrelParam.

getTID()

const LayerParameters& SiStripApvSimulationParameters::getTID ( layerid  wheel ) const

inline

Definition at line 48 of file SiStripApvSimulationParameters.h.

{ return m_endcapParam[wheel - 1]; }

References m_endcapParam, and makeMuonMisalignmentScenario::wheel.

getTOB()

const LayerParameters& SiStripApvSimulationParameters::getTOB ( layerid  layer ) const

inline

Definition at line 47 of file SiStripApvSimulationParameters.h.

{ return m_barrelParam[m_nTIB + layer - 1]; }

References m_barrelParam, and m_nTIB.

putTEC() [1/2]

bool SiStripApvSimulationParameters::putTEC ( layerid  wheel, const LayerParameters &  params  )

inline

Definition at line 43 of file SiStripApvSimulationParameters.h.

{ return putTEC(wheel, LayerParameters(params)); }

References CalibrationSummaryClient_cfi::params, putTEC(), and makeMuonMisalignmentScenario::wheel.

Referenced by putTEC().

putTEC() [2/2]

bool SiStripApvSimulationParameters::putTEC	(	SiStripApvSimulationParameters::layerid	wheel,
		SiStripApvSimulationParameters::LayerParameters &&	params
	)

Definition at line 92 of file SiStripApvSimulationParameters.cc.

                                                                                                    {
  if ((wheel > m_nTEC) || (wheel < 1)) {
    edm::LogError("SiStripApvSimulationParameters")
        << "[" << __PRETTY_FUNCTION__ << "] wheel index " << wheel << " out of range [1," << m_nTEC << ")";
    return false;
  }
  m_endcapParam[m_nTID + wheel - 1] = params;
  m_endcapParam_xInt[m_nTID + wheel - 1] = calculateXInt(params);
  return true;
}

References m_endcapParam, m_endcapParam_xInt, m_nTEC, m_nTID, CalibrationSummaryClient_cfi::params, and makeMuonMisalignmentScenario::wheel.

putTIB() [1/2]

bool SiStripApvSimulationParameters::putTIB ( layerid  layer, const LayerParameters &  params  )

inline

Definition at line 34 of file SiStripApvSimulationParameters.h.

{ return putTIB(layer, LayerParameters(params)); }

References CalibrationSummaryClient_cfi::params, and putTIB().

Referenced by putTIB().

putTIB() [2/2]

bool SiStripApvSimulationParameters::putTIB	(	SiStripApvSimulationParameters::layerid	layer,
		SiStripApvSimulationParameters::LayerParameters &&	params
	)

Definition at line 56 of file SiStripApvSimulationParameters.cc.

                                                                                                    {
  if ((layer > m_nTIB) || (layer < 1)) {
    edm::LogError("SiStripApvSimulationParameters")
        << "[" << __PRETTY_FUNCTION__ << "] layer index " << layer << " out of range [1," << m_nTIB << "]";
    return false;
  }
  m_barrelParam[layer - 1] = params;
  m_barrelParam_xInt[layer - 1] = calculateXInt(params);
  return true;
}

References m_barrelParam, m_barrelParam_xInt, m_nTIB, and CalibrationSummaryClient_cfi::params.

putTID() [1/2]

bool SiStripApvSimulationParameters::putTID ( layerid  wheel, const LayerParameters &  params  )

inline

Definition at line 40 of file SiStripApvSimulationParameters.h.

{ return putTID(wheel, LayerParameters(params)); }

References CalibrationSummaryClient_cfi::params, putTID(), and makeMuonMisalignmentScenario::wheel.

Referenced by putTID().

putTID() [2/2]

bool SiStripApvSimulationParameters::putTID	(	SiStripApvSimulationParameters::layerid	wheel,
		SiStripApvSimulationParameters::LayerParameters &&	params
	)

Definition at line 80 of file SiStripApvSimulationParameters.cc.

                                                                                                    {
  if ((wheel > m_nTID) || (wheel < 1)) {
    edm::LogError("SiStripApvSimulationParameters")
        << "[" << __PRETTY_FUNCTION__ << "] wheel index " << wheel << " out of range [1," << m_nTID << "]";
    return false;
  }
  m_endcapParam[wheel - 1] = params;
  m_endcapParam_xInt[wheel - 1] = calculateXInt(params);
  return true;
}

References m_endcapParam, m_endcapParam_xInt, m_nTID, CalibrationSummaryClient_cfi::params, and makeMuonMisalignmentScenario::wheel.

putTOB() [1/2]

bool SiStripApvSimulationParameters::putTOB ( layerid  layer, const LayerParameters &  params  )

inline

Definition at line 37 of file SiStripApvSimulationParameters.h.

{ return putTOB(layer, LayerParameters(params)); }

References CalibrationSummaryClient_cfi::params, and putTOB().

Referenced by putTOB().

putTOB() [2/2]

bool SiStripApvSimulationParameters::putTOB	(	SiStripApvSimulationParameters::layerid	layer,
		SiStripApvSimulationParameters::LayerParameters &&	params
	)

Definition at line 68 of file SiStripApvSimulationParameters.cc.

                                                                                                    {
  if ((layer > m_nTOB) || (layer < 1)) {
    edm::LogError("SiStripApvSimulationParameters")
        << "[" << __PRETTY_FUNCTION__ << "] layer index " << layer << " out of range [1," << m_nTOB << ")";
    return false;
  }
  m_barrelParam[m_nTIB + layer - 1] = params;
  m_barrelParam_xInt[m_nTIB + layer - 1] = calculateXInt(params);
  return true;
}

References m_barrelParam, m_barrelParam_xInt, m_nTIB, m_nTOB, and CalibrationSummaryClient_cfi::params.

sampleBarrel()

float SiStripApvSimulationParameters::sampleBarrel ( layerid  layerIdx, float  z, float  pu, CLHEP::HepRandomEngine *  engine  ) const

private

Definition at line 104 of file SiStripApvSimulationParameters.cc.

                                                                                       {
  if (m_barrelParam.size() != m_barrelParam_xInt.size()) {
    throw cms::Exception("LogicError") << "x-integrals of 3D histograms have not been calculated";
  }
  const auto layerParam = m_barrelParam[layerIdx];
  const int ip = layerParam.findBinY(pu);
  const int iz = layerParam.findBinZ(z);
  const float norm = m_barrelParam_xInt[layerIdx].binContent(ip, iz);
  const auto val = CLHEP::RandFlat::shoot(engine) * norm;
  if (val < layerParam.binContent(0, ip, iz)) {  // underflow
    return layerParam.rangeX().min;
  } else if (norm - val < layerParam.binContent(layerParam.numberOfBinsX() + 1, ip, iz)) {  // overflow
    return layerParam.rangeX().max;
  } else {  // loop over bins, return center of our bin
    float sum = layerParam.binContent(0, ip, iz);
    for (int i{1}; i != layerParam.numberOfBinsX() + 1; ++i) {
      sum += layerParam.binContent(i, ip, iz);
      if (sum > val) {
        return xBinPos(layerParam, i, (sum - val) / layerParam.binContent(i, ip, iz));
      }
    }
  }
  throw cms::Exception("LogicError") << "Problem drawing a random number from the distribution";
}

References Exception, mps_fire::i, m_barrelParam, m_barrelParam_xInt, muons2muons_cfi::pu, heppy_batch::val, and z.

Referenced by sampleTIB(), and sampleTOB().

sampleEndcap()

float SiStripApvSimulationParameters::sampleEndcap ( layerid  wheelIdx, float  r, float  pu, CLHEP::HepRandomEngine *  engine  ) const

private

Definition at line 132 of file SiStripApvSimulationParameters.cc.

                                                                                       {
  if (m_endcapParam.size() != m_endcapParam_xInt.size()) {
    throw cms::Exception("LogicError") << "x-integrals of 3D histograms have not been calculated";
  }
  const auto layerParam = m_endcapParam[wheelIdx];
  const int ip = layerParam.findBinY(pu);
  const int ir = layerParam.findBinZ(r);
  const float norm = m_endcapParam_xInt[wheelIdx].binContent(ip, ir);
  const auto val = CLHEP::RandFlat::shoot(engine) * norm;
  if (val < layerParam.binContent(0, ip, ir)) {  // underflow
    return layerParam.rangeX().min;
  } else if (norm - val < layerParam.binContent(layerParam.numberOfBinsX() + 1, ip, ir)) {  // overflow
    return layerParam.rangeX().max;
  } else {  // loop over bins, return center of our bin
    float sum = layerParam.binContent(0, ip, ir);
    for (int i{1}; i != layerParam.numberOfBinsX() + 1; ++i) {
      sum += layerParam.binContent(i, ip, ir);
      if (sum > val) {
        return xBinPos(layerParam, i, (sum - val) / layerParam.binContent(i, ip, ir));
      }
    }
  }
  throw cms::Exception("LogicError") << "Problem drawing a random number from the distribution";
}

References Exception, mps_fire::i, m_endcapParam, m_endcapParam_xInt, muons2muons_cfi::pu, alignCSCRings::r, and heppy_batch::val.

Referenced by sampleTEC(), and sampleTID().

sampleTEC()

float SiStripApvSimulationParameters::sampleTEC ( layerid  wheel, float  r, float  pu, CLHEP::HepRandomEngine *  engine  ) const

inline

Definition at line 60 of file SiStripApvSimulationParameters.h.

                                                                                        {
    return sampleEndcap(m_nTID + wheel - 1, r, pu, engine);
  }

References m_nTID, muons2muons_cfi::pu, alignCSCRings::r, sampleEndcap(), and makeMuonMisalignmentScenario::wheel.

Referenced by SiStripDigitizerAlgorithm::digitize(), and PreMixingSiStripWorker::put().

sampleTIB()

float SiStripApvSimulationParameters::sampleTIB ( layerid  layer, float  z, float  pu, CLHEP::HepRandomEngine *  engine  ) const

inline

Definition at line 51 of file SiStripApvSimulationParameters.h.

                                                                                        {
    return sampleBarrel(layer - 1, z, pu, engine);
  }

References muons2muons_cfi::pu, sampleBarrel(), and z.

Referenced by SiStripDigitizerAlgorithm::digitize(), and PreMixingSiStripWorker::put().

sampleTID()

float SiStripApvSimulationParameters::sampleTID ( layerid  wheel, float  r, float  pu, CLHEP::HepRandomEngine *  engine  ) const

inline

Definition at line 57 of file SiStripApvSimulationParameters.h.

                                                                                        {
    return sampleEndcap(wheel - 1, r, pu, engine);
  }

References muons2muons_cfi::pu, alignCSCRings::r, sampleEndcap(), and makeMuonMisalignmentScenario::wheel.

Referenced by SiStripDigitizerAlgorithm::digitize(), and PreMixingSiStripWorker::put().

sampleTOB()

float SiStripApvSimulationParameters::sampleTOB ( layerid  layer, float  z, float  pu, CLHEP::HepRandomEngine *  engine  ) const

inline

Definition at line 54 of file SiStripApvSimulationParameters.h.

                                                                                        {
    return sampleBarrel(m_nTIB + layer - 1, z, pu, engine);
  };

References m_nTIB, muons2muons_cfi::pu, sampleBarrel(), and z.

Referenced by SiStripDigitizerAlgorithm::digitize(), and PreMixingSiStripWorker::put().

serialize()

template<class Archive >

void SiStripApvSimulationParameters::serialize ( Archive &  ar, const unsigned int  version  )

private

Friends And Related Function Documentation

boost::serialization::access

friend class boost::serialization::access

friend

Definition at line 74 of file SiStripApvSimulationParameters.h.

cond::serialization::access

template<typename CondSerializationT , typename Enabled >

friend struct cond::serialization::access

friend

Definition at line 74 of file SiStripApvSimulationParameters.h.

Member Data Documentation

m_barrelParam

std::vector<PhysicsTools::Calibration::HistogramF3D> SiStripApvSimulationParameters::m_barrelParam

private

Definition at line 66 of file SiStripApvSimulationParameters.h.

Referenced by calculateIntegrals(), getTIB(), getTOB(), putTIB(), putTOB(), sampleBarrel(), and SiStripApvSimulationParameters().

m_barrelParam_xInt

std::vector<PhysicsTools::Calibration::HistogramF2D> SiStripApvSimulationParameters::m_barrelParam_xInt

private

Definition at line 67 of file SiStripApvSimulationParameters.h.

Referenced by calculateIntegrals(), putTIB(), putTOB(), sampleBarrel(), and SiStripApvSimulationParameters().

m_endcapParam

std::vector<PhysicsTools::Calibration::HistogramF3D> SiStripApvSimulationParameters::m_endcapParam

private

Definition at line 68 of file SiStripApvSimulationParameters.h.

Referenced by calculateIntegrals(), getTEC(), getTID(), putTEC(), putTID(), sampleEndcap(), and SiStripApvSimulationParameters().

m_endcapParam_xInt

std::vector<PhysicsTools::Calibration::HistogramF2D> SiStripApvSimulationParameters::m_endcapParam_xInt

private

Definition at line 69 of file SiStripApvSimulationParameters.h.

Referenced by calculateIntegrals(), putTEC(), putTID(), sampleEndcap(), and SiStripApvSimulationParameters().

m_nTEC

layerid SiStripApvSimulationParameters::m_nTEC

private

Definition at line 65 of file SiStripApvSimulationParameters.h.

Referenced by putTEC(), and SiStripApvSimulationParameters().

m_nTIB

layerid SiStripApvSimulationParameters::m_nTIB

private

Definition at line 65 of file SiStripApvSimulationParameters.h.

Referenced by getTOB(), putTIB(), putTOB(), sampleTOB(), and SiStripApvSimulationParameters().

m_nTID

layerid SiStripApvSimulationParameters::m_nTID

private

Definition at line 65 of file SiStripApvSimulationParameters.h.

Referenced by getTEC(), putTEC(), putTID(), sampleTEC(), and SiStripApvSimulationParameters().

m_nTOB

layerid SiStripApvSimulationParameters::m_nTOB

private

Definition at line 65 of file SiStripApvSimulationParameters.h.

Referenced by putTOB(), and SiStripApvSimulationParameters().