HcalHitCorrection Class Reference

#include <HcalHitCorrection.h>

Inheritance diagram for HcalHitCorrection:

Public Types
typedef std::map< DetId, double >	ChargeSumsByChannel

Public Member Functions
double	charge (const PCaloHit &hit) const
	how much charge we expect from this hit More...
void	clear ()
virtual double	delay (const PCaloHit &hit, CLHEP::HepRandomEngine *) const
	how much delay this hit will get More...
void	fillChargeSums (MixCollection< PCaloHit > &hits)
void	fillChargeSums (std::vector< PCaloHit > &hits)
	HcalHitCorrection (const CaloVSimParameterMap *parameterMap)
int	timeBin (const PCaloHit &hit) const
	which time bin the peak of the signal will fall in More...
double	timeOfFlight (const DetId &id) const
	simple average approximation More...
virtual	~HcalHitCorrection ()

Private Attributes
ChargeSumsByChannel	theChargeSumsForTimeBin [10]
const CaloVSimParameterMap *	theParameterMap

Detailed Description

Definition at line 22 of file HcalHitCorrection.h.

Member Typedef Documentation

typedef std::map<DetId, double> HcalHitCorrection::ChargeSumsByChannel

Definition at line 25 of file HcalHitCorrection.h.

Constructor & Destructor Documentation

HcalHitCorrection::HcalHitCorrection

(

const CaloVSimParameterMap *

parameterMap

)

Definition at line 13 of file HcalHitCorrection.cc.

  : theParameterMap(parameterMap)
{
}

virtual HcalHitCorrection::~HcalHitCorrection ( )

inlinevirtual

Definition at line 28 of file HcalHitCorrection.h.

{}

Member Function Documentation

double HcalHitCorrection::charge

(

const PCaloHit &

hit

)

const

how much charge we expect from this hit

Definition at line 69 of file HcalHitCorrection.cc.

References PCaloHit::energy(), PCaloHit::id(), Parameters::parameters, CaloSimParameters::simHitToPhotoelectrons(), CaloVSimParameterMap::simParameters(), and theParameterMap.

Referenced by delay(), and fillChargeSums().

{
  HcalDetId detId(hit.id());
  const CaloSimParameters & parameters = theParameterMap->simParameters(detId);
  double simHitToCharge = parameters.simHitToPhotoelectrons(detId)
                        * parameters.photoelectronsToAnalog(detId);
  return hit.energy() * simHitToCharge;
}

void HcalHitCorrection::clear

(

void

)

Definition at line 61 of file HcalHitCorrection.cc.

References i, and theChargeSumsForTimeBin.

Referenced by HcalDigitizer::finalizeEvent(), and HcalDigitizer::initializeEvent().

{
  for(int i = 0; i < 10; ++i)
  {
    theChargeSumsForTimeBin[i].clear();
  }
}

double HcalHitCorrection::delay ( const PCaloHit &  hit, CLHEP::HepRandomEngine *  engine  ) const

virtual

how much delay this hit will get

Implements CaloVHitCorrection.

Definition at line 79 of file HcalHitCorrection.cc.

References DetId::Calo, charge(), HcalTimeSlew::delay(), HcalSimParameters::doTimeSmear(), end, Exception, HcalBarrel, HcalEndcap, HcalOuter, PCaloHit::id(), LogDebug, HcalTimeSlew::Medium, plotscripts::rms(), CaloVSimParameterMap::simParameters(), HcalTimeSlew::Slow, HcalZDCDetId::SubdetectorId, theChargeSumsForTimeBin, theParameterMap, timeBin(), and HcalSimParameters::timeSmearRMS().

{
  // HO goes slow, HF shouldn't be used at all
  //ZDC not used for the moment

  DetId detId(hit.id());
  if(detId.det()==DetId::Calo && detId.subdetId()==HcalZDCDetId::SubdetectorId) return 0.0;
  HcalDetId hcalDetId(hit.id());
  double delay = 0.;

  if(hcalDetId.subdet() == HcalBarrel || hcalDetId.subdet() == HcalEndcap || hcalDetId.subdet() == HcalOuter ) {

    HcalTimeSlew::BiasSetting biasSetting = (hcalDetId.subdet() == HcalOuter) ?
      HcalTimeSlew::Slow :
      HcalTimeSlew::Medium;

    int tbin = timeBin(hit);
    double totalCharge=0;
    if(tbin >= 0 && tbin < 10)
      {
    ChargeSumsByChannel::const_iterator totalChargeItr = theChargeSumsForTimeBin[tbin].find(detId);
    if(totalChargeItr == theChargeSumsForTimeBin[tbin].end())
      {
        throw cms::Exception("HcalHitCorrection") << "Cannot find HCAL charge sum for hit " << hit;
      }
    totalCharge = totalChargeItr->second;
    delay = HcalTimeSlew::delay(totalCharge, biasSetting);
    LogDebug("HcalHitCorrection") << "TIMESLEWcharge " << charge(hit) 
                      << "  totalcharge " << totalCharge 
                      << " olddelay "  << HcalTimeSlew::delay(charge(hit), biasSetting) 
                      << " newdelay " << delay;
      }
    // now, the smearing
    const HcalSimParameters& params=static_cast<const HcalSimParameters&>(theParameterMap->simParameters(detId));
    if (params.doTimeSmear()) {
      double rms=params.timeSmearRMS(charge(hit));

      double smearns = CLHEP::RandGaussQ::shoot(engine) * rms;

      LogDebug("HcalHitCorrection") << "TimeSmear charge " << charge(hit) << " rms " << rms << " delay " << delay << " smearns " << smearns;
      delay+=smearns;
    }
  }

  return delay;
}

void HcalHitCorrection::fillChargeSums

(

MixCollection< PCaloHit > &

hits

)

Definition at line 19 of file HcalHitCorrection.cc.

References MixCollection< T >::begin(), charge(), TauDecayModes::dec, MixCollection< T >::end(), HcalBarrel, HcalEndcap, HcalOuter, LogDebug, theChargeSumsForTimeBin, and timeBin().

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

{
  for(MixCollection<PCaloHit>::MixItr hitItr = hits.begin();
      hitItr != hits.end(); ++hitItr)
  {
    HcalDetId hcalDetId(hitItr->id());
    if(hcalDetId.subdet() == HcalBarrel || hcalDetId.subdet() == HcalEndcap || hcalDetId.subdet() == HcalOuter )
    {
      LogDebug("HcalHitCorrection") << "HcalHitCorrection::Hit 0x" << std::hex
                          << hitItr->id() << std::dec;
      int tbin = timeBin(*hitItr);
      LogDebug("HcalHitCorrection") << "HcalHitCorrection::Hit tbin" << tbin;
      if(tbin >= 0 && tbin < 10) 
      {  
        theChargeSumsForTimeBin[tbin][DetId(hitItr->id())] += charge(*hitItr);
      }
    }
  }
}

void HcalHitCorrection::fillChargeSums

(

std::vector< PCaloHit > &

hits

)

Definition at line 40 of file HcalHitCorrection.cc.

References charge(), TauDecayModes::dec, HcalBarrel, HcalEndcap, HcalOuter, LogDebug, theChargeSumsForTimeBin, and timeBin().

                                                                 {

  for(std::vector<PCaloHit>::const_iterator hitItr = hits.begin();
      hitItr != hits.end(); ++hitItr)
  {
    HcalDetId hcalDetId(hitItr->id());
    if(hcalDetId.subdet() == HcalBarrel || hcalDetId.subdet() == HcalEndcap || hcalDetId.subdet() == HcalOuter )
    {
      LogDebug("HcalHitCorrection") << "HcalHitCorrection::Hit 0x" << std::hex
                          << hitItr->id() << std::dec;
      int tbin = timeBin(*hitItr);
      LogDebug("HcalHitCorrection") << "HcalHitCorrection::Hit tbin" << tbin;
      if(tbin >= 0 && tbin < 10) 
      {  
        theChargeSumsForTimeBin[tbin][DetId(hitItr->id())] += charge(*hitItr);
      }
    }
  }
}

int HcalHitCorrection::timeBin

(

const PCaloHit &

hit

)

const

which time bin the peak of the signal will fall in

Definition at line 126 of file HcalHitCorrection.cc.

References CaloSimParameters::binOfMaximum(), PCaloHit::id(), Parameters::parameters, CaloVSimParameterMap::simParameters(), tree::t, theParameterMap, PCaloHit::time(), timeOfFlight(), and CaloSimParameters::timePhase().

Referenced by delay(), and fillChargeSums().

{
  const CaloSimParameters & parameters = theParameterMap->simParameters(DetId(hit.id()));
  double t = hit.time() - timeOfFlight(DetId(hit.id())) + parameters.timePhase();
  return static_cast<int> (t / 25) + parameters.binOfMaximum() - 1;
}

double HcalHitCorrection::timeOfFlight

(

const DetId &

id

)

const

simple average approximation

Definition at line 134 of file HcalHitCorrection.cc.

References DetId::Calo, DetId::det(), Exception, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, HcalZDCDetId::SubdetectorId, and DetId::subdetId().

Referenced by timeBin().

{
  if(detId.det()==DetId::Calo && detId.subdetId()==HcalZDCDetId::SubdetectorId)
    return 37.666; 
  switch(detId.subdetId())
    {
    case HcalBarrel:
      return 8.4;
      break;
    case HcalEndcap:
      return 13.;
      break;
    case HcalOuter:
      return 18.7;
      break;
    case HcalForward:
      return 37.;
      break;
    default:
      throw cms::Exception("HcalHitCorrection") << "Bad Hcal subdetector " << detId.subdetId();
      break;
    }
}

Member Data Documentation

ChargeSumsByChannel HcalHitCorrection::theChargeSumsForTimeBin[10]

private

Definition at line 51 of file HcalHitCorrection.h.

Referenced by clear(), delay(), and fillChargeSums().

const CaloVSimParameterMap* HcalHitCorrection::theParameterMap

private

Definition at line 49 of file HcalHitCorrection.h.

Referenced by charge(), delay(), and timeBin().