HFFlexibleTimeCheck Class Reference

#include <HFFlexibleTimeCheck.h>

Inheritance diagram for HFFlexibleTimeCheck:

Public Member Functions
virtual bool	configure (const AbsHcalAlgoData *config) override
virtual bool	isConfigurable () const override
virtual HFRecHit	reconstruct (const HFPreRecHit &prehit, const HcalCalibrations &calibs, const bool flaggedBadInDB[2], bool expectSingleAnodePMT) override
virtual	~HFFlexibleTimeCheck ()
Public Member Functions inherited from 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, bool rejectAllFailures=true)
bool	rejectingAllFailures () const
unsigned	soiPhase () const
float	tfallIfNoTDC () const
float	timeShift () const
float	triseIfNoTDC () const
virtual	~HFSimpleTimeCheck ()
Public Member Functions inherited from AbsHFPhase1Algo
virtual	~AbsHFPhase1Algo ()

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

Private Attributes
const HFPhase1PMTParams *	algoConf_ = 0
const HFPhase1PMTData *	pmtInfo_ = 0

Detailed Description

Definition at line 7 of file HFFlexibleTimeCheck.h.

Constructor & Destructor Documentation

virtual HFFlexibleTimeCheck::~HFFlexibleTimeCheck ( )

inlinevirtual

Definition at line 12 of file HFFlexibleTimeCheck.h.

{}

Member Function Documentation

virtual bool HFFlexibleTimeCheck::configure ( const AbsHcalAlgoData *  config )

inlineoverridevirtual

Reimplemented from AbsHFPhase1Algo.

Definition at line 16 of file HFFlexibleTimeCheck.h.

References algoConf_, mps_alisetup::config, determineAnodeStatus(), and reconstruct().

    {
        algoConf_ = dynamic_cast<const HFPhase1PMTParams*>(config);
        return algoConf_;
    }

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

overrideprotectedvirtual

Reimplemented from HFSimpleTimeCheck.

Definition at line 10 of file HFFlexibleTimeCheck.cc.

References HFQIE10Info::charge(), ALCARECOTkAlJpsiMuMu_cff::charge, HFPhase1PMTData::cut(), HFAnodeStatus::FAILED_TIMING, HFAnodeStatus::HARDWARE_ERROR, HFQIE10Info::isDataframeOK(), HcalSpecialTimes::isSpecial(), HFPhase1PMTData::minCharge0(), HFPhase1PMTData::minCharge1(), HFAnodeStatus::OK, pmtInfo_, HFPhase1PMTData::T_0_MAX, HFPhase1PMTData::T_0_MIN, HFPhase1PMTData::T_1_MAX, HFPhase1PMTData::T_1_MIN, HFQIE10Info::timeRising(), and HFSimpleTimeCheck::timeShift().

Referenced by configure().

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

    // Require minimum charge for reliable timing measurement
    const float charge = anode.charge();
    const float minCharge = ianode ? pmtInfo_->minCharge1() :
                                     pmtInfo_->minCharge0();
    if (charge < minCharge)
    {
        *isTimingReliable = false;
        return HFAnodeStatus::OK;
    }

    // Check if the rise time information is meaningful
    float trise = anode.timeRising();
    if (HcalSpecialTimes::isSpecial(trise))
        return HFAnodeStatus::FAILED_TIMING;

    // Figure out the timing cuts for this PMT
    const AbsHcalFunctor& minTimeShape = pmtInfo_->cut(
        ianode ? HFPhase1PMTData::T_1_MIN : HFPhase1PMTData::T_0_MIN);
    const AbsHcalFunctor& maxTimeShape = pmtInfo_->cut(
        ianode ? HFPhase1PMTData::T_1_MAX : HFPhase1PMTData::T_0_MAX);

    // Apply the timing cuts
    trise += timeShift();
    if (minTimeShape(charge) <= trise && trise <= maxTimeShape(charge))
        return HFAnodeStatus::OK;
    else
        return HFAnodeStatus::FAILED_TIMING;
}

virtual bool HFFlexibleTimeCheck::isConfigurable ( ) const

inlineoverridevirtual

Reimplemented from HFSimpleTimeCheck.

Definition at line 15 of file HFFlexibleTimeCheck.h.

{return true;}

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

overridevirtual

Reimplemented from HFSimpleTimeCheck.

Definition at line 47 of file HFFlexibleTimeCheck.cc.

References algoConf_, HFPhase1PMTData::ASYMM_MAX, HFPhase1PMTData::ASYMM_MIN, HcalItemCollById< Item >::at(), HFPreRecHit::charge(), HFPreRecHit::chargeAsymmetry(), HFPhase1PMTData::cut(), Exception, HcalPhase1FlagLabels::HFSignalAsymmetry, HFRecHit::id(), HFPreRecHit::id(), HFPhase1PMTData::minChargeAsymm(), pmtInfo_, lumiQueryAPI::q, DetId::rawId(), HFSimpleTimeCheck::reconstruct(), CaloRecHit::setFlagField(), and mitigatedMETSequence_cff::U.

Referenced by configure().

{
    // The algorithm must be configured by now
    if (!algoConf_)
        throw cms::Exception("HFPhase1BadConfig")
            << "In HFFlexibleTimeCheck::reconstruct: algorithm is not configured";

    // Fetch the algorithm configuration data for this PMT
    pmtInfo_ = &algoConf_->at(prehit.id());

    // Run the reconstruction algorithm from the base class
    HFRecHit rh = HFSimpleTimeCheck::reconstruct(
        prehit, calibs, flaggedBadInDB, expectSingleAnodePMT);

    if (rh.id().rawId())
    {
        // Check the charge asymmetry between the two anodes
        bool passesAsymmetryCut = true;
        const std::pair<float,bool> qAsymm =
            prehit.chargeAsymmetry(pmtInfo_->minChargeAsymm());
        if (qAsymm.second)
        {
            const float q = prehit.charge();
            const float minAsymm = (pmtInfo_->cut(HFPhase1PMTData::ASYMM_MIN))(q);
            const float maxAsymm = (pmtInfo_->cut(HFPhase1PMTData::ASYMM_MAX))(q);
            passesAsymmetryCut = minAsymm <= qAsymm.first && qAsymm.first <= maxAsymm;
        }
        if (!passesAsymmetryCut)
            rh.setFlagField(1U, HcalPhase1FlagLabels::HFSignalAsymmetry);
    }

    return rh;
}

Member Data Documentation

const HFPhase1PMTParams* HFFlexibleTimeCheck::algoConf_ = 0

private

Definition at line 32 of file HFFlexibleTimeCheck.h.

Referenced by configure(), and reconstruct().

const HFPhase1PMTData* HFFlexibleTimeCheck::pmtInfo_ = 0

private

Definition at line 33 of file HFFlexibleTimeCheck.h.

Referenced by determineAnodeStatus(), and reconstruct().