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, const std::vector< double > &photoelectronsToAnalog, double samplingFactor, double timePhase, int readoutFrameSize, int binOfMaximum, bool doPhotostatistics, bool syncPhase, int firstRing, const std::vector< double > &samplingFactors)
	HcalSimParameters (const edm::ParameterSet &p)
virtual double	photoelectronsToAnalog (const DetId &detId) const
int	pixels () const
virtual double	samplingFactor (const DetId &detId) const
void	setDbService (const HcalDbService *service)
virtual double	simHitToPhotoelectrons (const DetId &detId) 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	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 >	thePE2fCByRing
int	thePixels
std::vector< double >	theSamplingFactors
bool	theSiPMSmearing
HcalTimeSmearSettings	theSmearSettings

Friends
class	HcalSimParameterMap

Detailed Description

Definition at line 10 of file HcalSimParameters.h.

Constructor & Destructor Documentation

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

Definition at line 11 of file HcalSimParameters.cc.

References defaultTimeSmearing().

: CaloSimParameters(simHitToPhotoelectrons,  photoelectronsToAnalog[0], samplingFactor, timePhase,
                    readoutFrameSize, binOfMaximum, doPhotostatistics, syncPhase),
  theDbService(0),
  theFirstRing(firstRing),
  theSamplingFactors(samplingFactors),
  thePE2fCByRing(photoelectronsToAnalog),
  thePixels(0),
  theSiPMSmearing(false),
  doTimeSmear_(true)
{
  defaultTimeSmearing();
}

HcalSimParameters::HcalSimParameters

(

const edm::ParameterSet &

p

)

Definition at line 29 of file HcalSimParameters.cc.

References defaultTimeSmearing(), edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), thePixels, and theSiPMSmearing.

:  CaloSimParameters(p),
   theDbService(0),
   theFirstRing( p.getParameter<int>("firstRing") ),
   theSamplingFactors( p.getParameter<std::vector<double> >("samplingFactors") ),
   thePE2fCByRing( p.getParameter<std::vector<double> >("photoelectronsToAnalog") ),
   thePixels(0), theSiPMSmearing(false),
   doTimeSmear_( p.getParameter<bool>("timeSmearing"))
{
  if(p.exists("pixels"))
  {
    thePixels = p.getParameter<int>("pixels");
  }
  if (p.exists("doSiPMSmearing"))
    theSiPMSmearing = p.getParameter<bool>("doSiPMSmearing");
  defaultTimeSmearing();
}

virtual HcalSimParameters::~HcalSimParameters ( )

inlinevirtual

Definition at line 20 of file HcalSimParameters.h.

{}

Member Function Documentation

void HcalSimParameters::defaultTimeSmearing ( )

private

Definition at line 97 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 38 of file HcalSimParameters.h.

References theSiPMSmearing.

Referenced by HcalSiPMHitResponse::makeSiPMSignal().

{ return theSiPMSmearing; }

bool HcalSimParameters::doTimeSmear ( ) const

inline

Definition at line 33 of file HcalSimParameters.h.

References doTimeSmear_.

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

{ return doTimeSmear_; }

double HcalSimParameters::fCtoGeV

(

const DetId &

detId

)

const

{

Definition at line 61 of file HcalSimParameters.cc.

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

Referenced by simHitToPhotoelectrons().

{
  assert(theDbService != 0);
  HcalGenericDetId hcalGenDetId(detId);
  const HcalGain* gains = theDbService->getGain(hcalGenDetId);
  const HcalGainWidth* gwidths = theDbService->getGainWidth(hcalGenDetId);
  if (!gains || !gwidths )
  {
    edm::LogError("HcalAmplifier") << "Could not fetch HCAL conditions for channel " << hcalGenDetId;
  }
  // only one gain will be recorded per channel, so just use capID 0 for now
  double 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 87 of file HcalSimParameters.cc.

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

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

int HcalSimParameters::pixels ( ) const

inline

Definition at line 37 of file HcalSimParameters.h.

References thePixels.

Referenced by HcalSiPMHitResponse::makeSiPMSignal().

{return thePixels;}

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

virtual

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

Definition at line 80 of file HcalSimParameters.cc.

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

Referenced by simHitToPhotoelectrons().

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

void HcalSimParameters::setDbService ( const HcalDbService *  service )

inline

Definition at line 22 of file HcalSimParameters.h.

References theDbService.

Referenced by HcalSimParameterMap::setDbService().

{theDbService = service;}

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

virtual

Reimplemented from CaloSimParameters.

Definition at line 47 of file HcalSimParameters.cc.

References fCtoGeV(), HcalGenericDetId::HcalGenForward, HcalGenericDetId::HcalGenZDC, CaloSimParameters::photoelectronsToAnalog(), query::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::timeSmearRMS

(

double

ampl

)

const

Definition at line 115 of file HcalSimParameters.cc.

References plotBeamSpotDB::first, i, and theSmearSettings.

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

                                                        {
  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 40 of file HcalSimParameters.h.

Member Data Documentation

bool HcalSimParameters::doTimeSmear_

private

Definition at line 50 of file HcalSimParameters.h.

Referenced by doTimeSmear().

const HcalDbService* HcalSimParameters::theDbService

private

Definition at line 44 of file HcalSimParameters.h.

Referenced by fCtoGeV(), and setDbService().

int HcalSimParameters::theFirstRing

private

Definition at line 45 of file HcalSimParameters.h.

Referenced by photoelectronsToAnalog(), and samplingFactor().

std::vector<double> HcalSimParameters::thePE2fCByRing

private

Definition at line 47 of file HcalSimParameters.h.

Referenced by photoelectronsToAnalog().

int HcalSimParameters::thePixels

private

Definition at line 48 of file HcalSimParameters.h.

Referenced by HcalSimParameterMap::HcalSimParameterMap(), HcalSimParameters(), and pixels().

std::vector<double> HcalSimParameters::theSamplingFactors

private

Definition at line 46 of file HcalSimParameters.h.

Referenced by samplingFactor().

bool HcalSimParameters::theSiPMSmearing

private

Definition at line 49 of file HcalSimParameters.h.

Referenced by doSiPMSmearing(), and HcalSimParameters().

HcalTimeSmearSettings HcalSimParameters::theSmearSettings

private

Definition at line 51 of file HcalSimParameters.h.

Referenced by defaultTimeSmearing(), and timeSmearRMS().