CaloHitResponse Class Reference

Creates electronics signals from hits. More...

#include <CaloHitResponse.h>

Inheritance diagram for CaloHitResponse:

Public Types
enum	{ BUNCHSPACE =25 }
typedef std::map< DetId, CaloSamples >	AnalogSignalMap

Public Member Functions
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 DetId &id, float energy, const CaloSimParameters &parameters) const
	CaloHitResponse (const CaloVSimParameterMap *parameterMap, const CaloVShape *shape)
	CaloHitResponse (const CaloVSimParameterMap *parameterMap, const CaloShapes *shapes)
void	clear ()
	frees up memory More...
virtual void	finalizeHits ()
	Finalize hits. More...
CaloSamples *	findSignal (const DetId &detId)
	users can look for the signal for a given cell More...
void	initHBHEScale ()
virtual void	initializeHits ()
	Initialize hits. 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...
virtual void	run (MixCollection< PCaloHit > &hits)
	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	setHBHEScale (std::string &)
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...
virtual void	setRandomEngine (CLHEP::HepRandomEngine &engine)
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 Attributes
bool	changeScale
float	hcal_en_scale [100][72][4]
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

Detailed Description

Creates electronics signals from hits.

Definition at line 31 of file CaloHitResponse.h.

Member Typedef Documentation

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

Definition at line 34 of file CaloHitResponse.h.

Member Enumeration Documentation

anonymous enum

Enumerator
BUNCHSPACE

Definition at line 36 of file CaloHitResponse.h.

{BUNCHSPACE=25};

Constructor & Destructor Documentation

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

Definition at line 25 of file CaloHitResponse.cc.

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

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

Definition at line 41 of file CaloHitResponse.cc.

: theAnalogSignalMap(),
  theParameterMap(parametersMap),
  theShapes(shapes),
  theShape(0),
  theHitCorrection(0),
  thePECorrection(0),
  theHitFilter(0),
  theGeometry(0),
  theRandPoisson(0),
  theMinBunch(-10),
  theMaxBunch(10),
  thePhaseShift_(1.),
  changeScale(false) {}

CaloHitResponse::~CaloHitResponse ( )

virtual

doesn't delete the pointers passed in

Definition at line 57 of file CaloHitResponse.cc.

References theRandPoisson.

                                  {
  delete theRandPoisson;
}

Member Function Documentation

void CaloHitResponse::add ( const PCaloHit &  hit )

virtual

process a single SimHit

Reimplemented in HcalSiPMHitResponse.

Definition at line 116 of file CaloHitResponse.cc.

References CaloVHitFilter::accepts(), i, edm::isNotFinite(), keep, keepBlank(), LogDebug, makeAnalogSignal(), CaloSamples::size(), findQualityFiles::size, theHitFilter, and PCaloHit::time().

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

                                               {
  // check the hit time makes sense
  if ( edm::isNotFinite(hit.time()) ) { return; }

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

    bool keep ( keepBlank() ) ;  // here we  check for blank signal if not keeping them
    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("CaloHitResponse") << signal;
    if( keep ) add(signal);
  }
}

void CaloHitResponse::add

(

const CaloSamples &

signal

)

add a signal, in units of pe

Definition at line 143 of file CaloHitResponse.cc.

References findSignal(), i, CaloSamples::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::addHit	(	const PCaloHit *	hit,
		CaloSamples &	frame
	)		const

adds the amplitude for a single hit to the frame

double CaloHitResponse::analogSignalAmplitude	(	const DetId &	id,
		float	energy,
		const CaloSimParameters &	parameters
	)		const

finds the amplitude contribution from this hit, applying photostatistics, if needed. Results are in photoelectrons

Definition at line 197 of file CaloHitResponse.cc.

References changeScale, CaloVPECorrection::correctPE(), HcalDetId::depth(), DetId::det(), CaloSimParameters::doPhotostatistics(), edm::hlt::Exception, edm::RandomNumberGenerator::getEngine(), DetId::Hcal, hcal_en_scale, HcalBarrel, HcalEndcap, HcalDetId::ieta(), HcalDetId::iphi(), edm::Service< T >::isAvailable(), LogDebug, CaloSimParameters::simHitToPhotoelectrons(), HcalDetId::subdet(), thePECorrection, and theRandPoisson.

Referenced by HcalSiPMHitResponse::add(), makeAnalogSignal(), and HcalSiPMHitResponse::makeSiPMSignal().

                                                                                                                           {

  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
  //GMA Smeared in photon production it self  
  double scl =1.0;
#ifdef ChangeHcalEnergyScale
  if (changeScale) {
    if (detId.det()==DetId::Hcal ) { 
      HcalDetId dId = HcalDetId(detId); 
      if (dId.subdet()==HcalBarrel || dId.subdet()==HcalEndcap) { 
    int ieta = dId.ieta()+50;
    int iphi = dId.iphi()-1;
    int idep = dId.depth()-1;
    scl = hcal_en_scale[ieta][iphi][idep];
    LogDebug("CaloHitResponse") << " ID " << dId << " Scale " << scl;
      }
    }
  } 
#endif
  double npe = scl * energy * parameters.simHitToPhotoelectrons(detId);
  // do we need to doPoisson statistics for the photoelectrons?
  if(parameters.doPhotostatistics()) {
    npe = theRandPoisson->fire(npe);
  }
  if(thePECorrection) npe = thePECorrection->correctPE(detId, npe);
  return npe;
}

void CaloHitResponse::clear ( void  )

inline

frees up memory

Definition at line 89 of file CaloHitResponse.h.

References theAnalogSignalMap.

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

{theAnalogSignalMap.clear();}

virtual void CaloHitResponse::finalizeHits ( )

inlinevirtual

Finalize hits.

Reimplemented in HcalSiPMHitResponse.

Definition at line 60 of file CaloHitResponse.h.

Referenced by CaloTDigitizer< Traits >::run().

{}

CaloSamples * CaloHitResponse::findSignal

(

const DetId &

detId

)

users can look for the signal for a given cell

Definition at line 239 of file CaloHitResponse.cc.

References query::result, and theAnalogSignalMap.

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

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

void CaloHitResponse::initHBHEScale

(

)

Definition at line 61 of file CaloHitResponse.cc.

References hcal_en_scale.

Referenced by HcalDigitizer::HcalDigitizer().

                                    {
#ifdef ChangeHcalEnergyScale
  for (int ij=0; ij<100; ij++) {
    for (int jk=0; jk<72; jk++) {   
      for (int kl=0; kl<4; kl++) {
    hcal_en_scale[ij][jk][kl] = 1.0;
      }
    }
  }
#endif
}

virtual void CaloHitResponse::initializeHits ( )

inlinevirtual

Initialize hits.

Reimplemented in HcalSiPMHitResponse.

Definition at line 57 of file CaloHitResponse.h.

Referenced by CaloTDigitizer< Traits >::initializeHits().

{}

virtual bool CaloHitResponse::keepBlank ( ) const

inlinevirtual

Definition at line 54 of file CaloHitResponse.h.

References funct::true.

Referenced by add(), and HcalSiPMHitResponse::finalizeHits().

{ return true ; }

CaloSamples CaloHitResponse::makeAnalogSignal ( const PCaloHit &  inputHit ) const

virtual

creates the signal corresponding to this hit

Definition at line 162 of file CaloHitResponse.cc.

References analogSignalAmplitude(), newFWLiteAna::bin, BUNCHSPACE, CaloVHitCorrection::delay(), PCaloHit::energy(), PCaloHit::id(), makeBlankSignal(), Parameters::parameters, query::result, CaloShapes::shape(), CaloVSimParameterMap::simParameters(), CaloSamples::size(), theHitCorrection, theParameterMap, thePhaseShift_, theShape, theShapes, PCaloHit::time(), cond::rpcobgas::time, timeOfFlight(), CaloVShape::timeToRise(), and tzero.

Referenced by add().

                                                                        {

  DetId detId(hit.id());
  const CaloSimParameters & parameters = theParameterMap->simParameters(detId);
  
  double signal = analogSignalAmplitude(detId, hit.energy(), parameters);

  double time = hit.time();
  if(theHitCorrection != 0) {
    time += theHitCorrection->delay(hit);
  }
  double jitter = hit.time() - timeOfFlight(detId);

  const CaloVShape * shape = theShape;
  if(!shape) {
    shape = theShapes->shape(detId);
  }
  // assume bins count from zero, go for center of bin
  const double tzero = ( shape->timeToRise()
             + parameters.timePhase() 
             - jitter 
             - BUNCHSPACE*( parameters.binOfMaximum()
                    - thePhaseShift_          ) ) ;
  double binTime = tzero;

  CaloSamples result(makeBlankSignal(detId));

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

CaloSamples CaloHitResponse::makeBlankSignal

(

const DetId &

detId

)

const

creates an empty signal for this DetId

Definition at line 251 of file CaloHitResponse.cc.

References CaloSimParameters::binOfMaximum(), Parameters::parameters, CaloSimParameters::readoutFrameSize(), query::result, CaloVSimParameterMap::simParameters(), and theParameterMap.

Referenced by makeAnalogSignal(), and CaloTDigitizer< 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 105 of file CaloHitResponse.h.

References theAnalogSignalMap.

Referenced by CaloTDigitizer< Traits >::run().

{return theAnalogSignalMap.size();}

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

virtual

Complete cell digitization.

Reimplemented in HcalSiPMHitResponse.

Definition at line 106 of file CaloHitResponse.cc.

References add(), MixCollection< T >::begin(), MixCollection< T >::end(), and withinBunchRange().

Referenced by EcalMixingModuleValidation::computeSDBunchDigi().

                                                        {

  for(MixCollection<PCaloHit>::MixItr hitItr = hits.begin();
      hitItr != hits.end(); ++hitItr) {
    if(withinBunchRange(hitItr.bunch())) {
      add(*hitItr);
    } // loop over hits
  }
}

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

tells it which pileup bunches to do

Definition at line 95 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 52 of file CaloHitResponse.h.

References geometry, and theGeometry.

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

{ theGeometry = geometry; }

void CaloHitResponse::setHBHEScale

(

std::string &

fileIn

)

Definition at line 73 of file CaloHitResponse.cc.

References changeScale, eta(), hcal_en_scale, EdgesToViz::infile, LogDebug, and phi.

Referenced by HcalDigitizer::HcalDigitizer().

                                                     {
  
  ifstream infile(fileIn.c_str());
  LogDebug("CaloHitResponse") << "Reading from " << fileIn;
#ifdef ChangeHcalEnergyScale
  if (!infile.is_open()) {
    edm::LogError("CaloHitResponse") << "** ERROR: Can't open '" << fileIn << "' for the input file";
  } else {
    int     eta, phi, depth;
    double  cFactor;
    while(1) {
      infile >> eta >> phi >> depth >> cFactor;
      if (!infile.good()) break;
      hcal_en_scale[eta][phi][depth] = cFactor;
      //      LogDebug("CaloHitResponse") << "hcal_en_scale[" << eta << "][" << phi << "][" << depth << "] = " << hcal_en_scale[eta][phi][depth];
    }
    infile.close();
  }
  changeScale = true;
#endif
}

void CaloHitResponse::setHitCorrection ( const CaloVHitCorrection *  hitCorrection )

inline

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

Definition at line 77 of file CaloHitResponse.h.

References theHitCorrection.

Referenced by CastorDigiProducer::CastorDigiProducer(), 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 72 of file CaloHitResponse.h.

References alcazmumu_cfi::filter, and theHitFilter.

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

                                                   {
    theHitFilter = filter;
  }

void CaloHitResponse::setPECorrection ( const CaloVPECorrection *  peCorrection )

inline

if you want to correct the photoelectrons

Definition at line 82 of file CaloHitResponse.h.

References thePECorrection.

Referenced by HcalDigitizer::HcalDigitizer().

                                                               {
    thePECorrection = peCorrection;
  }

void CaloHitResponse::setPhaseShift ( const double &  thePhaseShift )

inline

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

Definition at line 116 of file CaloHitResponse.h.

References thePhaseShift_.

Referenced by HcalTBDigiProducer::setPhaseShift().

{ thePhaseShift_ = thePhaseShift; }

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

virtual

Reimplemented in HcalSiPMHitResponse.

Definition at line 101 of file CaloHitResponse.cc.

References theRandPoisson.

Referenced by HcalSiPMHitResponse::setRandomEngine(), and CaloTDigitizer< Traits >::setRandomEngine().

                                                                   {
  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 259 of file CaloHitResponse.cc.

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

Referenced by HcalSiPMHitResponse::add(), makeAnalogSignal(), and HcalSiPMHitResponse::makeSiPMSignal().

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

bool CaloHitResponse::withinBunchRange ( int  bunchCrossing ) const

inline

check if crossing is within bunch range:

Definition at line 120 of file CaloHitResponse.h.

References theMaxBunch, and theMinBunch.

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

                                                 {
    return(bunchCrossing >= theMinBunch && bunchCrossing <= theMaxBunch);
  }

Member Data Documentation

bool CaloHitResponse::changeScale

protected

Definition at line 145 of file CaloHitResponse.h.

Referenced by analogSignalAmplitude(), and setHBHEScale().

float CaloHitResponse::hcal_en_scale[100][72][4]

protected

Definition at line 147 of file CaloHitResponse.h.

Referenced by analogSignalAmplitude(), initHBHEScale(), and setHBHEScale().

AnalogSignalMap CaloHitResponse::theAnalogSignalMap

protected

Definition at line 126 of file CaloHitResponse.h.

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

const CaloGeometry* CaloHitResponse::theGeometry

protected

Definition at line 135 of file CaloHitResponse.h.

Referenced by setGeometry(), and timeOfFlight().

const CaloVHitCorrection* CaloHitResponse::theHitCorrection

protected

Definition at line 131 of file CaloHitResponse.h.

Referenced by HcalSiPMHitResponse::add(), makeAnalogSignal(), HcalSiPMHitResponse::makeSiPMSignal(), and setHitCorrection().

const CaloVHitFilter* CaloHitResponse::theHitFilter

protected

Definition at line 133 of file CaloHitResponse.h.

Referenced by HcalSiPMHitResponse::add(), add(), HcalSiPMHitResponse::run(), and setHitFilter().

int CaloHitResponse::theMaxBunch

protected

Definition at line 140 of file CaloHitResponse.h.

Referenced by HcalSiPMHitResponse::run(), setBunchRange(), and withinBunchRange().

int CaloHitResponse::theMinBunch

protected

Definition at line 139 of file CaloHitResponse.h.

Referenced by HcalSiPMHitResponse::run(), setBunchRange(), and withinBunchRange().

const CaloVSimParameterMap* CaloHitResponse::theParameterMap

protected

Definition at line 128 of file CaloHitResponse.h.

Referenced by HcalSiPMHitResponse::add(), makeAnalogSignal(), HcalSiPMHitResponse::makeBlankSignal(), makeBlankSignal(), and HcalSiPMHitResponse::makeSiPMSignal().

const CaloVPECorrection* CaloHitResponse::thePECorrection

protected

Definition at line 132 of file CaloHitResponse.h.

Referenced by analogSignalAmplitude(), and setPECorrection().

double CaloHitResponse::thePhaseShift_

protected

Definition at line 142 of file CaloHitResponse.h.

Referenced by HcalSiPMHitResponse::add(), makeAnalogSignal(), HcalSiPMHitResponse::makeSiPMSignal(), and setPhaseShift().

CLHEP::RandPoissonQ* CaloHitResponse::theRandPoisson

mutableprotected

Definition at line 137 of file CaloHitResponse.h.

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

const CaloVShape* CaloHitResponse::theShape

protected

Definition at line 130 of file CaloHitResponse.h.

Referenced by makeAnalogSignal().

const CaloShapes* CaloHitResponse::theShapes

protected

Definition at line 129 of file CaloHitResponse.h.

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