SimpleHBHEPhase1Algo Class Reference

#include <SimpleHBHEPhase1Algo.h>

Inheritance diagram for SimpleHBHEPhase1Algo:

Public Member Functions
void	beginRun (const edm::Run &, const edm::EventSetup &) override
void	endRun () override
int	getFirstSampleShift () const
float	getPhaseNS () const
int	getRunNumber () const
int	getSamplesToAdd () const
float	getTimeShift () const
bool	isConfigurable () const override
bool	isCorrectingForPhaseContainment () const
HBHERecHit	reconstruct (const HBHEChannelInfo &info, const HcalRecoParam *params, const HcalCalibrations &calibs, bool isRealData) override
	SimpleHBHEPhase1Algo (int firstSampleShift, int samplesToAdd, float phaseNS, float timeShift, bool correctForPhaseContainment, bool applyLegacyHBMCorrection, bool applyFixPCC, std::unique_ptr< PulseShapeFitOOTPileupCorrection > m2, std::unique_ptr< HcalDeterministicFit > detFit, std::unique_ptr< MahiFit > mahi, edm::ConsumesCollector iC)
	~SimpleHBHEPhase1Algo () override
Public Member Functions inherited from AbsHBHEPhase1Algo
virtual bool	configure (const AbsHcalAlgoData *)
virtual	~AbsHBHEPhase1Algo ()

Public Attributes
edm::ESGetToken< HcalTimeSlew, HcalTimeSlewRecord >	delayToken_
const HcalTimeSlew *	hcalTimeSlew_delay_

Protected Member Functions
float	hbminusCorrectionFactor (const HcalDetId &cell, float energy, bool isRealData) const
float	m0Energy (const HBHEChannelInfo &info, double reconstructedCharge, bool applyContainmentCorrection, double phaseNS, int nSamplesToAdd)
float	m0Time (const HBHEChannelInfo &info, double reconstructedCharge, int nSamplesToExamine) const

Private Attributes
bool	applyLegacyHBMCorrection_
bool	corrFPC_
int	firstSampleShift_
std::unique_ptr < HcalDeterministicFit >	hltOOTpuCorr_
std::unique_ptr< MahiFit >	mahiOOTpuCorr_
float	phaseNS_
std::unique_ptr < PulseShapeFitOOTPileupCorrection >	psFitOOTpuCorr_
HcalPulseContainmentManager	pulseCorr_
int	runnum_
int	samplesToAdd_
HcalPulseShapes	theHcalPulseShapes_
float	timeShift_

Detailed Description

Definition at line 18 of file SimpleHBHEPhase1Algo.h.

Constructor & Destructor Documentation

SimpleHBHEPhase1Algo::SimpleHBHEPhase1Algo	(	int	firstSampleShift,
		int	samplesToAdd,
		float	phaseNS,
		float	timeShift,
		bool	correctForPhaseContainment,
		bool	applyLegacyHBMCorrection,
		bool	applyFixPCC,
		std::unique_ptr< PulseShapeFitOOTPileupCorrection >	m2,
		std::unique_ptr< HcalDeterministicFit >	detFit,
		std::unique_ptr< MahiFit >	mahi,
		edm::ConsumesCollector	iC
	)

Definition at line 22 of file SimpleHBHEPhase1Algo.cc.

References hcalTimeSlew_delay_.

    : delayToken_(iC.esConsumes<edm::Transition::BeginRun>(edm::ESInputTag("", "HBHE"))),
      pulseCorr_(PulseContainmentFractionalError, applyFixPCC, iC),
      firstSampleShift_(firstSampleShift),
      samplesToAdd_(samplesToAdd),
      phaseNS_(phaseNS),
      timeShift_(timeShift),
      runnum_(0),
      corrFPC_(correctForPhaseContainment),
      applyLegacyHBMCorrection_(applyLegacyHBMCorrection),
      psFitOOTpuCorr_(std::move(m2)),
      hltOOTpuCorr_(std::move(detFit)),
      mahiOOTpuCorr_(std::move(mahi)) {
  hcalTimeSlew_delay_ = nullptr;
}

SimpleHBHEPhase1Algo::~SimpleHBHEPhase1Algo ( )

inlineoverride

Definition at line 54 of file SimpleHBHEPhase1Algo.h.

{}

Member Function Documentation

void SimpleHBHEPhase1Algo::beginRun ( const edm::Run &  r, const edm::EventSetup &  es  )

overridevirtual

Reimplemented from AbsHBHEPhase1Algo.

Definition at line 48 of file SimpleHBHEPhase1Algo.cc.

References HcalPulseContainmentManager::beginRun(), delayToken_, edm::EventSetup::getData(), hcalTimeSlew_delay_, pulseCorr_, edm::RunBase::run(), and runnum_.

                                                                            {
  hcalTimeSlew_delay_ = &es.getData(delayToken_);

  runnum_ = r.run();
  pulseCorr_.beginRun(es);
}

void SimpleHBHEPhase1Algo::endRun ( )

overridevirtual

Reimplemented from AbsHBHEPhase1Algo.

Definition at line 55 of file SimpleHBHEPhase1Algo.cc.

References runnum_.

{ runnum_ = 0; }

int SimpleHBHEPhase1Algo::getFirstSampleShift ( ) const

inline

Definition at line 67 of file SimpleHBHEPhase1Algo.h.

References firstSampleShift_.

{ return firstSampleShift_; }

float SimpleHBHEPhase1Algo::getPhaseNS ( ) const

inline

Definition at line 69 of file SimpleHBHEPhase1Algo.h.

References phaseNS_.

{ return phaseNS_; }

int SimpleHBHEPhase1Algo::getRunNumber ( ) const

inline

Definition at line 72 of file SimpleHBHEPhase1Algo.h.

References runnum_.

{ return runnum_; }

int SimpleHBHEPhase1Algo::getSamplesToAdd ( ) const

inline

Definition at line 68 of file SimpleHBHEPhase1Algo.h.

References samplesToAdd_.

{ return samplesToAdd_; }

float SimpleHBHEPhase1Algo::getTimeShift ( ) const

inline

Definition at line 70 of file SimpleHBHEPhase1Algo.h.

References timeShift_.

{ return timeShift_; }

float SimpleHBHEPhase1Algo::hbminusCorrectionFactor ( const HcalDetId &  cell, float  energy, bool  isRealData  ) const

protected

Definition at line 155 of file SimpleHBHEPhase1Algo.cc.

References applyLegacyHBMCorrection_, alignCSCRings::corr, hbminus_special_ecorr(), HcalBarrel, HcalDetId::ieta(), HcalDetId::iphi(), runnum_, and HcalDetId::subdet().

Referenced by reconstruct().

                                                                                 {
  float corr = 1.f;
  if (applyLegacyHBMCorrection_ && isRealData && runnum_ > 0)
    if (cell.subdet() == HcalBarrel) {
      const int ieta = cell.ieta();
      const int iphi = cell.iphi();
      corr = hbminus_special_ecorr(ieta, iphi, energy, runnum_);
    }
  return corr;
}

bool SimpleHBHEPhase1Algo::isConfigurable ( ) const

inlineoverridevirtual

Implements AbsHBHEPhase1Algo.

Definition at line 60 of file SimpleHBHEPhase1Algo.h.

{ return false; }

bool SimpleHBHEPhase1Algo::isCorrectingForPhaseContainment ( ) const

inline

Definition at line 71 of file SimpleHBHEPhase1Algo.h.

References corrFPC_.

{ return corrFPC_; }

float SimpleHBHEPhase1Algo::m0Energy ( const HBHEChannelInfo &  info, double  reconstructedCharge, bool  applyContainmentCorrection, double  phaseNS, int  nSamplesToAdd  )

protected

Definition at line 168 of file SimpleHBHEPhase1Algo.cc.

References alignCSCRings::e, HBHEChannelInfo::energyInWindow(), firstSampleShift_, HcalPulseContainmentManager::get(), HBHEChannelInfo::id(), pulseCorr_, and HBHEChannelInfo::soi().

Referenced by reconstruct().

                                                              {
  int ibeg = static_cast<int>(info.soi()) + firstSampleShift_;
  if (ibeg < 0)
    ibeg = 0;
  double e = info.energyInWindow(ibeg, ibeg + nSamplesToAdd);

  // Pulse containment correction
  {
    double corrFactor = 1.0;
    if (applyContainmentCorrection)
      corrFactor = pulseCorr_.get(info.id(), nSamplesToAdd, phaseNs)->getCorrection(fc_ampl);
    e *= corrFactor;
  }

  return e;
}

float SimpleHBHEPhase1Algo::m0Time ( const HBHEChannelInfo &  info, double  reconstructedCharge, int  nSamplesToExamine  ) const

protected

Definition at line 189 of file SimpleHBHEPhase1Algo.cc.

References HcalTimeSlew::delay(), validate-o2o-wbm::f, firstSampleShift_, hcalTimeSlew_delay_, SiStripPI::max, findQualityFiles::maxI, HBHEChannelInfo::MAXSAMPLES, HcalTimeSlew::Medium, min(), hgc_digi::nSamples, HBHEChannelInfo::nSamples(), HBHEChannelInfo::peakEnergyTS(), position, HBHEChannelInfo::soi(), and HBHEChannelInfo::tsEnergy().

Referenced by reconstruct().

                                                                      {
  float time = -9999.f;  // historic value

  const unsigned nSamples = info.nSamples();
  if (nSamples > 2U) {
    const int soi = info.soi();
    int ibeg = soi + firstSampleShift_;
    if (ibeg < 0)
      ibeg = 0;
    const int iend = std::min(ibeg + nSamplesToExamine, (int)nSamples - 1);  // actual array

    unsigned maxI = info.peakEnergyTS((unsigned)ibeg, (unsigned)iend);  // requires unsigned params
    if (maxI < HBHEChannelInfo::MAXSAMPLES) {
      if (maxI >= nSamples)
        maxI = nSamples - 1U;  // just in case

      // Simplified evaluation for Phase1
      float emax0 = info.tsEnergy(maxI);
      float emax1 = 0.f;
      if (maxI < (nSamples - 1U))
        emax1 = info.tsEnergy(maxI + 1U);

      // consider soi reference for collisions
      int position = (int)maxI;
      if (nSamplesToExamine < (int)nSamples)
        position -= soi;

      time = 25.f * (float)position;
      if (emax0 > 0.f && emax1 > 0.f)
        time += 25.f * emax1 / (emax0 + emax1);  // 1st order corr.

      // TimeSlew correction
      time -= hcalTimeSlew_delay_->delay(std::max(1.0, fc_ampl), HcalTimeSlew::Medium);
    }
  }
  return time;
}

HBHERecHit SimpleHBHEPhase1Algo::reconstruct ( const HBHEChannelInfo &  info, const HcalRecoParam *  params, const HcalCalibrations &  calibs, bool  isRealData  )

overridevirtual

Implements AbsHBHEPhase1Algo.

Definition at line 57 of file SimpleHBHEPhase1Algo.cc.

References HBHEChannelInfo::chargeInWindow(), HLT_FULL_cff::chi2, HcalRecoParam::correctForPhaseContainment(), HcalRecoParam::correctionPhaseNS(), corrFPC_, firstSampleShift_, HcalPulseShapes::getShape(), HBHEChannelInfo::hasTimeInfo(), HcalPhase1FlagLabels::HBHEPulseFitBit, hbminusCorrectionFactor(), hcalTimeSlew_delay_, hltOOTpuCorr_, HBHEChannelInfo::id(), HcalSpecialTimes::isSpecial(), m0Energy(), m0Time(), mahiOOTpuCorr_, HBHEChannelInfo::nSamples(), HcalDeterministicFit::phase1Apply(), PulseShapeFitOOTPileupCorrection::phase1Apply(), MahiFit::phase1Apply(), phaseNS_, psFitOOTpuCorr_, HBHEChannelInfo::recoShape(), HcalRecoParam::samplesToAdd(), samplesToAdd_, HBHERecHitAuxSetter::setAux(), HBHERecHit::setAuxEnergy(), HBHERecHit::setChiSquared(), CaloRecHit::setFlagField(), HBHERecHit::setRawEnergy(), HBHEChannelInfo::soi(), HBHEChannelInfo::soiRiseTime(), theHcalPulseShapes_, and timeShift_.

                                                                {
  const HcalDetId channelId(info.id());

  // Calculate "Method 0" quantities
  float m0t = 0.f, m0E = 0.f;
  {
    int ibeg = static_cast<int>(info.soi()) + firstSampleShift_;
    if (ibeg < 0)
      ibeg = 0;
    const int nSamplesToAdd = params ? params->samplesToAdd() : samplesToAdd_;
    const double fc_ampl = info.chargeInWindow(ibeg, ibeg + nSamplesToAdd);
    const bool applyContainment = params ? params->correctForPhaseContainment() : corrFPC_;
    const float phasens = params ? params->correctionPhaseNS() : phaseNS_;
    m0E = m0Energy(info, fc_ampl, applyContainment, phasens, nSamplesToAdd);
    m0E *= hbminusCorrectionFactor(channelId, m0E, isData);
    m0t = m0Time(info, fc_ampl, nSamplesToAdd);
  }

  // Run "Method 2"
  float m2t = 0.f, m2E = 0.f, chi2 = -1.f;
  bool useTriple = false;
  const PulseShapeFitOOTPileupCorrection* method2 = psFitOOTpuCorr_.get();
  if (method2) {
    psFitOOTpuCorr_->setPulseShapeTemplate(
        theHcalPulseShapes_.getShape(info.recoShape()), !info.hasTimeInfo(), info.nSamples(), hcalTimeSlew_delay_);
    // "phase1Apply" call below sets m2E, m2t, useTriple, and chi2.
    // These parameters are pased by non-const reference.
    method2->phase1Apply(info, m2E, m2t, useTriple, chi2);
    m2E *= hbminusCorrectionFactor(channelId, m2E, isData);
  }

  // Run "Method 3"
  float m3t = 0.f, m3E = 0.f;
  const HcalDeterministicFit* method3 = hltOOTpuCorr_.get();
  if (method3) {
    // "phase1Apply" sets m3E and m3t (pased by non-const reference)
    method3->phase1Apply(info, m3E, m3t, hcalTimeSlew_delay_);
    m3E *= hbminusCorrectionFactor(channelId, m3E, isData);
  }

  // Run Mahi
  float m4E = 0.f, m4ESOIPlusOne = 0.f, m4chi2 = -1.f;
  float m4T = 0.f;
  bool m4UseTriple = false;

  const MahiFit* mahi = mahiOOTpuCorr_.get();

  if (mahi) {
    mahiOOTpuCorr_->setPulseShapeTemplate(
        info.recoShape(), theHcalPulseShapes_, info.hasTimeInfo(), hcalTimeSlew_delay_, info.nSamples());
    mahi->phase1Apply(info, m4E, m4ESOIPlusOne, m4T, m4UseTriple, m4chi2);
    m4E *= hbminusCorrectionFactor(channelId, m4E, isData);
    m4ESOIPlusOne *= hbminusCorrectionFactor(channelId, m4ESOIPlusOne, isData);
  }

  // Finally, construct the rechit
  HBHERecHit rh;

  float rhE = m0E;
  float rht = m0t;
  float rhX = -1.f;
  if (mahi) {
    rhE = m4E;
    rht = m4T;
    rhX = m4chi2;
  } else if (method2) {
    rhE = m2E;
    rht = m2t;
    rhX = chi2;
  } else if (method3) {
    rhE = m3E;
    rht = m3t;
  }
  float tdcTime = info.soiRiseTime();
  if (!HcalSpecialTimes::isSpecial(tdcTime))
    tdcTime += timeShift_;
  rh = HBHERecHit(channelId, rhE, rht, tdcTime);
  rh.setRawEnergy(m0E);
  if (method3) {
    rh.setAuxEnergy(m3E);
  } else {
    rh.setAuxEnergy(m4ESOIPlusOne);
  }
  rh.setChiSquared(rhX);

  // Set rechit aux words
  HBHERecHitAuxSetter::setAux(info, &rh);

  // Set some rechit flags (here, for Method 2/Mahi)
  if (useTriple || m4UseTriple)
    rh.setFlagField(1, HcalPhase1FlagLabels::HBHEPulseFitBit);

  return rh;
}

Member Data Documentation

bool SimpleHBHEPhase1Algo::applyLegacyHBMCorrection_

private

Definition at line 102 of file SimpleHBHEPhase1Algo.h.

Referenced by hbminusCorrectionFactor().

bool SimpleHBHEPhase1Algo::corrFPC_

private

Definition at line 101 of file SimpleHBHEPhase1Algo.h.

Referenced by isCorrectingForPhaseContainment(), and reconstruct().

edm::ESGetToken<HcalTimeSlew, HcalTimeSlewRecord> SimpleHBHEPhase1Algo::delayToken_

Definition at line 74 of file SimpleHBHEPhase1Algo.h.

Referenced by beginRun().

int SimpleHBHEPhase1Algo::firstSampleShift_

private

Definition at line 96 of file SimpleHBHEPhase1Algo.h.

Referenced by getFirstSampleShift(), m0Energy(), m0Time(), and reconstruct().

const HcalTimeSlew* SimpleHBHEPhase1Algo::hcalTimeSlew_delay_

Definition at line 75 of file SimpleHBHEPhase1Algo.h.

Referenced by beginRun(), m0Time(), reconstruct(), and SimpleHBHEPhase1Algo().

std::unique_ptr<HcalDeterministicFit> SimpleHBHEPhase1Algo::hltOOTpuCorr_

private

Definition at line 108 of file SimpleHBHEPhase1Algo.h.

Referenced by reconstruct().

std::unique_ptr<MahiFit> SimpleHBHEPhase1Algo::mahiOOTpuCorr_

private

Definition at line 111 of file SimpleHBHEPhase1Algo.h.

Referenced by reconstruct().

float SimpleHBHEPhase1Algo::phaseNS_

private

Definition at line 98 of file SimpleHBHEPhase1Algo.h.

Referenced by getPhaseNS(), and reconstruct().

std::unique_ptr<PulseShapeFitOOTPileupCorrection> SimpleHBHEPhase1Algo::psFitOOTpuCorr_

private

Definition at line 105 of file SimpleHBHEPhase1Algo.h.

Referenced by reconstruct().

HcalPulseContainmentManager SimpleHBHEPhase1Algo::pulseCorr_

private

Definition at line 94 of file SimpleHBHEPhase1Algo.h.

Referenced by beginRun(), and m0Energy().

int SimpleHBHEPhase1Algo::runnum_

private

Definition at line 100 of file SimpleHBHEPhase1Algo.h.

Referenced by beginRun(), endRun(), getRunNumber(), and hbminusCorrectionFactor().

int SimpleHBHEPhase1Algo::samplesToAdd_

private

Definition at line 97 of file SimpleHBHEPhase1Algo.h.

Referenced by getSamplesToAdd(), and reconstruct().

HcalPulseShapes SimpleHBHEPhase1Algo::theHcalPulseShapes_

private

Definition at line 113 of file SimpleHBHEPhase1Algo.h.

Referenced by reconstruct().

float SimpleHBHEPhase1Algo::timeShift_

private

Definition at line 99 of file SimpleHBHEPhase1Algo.h.

Referenced by getTimeShift(), and reconstruct().