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 (layerid nTIB, layerid nTOB, layerid nTID, layerid nTEC)
	SiStripApvSimulationParameters ()
	~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

using SiStripApvSimulationParameters::layerid = unsigned int

Definition at line 19 of file SiStripApvSimulationParameters.h.

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

Definition at line 20 of file SiStripApvSimulationParameters.h.

Constructor & Destructor Documentation

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

SiStripApvSimulationParameters::SiStripApvSimulationParameters ( )

inline

Definition at line 29 of file SiStripApvSimulationParameters.h.

{}

SiStripApvSimulationParameters::~SiStripApvSimulationParameters ( )

inline

Definition at line 30 of file SiStripApvSimulationParameters.h.

{}

Member Function Documentation

void SiStripApvSimulationParameters::calculateIntegrals

(

)

Definition at line 40 of file SiStripApvSimulationParameters.cc.

References mps_fire::i.

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

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

inline

Definition at line 49 of file SiStripApvSimulationParameters.h.

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

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

inline

Definition at line 46 of file SiStripApvSimulationParameters.h.

{ return m_barrelParam[layer - 1]; }

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

inline

Definition at line 48 of file SiStripApvSimulationParameters.h.

{ return m_endcapParam[wheel - 1]; }

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

inline

Definition at line 47 of file SiStripApvSimulationParameters.h.

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

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

inline

Definition at line 43 of file SiStripApvSimulationParameters.h.

References putTEC().

Referenced by putTEC().

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

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

Definition at line 91 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;
}

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

inline

Definition at line 34 of file SiStripApvSimulationParameters.h.

References putTIB().

Referenced by putTIB().

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

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

Definition at line 55 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;
}

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

inline

Definition at line 40 of file SiStripApvSimulationParameters.h.

References putTID(), and makeMuonMisalignmentScenario::wheel.

Referenced by putTID().

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

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

Definition at line 79 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;
}

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

inline

Definition at line 37 of file SiStripApvSimulationParameters.h.

References putTOB().

Referenced by putTOB().

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

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

Definition at line 67 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;
}

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

private

Definition at line 103 of file SiStripApvSimulationParameters.cc.

References Exception, mps_fire::i, and heppy_batch::val.

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

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

private

Definition at line 131 of file SiStripApvSimulationParameters.cc.

References Exception, mps_fire::i, and heppy_batch::val.

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

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

inline

Definition at line 60 of file SiStripApvSimulationParameters.h.

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

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

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

inline

Definition at line 51 of file SiStripApvSimulationParameters.h.

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

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

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

inline

Definition at line 57 of file SiStripApvSimulationParameters.h.

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

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

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

inline

Definition at line 54 of file SiStripApvSimulationParameters.h.

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

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

template<class Archive >

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

private

Friends And Related Function Documentation

friend class boost::serialization::access

friend

Definition at line 74 of file SiStripApvSimulationParameters.h.

template<typename CondSerializationT , typename Enabled >

friend struct cond::serialization::access

friend

Definition at line 74 of file SiStripApvSimulationParameters.h.

Member Data Documentation

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

private

Definition at line 66 of file SiStripApvSimulationParameters.h.

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

private

Definition at line 67 of file SiStripApvSimulationParameters.h.

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

private

Definition at line 68 of file SiStripApvSimulationParameters.h.

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

private

Definition at line 69 of file SiStripApvSimulationParameters.h.

layerid SiStripApvSimulationParameters::m_nTEC

private

Definition at line 65 of file SiStripApvSimulationParameters.h.

layerid SiStripApvSimulationParameters::m_nTIB

private

Definition at line 65 of file SiStripApvSimulationParameters.h.

layerid SiStripApvSimulationParameters::m_nTID

private

Definition at line 65 of file SiStripApvSimulationParameters.h.

layerid SiStripApvSimulationParameters::m_nTOB

private

Definition at line 65 of file SiStripApvSimulationParameters.h.