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 (double simHitToPhotoelectrons, double samplingFactor, double timePhase, int readoutFrameSize, int binOfMaximum, bool doPhotostatistics, bool syncPhase, int firstRing, const std::vector< double > &samplingFactors, double sipmTau)
	HcalSimParameters (const edm::ParameterSet &p)
virtual double	photoelectronsToAnalog (const DetId &detId) const
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
virtual double	simHitToPhotoelectrons (const DetId &detId) const
double	sipmCrossTalk (const DetId &detId) const
double	sipmDarkCurrentuA (const DetId &detId) const
std::vector< float >	sipmNonlinearity (const DetId &detId) const
double	sipmTau () const
double	timeSmearRMS (double ampl) const
virtual	~HcalSimParameters ()
Public Member Functions inherited from CaloSimParameters
int	binOfMaximum () const
	CaloSimParameters (double simHitToPhotoelectrons, double photoelectronsToAnalog, double samplingFactor, double timePhase, int readoutFrameSize, int binOfMaximum, bool doPhotostatistics, bool syncPhase=true)
	CaloSimParameters (const edm::ParameterSet &p, bool skipPe2Fc=false)
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
	choice of the ADC time alignment (synchronous for LHC, asynchronous for test beams) More...
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

Friends
class	HcalSimParameterMap

Detailed Description

Definition at line 10 of file HcalSimParameters.h.

Constructor & Destructor Documentation

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.

References defaultTimeSmearing(), and theSiPMSmearing.

: CaloSimParameters(simHitToPhotoelectrons,  0.0, samplingFactor, timePhase,
                    readoutFrameSize, binOfMaximum, doPhotostatistics, syncPhase),
  theDbService(0),
  theSiPMcharacteristics(0),
  theFirstRing(firstRing),
  theSamplingFactors(samplingFactors),
  theSiPMSmearing(false),
  doTimeSmear_(true),
  theSiPMTau(sipmTau)
{
  defaultTimeSmearing();

  edm::LogInfo("HcalSimParameters:") << " doSiPMsmearing    = " << theSiPMSmearing;
}

HcalSimParameters::HcalSimParameters

(

const edm::ParameterSet &

p

)

Definition at line 34 of file HcalSimParameters.cc.

References defaultTimeSmearing(), and theSiPMSmearing.

:  CaloSimParameters(p,true),
   theDbService(0),
   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") )
{
  defaultTimeSmearing();

  edm::LogInfo("HcalSimParameters:") << " doSiPMsmearing    = " << theSiPMSmearing;
}

virtual HcalSimParameters::~HcalSimParameters ( )

inlinevirtual

Definition at line 22 of file HcalSimParameters.h.

References fCtoGeV(), CaloSimParameters::photoelectronsToAnalog(), samplingFactor(), setDbService(), and CaloSimParameters::simHitToPhotoelectrons().

{}

Member Function Documentation

void HcalSimParameters::defaultTimeSmearing ( )

private

Definition at line 111 of file HcalSimParameters.cc.

References GeV2fC, and theSmearSettings.

Referenced by HcalSimParameters().

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

bool HcalSimParameters::doSiPMSmearing ( ) const

inline

Definition at line 40 of file HcalSimParameters.h.

References theSiPMSmearing.

Referenced by HcalSiPMHitResponse::makeSiPMSignal().

{ return theSiPMSmearing; }

bool HcalSimParameters::doTimeSmear ( ) const

inline

Definition at line 35 of file HcalSimParameters.h.

References doTimeSmear_, pixels(), and timeSmearRMS().

Referenced by HcalTimeSlewSim::delay().

{ return doTimeSmear_; }

double HcalSimParameters::fCtoGeV

(

const DetId &

detId

)

const

{

Definition at line 70 of file HcalSimParameters.cc.

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

Referenced by simHitToPhotoelectrons(), and ~HcalSimParameters().

{
  assert(theDbService != 0);
  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;
}

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

virtual

Reimplemented from CaloSimParameters.

Definition at line 100 of file HcalSimParameters.cc.

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

Referenced by HcalSiPMHitResponse::addPEnoise().

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

int HcalSimParameters::pixels

(

const DetId &

detId

)

const

Definition at line 154 of file HcalSimParameters.cc.

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

Referenced by doTimeSmear(), and HcalSiPMHitResponse::makeSiPMSignal().

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

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.

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

Referenced by simHitToPhotoelectrons(), and ~HcalSimParameters().

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

void HcalSimParameters::setDbService

(

const HcalDbService *

service

)

Definition at line 48 of file HcalSimParameters.cc.

References HcalDbService::getHcalSiPMCharacteristics(), theDbService, and theSiPMcharacteristics.

Referenced by HcalSimParameterMap::setDbService(), and ~HcalSimParameters().

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

unsigned int HcalSimParameters::signalShape

(

const DetId &

detId

)

const

Definition at line 179 of file HcalSimParameters.cc.

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

Referenced by HcalSiPMHitResponse::makeSiPMSignal(), and sipmTau().

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

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

virtual

Reimplemented from CaloSimParameters.

Definition at line 56 of file HcalSimParameters.cc.

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

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

double HcalSimParameters::sipmCrossTalk

(

const DetId &

detId

)

const

Definition at line 166 of file HcalSimParameters.cc.

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

Referenced by HcalSiPMHitResponse::makeSiPMSignal(), and sipmTau().

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

double HcalSimParameters::sipmDarkCurrentuA

(

const DetId &

detId

)

const

Definition at line 160 of file HcalSimParameters.cc.

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

Referenced by HcalSiPMHitResponse::addPEnoise(), and sipmTau().

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

std::vector< float > HcalSimParameters::sipmNonlinearity

(

const DetId &

detId

)

const

Definition at line 172 of file HcalSimParameters.cc.

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

Referenced by HcalSiPMHitResponse::makeSiPMSignal(), and sipmTau().

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

double HcalSimParameters::sipmTau ( ) const

inline

Definition at line 42 of file HcalSimParameters.h.

References signalShape(), sipmCrossTalk(), sipmDarkCurrentuA(), sipmNonlinearity(), and theSiPMTau.

Referenced by HcalSiPMHitResponse::makeSiPMSignal().

{ return theSiPMTau; }

double HcalSimParameters::timeSmearRMS

(

double

ampl

)

const

Definition at line 129 of file HcalSimParameters.cc.

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

Referenced by HcalTimeSlewSim::delay(), and doTimeSmear().

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

}

Friends And Related Function Documentation

friend class HcalSimParameterMap

friend

Definition at line 48 of file HcalSimParameters.h.

Member Data Documentation

bool HcalSimParameters::doTimeSmear_

private

Definition at line 57 of file HcalSimParameters.h.

Referenced by doTimeSmear().

const HcalDbService* HcalSimParameters::theDbService

private

Definition at line 52 of file HcalSimParameters.h.

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

int HcalSimParameters::theFirstRing

private

Definition at line 54 of file HcalSimParameters.h.

Referenced by samplingFactor().

std::vector<double> HcalSimParameters::theSamplingFactors

private

Definition at line 55 of file HcalSimParameters.h.

Referenced by samplingFactor().

const HcalSiPMCharacteristics* HcalSimParameters::theSiPMcharacteristics

private

Definition at line 53 of file HcalSimParameters.h.

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

bool HcalSimParameters::theSiPMSmearing

private

Definition at line 56 of file HcalSimParameters.h.

Referenced by doSiPMSmearing(), and HcalSimParameters().

double HcalSimParameters::theSiPMTau

private

Definition at line 59 of file HcalSimParameters.h.

Referenced by sipmTau().

HcalTimeSmearSettings HcalSimParameters::theSmearSettings

private

Definition at line 58 of file HcalSimParameters.h.

Referenced by defaultTimeSmearing(), and timeSmearRMS().