---

CaloHitResponse Class Reference

Creates electronics signals from hits. More...

#include <SimCalorimetry/CaloSimAlgos/interface/CaloHitResponse.h>

List of all members.

Public Types
enum	{ BUNCHSPACE = 25 }
typedef std::map< DetId, CaloSamples >	AnalogSignalMap
Public Member Functions
void	add (const CaloSamples &signal)
	add a signal, in units of pe
void	add (const PCaloHit &hit)
	process a single SimHit
void	addHit (const PCaloHit *hit, CaloSamples &frame) const
	adds the amplitude for a single hit to the frame
double	analogSignalAmplitude (const PCaloHit &hit, const CaloSimParameters &parameters) const
	finds the amplitude contribution from this hit, applying photostatistics, if needed.
	CaloHitResponse (const CaloVSimParameterMap *parameterMap, const CaloVShape *shape)
void	clear ()
	frees up memory
CaloSamples *	findSignal (const DetId &detId)
	users can look for the signal for a given cell
CaloSamples	makeAnalogSignal (const PCaloHit &inputHit) const
	creates the signal corresponding to this hit
CaloSamples	makeBlankSignal (const DetId &detId) const
	creates an empty signal for this DetId
int	nSignals () const
	number of signals in the current cache
void	run (MixCollection< PCaloHit > &hits)
	Complete cell digitization.
void	setBunchRange (int minBunch, int maxBunch)
	tells it which pileup bunches to do
void	setGeometry (const CaloGeometry *geometry)
	geometry needed for time-of-flight
void	setHitCorrection (const CaloVHitCorrection *hitCorrection)
	If you want to correct hits, for attenuation or delay, set this.
void	setHitFilter (const CaloVHitFilter *filter)
	if you want to reject hits, for example, from a certain subdetector, set this
void	setPhaseShift (const double &thePhaseShift)
	setting the phase shift for asynchronous trigger (e.g. test beams)
void	setRandomEngine (CLHEP::HepRandomEngine &engine)
double	timeOfFlight (const DetId &detId) const
	time-of-flight, in ns, to get to this cell returns 0 if no geometry has been set
	~CaloHitResponse ()
	doesn't delete the pointers passed in
Protected Attributes
AnalogSignalMap	theAnalogSignalMap
const CaloGeometry *	theGeometry
const CaloVHitCorrection *	theHitCorrection
const CaloVHitFilter *	theHitFilter
int	theMaxBunch
int	theMinBunch
const CaloVSimParameterMap *	theParameterMap
double	thePhaseShift_
CLHEP::RandPoissonQ *	theRandPoisson
const CaloVShape *	theShape

---

Detailed Description

Creates electronics signals from hits.

Definition at line 29 of file CaloHitResponse.h.

---

Member Typedef Documentation

typedef std::map<DetId, CaloSamples> CaloHitResponse::AnalogSignalMap

Definition at line 32 of file CaloHitResponse.h.

---

Member Enumeration Documentation

anonymous enum

Enumerator:

BUNCHSPACE

Definition at line 34 of file CaloHitResponse.h.

{BUNCHSPACE=25};

---

Constructor & Destructor Documentation

CaloHitResponse::CaloHitResponse	(	const CaloVSimParameterMap *	parameterMap,
		const CaloVShape *	shape
	)

Definition at line 21 of file CaloHitResponse.cc.

: theAnalogSignalMap(),
  theParameterMap(parametersMap), 
  theShape(shape),  
  theHitCorrection(0),
  theHitFilter(0),
  theGeometry(0),
  theRandPoisson(0),
  theMinBunch(-10), 
  theMaxBunch(10),
  thePhaseShift_(1.)
{
}

CaloHitResponse::~CaloHitResponse

(

)

doesn't delete the pointers passed in

Definition at line 37 of file CaloHitResponse.cc.

References theRandPoisson.

{
  delete theRandPoisson;
}

---

Member Function Documentation

void CaloHitResponse::add

(

const CaloSamples &

signal

)

add a signal, in units of pe

Definition at line 84 of file CaloHitResponse.cc.

References findSignal(), i, CaloSamples::id(), id, CaloSamples::presamples(), CaloSamples::size(), and theAnalogSignalMap.

{
  DetId id(signal.id());
  CaloSamples * oldSignal = findSignal(id);
  if (oldSignal == 0) {
    theAnalogSignalMap[id] = signal;
  } else  {
    // need a "+=" to CaloSamples
    int sampleSize =  oldSignal->size();
    assert(sampleSize == signal.size());
    assert(signal.presamples() == oldSignal->presamples());

    for(int i = 0; i < sampleSize; ++i) {
      (*oldSignal)[i] += signal[i];
    }
  }
}

void CaloHitResponse::add

(

const PCaloHit &

hit

)

process a single SimHit

Definition at line 69 of file CaloHitResponse.cc.

References CaloVHitFilter::accepts(), edm::isnan(), LogDebug, makeAnalogSignal(), signal, theHitFilter, and PCaloHit::time().

Referenced by CaloTDigitizer< Traits >::addNoiseHits(), ESFastTDigitizer::addNoiseHits(), and CaloTDigitizer< Traits >::addNoiseSignals().

{
  // check the hit time makes sense
  if ( isnan(hit.time()) ) { return; }

  // maybe it's not from this subdetector
  if(theHitFilter == 0 || theHitFilter->accepts(hit)) {
    LogDebug("CaloHitResponse") << hit;
    CaloSamples signal(makeAnalogSignal(hit));
    LogDebug("CaloHitResponse") << signal;
    add(signal);
  }
}

void CaloHitResponse::addHit	(	const PCaloHit *	hit,
		CaloSamples &	frame
	)			const

adds the amplitude for a single hit to the frame

double CaloHitResponse::analogSignalAmplitude	(	const PCaloHit &	hit,
		const CaloSimParameters &	parameters
	)			const

finds the amplitude contribution from this hit, applying photostatistics, if needed.

Results are in photoelectrons

Definition at line 135 of file CaloHitResponse.cc.

References CaloSimParameters::doPhotostatistics(), PCaloHit::energy(), Exception, PCaloHit::id(), edm::Service< T >::isAvailable(), CaloSimParameters::simHitToPhotoelectrons(), and theRandPoisson.

Referenced by makeAnalogSignal().

                                                                                                              {

  if(!theRandPoisson)
  {
    edm::Service<edm::RandomNumberGenerator> rng;
    if ( ! rng.isAvailable()) {
      throw cms::Exception("Configuration")
        << "CaloHitResponse requires the RandomNumberGeneratorService\n"
        "which is not present in the configuration file.  You must add the service\n"
        "in the configuration file or remove the modules that require it.";
    }
    theRandPoisson = new CLHEP::RandPoissonQ(rng->getEngine());
  }

  // OK, the "energy" in the hit could be a real energy, deposited energy,
  // or pe count.  This factor converts to photoelectrons
  double npe = hit.energy() * parameters.simHitToPhotoelectrons( DetId(hit.id()) );
  // do we need to doPoisson statistics for the photoelectrons?
  if(parameters.doPhotostatistics()) {
    npe = theRandPoisson->fire(npe);
  }
  return npe;
}

void CaloHitResponse::clear

(

void

)

[inline]

frees up memory

Definition at line 69 of file CaloHitResponse.h.

References theAnalogSignalMap.

Referenced by EcalMixingModuleValidation::computeSDBunchDigi(), ESFastTDigitizer::run(), CaloTDigitizer< Traits >::run(), and EcalTDigitizer< Traits >::run().

{theAnalogSignalMap.clear();}

CaloSamples * CaloHitResponse::findSignal

(

const DetId &

detId

)

users can look for the signal for a given cell

Definition at line 160 of file CaloHitResponse.cc.

References HLT_VtxMuL3::result, and theAnalogSignalMap.

Referenced by add(), EcalMixingModuleValidation::computeSDBunchDigi(), ESFastTDigitizer::run(), CaloTDigitizer< Traits >::run(), and EcalTDigitizer< Traits >::run().

                                                             {
  CaloSamples * result = 0;
  AnalogSignalMap::iterator signalItr = theAnalogSignalMap.find(detId);
  if(signalItr == theAnalogSignalMap.end()) {
     result = 0;
  } 
  else {
     result = &(signalItr->second);
  }
  return result;
}

CaloSamples CaloHitResponse::makeAnalogSignal

(

const PCaloHit &

inputHit

)

const

creates the signal corresponding to this hit

Definition at line 103 of file CaloHitResponse.cc.

References analogSignalAmplitude(), BUNCHSPACE, CaloVHitCorrection::correct(), detId, PCaloHit::id(), makeBlankSignal(), python::trackProbabilityAnalysis_cff::parameters, HLT_VtxMuL3::result, signal, CaloVSimParameterMap::simParameters(), CaloSamples::size(), theHitCorrection, theParameterMap, thePhaseShift_, PCaloHit::time(), timeOfFlight(), and tzero.

Referenced by add().

                                                                             {

  // see if we need to correct the hit 
  PCaloHit hit = inputHit;

  if(theHitCorrection != 0) {
    theHitCorrection->correct(hit);
  }

  DetId detId(hit.id());
  const CaloSimParameters & parameters = theParameterMap->simParameters(detId);

  double signal = analogSignalAmplitude(hit, parameters);

  double jitter = hit.time() - timeOfFlight(detId);

  // assume bins count from zero, go for center of bin
  const double tzero = parameters.timePhase() -jitter -
     BUNCHSPACE*(parameters.binOfMaximum()-thePhaseShift_);
  double binTime = tzero;

  CaloSamples result(makeBlankSignal(detId));

  for(int bin = 0; bin < result.size(); bin++) {
    result[bin] += (*theShape)(binTime)* signal;
    binTime += BUNCHSPACE;
  }

  return result;
} 

CaloSamples CaloHitResponse::makeBlankSignal

(

const DetId &

detId

)

const

creates an empty signal for this DetId

Definition at line 173 of file CaloHitResponse.cc.

References CaloSimParameters::binOfMaximum(), python::trackProbabilityAnalysis_cff::parameters, CaloSimParameters::readoutFrameSize(), HLT_VtxMuL3::result, CaloVSimParameterMap::simParameters(), and theParameterMap.

Referenced by makeAnalogSignal(), ESFastTDigitizer::run(), CaloTDigitizer< Traits >::run(), and EcalTDigitizer< Traits >::run().

                                                                      {
  const CaloSimParameters & parameters = theParameterMap->simParameters(detId);
  CaloSamples result(detId, parameters.readoutFrameSize());
  result.setPresamples(parameters.binOfMaximum()-1);
  return result;
}

int CaloHitResponse::nSignals

(

)

const [inline]

number of signals in the current cache

Definition at line 85 of file CaloHitResponse.h.

References theAnalogSignalMap.

Referenced by ESFastTDigitizer::run(), CaloTDigitizer< Traits >::run(), and EcalTDigitizer< Traits >::run().

{return theAnalogSignalMap.size();}

void CaloHitResponse::run

(

MixCollection< PCaloHit > &

hits

)

Complete cell digitization.

Definition at line 55 of file CaloHitResponse.cc.

References MixCollection< T >::begin(), MixCollection< T >::end(), theMaxBunch, and theMinBunch.

Referenced by EcalMixingModuleValidation::computeSDBunchDigi(), ESFastTDigitizer::run(), CaloTDigitizer< Traits >::run(), and EcalTDigitizer< Traits >::run().

                                                        {

  for(MixCollection<PCaloHit>::MixItr hitItr = hits.begin();
      hitItr != hits.end(); ++hitItr)
  {
    // check the bunch crossing range
    if ( hitItr.bunch() < theMinBunch || hitItr.bunch() > theMaxBunch ) 
      { continue; }
  
    add(*hitItr);
  } // loop over hits
}

void CaloHitResponse::setBunchRange	(	int	minBunch,
		int	maxBunch
	)

tells it which pileup bunches to do

Definition at line 43 of file CaloHitResponse.cc.

References theMaxBunch, and theMinBunch.

Referenced by EcalMixingModuleValidation::computeSDBunchDigi().

                                                              {
  theMinBunch = minBunch;
  theMaxBunch = maxBunch;
}

void CaloHitResponse::setGeometry

(

const CaloGeometry *

geometry

)

[inline]

geometry needed for time-of-flight

Definition at line 45 of file CaloHitResponse.h.

References theGeometry.

Referenced by CastorDigiProducer::checkGeometry(), HcalDigiProducer::checkGeometry(), EcalMixingModuleValidation::computeSDBunchDigi(), HcalTBDigiProducer::updateGeometry(), EcalTBDigiProducer::updateGeometry(), and EcalDigiProducer::updateGeometry().

{ theGeometry = geometry; }

void CaloHitResponse::setHitCorrection

(

const CaloVHitCorrection *

hitCorrection

)

[inline]

If you want to correct hits, for attenuation or delay, set this.

Definition at line 62 of file CaloHitResponse.h.

References theHitCorrection.

Referenced by CastorDigiProducer::CastorDigiProducer(), HcalDigiProducer::HcalDigiProducer(), and HcalTBDigiProducer::HcalTBDigiProducer().

                                                                  {
    theHitCorrection = hitCorrection;
  }

void CaloHitResponse::setHitFilter

(

const CaloVHitFilter *

filter

)

[inline]

if you want to reject hits, for example, from a certain subdetector, set this

Definition at line 57 of file CaloHitResponse.h.

References theHitFilter.

Referenced by CastorDigiProducer::CastorDigiProducer(), HcalDigiProducer::HcalDigiProducer(), and HcalTBDigiProducer::HcalTBDigiProducer().

                                                   {
    theHitFilter = filter;
  }

void CaloHitResponse::setPhaseShift

(

const double &

thePhaseShift

)

[inline]

setting the phase shift for asynchronous trigger (e.g. test beams)

Definition at line 96 of file CaloHitResponse.h.

References thePhaseShift_.

Referenced by EcalDigiProducer::EcalDigiProducer(), HcalTBDigiProducer::setPhaseShift(), and EcalTBDigiProducer::setPhaseShift().

{ thePhaseShift_ = thePhaseShift; }

void CaloHitResponse::setRandomEngine

(

CLHEP::HepRandomEngine &

engine

)

Definition at line 49 of file CaloHitResponse.cc.

References theRandPoisson.

{
  theRandPoisson = new CLHEP::RandPoissonQ(engine);
}

double CaloHitResponse::timeOfFlight

(

const DetId &

detId

)

const

time-of-flight, in ns, to get to this cell returns 0 if no geometry has been set

Definition at line 181 of file CaloHitResponse.cc.

References CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), PV3DBase< T, PVType, FrameType >::mag(), DetId::rawId(), HLT_VtxMuL3::result, and theGeometry.

Referenced by makeAnalogSignal().

                                                              {
  // not going to assume there's one of these per subdetector.
  // Take the whole CaloGeometry and find the right subdet
  double result = 0.;
  if(theGeometry == 0) {
    edm::LogWarning("CaloHitResponse") << "No Calo Geometry set, so no time of flight correction";
  } 
  else {
    const CaloCellGeometry* cellGeometry = theGeometry->getSubdetectorGeometry(detId)->getGeometry(detId);
    if(cellGeometry == 0) {
       edm::LogWarning("CaloHitResponse") << "No Calo cell found for ID"
         << detId.rawId() << " so no time-of-flight subtraction will be done";
    }
    else {
      double distance = cellGeometry->getPosition().mag();
      result =  distance * cm / c_light; // Units of c_light: mm/ns
    }
  }
  return result;
}

---

Member Data Documentation

AnalogSignalMap CaloHitResponse::theAnalogSignalMap [protected]

Definition at line 100 of file CaloHitResponse.h.

Referenced by add(), clear(), findSignal(), and nSignals().

const CaloGeometry* CaloHitResponse::theGeometry [protected]

Definition at line 107 of file CaloHitResponse.h.

Referenced by setGeometry(), and timeOfFlight().

const CaloVHitCorrection* CaloHitResponse::theHitCorrection [protected]

Definition at line 104 of file CaloHitResponse.h.

Referenced by makeAnalogSignal(), and setHitCorrection().

const CaloVHitFilter* CaloHitResponse::theHitFilter [protected]

Definition at line 105 of file CaloHitResponse.h.

Referenced by add(), and setHitFilter().

int CaloHitResponse::theMaxBunch [protected]

Definition at line 112 of file CaloHitResponse.h.

Referenced by run(), and setBunchRange().

int CaloHitResponse::theMinBunch [protected]

Definition at line 111 of file CaloHitResponse.h.

Referenced by run(), and setBunchRange().

const CaloVSimParameterMap* CaloHitResponse::theParameterMap [protected]

Definition at line 102 of file CaloHitResponse.h.

Referenced by makeAnalogSignal(), and makeBlankSignal().

double CaloHitResponse::thePhaseShift_ [protected]

Definition at line 114 of file CaloHitResponse.h.

Referenced by makeAnalogSignal(), and setPhaseShift().

CLHEP::RandPoissonQ* CaloHitResponse::theRandPoisson [mutable, protected]

Definition at line 109 of file CaloHitResponse.h.

Referenced by analogSignalAmplitude(), setRandomEngine(), and ~CaloHitResponse().

const CaloVShape* CaloHitResponse::theShape [protected]

Definition at line 103 of file CaloHitResponse.h.

---

The documentation for this class was generated from the following files:

• SimCalorimetry/CaloSimAlgos/interface/CaloHitResponse.h
• SimCalorimetry/CaloSimAlgos/src/CaloHitResponse.cc

---