CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
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)
 
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 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< AbsOOTPileupCorrectionhbhePileupCorr_
 
const HcalTimeSlewhcalTimeSlew_delay_
 
boost::shared_ptr< AbsOOTPileupCorrectionhfPileupCorr_
 
boost::shared_ptr< AbsOOTPileupCorrectionhoPileupCorr_
 
unsigned lenBunchCrossingInfo_
 
float phaseNS_
 
int pileupCleaningID_
 
int puCorrMethod_
 
std::unique_ptr< HcalPulseContainmentManagerpulseCorr_
 
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 35 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 20 of file HcalSimpleRecAlgo.cc.

References hcalTimeSlew_delay_, MaximumFractionalError, and pulseCorr_.

20  :
22  correctForPulse_(correctForPulse),
23  phaseNS_(phaseNS), runnum_(0), setLeakCorrection_(false), puCorrMethod_(0)
24 {
25  hcalTimeSlew_delay_ = nullptr;
26  pulseCorr_ = std::make_unique<HcalPulseContainmentManager>(MaximumFractionalError);
27 }
std::unique_ptr< HcalPulseContainmentManager > pulseCorr_
double MaximumFractionalError
const HcalTimeSlew * hcalTimeSlew_delay_

Member Function Documentation

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

Definition at line 30 of file HcalSimpleRecAlgo.cc.

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

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

31 {
33  es.get<HcalTimeSlewRecord>().get("HBHE", delay);
34  hcalTimeSlew_delay_ = &*delay;
35 
36  pulseCorr_->beginRun(es);
37 }
std::unique_ptr< HcalPulseContainmentManager > pulseCorr_
const HcalTimeSlew * hcalTimeSlew_delay_
void HcalSimpleRecAlgo::endRun ( )

Definition at line 40 of file HcalSimpleRecAlgo.cc.

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

41 {
42 }
void HcalSimpleRecAlgo::initPulseCorr ( int  toadd)

Definition at line 45 of file HcalSimpleRecAlgo.cc.

45  {
46 }
HFRecHit HcalSimpleRecAlgo::reconstruct ( const HFDataFrame digi,
int  first,
int  toadd,
const HcalCoder coder,
const HcalCalibrations calibs 
) const

Definition at line 352 of file HcalSimpleRecAlgo.cc.

References bunchCrossingInfo_, corr, correctForPulse_, correctForTimeslew_, hfPileupCorr_, HFDataFrame::id(), lenBunchCrossingInfo_, findQualityFiles::maxI, phaseNS_, pulseCorr_, HcalSimpleRecAlgoImpl::recoHFTime(), HcalSimpleRecAlgoImpl::removePileup(), setRawEnergy(), genVertex_cff::t0, protons_cff::time, and HcalCalibrations::timecorr().

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

357 {
358  const HcalPulseContainmentCorrection* corr = pulseCorr_->get(digi.id(), toadd, phaseNS_);
359 
360  double amp_fC, ampl, uncorr_ampl, maxA;
361  int nRead, maxI;
362  bool leakCorrApplied;
363  float t0, t2;
364 
365  HcalSimpleRecAlgoImpl::removePileup(digi, coder, calibs, first, toadd,
366  correctForPulse_, corr, hfPileupCorr_.get(),
368  &maxA, &ampl, &uncorr_ampl, &amp_fC, &nRead,
369  &maxI, &leakCorrApplied, &t0, &t2);
370 
371  float time=-9999.f;
372  if (maxI > 0 && maxI < (nRead - 1))
373  time = HcalSimpleRecAlgoImpl::recoHFTime(digi,maxI,amp_fC,correctForTimeslew_,maxA,t0,t2) -
374  calibs.timecorr();
375 
376  HFRecHit rh(digi.id(),ampl,time);
377  setRawEnergy(rh, static_cast<float>(uncorr_ampl));
378  return rh;
379 }
const BunchXParameter * bunchCrossingInfo_
std::unique_ptr< HcalPulseContainmentManager > pulseCorr_
void setRawEnergy(HcalRecHit &h, float e)
Definition: rawEnergy.h:215
HcalDetId const & id() const
Definition: HFDataFrame.h:26
unsigned lenBunchCrossingInfo_
void removePileup(const Digi &digi, const HcalCoder &coder, const HcalCalibrations &calibs, const int ifirst, const int n, const bool pulseCorrect, const HcalPulseContainmentCorrection *corr, const AbsOOTPileupCorrection *pileupCorrection, const BunchXParameter *bxInfo, const unsigned lenInfo, double *p_maxA, double *p_ampl, double *p_uncorr_ampl, double *p_fc_ampl, int *p_nRead, int *p_maxI, bool *leakCorrApplied, float *p_t0, float *p_t2)
JetCorrectorParameters corr
Definition: classes.h:5
constexpr double timecorr() const
get time correction factor
boost::shared_ptr< AbsOOTPileupCorrection > hfPileupCorr_
float recoHFTime(const Digi &digi, const int maxI, const double amp_fC, const bool slewCorrect, double maxA, float t0, float t2)
HORecHit HcalSimpleRecAlgo::reconstruct ( const HODataFrame digi,
int  first,
int  toadd,
const HcalCoder coder,
const HcalCalibrations calibs 
) const

Definition at line 328 of file HcalSimpleRecAlgo.cc.

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

328  {
329  return HcalSimpleRecAlgoImpl::reco<HODataFrame,HORecHit>(digi,coder,calibs,
331  pulseCorr_->get(digi.id(), toadd, phaseNS_),
333  runnum_, false, hoPileupCorr_.get(),
337 }
const BunchXParameter * bunchCrossingInfo_
boost::shared_ptr< AbsOOTPileupCorrection > hoPileupCorr_
std::unique_ptr< HcalPulseContainmentManager > pulseCorr_
unsigned lenBunchCrossingInfo_
const HcalTimeSlew * hcalTimeSlew_delay_
HcalCalibRecHit HcalSimpleRecAlgo::reconstruct ( const HcalCalibDataFrame digi,
int  first,
int  toadd,
const HcalCoder coder,
const HcalCalibrations calibs 
) const

Definition at line 340 of file HcalSimpleRecAlgo.cc.

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

340  {
341  return HcalSimpleRecAlgoImpl::reco<HcalCalibDataFrame,HcalCalibRecHit>(digi,coder,calibs,
343  pulseCorr_->get(digi.id(), toadd, phaseNS_),
345  runnum_, false, nullptr,
349 }
const BunchXParameter * bunchCrossingInfo_
std::unique_ptr< HcalPulseContainmentManager > pulseCorr_
unsigned lenBunchCrossingInfo_
const HcalTimeSlew * hcalTimeSlew_delay_
HFRecHit HcalSimpleRecAlgo::reconstructQIE10 ( const QIE10DataFrame digi,
int  first,
int  toadd,
const HcalCoder coder,
const HcalCalibrations calibs 
) const

Definition at line 381 of file HcalSimpleRecAlgo.cc.

References bunchCrossingInfo_, corr, correctForPulse_, correctForTimeslew_, hfPileupCorr_, QIE10DataFrame::id(), lenBunchCrossingInfo_, findQualityFiles::maxI, phaseNS_, pulseCorr_, HcalSimpleRecAlgoImpl::recoHFTime(), HcalSimpleRecAlgoImpl::removePileup(), setRawEnergy(), genVertex_cff::t0, protons_cff::time, and HcalCalibrations::timecorr().

386 {
387  const HcalPulseContainmentCorrection* corr = pulseCorr_->get(digi.id(), toadd, phaseNS_);
388 
389  double amp_fC, ampl, uncorr_ampl, maxA;
390  int nRead, maxI;
391  bool leakCorrApplied;
392  float t0, t2;
393 
394  HcalSimpleRecAlgoImpl::removePileup(digi, coder, calibs, first, toadd,
395  correctForPulse_, corr, hfPileupCorr_.get(),
397  &maxA, &ampl, &uncorr_ampl, &amp_fC, &nRead,
398  &maxI, &leakCorrApplied, &t0, &t2);
399 
400  float time=-9999.f;
401  if (maxI > 0 && maxI < (nRead - 1))
402  time = HcalSimpleRecAlgoImpl::recoHFTime(digi,maxI,amp_fC,correctForTimeslew_,maxA,t0,t2) -
403  calibs.timecorr();
404 
405  HFRecHit rh(digi.id(),ampl,time);
406  setRawEnergy(rh, static_cast<float>(uncorr_ampl));
407  return rh;
408 }
edm::DataFrame::id_type id() const
const BunchXParameter * bunchCrossingInfo_
std::unique_ptr< HcalPulseContainmentManager > pulseCorr_
void setRawEnergy(HcalRecHit &h, float e)
Definition: rawEnergy.h:215
unsigned lenBunchCrossingInfo_
void removePileup(const Digi &digi, const HcalCoder &coder, const HcalCalibrations &calibs, const int ifirst, const int n, const bool pulseCorrect, const HcalPulseContainmentCorrection *corr, const AbsOOTPileupCorrection *pileupCorrection, const BunchXParameter *bxInfo, const unsigned lenInfo, double *p_maxA, double *p_ampl, double *p_uncorr_ampl, double *p_fc_ampl, int *p_nRead, int *p_maxI, bool *leakCorrApplied, float *p_t0, float *p_t2)
JetCorrectorParameters corr
Definition: classes.h:5
constexpr double timecorr() const
get time correction factor
boost::shared_ptr< AbsOOTPileupCorrection > hfPileupCorr_
float recoHFTime(const Digi &digi, const int maxI, const double amp_fC, const bool slewCorrect, double maxA, float t0, float t2)
void HcalSimpleRecAlgo::setBXInfo ( const BunchXParameter *  info,
unsigned  lenInfo 
)

Definition at line 70 of file HcalSimpleRecAlgo.cc.

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

72 {
74  lenBunchCrossingInfo_ = lenInfo;
75 }
static const TGPicture * info(bool iBackgroundIsBlack)
const BunchXParameter * bunchCrossingInfo_
unsigned lenBunchCrossingInfo_
void HcalSimpleRecAlgo::setHFPileupCorrection ( boost::shared_ptr< AbsOOTPileupCorrection corr)

Definition at line 58 of file HcalSimpleRecAlgo.cc.

References corr, and hfPileupCorr_.

60 {
62 }
JetCorrectorParameters corr
Definition: classes.h:5
boost::shared_ptr< AbsOOTPileupCorrection > hfPileupCorr_
void HcalSimpleRecAlgo::setHOPileupCorrection ( boost::shared_ptr< AbsOOTPileupCorrection corr)

Definition at line 64 of file HcalSimpleRecAlgo.cc.

References corr, and hoPileupCorr_.

66 {
68 }
boost::shared_ptr< AbsOOTPileupCorrection > hoPileupCorr_
JetCorrectorParameters corr
Definition: classes.h:5
void HcalSimpleRecAlgo::setLeakCorrection ( )

Definition at line 56 of file HcalSimpleRecAlgo.cc.

References setLeakCorrection_.

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

56 { setLeakCorrection_ = true;}
void HcalSimpleRecAlgo::setRecoParams ( bool  correctForTimeslew,
bool  correctForPulse,
bool  setLeakCorrection,
int  pileupCleaningID,
float  phaseNS 
)

Member Data Documentation

const BunchXParameter* HcalSimpleRecAlgo::bunchCrossingInfo_
private

Definition at line 75 of file HcalSimpleRecAlgo.h.

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

bool HcalSimpleRecAlgo::correctForPulse_
private

Definition at line 69 of file HcalSimpleRecAlgo.h.

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

bool HcalSimpleRecAlgo::correctForTimeslew_
private

Definition at line 68 of file HcalSimpleRecAlgo.h.

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

boost::shared_ptr<AbsOOTPileupCorrection> HcalSimpleRecAlgo::hbhePileupCorr_
private

Definition at line 77 of file HcalSimpleRecAlgo.h.

const HcalTimeSlew* HcalSimpleRecAlgo::hcalTimeSlew_delay_
private

Definition at line 85 of file HcalSimpleRecAlgo.h.

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

boost::shared_ptr<AbsOOTPileupCorrection> HcalSimpleRecAlgo::hfPileupCorr_
private

Definition at line 78 of file HcalSimpleRecAlgo.h.

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

boost::shared_ptr<AbsOOTPileupCorrection> HcalSimpleRecAlgo::hoPileupCorr_
private

Definition at line 79 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct(), and setHOPileupCorrection().

unsigned HcalSimpleRecAlgo::lenBunchCrossingInfo_
private

Definition at line 76 of file HcalSimpleRecAlgo.h.

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

float HcalSimpleRecAlgo::phaseNS_
private

Definition at line 70 of file HcalSimpleRecAlgo.h.

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

int HcalSimpleRecAlgo::pileupCleaningID_
private

Definition at line 74 of file HcalSimpleRecAlgo.h.

Referenced by setRecoParams().

int HcalSimpleRecAlgo::puCorrMethod_
private

Definition at line 83 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct().

std::unique_ptr<HcalPulseContainmentManager> HcalSimpleRecAlgo::pulseCorr_
private

Definition at line 71 of file HcalSimpleRecAlgo.h.

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

int HcalSimpleRecAlgo::runnum_
private

Definition at line 72 of file HcalSimpleRecAlgo.h.

Referenced by reconstruct().

bool HcalSimpleRecAlgo::setLeakCorrection_
private

Definition at line 73 of file HcalSimpleRecAlgo.h.

Referenced by setLeakCorrection(), and setRecoParams().

HcalPulseShapes HcalSimpleRecAlgo::theHcalPulseShapes_
private

Definition at line 81 of file HcalSimpleRecAlgo.h.