HBHENegativeEFilter Class Reference

#include <HBHENegativeEFilter.h>

Public Member Functions
bool	checkPassFilter (const HcalDetId &id, const double *ts, unsigned lenTS) const
const PiecewiseScalingPolynomial &	getA1 (const HcalDetId &id) const
const PiecewiseScalingPolynomial &	getA2 (const HcalDetId &id) const
const std::vector< std::pair< double, double > > &	getCut () const
const std::vector< uint32_t > &	getEtaLimits () const
unsigned	getFirstTimeSlice () const
unsigned	getLastTimeSlice () const
double	getMinCharge () const
	HBHENegativeEFilter ()
	HBHENegativeEFilter (const std::vector< PiecewiseScalingPolynomial > &a1vec, const std::vector< PiecewiseScalingPolynomial > &a2vec, const std::vector< uint32_t > &iEtaLimits, const std::vector< std::pair< double, double > > &cut, double minCharge, unsigned firstTimeSlice, unsigned lastTimeSlice)
bool	isEnabled () const
bool	operator!= (const HBHENegativeEFilter &r) const
bool	operator== (const HBHENegativeEFilter &r) const

Private Member Functions
unsigned	getEtaIndex (const HcalDetId &id) const
template<class Archive >
void	load (Archive &ar, const unsigned)
template<class Archive >
void	save (Archive &ar, const unsigned) const
bool	validate () const

Private Attributes
std::vector< PiecewiseScalingPolynomial >	a1v_
std::vector< PiecewiseScalingPolynomial >	a2v_
std::vector< std::pair< double, double > >	cut_
std::vector< uint32_t >	iEtaLimits_
double	minCharge_
uint32_t	tFirst_
uint32_t	tLast_

Friends
class	boost::serialization::access

Detailed Description

Definition at line 16 of file HBHENegativeEFilter.h.

Constructor & Destructor Documentation

HBHENegativeEFilter::HBHENegativeEFilter ( )

inline

Definition at line 19 of file HBHENegativeEFilter.h.

References checkPassFilter(), and TkAlMuonSelectors_cfi::cut.

: minCharge_(0.), tFirst_(0), tLast_(0) {}

HBHENegativeEFilter::HBHENegativeEFilter	(	const std::vector< PiecewiseScalingPolynomial > &	a1vec,
		const std::vector< PiecewiseScalingPolynomial > &	a2vec,
		const std::vector< uint32_t > &	iEtaLimits,
		const std::vector< std::pair< double, double > > &	cut,
		double	minCharge,
		unsigned	firstTimeSlice,
		unsigned	lastTimeSlice
	)

Definition at line 6 of file HBHENegativeEFilter.cc.

References Exception, and validate().

    : a1v_(a1vec),
      a2v_(a2vec),
      iEtaLimits_(iEtaLimits),
      cut_(cut),
      minCharge_(minCharge),
      tFirst_(firstTimeSlice),
      tLast_(lastTimeSlice)
{
    if (!validate()) throw cms::Exception(
        "Invalid HBHENegativeEFilter constructor arguments");
}

Member Function Documentation

bool HBHENegativeEFilter::checkPassFilter	(	const HcalDetId &	id,
		const double *	ts,
		unsigned	lenTS
	)		const

Definition at line 87 of file HBHENegativeEFilter.cc.

References a1v_, a2v_, TkAlMuonSelectors_cfi::cut, cut_, plotBeamSpotDB::first, getEtaIndex(), mps_fire::i, minCharge_, slope, tFirst_, tLast_, and mitigatedMETSequence_cff::U.

Referenced by HBHENegativeEFilter(), HBHEPhase1Reconstructor::runHBHENegativeEFilter(), and HBHENegativeFlagSetter::setPulseShapeFlags().

{
    bool passes = true;
    const unsigned sz = cut_.size();
    if (sz)
    {
        double chargeInWindow = 0.0;
        for (unsigned i=tFirst_; i<=tLast_ && i<lenTS; ++i)
            chargeInWindow += ts[i];
        if (chargeInWindow >= minCharge_)
        {
            // Figure out the cut value for this charge
            const std::pair<double,double>* cut = &cut_[0];
            double cutValue = cut[0].second;
            if (sz > 1U)
            {
                // First point larger than charge
                unsigned largerPoint = 0;
                for (; cut[largerPoint].first <= chargeInWindow; ++largerPoint) {}

                // Constant extrapolation beyond min and max coords
                if (largerPoint >= sz)
                    cutValue = cut[sz - 1U].second;
                else if (largerPoint)
                {
                    const double slope = (cut[largerPoint].second - cut[largerPoint-1U].second)/
                                         (cut[largerPoint].first - cut[largerPoint-1U].first);
                    cutValue = cut[largerPoint-1U].second + slope*
                        (chargeInWindow - cut[largerPoint-1U].first);
                }
            }

            // Compare the modified time slices with the cut
            const unsigned itaIdx = getEtaIndex(id);
            const PiecewiseScalingPolynomial& a1(a1v_[itaIdx]);
            const PiecewiseScalingPolynomial& a2(a2v_[itaIdx]);

            for (unsigned i=tFirst_; i<=tLast_ && i<lenTS && passes; ++i)
            {
                const double ecorr = ts[i] - a1(ts[i-1U]) - a2(ts[i-2U]);
                passes = ecorr >= cutValue;
            }
        }
    }
    return passes;
}

const PiecewiseScalingPolynomial& HBHENegativeEFilter::getA1 ( const HcalDetId &  id ) const

inline

Definition at line 34 of file HBHENegativeEFilter.h.

References a1v_, and getEtaIndex().

        {return a1v_.at(getEtaIndex(id));}

const PiecewiseScalingPolynomial& HBHENegativeEFilter::getA2 ( const HcalDetId &  id ) const

inline

Definition at line 36 of file HBHENegativeEFilter.h.

References a2v_, and getEtaIndex().

        {return a2v_.at(getEtaIndex(id));}

const std::vector<std::pair<double,double> >& HBHENegativeEFilter::getCut ( ) const

inline

Definition at line 40 of file HBHENegativeEFilter.h.

References cut_.

        {return cut_;}

unsigned HBHENegativeEFilter::getEtaIndex ( const HcalDetId &  id ) const

private

Definition at line 72 of file HBHENegativeEFilter.cc.

References funct::abs(), iEtaLimits_, mitigatedMETSequence_cff::U, and eostools::which().

Referenced by checkPassFilter(), getA1(), getA2(), and operator!=().

{
    const unsigned nLimits = iEtaLimits_.size();
    unsigned which(0U);
    if (nLimits)
    {
        const uint32_t uEta = std::abs(id.ieta());
        const uint32_t* limits(&iEtaLimits_[0]);
        for (; which<nLimits; ++which)
            if (uEta < limits[which])
                break;
    }
    return which;
}

const std::vector<uint32_t>& HBHENegativeEFilter::getEtaLimits ( ) const

inline

Definition at line 38 of file HBHENegativeEFilter.h.

References iEtaLimits_.

        {return iEtaLimits_;}

unsigned HBHENegativeEFilter::getFirstTimeSlice ( ) const

inline

Definition at line 43 of file HBHENegativeEFilter.h.

References tFirst_.

{return tFirst_;}

unsigned HBHENegativeEFilter::getLastTimeSlice ( ) const

inline

Definition at line 44 of file HBHENegativeEFilter.h.

References tLast_.

{return tLast_;}

double HBHENegativeEFilter::getMinCharge ( ) const

inline

Definition at line 42 of file HBHENegativeEFilter.h.

References minCharge_.

{return minCharge_;}

bool HBHENegativeEFilter::isEnabled ( ) const

inline

Definition at line 45 of file HBHENegativeEFilter.h.

References cut_, operator==(), and alignCSCRings::r.

{return !cut_.empty();}

template<class Archive >

void HBHENegativeEFilter::load ( Archive &  ar, const unsigned    )

inlineprivate

Definition at line 75 of file HBHENegativeEFilter.h.

References Exception, tLast_, and validate().

Referenced by MatrixToProcess.MatrixToProcess::getProcess(), MatrixToProcess.MatrixToProcess::listAll(), and ConfigBuilder.ConfigBuilder::prepare_FILTER().

    {
        ar & a1v_ & a2v_ & iEtaLimits_ & cut_ & minCharge_ & tFirst_ & tLast_;
        if (!validate()) throw cms::Exception(
            "In HBHENegativeEFilter::load: invalid data");
    }

bool HBHENegativeEFilter::operator!= ( const HBHENegativeEFilter &  r ) const

inline

Definition at line 49 of file HBHENegativeEFilter.h.

References getEtaIndex(), alignCSCRings::r, and validate().

        {return !(*this == r);}

bool HBHENegativeEFilter::operator==

(

const HBHENegativeEFilter &

r

)

const

Definition at line 58 of file HBHENegativeEFilter.cc.

References a1v_, a2v_, cut_, iEtaLimits_, minCharge_, tFirst_, and tLast_.

Referenced by isEnabled().

{
    if (cut_.empty() && r.cut_.empty())
        return true;
    else
        return a1v_ == r.a1v_ &&
               a2v_ == r.a2v_ &&
               iEtaLimits_ == r.iEtaLimits_ &&
               cut_ == r.cut_ &&
               minCharge_ == r.minCharge_ &&
               tFirst_ == r.tFirst_ &&
               tLast_ == r.tLast_;
}

template<class Archive >

void HBHENegativeEFilter::save ( Archive &  ar, const unsigned    ) const

inlineprivate

Definition at line 67 of file HBHENegativeEFilter.h.

References Exception, tLast_, and validate().

Referenced by Vispa.Main.TabController.TabController::allowClose(), Vispa.Main.TabController.TabController::checkModificationTimestamp(), and SpecificationBuilder_cfi.Specification::saveAll().

    {
        if (!validate()) throw cms::Exception(
            "In HBHENegativeEFilter::save: invalid data");
        ar & a1v_ & a2v_ & iEtaLimits_ & cut_ & minCharge_ & tFirst_ & tLast_;
    }

bool HBHENegativeEFilter::validate ( ) const

private

Definition at line 26 of file HBHENegativeEFilter.cc.

References a1v_, a2v_, cut_, plotBeamSpotDB::first, mps_fire::i, iEtaLimits_, tFirst_, tLast_, and mitigatedMETSequence_cff::U.

Referenced by HBHENegativeEFilter(), load(), operator!=(), and save().

{
    if (cut_.empty())
        return true;

    const std::size_t nLimits(iEtaLimits_.size());
    if (nLimits >= static_cast<std::size_t>(UINT_MAX - 1U))
        return false;
    for (std::size_t i=1; i<nLimits; ++i)
        if (!(iEtaLimits_[i-1] < iEtaLimits_[i]))
            return false;

    if (a1v_.size() != nLimits + 1)
        return false;
    if (a2v_.size() != nLimits + 1)
        return false;

    const std::size_t sz = cut_.size();
    if (sz >= static_cast<std::size_t>(UINT_MAX - 1U))
        return false;
    for (std::size_t i=1; i<sz; ++i)
        if (!(cut_[i-1U].first < cut_[i].first))
            return false;

    if (tFirst_ < 2U)
        return false;
    if (!(tFirst_ <= tLast_))
        return false;

    return true;
}

Friends And Related Function Documentation

friend class boost::serialization::access

friend

Definition at line 64 of file HBHENegativeEFilter.h.

Member Data Documentation

std::vector<PiecewiseScalingPolynomial> HBHENegativeEFilter::a1v_

private

Definition at line 56 of file HBHENegativeEFilter.h.

Referenced by checkPassFilter(), getA1(), operator==(), and validate().

std::vector<PiecewiseScalingPolynomial> HBHENegativeEFilter::a2v_

private

Definition at line 57 of file HBHENegativeEFilter.h.

Referenced by checkPassFilter(), getA2(), operator==(), and validate().

std::vector<std::pair<double,double> > HBHENegativeEFilter::cut_

private

Definition at line 59 of file HBHENegativeEFilter.h.

Referenced by checkPassFilter(), getCut(), isEnabled(), operator==(), and validate().

std::vector<uint32_t> HBHENegativeEFilter::iEtaLimits_

private

Definition at line 58 of file HBHENegativeEFilter.h.

Referenced by getEtaIndex(), getEtaLimits(), operator==(), and validate().

double HBHENegativeEFilter::minCharge_

private

Definition at line 60 of file HBHENegativeEFilter.h.

Referenced by checkPassFilter(), getMinCharge(), and operator==().

uint32_t HBHENegativeEFilter::tFirst_

private

Definition at line 61 of file HBHENegativeEFilter.h.

Referenced by checkPassFilter(), getFirstTimeSlice(), operator==(), and validate().

uint32_t HBHENegativeEFilter::tLast_

private

Definition at line 62 of file HBHENegativeEFilter.h.

Referenced by checkPassFilter(), getLastTimeSlice(), load(), operator==(), save(), and validate().