HBHELinearMap Class Reference

#include <HBHELinearMap.h>

Classes
class	HBHEChannelId

Public Types
enum	{ ChannelCount = 5184U }

Public Member Functions
void	getChannelTriple (unsigned index, unsigned *depth, int *ieta, unsigned *iphi) const
	HBHELinearMap ()
bool	isValidTriple (unsigned depth, int ieta, unsigned iphi) const
unsigned	linearIndex (unsigned depth, int ieta, unsigned iphi) const

Static Public Member Functions
static HcalSubdetector	getSubdetector (unsigned depth, int ieta)

Private Types
typedef std::vector< MapPair >	ChannelMap
typedef std::pair< HBHEChannelId, unsigned >	MapPair

Private Member Functions
unsigned	find (unsigned depth, int ieta, unsigned iphi) const

Private Attributes
ChannelMap	inverse_
HBHEChannelId	lookup_ [ChannelCount]

Detailed Description

Definition at line 16 of file HBHELinearMap.h.

Member Typedef Documentation

typedef std::vector<MapPair> HBHELinearMap::ChannelMap

private

Definition at line 80 of file HBHELinearMap.h.

typedef std::pair<HBHEChannelId,unsigned> HBHELinearMap::MapPair

private

Definition at line 79 of file HBHELinearMap.h.

Member Enumeration Documentation

anonymous enum

Enumerator
ChannelCount

Definition at line 19 of file HBHELinearMap.h.

{ChannelCount = 5184U};

Constructor & Destructor Documentation

HBHELinearMap::HBHELinearMap

(

)

Definition at line 51 of file HBHELinearMap.cc.

References ChannelCount, egammaForCoreTracking_cff::depth, mps_fire::i, inverse_, checklumidiff::l, lookup_, and jetUpdater_cfi::sort.

{
    unsigned l = 0;
    unsigned depth = 1;

    for (int ieta = -29; ieta <= -21; ++ieta)
        for (unsigned iphi=1; iphi<72; iphi+=2)
            lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

    for (int ieta = -20; ieta <= 20; ++ieta)
        if (ieta)
            for (unsigned iphi=1; iphi<=72; ++iphi)
                lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

    for (int ieta = 21; ieta <= 29; ++ieta)
        for (unsigned iphi=1; iphi<72; iphi+=2)
            lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

    depth = 2;

    for (int ieta = -29; ieta <= -21; ++ieta)
        for (unsigned iphi=1; iphi<72; iphi+=2)
            lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

    for (int ieta = -20; ieta <= -18; ++ieta)
        for (unsigned iphi=1; iphi<=72; ++iphi)
            lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

    for (int ieta = -16; ieta <= -15; ++ieta)
        for (unsigned iphi=1; iphi<=72; ++iphi)
            lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

    for (int ieta = 15; ieta <= 16; ++ieta)
        for (unsigned iphi=1; iphi<=72; ++iphi)
            lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

    for (int ieta = 18; ieta <= 20; ++ieta)
        for (unsigned iphi=1; iphi<=72; ++iphi)
            lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

    for (int ieta = 21; ieta <= 29; ++ieta)
        for (unsigned iphi=1; iphi<72; iphi+=2)
            lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

    depth = 3;

    for (int ieta = -28; ieta <= -27; ++ieta)
        for (unsigned iphi=1; iphi<72; iphi+=2)
            lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

    for (int ieta = -16; ieta <= -16; ++ieta)
        for (unsigned iphi=1; iphi<=72; ++iphi)
             lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

    for (int ieta = 16; ieta <= 16; ++ieta)
        for (unsigned iphi=1; iphi<=72; ++iphi)
             lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

    for (int ieta = 27; ieta <= 28; ++ieta)
        for (unsigned iphi=1; iphi<72; iphi+=2)
            lookup_[l++] = HBHEChannelId(depth, ieta, iphi);

    assert(l == ChannelCount);

    inverse_.reserve(ChannelCount);
    for (unsigned i=0; i<ChannelCount; ++i)
        inverse_.push_back(MapPair(lookup_[i], i));
    std::sort(inverse_.begin(), inverse_.end());
}

Member Function Documentation

unsigned HBHELinearMap::find ( unsigned  depth, int  ieta, unsigned  iphi  ) const

private

Definition at line 22 of file HBHELinearMap.cc.

References ChannelCount, plotBeamSpotDB::first, triggerObjects_cff::id, inverse_, create_public_lumi_plots::loc, pfDeepBoostedJetPreprocessParams_cfi::lower_bound, and mitigatedMETSequence_cff::U.

Referenced by BeautifulSoup.Tag::__getattr__(), BeautifulSoup.Tag::firstText(), isValidTriple(), and linearIndex().

{
    const HBHEChannelId id(depth, ieta, iphi);
    const unsigned loc = std::lower_bound(
        inverse_.begin(), inverse_.end(), MapPair(id, 0U)) - inverse_.begin();
    if (loc < ChannelCount)
        if (inverse_[loc].first == id)
            return inverse_[loc].second;
    return ChannelCount;
}

void HBHELinearMap::getChannelTriple	(	unsigned	index,
		unsigned *	depth,
		int *	ieta,
		unsigned *	iphi
	)		const

Definition at line 7 of file HBHELinearMap.cc.

References ChannelCount, HBHELinearMap::HBHEChannelId::depth(), Exception, and lookup_.

{
    if (index >= ChannelCount)
        throw cms::Exception("In HBHELinearMap::getChannelTriple: "
                                    "input index out of range");
    const HBHEChannelId& id = lookup_[index];
    if (depth)
        *depth = id.depth();
    if (ieta)
        *ieta = id.ieta();
    if (iphi)
        *iphi = id.iphi();
}

HcalSubdetector HBHELinearMap::getSubdetector ( unsigned  depth, int  ieta  )

static

Definition at line 121 of file HBHELinearMap.cc.

References funct::abs(), Exception, HcalBarrel, HcalEndcap, and mitigatedMETSequence_cff::U.

{
    const unsigned abseta = std::abs(ieta);

    // Make sure the arguments are in range
    if (!(abseta <= 29U))
        throw cms::Exception("In HBHELinearMap::getSubdetector: "
                                    "eta argument out of range");
    if (!(depth > 0U && depth < 4U))
        throw cms::Exception("In HBHELinearMap::getSubdetector: "
                                    "depth argument out of range");
    if (abseta == 29U)
        if (!(depth <= 2U))
            throw cms::Exception("In HBHELinearMap::getSubdetector: "
                                        "depth argument out of range "
                                        "for |ieta| = 29");
    if (abseta <= 15U)
        return HcalBarrel;
    else if (abseta == 16U)
        return depth <= 2U ? HcalBarrel : HcalEndcap;
    else
        return HcalEndcap;
}

bool HBHELinearMap::isValidTriple	(	unsigned	depth,
		int	ieta,
		unsigned	iphi
	)		const

Definition at line 34 of file HBHELinearMap.cc.

References ChannelCount, and find().

{
    const unsigned ind = find(depth, ieta, iphi);
    return ind < ChannelCount;
}

unsigned HBHELinearMap::linearIndex	(	unsigned	depth,
		int	ieta,
		unsigned	iphi
	)		const

Definition at line 41 of file HBHELinearMap.cc.

References ChannelCount, Exception, and find().

Referenced by HcalInterpolatedPulseColl::getChannelPulse().

{
    const unsigned ind = find(depth, ieta, iphi);
    if (ind >= ChannelCount)
        throw cms::Exception("In HBHELinearMap::linearIndex: "
                                    "invalid channel triple");
    return ind;
}

Member Data Documentation

ChannelMap HBHELinearMap::inverse_

private

Definition at line 85 of file HBHELinearMap.h.

Referenced by find(), and HBHELinearMap().

HBHEChannelId HBHELinearMap::lookup_[ChannelCount]

private

Definition at line 84 of file HBHELinearMap.h.

Referenced by getChannelTriple(), and HBHELinearMap().