#include <HcalSimpleRecAlgo.h>

Public Member Functions
void	beginRun (edm::EventSetup const &es)

void	endRun ()

	HcalSimpleRecAlgo (bool correctForTimeslew, bool correctForContainment, float fixedPhaseNs)

void	initPulseCorr (int toadd)

HcalCalibRecHit	reconstruct (const HcalCalibDataFrame &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

HFRecHit	reconstructQIE10 (const QIE10DataFrame &digi, int first, int toadd, const HcalCoder &coder, const HcalCalibrations &calibs) const

void	setBXInfo (const BunchXParameter *info, unsigned lenInfo)

void	setHFPileupCorrection (std::shared_ptr< AbsOOTPileupCorrection > corr)

void	setHOPileupCorrection (std::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_

std::shared_ptr< AbsOOTPileupCorrection >	hbhePileupCorr_

const HcalTimeSlew *	hcalTimeSlew_delay_

std::shared_ptr< AbsOOTPileupCorrection >	hfPileupCorr_

std::shared_ptr< AbsOOTPileupCorrection >	hoPileupCorr_

unsigned	lenBunchCrossingInfo_

float	phaseNS_

int	pileupCleaningID_

int	puCorrMethod_

std::unique_ptr< HcalPulseContainmentManager >	pulseCorr_

int	runnum_

bool	setLeakCorrection_

HcalPulseShapes	theHcalPulseShapes_

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::HcalSimpleRecAlgo	(	bool	correctForTimeslew,
		bool	correctForContainment,
		float	fixedPhaseNs
	)

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

Definition at line 20 of file HcalSimpleRecAlgo.cc.

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

References hcalTimeSlew_delay_, MaximumFractionalError, and pulseCorr_.

Member Function Documentation

◆ beginRun()

void HcalSimpleRecAlgo::beginRun ( edm::EventSetup const & es )

Definition at line 31 of file HcalSimpleRecAlgo.cc.

                                                         {
   edm::ESHandle<HcalTimeSlew> delay;
   es.get<HcalTimeSlewRecord>().get("HBHE", delay);
   hcalTimeSlew_delay_ = &*delay;
  
   pulseCorr_->beginRun(es);
 }

References phase2TrackerDigitizer_cfi::delay, edm::EventSetup::get(), get, hcalTimeSlew_delay_, and pulseCorr_.

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

◆ endRun()

void HcalSimpleRecAlgo::endRun ( )

Definition at line 39 of file HcalSimpleRecAlgo.cc.

39 {}

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

◆ initPulseCorr()

void HcalSimpleRecAlgo::initPulseCorr ( int toadd )

Definition at line 41 of file HcalSimpleRecAlgo.cc.

41 {}

◆ reconstruct() [1/3]

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

Definition at line 355 of file HcalSimpleRecAlgo.cc.

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

References bunchCrossingInfo_, correctForPulse_, correctForTimeslew_, HcalTimeSlew::Fast, dqmdumpme::first, hcalTimeSlew_delay_, HcalCalibDataFrame::id(), lenBunchCrossingInfo_, phaseNS_, puCorrMethod_, pulseCorr_, and runnum_.

◆ reconstruct() [2/3]

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

Definition at line 378 of file HcalSimpleRecAlgo.cc.

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

References bunchCrossingInfo_, alignCSCRings::corr, correctForPulse_, correctForTimeslew_, dqmdumpme::first, hfPileupCorr_, HFDataFrame::id(), lenBunchCrossingInfo_, findQualityFiles::maxI, phaseNS_, pulseCorr_, HcalSimpleRecAlgoImpl::recoHFTime(), HcalSimpleRecAlgoImpl::removePileup(), setRawEnergy(), FrontierCondition_GT_autoExpress_cfi::t0, RandomServiceHelper::t2, ntuplemaker::time, and HcalCalibrations::timecorr().

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

◆ reconstruct() [3/3]

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

Definition at line 335 of file HcalSimpleRecAlgo.cc.

                                                                                                                  {
   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_,
                                                             puCorrMethod_,
                                                             hcalTimeSlew_delay_);
 }

References bunchCrossingInfo_, correctForPulse_, correctForTimeslew_, dqmdumpme::first, hcalTimeSlew_delay_, hoPileupCorr_, HODataFrame::id(), lenBunchCrossingInfo_, phaseNS_, puCorrMethod_, pulseCorr_, runnum_, and HcalTimeSlew::Slow.

◆ reconstructQIE10()

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

Definition at line 419 of file HcalSimpleRecAlgo.cc.

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

References bunchCrossingInfo_, alignCSCRings::corr, correctForPulse_, correctForTimeslew_, dqmdumpme::first, hfPileupCorr_, QIE10DataFrame::id(), lenBunchCrossingInfo_, findQualityFiles::maxI, phaseNS_, pulseCorr_, HcalSimpleRecAlgoImpl::recoHFTime(), HcalSimpleRecAlgoImpl::removePileup(), setRawEnergy(), FrontierCondition_GT_autoExpress_cfi::t0, RandomServiceHelper::t2, ntuplemaker::time, and HcalCalibrations::timecorr().

◆ setBXInfo()

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

Definition at line 58 of file HcalSimpleRecAlgo.cc.

                                                                                      {
   bunchCrossingInfo_ = info;
   lenBunchCrossingInfo_ = lenInfo;
 }

References bunchCrossingInfo_, info(), and lenBunchCrossingInfo_.

◆ setHFPileupCorrection()

void HcalSimpleRecAlgo::setHFPileupCorrection ( std::shared_ptr< AbsOOTPileupCorrection > corr )

Definition at line 54 of file HcalSimpleRecAlgo.cc.

54 { hfPileupCorr_ = corr; }

References alignCSCRings::corr, and hfPileupCorr_.

◆ setHOPileupCorrection()

void HcalSimpleRecAlgo::setHOPileupCorrection ( std::shared_ptr< AbsOOTPileupCorrection > corr )

Definition at line 56 of file HcalSimpleRecAlgo.cc.

56 { hoPileupCorr_ = corr; }

References alignCSCRings::corr, and hoPileupCorr_.

◆ setLeakCorrection()

void HcalSimpleRecAlgo::setLeakCorrection ( )

Definition at line 52 of file HcalSimpleRecAlgo.cc.

52 { setLeakCorrection_ = true; }

References setLeakCorrection_.

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

◆ setRecoParams()

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

Definition at line 43 of file HcalSimpleRecAlgo.cc.

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

References correctForPulse_, castor_dqm_sourceclient-live_cfg::correctForTimeslew, correctForTimeslew_, phaseNS_, pileupCleaningID_, setLeakCorrection(), and setLeakCorrection_.

Referenced by HcalHitReconstructor::produce().

Member Data Documentation

◆ bunchCrossingInfo_

const BunchXParameter* HcalSimpleRecAlgo::bunchCrossingInfo_

private

Definition at line 80 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), reconstructQIE10(), and setBXInfo().

◆ correctForPulse_

bool HcalSimpleRecAlgo::correctForPulse_

private

Definition at line 74 of file HcalSimpleRecAlgo.h.

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

◆ correctForTimeslew_

bool HcalSimpleRecAlgo::correctForTimeslew_

private

Definition at line 73 of file HcalSimpleRecAlgo.h.

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

◆ hbhePileupCorr_

std::shared_ptr<AbsOOTPileupCorrection> HcalSimpleRecAlgo::hbhePileupCorr_

private

Definition at line 82 of file HcalSimpleRecAlgo.h.

◆ hcalTimeSlew_delay_

const HcalTimeSlew* HcalSimpleRecAlgo::hcalTimeSlew_delay_

private

Definition at line 90 of file HcalSimpleRecAlgo.h.

Referenced by beginRun(), HcalSimpleRecAlgo(), and reconstruct().

◆ hfPileupCorr_

std::shared_ptr<AbsOOTPileupCorrection> HcalSimpleRecAlgo::hfPileupCorr_

private

Definition at line 83 of file HcalSimpleRecAlgo.h.

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

◆ hoPileupCorr_

std::shared_ptr<AbsOOTPileupCorrection> HcalSimpleRecAlgo::hoPileupCorr_

private

Definition at line 84 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), and setHOPileupCorrection().

◆ lenBunchCrossingInfo_

unsigned HcalSimpleRecAlgo::lenBunchCrossingInfo_

private

Definition at line 81 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), reconstructQIE10(), and setBXInfo().

◆ phaseNS_

float HcalSimpleRecAlgo::phaseNS_

private

Definition at line 75 of file HcalSimpleRecAlgo.h.

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

◆ pileupCleaningID_

int HcalSimpleRecAlgo::pileupCleaningID_

private

Definition at line 79 of file HcalSimpleRecAlgo.h.

Referenced by setRecoParams().

◆ puCorrMethod_

int HcalSimpleRecAlgo::puCorrMethod_

private

Definition at line 88 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct().

◆ pulseCorr_

std::unique_ptr<HcalPulseContainmentManager> HcalSimpleRecAlgo::pulseCorr_

private

Definition at line 76 of file HcalSimpleRecAlgo.h.

Referenced by beginRun(), HcalSimpleRecAlgo(), reconstruct(), and reconstructQIE10().

◆ runnum_

int HcalSimpleRecAlgo::runnum_

private

Definition at line 77 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct().

◆ setLeakCorrection_

bool HcalSimpleRecAlgo::setLeakCorrection_

private

Definition at line 78 of file HcalSimpleRecAlgo.h.

Referenced by setLeakCorrection(), and setRecoParams().

◆ theHcalPulseShapes_

HcalPulseShapes HcalSimpleRecAlgo::theHcalPulseShapes_

private

Definition at line 86 of file HcalSimpleRecAlgo.h.

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ HcalSimpleRecAlgo()

Member Function Documentation

◆ beginRun()

◆ endRun()

◆ initPulseCorr()

◆ reconstruct() [1/3]

◆ reconstruct() [2/3]

◆ reconstruct() [3/3]

◆ reconstructQIE10()

◆ setBXInfo()

◆ setHFPileupCorrection()

◆ setHOPileupCorrection()

◆ setLeakCorrection()

◆ setRecoParams()

Member Data Documentation

◆ bunchCrossingInfo_

◆ correctForPulse_

◆ correctForTimeslew_

◆ hbhePileupCorr_

◆ hcalTimeSlew_delay_

◆ hfPileupCorr_

◆ hoPileupCorr_

◆ lenBunchCrossingInfo_

◆ phaseNS_

◆ pileupCleaningID_

◆ puCorrMethod_

◆ pulseCorr_

◆ runnum_

◆ setLeakCorrection_

◆ theHcalPulseShapes_