#include <HcalSiPMHitResponse.h>

Inheritance diagram for HcalSiPMHitResponse:

Public Types
typedef std::vector< unsigned int >	photonTimeHist

typedef std::map< DetId, photonTimeHist >	photonTimeMap

Public Types inherited from CaloHitResponse
enum	{ BUNCHSPACE =25 }

typedef std::map< DetId, CaloSamples >	AnalogSignalMap

Public Member Functions
virtual void	add (const PCaloHit &hit, CLHEP::HepRandomEngine *) override
	process a single SimHit More...

virtual void	add (const CaloSamples &signal)

virtual void	addPEnoise (CLHEP::HepRandomEngine *engine)

virtual void	finalizeHits (CLHEP::HepRandomEngine *) override
	Finalize hits. More...

	HcalSiPMHitResponse (const CaloVSimParameterMap parameterMap, const CaloShapes shapes, bool PreMix1=false)

virtual void	initializeHits () override
	Initialize hits. More...

virtual CaloSamples	makeBlankSignal (const DetId &detId) const

virtual void	setDetIds (const std::vector< DetId > &detIds)

virtual	~HcalSiPMHitResponse ()

Public Member Functions inherited from CaloHitResponse
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 DetId &id, float energy, const CaloSimParameters &parameters, CLHEP::HepRandomEngine *) 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, CLHEP::HepRandomEngine *) 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...

virtual void	run (MixCollection< PCaloHit > &hits, CLHEP::HepRandomEngine *)
	Complete cell digitization. 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	setIgnoreGeantTime (bool gt)

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...

void	setStorePrecise (bool sp)

double	timeOfFlight (const DetId &detId) const

bool	withinBunchRange (int bunchCrossing) const
	check if crossing is within bunch range: More...

virtual	~CaloHitResponse ()
	doesn't delete the pointers passed in More...

Protected Member Functions
virtual CaloSamples	makeSiPMSignal (DetId const &id, photonTimeHist const &photons, CLHEP::HepRandomEngine *)

Private Attributes
double	dt

double	invdt

int	nbins

photonTimeMap	precisionTimedPhotons

bool	PreMixDigis

const std::vector< DetId > *	theDetIds

HcalSiPM	theSiPM

Additional Inherited Members
Protected Attributes inherited from CaloHitResponse
bool	ignoreTime

bool	storePrecise

AnalogSignalMap	theAnalogSignalMap

const CaloGeometry *	theGeometry

const CaloVHitCorrection *	theHitCorrection

const CaloVHitFilter *	theHitFilter

int	theMaxBunch

int	theMinBunch

const CaloVSimParameterMap *	theParameterMap

const CaloVPECorrection *	thePECorrection

double	thePhaseShift_

const CaloVShape *	theShape

const CaloShapes *	theShapes

Detailed Description

Definition at line 24 of file HcalSiPMHitResponse.h.

Member Typedef Documentation

typedef std::vector<unsigned int> HcalSiPMHitResponse::photonTimeHist

Definition at line 32 of file HcalSiPMHitResponse.h.

typedef std::map< DetId, photonTimeHist > HcalSiPMHitResponse::photonTimeMap

Definition at line 33 of file HcalSiPMHitResponse.h.

Constructor & Destructor Documentation

HcalSiPMHitResponse::HcalSiPMHitResponse	(	const CaloVSimParameterMap *	parameterMap,
		const CaloShapes *	shapes,
		bool	PreMix1 = `false`
	)

Definition at line 20 of file HcalSiPMHitResponse.cc.

                                                               :
   CaloHitResponse(parameterMap, shapes), theSiPM(), PreMixDigis(PreMix1),
   nbins(BUNCHSPACE*HcalPulseShapes::invDeltaTSiPM_), dt(HcalPulseShapes::deltaTSiPM_), invdt(HcalPulseShapes::invDeltaTSiPM_) {}

HcalSiPMHitResponse::~HcalSiPMHitResponse ( )

virtual

Definition at line 25 of file HcalSiPMHitResponse.cc.

25 {}

Member Function Documentation

void HcalSiPMHitResponse::add	(	const PCaloHit &	hit,
		CLHEP::HepRandomEngine *	engine
	)

overridevirtual

process a single SimHit

Reimplemented from CaloHitResponse.

Definition at line 68 of file HcalSiPMHitResponse.cc.

References CaloVHitFilter::accepts(), CaloHitResponse::analogSignalAmplitude(), CaloHitResponse::BUNCHSPACE, dt, PCaloHit::energy(), HcalPulseShapes::generatePhotonTime(), PCaloHit::id(), CaloHitResponse::ignoreTime, invdt, edm::isNotFinite(), LogDebug, nbins, interactiveExample::photons, precisionTimedPhotons, CaloVSimParameterMap::simParameters(), findQualityFiles::size, CaloHitResponse::theHitFilter, CaloHitResponse::theParameterMap, CaloHitResponse::thePhaseShift_, PCaloHit::time(), CaloHitResponse::timeOfFlight(), and tzero.

Referenced by finalizeHits(), and counter.Counter::register().

                                                                                {
     if (!edm::isNotFinite(hit.time()) &&
     ((theHitFilter == 0) || (theHitFilter->accepts(hit)))) {
       HcalDetId id(hit.id());
       const HcalSimParameters& pars = dynamic_cast<const HcalSimParameters&>(theParameterMap->simParameters(id));
       //divide out mean of crosstalk distribution 1/(1-lambda) = multiply by (1-lambda)
       double signal(analogSignalAmplitude(id, hit.energy(), pars, engine)*(1-pars.sipmCrossTalk(id)));
       unsigned int photons(signal + 0.5);
       double tof( timeOfFlight(id) );
       double time( hit.time() );
       if(ignoreTime) time = tof;
 
       if (photons > 0)
         if (precisionTimedPhotons.find(id)==precisionTimedPhotons.end()) {
           precisionTimedPhotons.insert(
             std::pair<DetId, photonTimeHist >(id, 
               photonTimeHist(nbins * pars.readoutFrameSize(), 0)
             )
           );
         }
 
       LogDebug("HcalSiPMHitResponse") << id;
       LogDebug("HcalSiPMHitResponse") << " fCtoGeV: " << pars.fCtoGeV(id)
         << " samplingFactor: " << pars.samplingFactor(id)
         << " photoelectronsToAnalog: " << pars.photoelectronsToAnalog(id)
         << " simHitToPhotoelectrons: " << pars.simHitToPhotoelectrons(id);
       LogDebug("HcalSiPMHitResponse") << " energy: " << hit.energy()
         << " photons: " << photons 
         << " time: " << time;
       LogDebug("HcalSiPMHitResponse") << " timePhase: " << pars.timePhase()
         << " tof: " << tof
         << " binOfMaximum: " << pars.binOfMaximum()
         << " phaseShift: " << thePhaseShift_;
       double tzero(0.0 + pars.timePhase() - 
            (time - tof) - 
            BUNCHSPACE*( pars.binOfMaximum() - thePhaseShift_));
       LogDebug("HcalSiPMHitResponse") << " tzero: " << tzero;
       double tzero_bin(-tzero*invdt);
       LogDebug("HcalSiPMHitResponse") << " corrected tzero: " << tzero_bin << '\n';
       double t_pe(0.);
       int t_bin(0);
       for (unsigned int pe(0); pe<photons; ++pe) {
         t_pe = HcalPulseShapes::generatePhotonTime(engine);
         t_bin = int(t_pe*invdt + tzero_bin + 0.5);
         LogDebug("HcalSiPMHitResponse") << "t_pe: " << t_pe << " t_pe + tzero: " << (t_pe+tzero_bin*dt)
                   << " t_bin: " << t_bin << '\n';
         if ((t_bin >= 0) && 
             (static_cast<unsigned int>(t_bin) < precisionTimedPhotons[id].size()))
             precisionTimedPhotons[id][t_bin] += 1;
       }
     }
 }

void HcalSiPMHitResponse::add ( const CaloSamples & signal )

virtual

Definition at line 58 of file HcalSiPMHitResponse.cc.

References CaloHitResponse::findSignal(), CaloSamples::id(), and CaloHitResponse::theAnalogSignalMap.

Referenced by counter.Counter::register().

                                                        {
   DetId id(signal.id());
   CaloSamples * oldSignal = findSignal(id);
   if (oldSignal == 0) {
     theAnalogSignalMap[id] = signal;
   } else {
     (*oldSignal) += signal;
   }
 }

void HcalSiPMHitResponse::addPEnoise ( CLHEP::HepRandomEngine * engine )

virtual

Definition at line 121 of file HcalSiPMHitResponse.cc.

References dt, LogDebug, nbins, HcalSimParameters::photoelectronsToAnalog(), interactiveExample::photons, precisionTimedPhotons, CaloSimParameters::readoutFrameSize(), CaloVSimParameterMap::simParameters(), HcalSimParameters::sipmDarkCurrentuA(), theDetIds, and CaloHitResponse::theParameterMap.

Referenced by finalizeHits().

 {
   // Add SiPM dark current noise to all cells
   for(std::vector<DetId>::const_iterator idItr = theDetIds->begin();
       idItr != theDetIds->end(); ++idItr) {
     HcalDetId id(*idItr);
     const HcalSimParameters& pars =
       static_cast<const HcalSimParameters&>(theParameterMap->simParameters(id));
 
     // uA * ns / (fC/pe) = pe!
     double dc_pe_avg =
       pars.sipmDarkCurrentuA(id) * dt / 
       pars.photoelectronsToAnalog(id);
 
     if (dc_pe_avg <= 0.) continue;
 
     int nPreciseBins = nbins * pars.readoutFrameSize();
 
     unsigned int sumnoisePE(0);
     double  elapsedTime(0.);
     for (int tprecise(0); tprecise < nPreciseBins; ++tprecise) {
       int noisepe = CLHEP::RandPoissonQ::shoot(engine, dc_pe_avg); // add dark current noise
 
       if (noisepe > 0) {
     if (precisionTimedPhotons.find(id)==precisionTimedPhotons.end()) {
       photonTimeHist photons(nPreciseBins, 0);
       photons[tprecise] = noisepe;
       precisionTimedPhotons.insert
         (std::pair<DetId, photonTimeHist >(id, photons ) );
     } else {
       precisionTimedPhotons[id][tprecise] += noisepe;
     }
 
     sumnoisePE += noisepe;
       }
       elapsedTime += dt;
 
     } // precise time loop
 
     LogDebug("HcalSiPMHitResponse") << id;
     LogDebug("HcalSiPMHitResponse") << " total noise (PEs): " << sumnoisePE;
 
   } // detId loop
 }                                               // HcalSiPMHitResponse::addPEnoise()

void HcalSiPMHitResponse::finalizeHits ( CLHEP::HepRandomEngine * )

overridevirtual

Finalize hits.

Reimplemented from CaloHitResponse.

Definition at line 31 of file HcalSiPMHitResponse.cc.

References add(), addPEnoise(), i, CaloSamples::id(), keep, CaloHitResponse::keepBlank(), LogDebug, makeSiPMSignal(), precisionTimedPhotons, PreMixDigis, CaloSamples::size(), and findQualityFiles::size.

                                                                    {
   //do not add PE noise for initial premix
   if(!PreMixDigis) addPEnoise(engine);
 
   photonTimeMap::iterator channelPhotons;
   for (channelPhotons = precisionTimedPhotons.begin();
        channelPhotons != precisionTimedPhotons.end();
        ++channelPhotons) {
     CaloSamples signal(makeSiPMSignal(channelPhotons->first, 
                       channelPhotons->second,
                                       engine));
     bool keep( keepBlank() );
     if (!keep) {
       const unsigned int size ( signal.size() ) ;
       if( 0 != size ) {
     for( unsigned int i ( 0 ) ; i != size ; ++i ) {
       keep = keep || signal[i] > 1.e-7 ;
     }
       }
     }
 
     LogDebug("HcalSiPMHitResponse") << HcalDetId(signal.id()) << ' ' << signal;
 
     if (keep) add(signal);
   }
 }

void HcalSiPMHitResponse::initializeHits ( )

overridevirtual

Initialize hits.

Reimplemented from CaloHitResponse.

Definition at line 27 of file HcalSiPMHitResponse.cc.

References precisionTimedPhotons.

                                          {
   precisionTimedPhotons.clear();
 }

CaloSamples HcalSiPMHitResponse::makeBlankSignal ( const DetId & detId ) const

virtual

Definition at line 166 of file HcalSiPMHitResponse.cc.

References CaloSimParameters::binOfMaximum(), dt, nbins, HLT_FULL_cff::parameters, CaloSimParameters::readoutFrameSize(), mps_fire::result, CaloVSimParameterMap::simParameters(), and CaloHitResponse::theParameterMap.

Referenced by makeSiPMSignal().

                                                                          {
   const CaloSimParameters & parameters = theParameterMap->simParameters(detId);
   int preciseSize(parameters.readoutFrameSize() * nbins);
   CaloSamples result(detId, parameters.readoutFrameSize(), preciseSize);
   result.setPresamples(parameters.binOfMaximum()-1);
   result.setPrecise(result.presamples() * nbins, dt);
   return result;
 }

CaloSamples HcalSiPMHitResponse::makeSiPMSignal	(	DetId const &	id,
		photonTimeHist const &	photons,
		CLHEP::HepRandomEngine *	engine
	)

protectedvirtual

Definition at line 175 of file HcalSiPMHitResponse.cc.

References HcalSimParameters::doSiPMSmearing(), dt, alignCSCRings::e, HcalSiPM::hitCells(), invdt, LogDebug, makeBlankSignal(), nbins, HcalSimParameters::pixels(), CaloSamples::preciseAtMod(), CaloSamples::preciseSize(), CaloSamples::resetPrecise(), HcalSiPM::setCrossTalk(), HcalSiPM::setNCells(), HcalSiPM::setSaturationPars(), HcalSiPM::setTau(), HcalSimParameters::signalShape(), CaloVSimParameterMap::simParameters(), HcalSimParameters::sipmCrossTalk(), HcalSimParameters::sipmNonlinearity(), HcalSimParameters::sipmTau(), CaloHitResponse::theParameterMap, and theSiPM.

Referenced by finalizeHits().

                                                                               {
   const HcalSimParameters& pars = static_cast<const HcalSimParameters&>(theParameterMap->simParameters(id));  
   theSiPM.setNCells(pars.pixels(id));
   theSiPM.setTau(pars.sipmTau());
   theSiPM.setCrossTalk(pars.sipmCrossTalk(id));
   theSiPM.setSaturationPars(pars.sipmNonlinearity(id));
 
   //use to make signal
   CaloSamples signal( makeBlankSignal(id) );
   int sampleBin(0), preciseBin(0);
   signal.resetPrecise();
   unsigned int pe(0);
   double hitPixels(0.), elapsedTime(0.);
   unsigned int sumPE(0);
   double sumHits(0.);
 
   HcalSiPMShape sipmPulseShape(pars.signalShape(id));
 
   std::list< std::pair<double, double> > pulses;
   std::list< std::pair<double, double> >::iterator pulse;
   double timeDiff, pulseBit;
   LogDebug("HcalSiPMHitResponse") << "makeSiPMSignal for " << HcalDetId(id);
 
   for (unsigned int tbin(0); tbin < photonTimeBins.size(); ++tbin) {
     pe = photonTimeBins[tbin];
     sumPE += pe;
     preciseBin = tbin;
     sampleBin = preciseBin/nbins;
     if (pe > 0) {
       hitPixels = theSiPM.hitCells(engine, pe, 0., elapsedTime);
       sumHits += hitPixels;
       LogDebug("HcalSiPMHitResponse") << " elapsedTime: " << elapsedTime
                       << " sampleBin: " << sampleBin
                       << " preciseBin: " << preciseBin
                       << " pe: " << pe 
                       << " hitPixels: " << hitPixels ;
       if (pars.doSiPMSmearing()) {
     pulses.push_back( std::pair<double, double>(elapsedTime, hitPixels) );
       } else {
     signal[sampleBin] += hitPixels;
     hitPixels *= invdt;
     signal.preciseAtMod(preciseBin) += 0.6*hitPixels;
     if (preciseBin > 0)
       signal.preciseAtMod(preciseBin-1) += 0.2*hitPixels;
     if (preciseBin < signal.preciseSize() -1)
       signal.preciseAtMod(preciseBin+1) += 0.2*hitPixels;
       }
     }
     
     if (pars.doSiPMSmearing()) {
       pulse = pulses.begin();
       while (pulse != pulses.end()) {
     timeDiff = elapsedTime - pulse->first;
     pulseBit = sipmPulseShape(timeDiff)*pulse->second;
     LogDebug("HcalSiPMHitResponse") << " pulse t: " << pulse->first 
                     << " pulse A: " << pulse->second
                     << " timeDiff: " << timeDiff
                     << " pulseBit: " << pulseBit;
     signal[sampleBin] += pulseBit;
     signal.preciseAtMod(preciseBin) += pulseBit*invdt;
 
     if (timeDiff > 1 && sipmPulseShape(timeDiff) < 1e-7)
       pulse = pulses.erase(pulse);
     else
       ++pulse;
       }
     }
     elapsedTime += dt;
   }
 
   return signal;
 }

void HcalSiPMHitResponse::setDetIds ( const std::vector< DetId > & detIds )

virtual

Definition at line 250 of file HcalSiPMHitResponse.cc.

References theDetIds.

                                                                    {
   theDetIds = &detIds;
 }

Member Data Documentation

double HcalSiPMHitResponse::dt

private

Definition at line 56 of file HcalSiPMHitResponse.h.

Referenced by add(), addPEnoise(), makeBlankSignal(), and makeSiPMSignal().

double HcalSiPMHitResponse::invdt

private

Definition at line 56 of file HcalSiPMHitResponse.h.

Referenced by add(), and makeSiPMSignal().

int HcalSiPMHitResponse::nbins

private

Definition at line 55 of file HcalSiPMHitResponse.h.

Referenced by add(), addPEnoise(), makeBlankSignal(), and makeSiPMSignal().

photonTimeMap HcalSiPMHitResponse::precisionTimedPhotons

private

Definition at line 58 of file HcalSiPMHitResponse.h.

Referenced by add(), addPEnoise(), finalizeHits(), and initializeHits().

bool HcalSiPMHitResponse::PreMixDigis

private

Definition at line 54 of file HcalSiPMHitResponse.h.

Referenced by finalizeHits().

const std::vector<DetId>* HcalSiPMHitResponse::theDetIds

private

Definition at line 60 of file HcalSiPMHitResponse.h.

Referenced by addPEnoise(), and setDetIds().

HcalSiPM HcalSiPMHitResponse::theSiPM

private

Definition at line 53 of file HcalSiPMHitResponse.h.

Referenced by makeSiPMSignal().

Public Types

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation