HcalSimParameters Class Reference

#include <HcalSimParameters.h>

Inheritance diagram for HcalSimParameters:

Public Member Functions
bool	doSiPMSmearing () const
bool	doTimeSmear () const
double	fCtoGeV (const DetId &detId) const
	HcalSimParameters (const edm::ParameterSet &p)
	HcalSimParameters (double simHitToPhotoelectrons, double samplingFactor, double timePhase, int readoutFrameSize, int binOfMaximum, bool doPhotostatistics, bool syncPhase, int firstRing, const std::vector< double > &samplingFactors, double sipmTau)
double	photoelectronsToAnalog (const DetId &detId) const override
int	pixels (const DetId &detId) const
virtual double	samplingFactor (const DetId &detId) const
void	setDbService (const HcalDbService *service)
unsigned int	signalShape (const DetId &detId) const
double	simHitToPhotoelectrons (const DetId &detId) const override
double	sipmCrossTalk (const DetId &detId) const
double	sipmDarkCurrentuA (const DetId &detId) const
std::vector< float >	sipmNonlinearity (const DetId &detId) const
double	sipmTau () const
double	threshold_currentTDC () const
double	timeSmearRMS (double ampl) const
	~HcalSimParameters () override
Public Member Functions inherited from CaloSimParameters
int	binOfMaximum () const
	CaloSimParameters (const edm::ParameterSet &p, bool skipPe2Fc=false)
	CaloSimParameters (double simHitToPhotoelectrons, double photoelectronsToAnalog, double samplingFactor, double timePhase, int readoutFrameSize, int binOfMaximum, bool doPhotostatistics, bool syncPhase=true)
bool	doPhotostatistics () const
	whether or not to apply Poisson statistics to photoelectrons More...
double	photoelectronsToAnalog () const
	the factor which goes from photoelectrons to whatever gets read by ADCs More...
int	readoutFrameSize () const
	for now, the LinearFrames and trhe digis will be one-to-one. More...
void	setBinOfMaximum (int binOfMax)
void	setReadoutFrameSize (int frameSize)
	some datamixing apps need this to be set dynamically More...
double	simHitToPhotoelectrons () const
bool	syncPhase () const
double	timePhase () const
	the adjustment you need to apply to get the signal where you want it More...
virtual	~CaloSimParameters ()

Private Member Functions
void	defaultTimeSmearing ()

Private Attributes
bool	doTimeSmear_
const HcalDbService *	theDbService
int	theFirstRing
std::vector< double >	theSamplingFactors
const HcalSiPMCharacteristics *	theSiPMcharacteristics
bool	theSiPMSmearing
double	theSiPMTau
HcalTimeSmearSettings	theSmearSettings
double	threshold_currentTDC_

Friends
class	HcalSimParameterMap

Detailed Description

Definition at line 10 of file HcalSimParameters.h.

Constructor & Destructor Documentation

HcalSimParameters() [1/2]

HcalSimParameters::HcalSimParameters	(	double	simHitToPhotoelectrons,
		double	samplingFactor,
		double	timePhase,
		int	readoutFrameSize,
		int	binOfMaximum,
		bool	doPhotostatistics,
		bool	syncPhase,
		int	firstRing,
		const std::vector< double > &	samplingFactors,
		double	sipmTau
	)

Definition at line 12 of file HcalSimParameters.cc.

    : CaloSimParameters(simHitToPhotoelectrons,
                        0.0,
                        samplingFactor,
                        timePhase,
                        readoutFrameSize,
                        binOfMaximum,
                        doPhotostatistics,
                        syncPhase),
      theDbService(nullptr),
      theSiPMcharacteristics(nullptr),
      theFirstRing(firstRing),
      theSamplingFactors(samplingFactors),
      theSiPMSmearing(false),
      doTimeSmear_(true),
      theSiPMTau(sipmTau) {
  defaultTimeSmearing();
 
  edm::LogInfo("HcalSimParameters:") << " doSiPMsmearing    = " << theSiPMSmearing;
}

References defaultTimeSmearing(), and theSiPMSmearing.

HcalSimParameters() [2/2]

HcalSimParameters::HcalSimParameters

(

const edm::ParameterSet &

p

)

Definition at line 42 of file HcalSimParameters.cc.

    : CaloSimParameters(p, true),
      theDbService(nullptr),
      theFirstRing(p.getParameter<int>("firstRing")),
      theSamplingFactors(p.getParameter<std::vector<double> >("samplingFactors")),
      theSiPMSmearing(p.getParameter<bool>("doSiPMSmearing")),
      doTimeSmear_(p.getParameter<bool>("timeSmearing")),
      theSiPMTau(p.getParameter<double>("sipmTau")),
      threshold_currentTDC_(p.getParameter<double>("threshold_currentTDC")) {
  defaultTimeSmearing();
 
  edm::LogInfo("HcalSimParameters:") << " doSiPMsmearing    = " << theSiPMSmearing;
}

References defaultTimeSmearing(), and theSiPMSmearing.

~HcalSimParameters()

HcalSimParameters::~HcalSimParameters ( )

inlineoverride

Definition at line 24 of file HcalSimParameters.h.

{}

Member Function Documentation

defaultTimeSmearing()

void HcalSimParameters::defaultTimeSmearing ( )

private

Definition at line 107 of file HcalSimParameters.cc.

                                            {
  // GeV->ampl (fC), time (ns)
  theSmearSettings.emplace_back(4.00 * GeV2fC, 4.050);
  theSmearSettings.emplace_back(20.00 * GeV2fC, 3.300);
  theSmearSettings.emplace_back(25.00 * GeV2fC, 2.925);
  theSmearSettings.emplace_back(30.00 * GeV2fC, 2.714);
  theSmearSettings.emplace_back(37.00 * GeV2fC, 2.496);
  theSmearSettings.emplace_back(44.50 * GeV2fC, 2.278);
  theSmearSettings.emplace_back(56.00 * GeV2fC, 2.138);
  theSmearSettings.emplace_back(63.50 * GeV2fC, 2.022);
  theSmearSettings.emplace_back(81.00 * GeV2fC, 1.788);
  theSmearSettings.emplace_back(88.50 * GeV2fC, 1.695);
  theSmearSettings.emplace_back(114.50 * GeV2fC, 1.716);
  theSmearSettings.emplace_back(175.50 * GeV2fC, 1.070);
  theSmearSettings.emplace_back(350.00 * GeV2fC, 1.564);
  theSmearSettings.emplace_back(99999.00 * GeV2fC, 1.564);
}

References GeV2fC, and theSmearSettings.

Referenced by HcalSimParameters().

doSiPMSmearing()

bool HcalSimParameters::doSiPMSmearing ( ) const

inline

Definition at line 42 of file HcalSimParameters.h.

{ return theSiPMSmearing; }

References theSiPMSmearing.

Referenced by HcalSiPMHitResponse::makeSiPMSignal().

doTimeSmear()

bool HcalSimParameters::doTimeSmear ( ) const

inline

Definition at line 37 of file HcalSimParameters.h.

{ return doTimeSmear_; }

References doTimeSmear_.

fCtoGeV()

double HcalSimParameters::fCtoGeV

(

const DetId &

detId

)

const

{

Definition at line 74 of file HcalSimParameters.cc.

                                                          {
  assert(theDbService != nullptr);
  HcalGenericDetId hcalGenDetId(detId);
  const HcalGain* gains = theDbService->getGain(hcalGenDetId);
  const HcalGainWidth* gwidths = theDbService->getGainWidth(hcalGenDetId);
  double result = 0.0;
  if (!gains || !gwidths) {
    edm::LogError("HcalAmplifier") << "Could not fetch HCAL conditions for channel " << hcalGenDetId;
  } else {
    // only one gain will be recorded per channel, so just use capID 0 for now
    result = gains->getValue(0);
    //  if(doNoise_)
    //    result += CLHEP::RandGaussQ::shoot(0.,  gwidths->getValue(0));
    //  }
  }
  return result;
}

References cms::cuda::assert(), HcalDbService::getGain(), HcalDbService::getGainWidth(), HcalGain::getValue(), mps_fire::result, and theDbService.

Referenced by simHitToPhotoelectrons().

photoelectronsToAnalog()

double HcalSimParameters::photoelectronsToAnalog ( const DetId &  detId ) const

overridevirtual

Reimplemented from CaloSimParameters.

Definition at line 98 of file HcalSimParameters.cc.

                                                                         {
  //now always taken from database for HPDs or SiPMs (HB, HE, HO)
  assert(theDbService);
  return theDbService->getHcalSiPMParameter(detId)->getFCByPE();
}

References cms::cuda::assert(), HcalSiPMParameter::getFCByPE(), HcalDbService::getHcalSiPMParameter(), and theDbService.

Referenced by HcalSiPMHitResponse::addPEnoise().

pixels()

int HcalSimParameters::pixels

(

const DetId &

detId

)

const

Definition at line 149 of file HcalSimParameters.cc.

                                                      {
  assert(theDbService);
  int type = theDbService->getHcalSiPMParameter(detId)->getType();
  return theSiPMcharacteristics->getPixels(type);
}

References cms::cuda::assert(), HcalDbService::getHcalSiPMParameter(), HcalSiPMCharacteristics::getPixels(), HcalSiPMParameter::getType(), theDbService, and theSiPMcharacteristics.

Referenced by HcalSiPMHitResponse::makeSiPMSignal().

samplingFactor()

double HcalSimParameters::samplingFactor ( const DetId &  detId ) const

virtual

the ratio of actual incident energy to deposited energy in the SimHit

Definition at line 93 of file HcalSimParameters.cc.

                                                                 {
  HcalDetId hcalDetId(detId);
  return theSamplingFactors.at(hcalDetId.ietaAbs() - theFirstRing);
}

References HcalDetId::ietaAbs(), theFirstRing, and theSamplingFactors.

Referenced by simHitToPhotoelectrons().

setDbService()

void HcalSimParameters::setDbService

(

const HcalDbService *

service

)

Definition at line 56 of file HcalSimParameters.cc.

                                                                 {
  assert(service);
  theDbService = service;
  theSiPMcharacteristics = service->getHcalSiPMCharacteristics();
  assert(theSiPMcharacteristics);
}

References cms::cuda::assert(), theDbService, and theSiPMcharacteristics.

Referenced by HcalSimParameterMap::setDbService().

signalShape()

unsigned int HcalSimParameters::signalShape

(

const DetId &

detId

)

const

Definition at line 171 of file HcalSimParameters.cc.

                                                                    {
  assert(theDbService);
  return theDbService->getHcalMCParam(detId)->signalShape();
}

References cms::cuda::assert(), HcalDbService::getHcalMCParam(), HcalMCParam::signalShape(), and theDbService.

Referenced by HcalSiPMHitResponse::makeSiPMSignal().

simHitToPhotoelectrons()

double HcalSimParameters::simHitToPhotoelectrons ( const DetId &  detId ) const

overridevirtual

Reimplemented from CaloSimParameters.

Definition at line 63 of file HcalSimParameters.cc.

                                                                         {
  // the gain is in units of GeV/fC.  We want a constant with pe/dGeV
  // pe/dGeV = (GeV/dGeV) / (GeV/fC) / (fC/pe)
  double result = CaloSimParameters::simHitToPhotoelectrons(detId);
  if (HcalGenericDetId(detId).genericSubdet() != HcalGenericDetId::HcalGenForward ||
      HcalGenericDetId(detId).genericSubdet() != HcalGenericDetId::HcalGenZDC) {
    result = samplingFactor(detId) / fCtoGeV(detId) / photoelectronsToAnalog(detId);
  }
  return result;
}

References fCtoGeV(), HcalGenericDetId::HcalGenForward, HcalGenericDetId::HcalGenZDC, CaloSimParameters::photoelectronsToAnalog(), mps_fire::result, samplingFactor(), and CaloSimParameters::simHitToPhotoelectrons().

sipmCrossTalk()

double HcalSimParameters::sipmCrossTalk

(

const DetId &

detId

)

const

Definition at line 160 of file HcalSimParameters.cc.

                                                                {
  assert(theDbService);
  int type = theDbService->getHcalSiPMParameter(detId)->getType();
  return theSiPMcharacteristics->getCrossTalk(type);
}

References cms::cuda::assert(), HcalSiPMCharacteristics::getCrossTalk(), HcalDbService::getHcalSiPMParameter(), HcalSiPMParameter::getType(), theDbService, and theSiPMcharacteristics.

Referenced by HcalSiPMHitResponse::makeSiPMSignal().

sipmDarkCurrentuA()

double HcalSimParameters::sipmDarkCurrentuA

(

const DetId &

detId

)

const

Definition at line 155 of file HcalSimParameters.cc.

                                                                    {
  assert(theDbService);
  return theDbService->getHcalSiPMParameter(detId)->getDarkCurrent();
}

References cms::cuda::assert(), HcalSiPMParameter::getDarkCurrent(), HcalDbService::getHcalSiPMParameter(), and theDbService.

Referenced by HcalSiPMHitResponse::addPEnoise().

sipmNonlinearity()

std::vector< float > HcalSimParameters::sipmNonlinearity

(

const DetId &

detId

)

const

Definition at line 165 of file HcalSimParameters.cc.

                                                                             {
  assert(theDbService);
  int type = theDbService->getHcalSiPMParameter(detId)->getType();
  return theSiPMcharacteristics->getNonLinearities(type);
}

References cms::cuda::assert(), HcalDbService::getHcalSiPMParameter(), HcalSiPMCharacteristics::getNonLinearities(), HcalSiPMParameter::getType(), theDbService, and theSiPMcharacteristics.

Referenced by HcalSiPMHitResponse::makeSiPMSignal().

sipmTau()

double HcalSimParameters::sipmTau ( ) const

inline

Definition at line 45 of file HcalSimParameters.h.

{ return theSiPMTau; }

References theSiPMTau.

Referenced by HcalSiPMHitResponse::makeSiPMSignal().

threshold_currentTDC()

double HcalSimParameters::threshold_currentTDC ( ) const

inline

Definition at line 44 of file HcalSimParameters.h.

{ return threshold_currentTDC_; }

References threshold_currentTDC_.

Referenced by HcalElectronicsSim::analogToDigital().

timeSmearRMS()

double HcalSimParameters::timeSmearRMS

(

double

ampl

)

const

Definition at line 125 of file HcalSimParameters.cc.

                                                        {
  HcalTimeSmearSettings::size_type i;
  double smearsigma = 0;
 
  for (i = 0; i < theSmearSettings.size(); i++)
    if (theSmearSettings[i].first > ampl)
      break;
 
  // Smearing occurs only within the envelope definitions.
  if (i != 0 && (i < theSmearSettings.size())) {
    double energy1 = theSmearSettings[i - 1].first;
    double sigma1 = theSmearSettings[i - 1].second;
    double energy2 = theSmearSettings[i].first;
    double sigma2 = theSmearSettings[i].second;
 
    if (energy2 != energy1)
      smearsigma = sigma1 + ((sigma2 - sigma1) * (ampl - energy1) / (energy2 - energy1));
    else
      smearsigma = (sigma2 + sigma1) / 2.;
  }
 
  return smearsigma;
}

References dqmdumpme::first, mps_fire::i, and theSmearSettings.

Friends And Related Function Documentation

HcalSimParameterMap

friend class HcalSimParameterMap

friend

Definition at line 51 of file HcalSimParameters.h.

Member Data Documentation

doTimeSmear_

bool HcalSimParameters::doTimeSmear_

private

Definition at line 60 of file HcalSimParameters.h.

Referenced by doTimeSmear().

theDbService

const HcalDbService* HcalSimParameters::theDbService

private

Definition at line 55 of file HcalSimParameters.h.

Referenced by fCtoGeV(), photoelectronsToAnalog(), pixels(), setDbService(), signalShape(), sipmCrossTalk(), sipmDarkCurrentuA(), and sipmNonlinearity().

theFirstRing

int HcalSimParameters::theFirstRing

private

Definition at line 57 of file HcalSimParameters.h.

Referenced by samplingFactor().

theSamplingFactors

std::vector<double> HcalSimParameters::theSamplingFactors

private

Definition at line 58 of file HcalSimParameters.h.

Referenced by samplingFactor().

theSiPMcharacteristics

const HcalSiPMCharacteristics* HcalSimParameters::theSiPMcharacteristics

private

Definition at line 56 of file HcalSimParameters.h.

Referenced by pixels(), setDbService(), sipmCrossTalk(), and sipmNonlinearity().

theSiPMSmearing

bool HcalSimParameters::theSiPMSmearing

private

Definition at line 59 of file HcalSimParameters.h.

Referenced by doSiPMSmearing(), and HcalSimParameters().

theSiPMTau

double HcalSimParameters::theSiPMTau

private

Definition at line 62 of file HcalSimParameters.h.

Referenced by sipmTau().

theSmearSettings

HcalTimeSmearSettings HcalSimParameters::theSmearSettings

private

Definition at line 61 of file HcalSimParameters.h.

Referenced by defaultTimeSmearing(), and timeSmearRMS().

threshold_currentTDC_

double HcalSimParameters::threshold_currentTDC_

private

Definition at line 63 of file HcalSimParameters.h.

Referenced by threshold_currentTDC().