HFSimpleTimeCheck Class Reference

#include <HFSimpleTimeCheck.h>

Inheritance diagram for HFSimpleTimeCheck:

Public Member Functions
bool	alwaysCalculatingQAsym () const
	HFSimpleTimeCheck (const std::pair< float, float > tlimits[2], const float energyWeights[2 *HFAnodeStatus::N_POSSIBLE_STATES - 1][2], unsigned soiPhase, float timeShift, float triseIfNoTDC, float tfallIfNoTDC, float minChargeForUndershoot, float minChargeForOvershoot, bool rejectAllFailures=true, bool alwaysCalculateChargeAsymmetry=true)
bool	isConfigurable () const override
float	minChargeForOvershoot () const
float	minChargeForUndershoot () const
HFRecHit	reconstruct (const HFPreRecHit &prehit, const HcalCalibrations &calibs, const bool flaggedBadInDB[2], bool expectSingleAnodePMT) override
bool	rejectingAllFailures () const
unsigned	soiPhase () const
float	tfallIfNoTDC () const
float	timeShift () const
float	triseIfNoTDC () const
	~HFSimpleTimeCheck () override
Public Member Functions inherited from AbsHFPhase1Algo
virtual bool	configure (const AbsHcalAlgoData *)
virtual	~AbsHFPhase1Algo ()

Protected Member Functions
virtual unsigned	determineAnodeStatus (unsigned anodeNumber, const HFQIE10Info &anode, bool *isTimingReliable) const

Private Member Functions
unsigned	mapStatusIntoIndex (const unsigned states[2]) const

Private Attributes
bool	alwaysQAsym_
float	energyWeights_ [2 *HFAnodeStatus::N_POSSIBLE_STATES - 1][2]
float	minChargeForOvershoot_
float	minChargeForUndershoot_
bool	rejectAllFailures_
unsigned	soiPhase_
float	tfallIfNoTDC_
float	timeShift_
std::pair< float, float >	tlimits_ [2]
float	triseIfNoTDC_

Detailed Description

Definition at line 9 of file HFSimpleTimeCheck.h.

Constructor & Destructor Documentation

HFSimpleTimeCheck()

HFSimpleTimeCheck::HFSimpleTimeCheck	(	const std::pair< float, float >	tlimits[2],
		const float	energyWeights[2 *HFAnodeStatus::N_POSSIBLE_STATES - 1][2],
		unsigned	soiPhase,
		float	timeShift,
		float	triseIfNoTDC,
		float	tfallIfNoTDC,
		float	minChargeForUndershoot,
		float	minChargeForOvershoot,
		bool	rejectAllFailures = true,
		bool	alwaysCalculateChargeAsymmetry = true
	)

Definition at line 32 of file HFSimpleTimeCheck.cc.

References HLT_2023v12_cff::energyWeights, energyWeights_, HLT_2023v12_cff::tlimits, and tlimits_.

    : soiPhase_(i_soiPhase),
      timeShift_(i_timeShift),
      triseIfNoTDC_(i_triseIfNoTDC),
      tfallIfNoTDC_(i_tfallIfNoTDC),
      minChargeForUndershoot_(i_minChargeForUndershoot),
      minChargeForOvershoot_(i_minChargeForOvershoot),
      rejectAllFailures_(rejectAllFailures),
      alwaysQAsym_(alwaysCalculateQAsymmetry) {
  tlimits_[0] = tlimits[0];
  tlimits_[1] = tlimits[1];
  float* to = &energyWeights_[0][0];
  const float* from = &energyWeights[0][0];
  memcpy(to, from, sizeof(energyWeights_));
}

~HFSimpleTimeCheck()

HFSimpleTimeCheck::~HFSimpleTimeCheck ( )

inlineoverride

Definition at line 71 of file HFSimpleTimeCheck.h.

{}

Member Function Documentation

alwaysCalculatingQAsym()

bool HFSimpleTimeCheck::alwaysCalculatingQAsym ( ) const

inline

Definition at line 87 of file HFSimpleTimeCheck.h.

References alwaysQAsym_.

Referenced by HFFlexibleTimeCheck::reconstruct().

{ return alwaysQAsym_; }

determineAnodeStatus()

unsigned HFSimpleTimeCheck::determineAnodeStatus ( unsigned  anodeNumber, const HFQIE10Info &  anode, bool *  isTimingReliable  ) const

protectedvirtual

Reimplemented in HFFlexibleTimeCheck.

Definition at line 57 of file HFSimpleTimeCheck.cc.

References HFAnodeStatus::FAILED_TIMING, dqmdumpme::first, HFAnodeStatus::HARDWARE_ERROR, HFQIE10Info::isDataframeOK(), HcalSpecialTimes::isSpecial(), HFAnodeStatus::OK, edm::second(), HFQIE10Info::timeRising(), timeShift_, and tlimits_.

Referenced by reconstruct().

                                                                                                             {
  // Check if this anode has a dataframe error
  if (!anode.isDataframeOK())
    return HFAnodeStatus::HARDWARE_ERROR;

  // Check the time limits
  float trise = anode.timeRising();
  const bool timeIsKnown = !HcalSpecialTimes::isSpecial(trise);
  trise += timeShift_;
  if (timeIsKnown && tlimits_[ianode].first <= trise && trise <= tlimits_[ianode].second)
    return HFAnodeStatus::OK;
  else
    return HFAnodeStatus::FAILED_TIMING;
}

isConfigurable()

bool HFSimpleTimeCheck::isConfigurable ( ) const

inlineoverridevirtual

Implements AbsHFPhase1Algo.

Definition at line 73 of file HFSimpleTimeCheck.h.

{ return false; }

mapStatusIntoIndex()

unsigned HFSimpleTimeCheck::mapStatusIntoIndex ( const unsigned  states[2] ) const

private

Definition at line 72 of file HFSimpleTimeCheck.cc.

References HFAnodeStatus::FAILED_OTHER, HFAnodeStatus::FAILED_TIMING, mps_fire::i, HFAnodeStatus::N_POSSIBLE_STATES, HFAnodeStatus::OK, and rejectAllFailures_.

Referenced by reconstruct().

                                                                             {
  unsigned eStates[2];
  eStates[0] = states[0];
  eStates[1] = states[1];
  if (!rejectAllFailures_)
    for (unsigned i = 0; i < 2; ++i)
      if (eStates[i] == HFAnodeStatus::FAILED_TIMING || eStates[i] == HFAnodeStatus::FAILED_OTHER)
        eStates[i] = HFAnodeStatus::OK;
  if (eStates[0] == HFAnodeStatus::OK)
    return eStates[1];
  else if (eStates[1] == HFAnodeStatus::OK)
    return HFAnodeStatus::N_POSSIBLE_STATES + eStates[0] - 1;
  else
    return UINT_MAX;
}

minChargeForOvershoot()

float HFSimpleTimeCheck::minChargeForOvershoot ( ) const

inline

Definition at line 85 of file HFSimpleTimeCheck.h.

References minChargeForOvershoot_.

Referenced by HFFlexibleTimeCheck::determineAnodeStatus().

{ return minChargeForOvershoot_; }

minChargeForUndershoot()

float HFSimpleTimeCheck::minChargeForUndershoot ( ) const

inline

Definition at line 84 of file HFSimpleTimeCheck.h.

References minChargeForUndershoot_.

Referenced by HFFlexibleTimeCheck::determineAnodeStatus().

{ return minChargeForUndershoot_; }

reconstruct()

HFRecHit HFSimpleTimeCheck::reconstruct ( const HFPreRecHit &  prehit, const HcalCalibrations &  calibs, const bool  flaggedBadInDB[2], bool  expectSingleAnodePMT  )

overridevirtual

Implements AbsHFPhase1Algo.

Definition at line 88 of file HFSimpleTimeCheck.cc.

References determineAnodeStatus(), hcalRecHitTable_cff::energy, HFQIE10Info::energy(), energyWeights_, f, HFAnodeStatus::FAILED_TIMING, RemoveAddSevLevel::flag, HFAnodeStatus::FLAGGED_BAD, HFPreRecHit::getHFQIE10Info(), HFPreRecHit::id(), HcalSpecialTimes::isSpecial(), mapStatusIntoIndex(), HFAnodeStatus::NOT_DUAL, HFAnodeStatus::NOT_READ_OUT, HFRecHitAuxSetter::setAux(), CaloRecHit::setFlagField(), soiPhase_, tfallIfNoTDC_, HFQIE10Info::timeFalling(), HFQIE10Info::timeRising(), timeShift_, HcalPhase1FlagLabels::TimingFromTDC, triseIfNoTDC_, mitigatedMETSequence_cff::U, and hltDeepSecondaryVertexTagInfosPFPuppi_cfi::weights.

Referenced by HFFlexibleTimeCheck::reconstruct().

                                                                         {
  HFRecHit rh;

  // Determine the status of each anode
  unsigned states[2] = {HFAnodeStatus::NOT_READ_OUT, HFAnodeStatus::NOT_READ_OUT};
  if (expectSingleAnodePMT)
    states[1] = HFAnodeStatus::NOT_DUAL;

  bool isTimingReliable[2] = {true, true};
  for (unsigned ianode = 0; ianode < 2; ++ianode) {
    if (flaggedBadInDB[ianode])
      states[ianode] = HFAnodeStatus::FLAGGED_BAD;
    else {
      const HFQIE10Info* anodeInfo = prehit.getHFQIE10Info(ianode);
      if (anodeInfo)
        states[ianode] = determineAnodeStatus(ianode, *anodeInfo, &isTimingReliable[ianode]);
    }
  }

  // Reconstruct energy and time
  const unsigned lookupInd = mapStatusIntoIndex(states);
  if (lookupInd != UINT_MAX) {
    // In this scope, at least one of states[i] is HFAnodeStatus::OK
    // or was mapped into that status by "mapStatusIntoIndex" method
    //
    const float* weights = &energyWeights_[lookupInd][0];
    float energy = 0.f, tfallWeightedEnergySum = 0.f, triseWeightedEnergySum = 0.f;
    float tfallWeightedSum = 0.f, triseWeightedSum = 0.f;
    float tfallSum = 0.f, triseSum = 0.f;
    unsigned tfallCount = 0, triseCount = 0;

    for (unsigned ianode = 0; ianode < 2; ++ianode) {
      const HFQIE10Info* anodeInfo = prehit.getHFQIE10Info(ianode);
      if (anodeInfo && weights[ianode] > 0.f) {
        const float weightedEnergy = weights[ianode] * anodeInfo->energy();
        energy += weightedEnergy;

        if (isTimingReliable[ianode] && states[ianode] != HFAnodeStatus::FAILED_TIMING) {
          float trise = anodeInfo->timeRising();
          if (!HcalSpecialTimes::isSpecial(trise)) {
            trise += timeShift_;
            triseSum += trise;
            ++triseCount;
            if (weightedEnergy > 0.f) {
              triseWeightedSum += trise * weightedEnergy;
              triseWeightedEnergySum += weightedEnergy;
            }
          }

          float tfall = anodeInfo->timeFalling();
          if (!HcalSpecialTimes::isSpecial(tfall)) {
            tfall += timeShift_;
            tfallSum += tfall;
            ++tfallCount;
            if (weightedEnergy > 0.f) {
              tfallWeightedSum += tfall * weightedEnergy;
              tfallWeightedEnergySum += weightedEnergy;
            }
          }
        }
      }
    }

    bool triseFromTDC = false;
    const float trise =
        build_rechit_time(triseWeightedEnergySum, triseWeightedSum, triseSum, triseCount, triseIfNoTDC_, &triseFromTDC);

    bool tfallFromTDC = false;
    const float tfall =
        build_rechit_time(tfallWeightedEnergySum, tfallWeightedSum, tfallSum, tfallCount, tfallIfNoTDC_, &tfallFromTDC);

    rh = HFRecHit(prehit.id(), energy, trise, tfall);
    HFRecHitAuxSetter::setAux(prehit, states, soiPhase_, &rh);

    // Set the "timing from TDC" flag
    const uint32_t flag = triseFromTDC ? 1U : 0U;
    rh.setFlagField(flag, HcalPhase1FlagLabels::TimingFromTDC);
  }

  return rh;
}

rejectingAllFailures()

bool HFSimpleTimeCheck::rejectingAllFailures ( ) const

inline

Definition at line 86 of file HFSimpleTimeCheck.h.

References rejectAllFailures_.

{ return rejectAllFailures_; }

soiPhase()

unsigned HFSimpleTimeCheck::soiPhase ( ) const

inline

Definition at line 80 of file HFSimpleTimeCheck.h.

References soiPhase_.

{ return soiPhase_; }

tfallIfNoTDC()

float HFSimpleTimeCheck::tfallIfNoTDC ( ) const

inline

Definition at line 83 of file HFSimpleTimeCheck.h.

References tfallIfNoTDC_.

{ return tfallIfNoTDC_; }

timeShift()

float HFSimpleTimeCheck::timeShift ( ) const

inline

Definition at line 81 of file HFSimpleTimeCheck.h.

References timeShift_.

Referenced by HFFlexibleTimeCheck::determineAnodeStatus().

{ return timeShift_; }

triseIfNoTDC()

float HFSimpleTimeCheck::triseIfNoTDC ( ) const

inline

Definition at line 82 of file HFSimpleTimeCheck.h.

References triseIfNoTDC_.

{ return triseIfNoTDC_; }

Member Data Documentation

alwaysQAsym_

bool HFSimpleTimeCheck::alwaysQAsym_

private

Definition at line 105 of file HFSimpleTimeCheck.h.

Referenced by alwaysCalculatingQAsym().

energyWeights_

float HFSimpleTimeCheck::energyWeights_[2 *HFAnodeStatus::N_POSSIBLE_STATES - 1][2]

private

Definition at line 97 of file HFSimpleTimeCheck.h.

Referenced by HFSimpleTimeCheck(), and reconstruct().

minChargeForOvershoot_

float HFSimpleTimeCheck::minChargeForOvershoot_

private

Definition at line 103 of file HFSimpleTimeCheck.h.

Referenced by minChargeForOvershoot().

minChargeForUndershoot_

float HFSimpleTimeCheck::minChargeForUndershoot_

private

Definition at line 102 of file HFSimpleTimeCheck.h.

Referenced by minChargeForUndershoot().

rejectAllFailures_

bool HFSimpleTimeCheck::rejectAllFailures_

private

Definition at line 104 of file HFSimpleTimeCheck.h.

Referenced by mapStatusIntoIndex(), and rejectingAllFailures().

soiPhase_

unsigned HFSimpleTimeCheck::soiPhase_

private

Definition at line 98 of file HFSimpleTimeCheck.h.

Referenced by reconstruct(), and soiPhase().

tfallIfNoTDC_

float HFSimpleTimeCheck::tfallIfNoTDC_

private

Definition at line 101 of file HFSimpleTimeCheck.h.

Referenced by reconstruct(), and tfallIfNoTDC().

timeShift_

float HFSimpleTimeCheck::timeShift_

private

Definition at line 99 of file HFSimpleTimeCheck.h.

Referenced by determineAnodeStatus(), reconstruct(), and timeShift().

tlimits_

std::pair<float, float> HFSimpleTimeCheck::tlimits_[2]

private

Definition at line 96 of file HFSimpleTimeCheck.h.

Referenced by determineAnodeStatus(), and HFSimpleTimeCheck().

triseIfNoTDC_

float HFSimpleTimeCheck::triseIfNoTDC_

private

Definition at line 100 of file HFSimpleTimeCheck.h.

Referenced by reconstruct(), and triseIfNoTDC().