HcalSiPMHitResponse Class Reference

#include <HcalSiPMHitResponse.h>

Inheritance diagram for HcalSiPMHitResponse:

Public Member Functions
	HcalSiPMHitResponse (const CaloVSimParameterMap *parameterMap, const CaloShapes *shapes)
virtual void	run (MixCollection< PCaloHit > &hits)
	Complete cell digitization. More...
virtual void	setRandomEngine (CLHEP::HepRandomEngine &engine)
virtual	~HcalSiPMHitResponse ()
Public Member Functions inherited from CaloHitResponse
virtual void	add (const PCaloHit &hit)
	process a single SimHit More...
void	add (const CaloSamples &signal)
	add a signal, in units of pe More...
void	addHit (const PCaloHit *hit, CaloSamples &frame) const
	adds the amplitude for a single hit to the frame More...
double	analogSignalAmplitude (const PCaloHit &hit, const CaloSimParameters &parameters) const
	CaloHitResponse (const CaloVSimParameterMap *parameterMap, const CaloVShape *shape)
	CaloHitResponse (const CaloVSimParameterMap *parameterMap, const CaloShapes *shapes)
void	clear ()
	frees up memory More...
CaloSamples *	findSignal (const DetId &detId)
	users can look for the signal for a given cell More...
virtual bool	keepBlank () const
virtual CaloSamples	makeAnalogSignal (const PCaloHit &inputHit) const
	creates the signal corresponding to this hit More...
CaloSamples	makeBlankSignal (const DetId &detId) const
	creates an empty signal for this DetId More...
int	nSignals () const
	number of signals in the current cache More...
void	setBunchRange (int minBunch, int maxBunch)
	tells it which pileup bunches to do More...
void	setGeometry (const CaloGeometry *geometry)
	geometry needed for time-of-flight More...
void	setHitCorrection (const CaloVHitCorrection *hitCorrection)
	If you want to correct hits, for attenuation or delay, set this. More...
void	setHitFilter (const CaloVHitFilter *filter)
	if you want to reject hits, for example, from a certain subdetector, set this More...
void	setPECorrection (const CaloVPECorrection *peCorrection)
	if you want to correct the photoelectrons More...
void	setPhaseShift (const double &thePhaseShift)
	setting the phase shift for asynchronous trigger (e.g. test beams) More...
double	timeOfFlight (const DetId &detId) const
virtual	~CaloHitResponse ()
	doesn't delete the pointers passed in More...

Protected Member Functions
virtual CaloSamples	makeSiPMSignal (const PCaloHit &inHit, int &integral) const

Protected Attributes
double	theRecoveryTime
HcalSiPM *	theSiPM
Protected Attributes inherited from CaloHitResponse
AnalogSignalMap	theAnalogSignalMap
const CaloGeometry *	theGeometry
const CaloVHitCorrection *	theHitCorrection
const CaloVHitFilter *	theHitFilter
int	theMaxBunch
int	theMinBunch
const CaloVSimParameterMap *	theParameterMap
const CaloVPECorrection *	thePECorrection
double	thePhaseShift_
CLHEP::RandPoissonQ *	theRandPoisson
const CaloVShape *	theShape
const CaloShapes *	theShapes

Additional Inherited Members
Public Types inherited from CaloHitResponse
enum	{ BUNCHSPACE =25 }
typedef std::map< DetId, CaloSamples >	AnalogSignalMap

Detailed Description

Definition at line 9 of file HcalSiPMHitResponse.h.

Constructor & Destructor Documentation

HcalSiPMHitResponse::HcalSiPMHitResponse	(	const CaloVSimParameterMap *	parameterMap,
		const CaloShapes *	shapes
	)

Definition at line 25 of file HcalSiPMHitResponse.cc.

References theSiPM.

                                                :
  CaloHitResponse(parameterMap, shapes), theSiPM(0), theRecoveryTime(250.) {
  theSiPM = new HcalSiPM(14000);
}

HcalSiPMHitResponse::~HcalSiPMHitResponse ( )

virtual

Definition at line 31 of file HcalSiPMHitResponse.cc.

References theSiPM.

                                          {
  if (theSiPM)
    delete theSiPM;
}

Member Function Documentation

CaloSamples HcalSiPMHitResponse::makeSiPMSignal ( const PCaloHit &  inHit, int &  integral  ) const

protectedvirtual

Definition at line 78 of file HcalSiPMHitResponse.cc.

References CaloHitResponse::analogSignalAmplitude(), newFWLiteAna::bin, CaloHitResponse::BUNCHSPACE, CaloVHitCorrection::correct(), HcalSiPM::hitCells(), PCaloHit::id(), errorMatrix2Lands_multiChannel::id, CaloHitResponse::makeBlankSignal(), interactiveExample::photons, query::result, HcalSiPM::setNCells(), CaloShapes::shape(), CaloVSimParameterMap::simParameters(), CaloSamples::size(), CaloHitResponse::theHitCorrection, CaloHitResponse::theParameterMap, CaloHitResponse::thePhaseShift_, CaloHitResponse::theShapes, theSiPM, PCaloHit::time(), CaloHitResponse::timeOfFlight(), and tzero.

Referenced by run().

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

  DetId id(hit.id());
  const HcalSimParameters& pars = dynamic_cast<const HcalSimParameters&>(theParameterMap->simParameters(id));
  theSiPM->setNCells(pars.pixels());

  double signal = analogSignalAmplitude(hit, pars);
  int photons = static_cast<int>(signal + 0.5);
  int pixels = theSiPM->hitCells(photons, integral);
  integral += pixels;
  signal = double(pixels);

  CaloSamples result(makeBlankSignal(id));

  if(pixels > 0)
  {
    const CaloVShape * shape = theShapes->shape(id);
    double jitter = hit.time() - timeOfFlight(id);

    const double tzero = pars.timePhase() - jitter -
      BUNCHSPACE*(pars.binOfMaximum() - thePhaseShift_);
    double binTime = tzero;

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

  return result;
}

void HcalSiPMHitResponse::run ( MixCollection< PCaloHit > &  hits )

virtual

Complete cell digitization.

Reimplemented from CaloHitResponse.

Definition at line 36 of file HcalSiPMHitResponse.cc.

References CaloVHitFilter::accepts(), CaloHitResponse::add(), HcalSiPMRecovery::addToHistory(), MixCollection< T >::begin(), HcalSiPMRecovery::clearHistory(), MixCollection< T >::end(), HcalSiPMRecovery::getIntegral(), i, errorMatrix2Lands_multiChannel::id, edm::detail::isnan(), makeSiPMSignal(), CaloHitResponse::theHitFilter, CaloHitResponse::theMaxBunch, CaloHitResponse::theMinBunch, theRecoveryTime, and PCaloHit::time().

                                                            {

  typedef std::multiset <const PCaloHit *, PCaloHitCompareTimes> SortedHitSet;

  std::map< DetId, SortedHitSet > sortedhits;
  for (MixCollection<PCaloHit>::MixItr hitItr = hits.begin();
       hitItr != hits.end(); ++hitItr) {
    if (!((hitItr.bunch() < theMinBunch) || (hitItr.bunch() > theMaxBunch)) &&
    !(isnan(hitItr->time())) &&
    ((theHitFilter == 0) || (theHitFilter->accepts(*hitItr)))) {
      DetId id(hitItr->id());
      if (sortedhits.find(id)==sortedhits.end())
    sortedhits.insert(std::pair<DetId, SortedHitSet>(id, SortedHitSet()));
      sortedhits[id].insert(&(*hitItr));
    }
  }
  int pixelIntegral, oldIntegral;
  HcalSiPMRecovery pixelHistory(theRecoveryTime);
  const PCaloHit * hit;
  for (std::map<DetId, SortedHitSet>::iterator i = sortedhits.begin(); 
       i!=sortedhits.end(); ++i) {
    pixelHistory.clearHistory();
    for (SortedHitSet::iterator itr = i->second.begin(); 
     itr != i->second.end(); ++itr) {
      hit = *itr;
      pixelIntegral = pixelHistory.getIntegral(hit->time());
      oldIntegral = pixelIntegral;
      CaloSamples signal(makeSiPMSignal(*hit, pixelIntegral));
      pixelHistory.addToHistory(hit->time(), pixelIntegral-oldIntegral);
      add(signal);
    }
  }
}

void HcalSiPMHitResponse::setRandomEngine ( CLHEP::HepRandomEngine &  engine )

virtual

Reimplemented from CaloHitResponse.

Definition at line 71 of file HcalSiPMHitResponse.cc.

References HcalSiPM::initRandomEngine(), CaloHitResponse::setRandomEngine(), and theSiPM.

{
  theSiPM->initRandomEngine(engine);
  CaloHitResponse::setRandomEngine(engine);
}

Member Data Documentation

double HcalSiPMHitResponse::theRecoveryTime

protected

Definition at line 26 of file HcalSiPMHitResponse.h.

Referenced by run().

HcalSiPM* HcalSiPMHitResponse::theSiPM

protected

Definition at line 25 of file HcalSiPMHitResponse.h.

Referenced by HcalSiPMHitResponse(), makeSiPMSignal(), setRandomEngine(), and ~HcalSiPMHitResponse().