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 78 of file HBHELinearMap.h.

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

private

Definition at line 77 of file HBHELinearMap.h.

Member Enumeration Documentation

anonymous enum

Enumerator
ChannelCount

Definition at line 18 of file HBHELinearMap.h.

{ ChannelCount = 5184U };

Constructor & Destructor Documentation

HBHELinearMap::HBHELinearMap

(

)

Definition at line 44 of file HBHELinearMap.cc.

References cms::cuda::assert(), ChannelCount, HLT_FULL_cff::depth, mps_fire::i, inverse_, cmsLHEtoEOSManager::l, and lookup_.

                             {
  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 21 of file HBHELinearMap.cc.

References ChannelCount, first, gpuClustering::id, inverse_, and cuda_std::lower_bound().

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 113 of file HBHELinearMap.cc.

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

                                                                                  {
  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 30 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 35 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 83 of file HBHELinearMap.h.

Referenced by find(), and HBHELinearMap().

HBHEChannelId HBHELinearMap::lookup_[ChannelCount]

private

Definition at line 82 of file HBHELinearMap.h.

Referenced by getChannelTriple(), and HBHELinearMap().