HcalSimpleRecAlgo Class Reference

#include <HcalSimpleRecAlgo.h>

Public Member Functions
void	beginRun (edm::EventSetup const &es)
void	endRun ()
	HcalSimpleRecAlgo (bool correctForTimeslew, bool correctForContainment, float fixedPhaseNs)
	HcalSimpleRecAlgo ()
void	initPulseCorr (int toadd)
HBHERecHit	reconstruct (const HBHEDataFrame &digi, int first, int toadd, const HcalCoder &coder, const HcalCalibrations &calibs) const
HFRecHit	reconstruct (const HFDataFrame &digi, int first, int toadd, const HcalCoder &coder, const HcalCalibrations &calibs) const
HORecHit	reconstruct (const HODataFrame &digi, int first, int toadd, const HcalCoder &coder, const HcalCalibrations &calibs) const
HcalCalibRecHit	reconstruct (const HcalCalibDataFrame &digi, int first, int toadd, const HcalCoder &coder, const HcalCalibrations &calibs) const
HBHERecHit	reconstructHBHEUpgrade (const HcalUpgradeDataFrame &digi, int first, int toadd, const HcalCoder &coder, const HcalCalibrations &calibs) const
HFRecHit	reconstructHFUpgrade (const HcalUpgradeDataFrame &digi, int first, int toadd, const HcalCoder &coder, const HcalCalibrations &calibs) const
void	setBXInfo (const BunchXParameter *info, unsigned lenInfo)
void	setForData (int runnum)
void	setHBHEPileupCorrection (boost::shared_ptr< AbsOOTPileupCorrection > corr)
void	setHFPileupCorrection (boost::shared_ptr< AbsOOTPileupCorrection > corr)
void	setHOPileupCorrection (boost::shared_ptr< AbsOOTPileupCorrection > corr)
void	setLeakCorrection ()
void	setRecoParams (bool correctForTimeslew, bool correctForPulse, bool setLeakCorrection, int pileupCleaningID, float phaseNS)

Private Attributes
const BunchXParameter *	bunchCrossingInfo_
bool	correctForPulse_
bool	correctForTimeslew_
boost::shared_ptr < AbsOOTPileupCorrection >	hbhePileupCorr_
boost::shared_ptr < AbsOOTPileupCorrection >	hfPileupCorr_
boost::shared_ptr < AbsOOTPileupCorrection >	hoPileupCorr_
unsigned	lenBunchCrossingInfo_
float	phaseNS_
int	pileupCleaningID_
std::auto_ptr < HcalPulseContainmentManager >	pulseCorr_
int	runnum_
bool	setLeakCorrection_

Detailed Description

This class reconstructs RecHits from Digis for HBHE, HF, and HO by addition of selected time samples, pedestal subtraction, and gain application. The time of the hit is reconstructed using a weighted peak bin calculation supplemented by precise time lookup table. A consumer of this class also has the option of correcting the reconstructed time for energy-dependent time slew associated with the QIE.

Author

J. Mans - Minnesota

Definition at line 34 of file HcalSimpleRecAlgo.h.

Constructor & Destructor Documentation

HcalSimpleRecAlgo::HcalSimpleRecAlgo	(	bool	correctForTimeslew,
		bool	correctForContainment,
		float	fixedPhaseNs
	)

Full featured constructor for HB/HE and HO (HPD-based detectors)

Definition at line 15 of file HcalSimpleRecAlgo.cc.

References MaximumFractionalError, and pulseCorr_.

                                                                                                 : 
  correctForTimeslew_(correctForTimeslew),
  correctForPulse_(correctForPulse),
  phaseNS_(phaseNS), runnum_(0), setLeakCorrection_(false)
{ 
  
  pulseCorr_ = std::auto_ptr<HcalPulseContainmentManager>(
           new HcalPulseContainmentManager(MaximumFractionalError)
                              );
}

HcalSimpleRecAlgo::HcalSimpleRecAlgo

(

)

Simple constructor for PMT-based detectors

Definition at line 27 of file HcalSimpleRecAlgo.cc.

                                     : 
  correctForTimeslew_(false), runnum_(0) { }

Member Function Documentation

void HcalSimpleRecAlgo::beginRun

(

edm::EventSetup const &

es

)

Definition at line 31 of file HcalSimpleRecAlgo.cc.

References pulseCorr_.

Referenced by HcalSimpleReconstructor::beginRun(), and HcalHitReconstructor::beginRun().

{
  pulseCorr_->beginRun(es);
}

void HcalSimpleRecAlgo::endRun

(

void

)

Definition at line 37 of file HcalSimpleRecAlgo.cc.

References pulseCorr_.

Referenced by HcalSimpleReconstructor::endRun(), and HcalHitReconstructor::endRun().

{
  pulseCorr_->endRun();
}

void HcalSimpleRecAlgo::initPulseCorr

(

int

toadd

)

Definition at line 43 of file HcalSimpleRecAlgo.cc.

                                               {
}

HBHERecHit HcalSimpleRecAlgo::reconstruct	(	const HBHEDataFrame &	digi,
		int	first,
		int	toadd,
		const HcalCoder &	coder,
		const HcalCalibrations &	calibs
	)		const

Definition at line 348 of file HcalSimpleRecAlgo.cc.

References bunchCrossingInfo_, correctForPulse_, correctForTimeslew_, first, hbhePileupCorr_, lenBunchCrossingInfo_, HcalTimeSlew::Medium, phaseNS_, pulseCorr_, runnum_, and setLeakCorrection_.

Referenced by HcalSimpleReconstructor::process(), and HcalHitReconstructor::produce().

                                                                                                                                                       {
  return HcalSimpleRecAlgoImpl::reco<HBHEDataFrame,HBHERecHit>(digi,coder,calibs,
                                   first,toadd,correctForTimeslew_, correctForPulse_,
                                   pulseCorr_->get(digi.id(), toadd, phaseNS_),
                                   HcalTimeSlew::Medium,
                                                               runnum_, setLeakCorrection_,
                                                               hbhePileupCorr_.get(),
                                                               bunchCrossingInfo_, lenBunchCrossingInfo_);
}

HFRecHit HcalSimpleRecAlgo::reconstruct	(	const HFDataFrame &	digi,
		int	first,
		int	toadd,
		const HcalCoder &	coder,
		const HcalCalibrations &	calibs
	)		const

Definition at line 392 of file HcalSimpleRecAlgo.cc.

References bunchCrossingInfo_, corr, correctForPulse_, correctForTimeslew_, hfPileupCorr_, HFDataFrame::id(), lenBunchCrossingInfo_, findQualityFiles::maxI, phaseNS_, pulseCorr_, HcalSimpleRecAlgoImpl::recoHFTime(), HcalSimpleRecAlgoImpl::removePileup(), setRawEnergy(), reco::t2, cond::rpcobgas::time, and HcalCalibrations::timecorr().

{
  const HcalPulseContainmentCorrection* corr = pulseCorr_->get(digi.id(), toadd, phaseNS_);

  double amp_fC, ampl, uncorr_ampl, maxA;
  int nRead, maxI;
  bool leakCorrApplied;
  float t0, t2;

  HcalSimpleRecAlgoImpl::removePileup(digi, coder, calibs, first, toadd,
                                      correctForPulse_, corr, hfPileupCorr_.get(),
                                      bunchCrossingInfo_, lenBunchCrossingInfo_,
                                      &maxA, &ampl, &uncorr_ampl, &amp_fC, &nRead,
                                      &maxI, &leakCorrApplied, &t0, &t2);

  float time=-9999.f;
  if (maxI > 0 && maxI < (nRead - 1))
      time = HcalSimpleRecAlgoImpl::recoHFTime(digi,maxI,amp_fC,correctForTimeslew_,maxA,t0,t2) -
             calibs.timecorr();

  HFRecHit rh(digi.id(),ampl,time);
  setRawEnergy(rh, static_cast<float>(uncorr_ampl));
  return rh;
}

HORecHit HcalSimpleRecAlgo::reconstruct	(	const HODataFrame &	digi,
		int	first,
		int	toadd,
		const HcalCoder &	coder,
		const HcalCalibrations &	calibs
	)		const

Definition at line 359 of file HcalSimpleRecAlgo.cc.

References bunchCrossingInfo_, correctForPulse_, correctForTimeslew_, first, hoPileupCorr_, lenBunchCrossingInfo_, phaseNS_, pulseCorr_, runnum_, and HcalTimeSlew::Slow.

                                                                                                                                                   {
  return HcalSimpleRecAlgoImpl::reco<HODataFrame,HORecHit>(digi,coder,calibs,
                               first,toadd,correctForTimeslew_,correctForPulse_,
                               pulseCorr_->get(digi.id(), toadd, phaseNS_),
                               HcalTimeSlew::Slow,
                                                           runnum_, false, hoPileupCorr_.get(),
                                                           bunchCrossingInfo_, lenBunchCrossingInfo_);
}

HcalCalibRecHit HcalSimpleRecAlgo::reconstruct	(	const HcalCalibDataFrame &	digi,
		int	first,
		int	toadd,
		const HcalCoder &	coder,
		const HcalCalibrations &	calibs
	)		const

Definition at line 369 of file HcalSimpleRecAlgo.cc.

References bunchCrossingInfo_, correctForPulse_, correctForTimeslew_, HcalTimeSlew::Fast, first, lenBunchCrossingInfo_, phaseNS_, pulseCorr_, and runnum_.

                                                                                                                                                                 {
  return HcalSimpleRecAlgoImpl::reco<HcalCalibDataFrame,HcalCalibRecHit>(digi,coder,calibs,
                                     first,toadd,correctForTimeslew_,correctForPulse_,
                                     pulseCorr_->get(digi.id(), toadd, phaseNS_),
                                     HcalTimeSlew::Fast,
                                                                         runnum_, false, 0,
                                                                         bunchCrossingInfo_, lenBunchCrossingInfo_);
}

HBHERecHit HcalSimpleRecAlgo::reconstructHBHEUpgrade	(	const HcalUpgradeDataFrame &	digi,
		int	first,
		int	toadd,
		const HcalCoder &	coder,
		const HcalCalibrations &	calibs
	)		const

Definition at line 379 of file HcalSimpleRecAlgo.cc.

References bunchCrossingInfo_, correctForPulse_, correctForTimeslew_, first, hbhePileupCorr_, lenBunchCrossingInfo_, HcalTimeSlew::Medium, phaseNS_, pulseCorr_, query::result, and EcalTBTDCReconstructor_cfi::tdcReco.

Referenced by HcalSimpleReconstructor::processUpgrade().

                                                                                                                                                                         {
  HBHERecHit result = HcalSimpleRecAlgoImpl::reco<HcalUpgradeDataFrame,HBHERecHit>(digi, coder, calibs,
                                                                                   first, toadd, correctForTimeslew_, correctForPulse_,
                                                                                   pulseCorr_->get(digi.id(), toadd, phaseNS_),
                                                                                   HcalTimeSlew::Medium, 0, false,
                                                                                   hbhePileupCorr_.get(),
                                                                                   bunchCrossingInfo_, lenBunchCrossingInfo_);
  HcalTDCReco tdcReco;
  tdcReco.reconstruct(digi, result);
  return result;
}

HFRecHit HcalSimpleRecAlgo::reconstructHFUpgrade	(	const HcalUpgradeDataFrame &	digi,
		int	first,
		int	toadd,
		const HcalCoder &	coder,
		const HcalCalibrations &	calibs
	)		const

Definition at line 424 of file HcalSimpleRecAlgo.cc.

References bunchCrossingInfo_, corr, correctForPulse_, correctForTimeslew_, hfPileupCorr_, HcalUpgradeDataFrame::id(), lenBunchCrossingInfo_, findQualityFiles::maxI, phaseNS_, pulseCorr_, HcalSimpleRecAlgoImpl::recoHFTime(), HcalSimpleRecAlgoImpl::removePileup(), setRawEnergy(), reco::t2, cond::rpcobgas::time, and HcalCalibrations::timecorr().

Referenced by HcalSimpleReconstructor::processUpgrade().

{
  const HcalPulseContainmentCorrection* corr = pulseCorr_->get(digi.id(), toadd, phaseNS_);

  double amp_fC, ampl, uncorr_ampl, maxA;
  int nRead, maxI;
  bool leakCorrApplied;
  float t0, t2;

  HcalSimpleRecAlgoImpl::removePileup(digi, coder, calibs, first, toadd,
                                      correctForPulse_, corr, hfPileupCorr_.get(),
                                      bunchCrossingInfo_, lenBunchCrossingInfo_,
                                      &maxA, &ampl, &uncorr_ampl, &amp_fC, &nRead,
                                      &maxI, &leakCorrApplied, &t0, &t2);

  float time=-9999.f;
  if (maxI > 0 && maxI < (nRead - 1))
      time = HcalSimpleRecAlgoImpl::recoHFTime(digi,maxI,amp_fC,correctForTimeslew_,maxA,t0,t2) -
             calibs.timecorr();

  HFRecHit rh(digi.id(),ampl,time); // new RecHit gets second time = 0.
  setRawEnergy(rh, static_cast<float>(uncorr_ampl));
  return rh;
}

void HcalSimpleRecAlgo::setBXInfo	(	const BunchXParameter *	info,
		unsigned	lenInfo
	)

Definition at line 76 of file HcalSimpleRecAlgo.cc.

References bunchCrossingInfo_, info(), and lenBunchCrossingInfo_.

{
    bunchCrossingInfo_ = info;
    lenBunchCrossingInfo_ = lenInfo;
}

void HcalSimpleRecAlgo::setForData

(

int

runnum

)

Definition at line 54 of file HcalSimpleRecAlgo.cc.

References summaryLumi::runnum, and runnum_.

Referenced by HcalSimpleReconstructor::produce(), and HcalHitReconstructor::produce().

{ runnum_ = runnum;}

void HcalSimpleRecAlgo::setHBHEPileupCorrection

(

boost::shared_ptr< AbsOOTPileupCorrection >

corr

)

Definition at line 58 of file HcalSimpleRecAlgo.cc.

References corr, and hbhePileupCorr_.

Referenced by HcalHitReconstructor::HcalHitReconstructor().

{
    hbhePileupCorr_ = corr;
}

void HcalSimpleRecAlgo::setHFPileupCorrection

(

boost::shared_ptr< AbsOOTPileupCorrection >

corr

)

Definition at line 64 of file HcalSimpleRecAlgo.cc.

References corr, and hfPileupCorr_.

Referenced by HcalHitReconstructor::HcalHitReconstructor().

{
    hfPileupCorr_ = corr;
}

void HcalSimpleRecAlgo::setHOPileupCorrection

(

boost::shared_ptr< AbsOOTPileupCorrection >

corr

)

Definition at line 70 of file HcalSimpleRecAlgo.cc.

References corr, and hoPileupCorr_.

Referenced by HcalHitReconstructor::HcalHitReconstructor().

{
    hoPileupCorr_ = corr;
}

void HcalSimpleRecAlgo::setLeakCorrection

(

)

Definition at line 56 of file HcalSimpleRecAlgo.cc.

References setLeakCorrection_.

Referenced by HcalHitReconstructor::produce(), and setRecoParams().

{ setLeakCorrection_ = true;}

void HcalSimpleRecAlgo::setRecoParams	(	bool	correctForTimeslew,
		bool	correctForPulse,
		bool	setLeakCorrection,
		int	pileupCleaningID,
		float	phaseNS
	)

Definition at line 46 of file HcalSimpleRecAlgo.cc.

References correctForPulse_, correctForTimeslew_, phaseNS_, pileupCleaningID_, setLeakCorrection(), and setLeakCorrection_.

Referenced by HcalHitReconstructor::produce().

                                                                                                                                               {
   correctForTimeslew_=correctForTimeslew;
   correctForPulse_=correctForPulse;
   phaseNS_=phaseNS;
   setLeakCorrection_=setLeakCorrection;
   pileupCleaningID_=pileupCleaningID;
}

Member Data Documentation

const BunchXParameter* HcalSimpleRecAlgo::bunchCrossingInfo_

private

Definition at line 85 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), reconstructHBHEUpgrade(), reconstructHFUpgrade(), and setBXInfo().

bool HcalSimpleRecAlgo::correctForPulse_

private

Definition at line 79 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), reconstructHBHEUpgrade(), reconstructHFUpgrade(), and setRecoParams().

bool HcalSimpleRecAlgo::correctForTimeslew_

private

Definition at line 78 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), reconstructHBHEUpgrade(), reconstructHFUpgrade(), and setRecoParams().

boost::shared_ptr<AbsOOTPileupCorrection> HcalSimpleRecAlgo::hbhePileupCorr_

private

Definition at line 87 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), reconstructHBHEUpgrade(), and setHBHEPileupCorrection().

boost::shared_ptr<AbsOOTPileupCorrection> HcalSimpleRecAlgo::hfPileupCorr_

private

Definition at line 88 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), reconstructHFUpgrade(), and setHFPileupCorrection().

boost::shared_ptr<AbsOOTPileupCorrection> HcalSimpleRecAlgo::hoPileupCorr_

private

Definition at line 89 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), and setHOPileupCorrection().

unsigned HcalSimpleRecAlgo::lenBunchCrossingInfo_

private

Definition at line 86 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), reconstructHBHEUpgrade(), reconstructHFUpgrade(), and setBXInfo().

float HcalSimpleRecAlgo::phaseNS_

private

Definition at line 80 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), reconstructHBHEUpgrade(), reconstructHFUpgrade(), and setRecoParams().

int HcalSimpleRecAlgo::pileupCleaningID_

private

Definition at line 84 of file HcalSimpleRecAlgo.h.

Referenced by setRecoParams().

std::auto_ptr<HcalPulseContainmentManager> HcalSimpleRecAlgo::pulseCorr_

private

Definition at line 81 of file HcalSimpleRecAlgo.h.

Referenced by beginRun(), endRun(), HcalSimpleRecAlgo(), reconstruct(), reconstructHBHEUpgrade(), and reconstructHFUpgrade().

int HcalSimpleRecAlgo::runnum_

private

Definition at line 82 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), and setForData().

bool HcalSimpleRecAlgo::setLeakCorrection_

private

Definition at line 83 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), setLeakCorrection(), and setRecoParams().