HcalSimpleRecAlgo Class Reference

#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)
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
HFRecHit	reconstructQIE10 (const QIE10DataFrame &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	setMeth3Params (bool iApplyTimeSlew, float iPedSubThreshold, int iTimeSlewParsType, std::vector< double > iTimeSlewPars, double irespCorrM3)
void	setpuCorrMethod (int method)
void	setpuCorrParams (bool iPedestalConstraint, bool iTimeConstraint, bool iAddPulseJitter, bool iApplyTimeSlew, double iTS4Min, const std::vector< double > &iTS4Max, double iPulseJitter, double iTimeMean, double iTimeSig, double iTimeSigSiPM, double iPedMean, double iPedSig, double iPedSigSiPM, double iNoise, double iNoiseSiPM, double iTMin, double iTMax, const std::vector< double > &its4Chi2, int iFitTimes)
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_
std::unique_ptr< HcalDeterministicFit >	hltOOTpuCorr_
boost::shared_ptr< AbsOOTPileupCorrection >	hoPileupCorr_
unsigned	lenBunchCrossingInfo_
std::unique_ptr< PedestalSub >	pedSubFxn_
float	phaseNS_
int	pileupCleaningID_
std::unique_ptr< PulseShapeFitOOTPileupCorrection >	psFitOOTpuCorr_
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 40 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 19 of file HcalSimpleRecAlgo.cc.

References hltOOTpuCorr_, MaximumFractionalError, pedSubFxn_, and pulseCorr_.

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

Member Function Documentation

void HcalSimpleRecAlgo::beginRun

(

edm::EventSetup const &

es

)

Definition at line 30 of file HcalSimpleRecAlgo.cc.

References pulseCorr_.

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

{
  pulseCorr_->beginRun(es);
}

void HcalSimpleRecAlgo::endRun

(

)

Definition at line 36 of file HcalSimpleRecAlgo.cc.

References pulseCorr_.

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

{
  pulseCorr_->endRun();
}

void HcalSimpleRecAlgo::initPulseCorr

(

int

toadd

)

Definition at line 42 of file HcalSimpleRecAlgo.cc.

                                               {
}

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

Definition at line 520 of file HcalSimpleRecAlgo.cc.

References bunchCrossingInfo_, correctForPulse_, correctForTimeslew_, plotBeamSpotDB::first, hbhePileupCorr_, hltOOTpuCorr_, lenBunchCrossingInfo_, HcalTimeSlew::Medium, pedSubFxn_, phaseNS_, psFitOOTpuCorr_, puCorrMethod_, pulseCorr_, runnum_, and setLeakCorrection_.

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

                                                                                                                                                       {
  return HcalSimpleRecAlgoImpl::recoHBHE<HBHEDataFrame,HBHERecHit>(digi,coder,calibs,
                                   first,toadd,correctForTimeslew_, correctForPulse_,
                                   pulseCorr_->get(digi.id(), toadd, phaseNS_),
                                   HcalTimeSlew::Medium,
                                                               runnum_, setLeakCorrection_,
                                                               hbhePileupCorr_.get(),
                                                               bunchCrossingInfo_, lenBunchCrossingInfo_, puCorrMethod_, psFitOOTpuCorr_.get(),/*hlt*/hltOOTpuCorr_.get(),pedSubFxn_.get());
}

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

Definition at line 551 of file HcalSimpleRecAlgo.cc.

References bunchCrossingInfo_, corr, correctForPulse_, correctForTimeslew_, hfPileupCorr_, HFDataFrame::id(), lenBunchCrossingInfo_, findQualityFiles::maxI, phaseNS_, pulseCorr_, HcalSimpleRecAlgoImpl::recoHFTime(), HcalSimpleRecAlgoImpl::removePileup(), setRawEnergy(), cscNeutronWriter_cfi::t0, reco::t2, ntuplemaker::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 531 of file HcalSimpleRecAlgo.cc.

References bunchCrossingInfo_, correctForPulse_, correctForTimeslew_, plotBeamSpotDB::first, hltOOTpuCorr_, hoPileupCorr_, lenBunchCrossingInfo_, pedSubFxn_, phaseNS_, psFitOOTpuCorr_, puCorrMethod_, 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_, puCorrMethod_, psFitOOTpuCorr_.get(),/*hlt*/hltOOTpuCorr_.get(),pedSubFxn_.get());
}

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

Definition at line 541 of file HcalSimpleRecAlgo.cc.

References bunchCrossingInfo_, correctForPulse_, correctForTimeslew_, HcalTimeSlew::Fast, plotBeamSpotDB::first, hltOOTpuCorr_, lenBunchCrossingInfo_, pedSubFxn_, phaseNS_, psFitOOTpuCorr_, puCorrMethod_, 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_, puCorrMethod_, psFitOOTpuCorr_.get(),/*hlt*/hltOOTpuCorr_.get(),pedSubFxn_.get());
}

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

Definition at line 580 of file HcalSimpleRecAlgo.cc.

References bunchCrossingInfo_, corr, correctForPulse_, correctForTimeslew_, hfPileupCorr_, QIE10DataFrame::id(), lenBunchCrossingInfo_, findQualityFiles::maxI, phaseNS_, pulseCorr_, HcalSimpleRecAlgoImpl::recoHFTime(), HcalSimpleRecAlgoImpl::removePileup(), setRawEnergy(), cscNeutronWriter_cfi::t0, reco::t2, ntuplemaker::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;
}

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

Definition at line 112 of file HcalSimpleRecAlgo.cc.

References bunchCrossingInfo_, info(), leakCorr(), lenBunchCrossingInfo_, and timeshift_ns_hf().

{
    bunchCrossingInfo_ = info;
    lenBunchCrossingInfo_ = lenInfo;
}

void HcalSimpleRecAlgo::setForData

(

int

runnum

)

Definition at line 80 of file HcalSimpleRecAlgo.cc.

References HcalPulseShapes::getShape(), HPDShapev3DataNum, HPDShapev3MCNum, psFitOOTpuCorr_, puCorrMethod_, summaryLumi::runnum, runnum_, and theHcalPulseShapes_.

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

                                              { 
   runnum_ = runnum;
   if( puCorrMethod_ ==2 ){
      int shapeNum = HPDShapev3MCNum;
      if( runnum_ > 0 ){
         shapeNum = HPDShapev3DataNum;
      }
      bool isHPD=true;
      psFitOOTpuCorr_->setPulseShapeTemplate(theHcalPulseShapes_.getShape(shapeNum),isHPD);
   }
}

void HcalSimpleRecAlgo::setHBHEPileupCorrection

(

boost::shared_ptr< AbsOOTPileupCorrection >

corr

)

Definition at line 94 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 100 of file HcalSimpleRecAlgo.cc.

References corr, and hfPileupCorr_.

{
    hfPileupCorr_ = corr;
}

void HcalSimpleRecAlgo::setHOPileupCorrection

(

boost::shared_ptr< AbsOOTPileupCorrection >

corr

)

Definition at line 106 of file HcalSimpleRecAlgo.cc.

References corr, and hoPileupCorr_.

{
    hoPileupCorr_ = corr;
}

void HcalSimpleRecAlgo::setLeakCorrection

(

)

Definition at line 92 of file HcalSimpleRecAlgo.cc.

References setLeakCorrection_.

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

{ setLeakCorrection_ = true;}

void HcalSimpleRecAlgo::setMeth3Params	(	bool	iApplyTimeSlew,
		float	iPedSubThreshold,
		int	iTimeSlewParsType,
		std::vector< double >	iTimeSlewPars,
		double	irespCorrM3
	)

Definition at line 73 of file HcalSimpleRecAlgo.cc.

References hltOOTpuCorr_, HcalTimeSlew::Medium, and pedSubFxn_.

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

                                                                                                                                                               {

  pedSubFxn_->init(0, iPedSubThreshold, 0.0);
  hltOOTpuCorr_->init((HcalTimeSlew::ParaSource)iTimeSlewParsType, HcalTimeSlew::Medium, iApplyTimeSlew, *pedSubFxn_, iTimeSlewPars,irespCorrM3);

}

void HcalSimpleRecAlgo::setpuCorrMethod ( int  method )

inline

Definition at line 78 of file HcalSimpleRecAlgo.h.

References AlcaSiPixelAliHarvester0T_cff::method, psFitOOTpuCorr_, puCorrMethod_, setMeth3Params(), and setpuCorrParams().

Referenced by HcalHitReconstructor::HcalHitReconstructor().

                                  { 
    puCorrMethod_ = method;
    if( puCorrMethod_ == 2 )
        psFitOOTpuCorr_ = std::make_unique<PulseShapeFitOOTPileupCorrection>();
  }

void HcalSimpleRecAlgo::setpuCorrParams	(	bool	iPedestalConstraint,
		bool	iTimeConstraint,
		bool	iAddPulseJitter,
		bool	iApplyTimeSlew,
		double	iTS4Min,
		const std::vector< double > &	iTS4Max,
		double	iPulseJitter,
		double	iTimeMean,
		double	iTimeSig,
		double	iTimeSigSiPM,
		double	iPedMean,
		double	iPedSig,
		double	iPedSigSiPM,
		double	iNoise,
		double	iNoiseSiPM,
		double	iTMin,
		double	iTMax,
		const std::vector< double > &	its4Chi2,
		int	iFitTimes
	)

Definition at line 53 of file HcalSimpleRecAlgo.cc.

References HcalTimeSlew::Medium, and psFitOOTpuCorr_.

Referenced by HcalHitReconstructor::HcalHitReconstructor(), and setpuCorrMethod().

                                                                       {
  if( iPedestalConstraint ) assert ( iPedSig );
  if( iTimeConstraint ) assert( iTimeSig );
  psFitOOTpuCorr_->setPUParams(iPedestalConstraint,iTimeConstraint,iAddPulseJitter,iApplyTimeSlew,
                   iTS4Min, iTS4Max, iPulseJitter,iTimeMean,iTimeSig,iTimeSigSiPM,iPedMean,iPedSig,iPedSigSiPM,iNoise,iNoiseSiPM,iTMin,iTMax,its4Chi2,
                   HcalTimeSlew::Medium, iFitTimes);

  psFitOOTpuCorr_->setChi2Term(1); // isHPD all the time

//  int shapeNum = HPDShapev3MCNum;
//  psFitOOTpuCorr_->setPulseShapeTemplate(theHcalPulseShapes_.getShape(shapeNum));
}

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

Definition at line 45 of file HcalSimpleRecAlgo.cc.

References correctForPulse_, CastorSimpleReconstructor_cfi::correctForTimeslew, 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 101 of file HcalSimpleRecAlgo.h.

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

bool HcalSimpleRecAlgo::correctForPulse_

private

Definition at line 95 of file HcalSimpleRecAlgo.h.

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

bool HcalSimpleRecAlgo::correctForTimeslew_

private

Definition at line 94 of file HcalSimpleRecAlgo.h.

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

boost::shared_ptr<AbsOOTPileupCorrection> HcalSimpleRecAlgo::hbhePileupCorr_

private

Definition at line 103 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), and setHBHEPileupCorrection().

boost::shared_ptr<AbsOOTPileupCorrection> HcalSimpleRecAlgo::hfPileupCorr_

private

Definition at line 104 of file HcalSimpleRecAlgo.h.

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

std::unique_ptr<HcalDeterministicFit> HcalSimpleRecAlgo::hltOOTpuCorr_

private

Definition at line 116 of file HcalSimpleRecAlgo.h.

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

boost::shared_ptr<AbsOOTPileupCorrection> HcalSimpleRecAlgo::hoPileupCorr_

private

Definition at line 105 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), and setHOPileupCorrection().

unsigned HcalSimpleRecAlgo::lenBunchCrossingInfo_

private

Definition at line 102 of file HcalSimpleRecAlgo.h.

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

std::unique_ptr<PedestalSub> HcalSimpleRecAlgo::pedSubFxn_

private

Definition at line 113 of file HcalSimpleRecAlgo.h.

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

float HcalSimpleRecAlgo::phaseNS_

private

Definition at line 96 of file HcalSimpleRecAlgo.h.

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

int HcalSimpleRecAlgo::pileupCleaningID_

private

Definition at line 100 of file HcalSimpleRecAlgo.h.

Referenced by setRecoParams().

std::unique_ptr<PulseShapeFitOOTPileupCorrection> HcalSimpleRecAlgo::psFitOOTpuCorr_

private

Definition at line 111 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), setForData(), setpuCorrMethod(), and setpuCorrParams().

int HcalSimpleRecAlgo::puCorrMethod_

private

Definition at line 109 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), setForData(), and setpuCorrMethod().

std::unique_ptr<HcalPulseContainmentManager> HcalSimpleRecAlgo::pulseCorr_

private

Definition at line 97 of file HcalSimpleRecAlgo.h.

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

int HcalSimpleRecAlgo::runnum_

private

Definition at line 98 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), and setForData().

bool HcalSimpleRecAlgo::setLeakCorrection_

private

Definition at line 99 of file HcalSimpleRecAlgo.h.

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

HcalPulseShapes HcalSimpleRecAlgo::theHcalPulseShapes_

private

Definition at line 107 of file HcalSimpleRecAlgo.h.

Referenced by setForData().