InttDeadMap.cc

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file InttDeadMap.cc

#include "InttDeadMap.h"

#include <g4detectors/PHG4CellDefs.h>

#include <cassert>
#include <iostream>
#include <limits>

using namespace std;

int InttDeadMap::s_wildCardID = -1;

const InttDeadMap::Map&
InttDeadMap::getDeadChannels(void) const
{
  static Map tmp_map;
  return tmp_map;
}

InttDeadMap::Map&
InttDeadMap::getDeadChannels(void)
{
  static Map tmp_map;
  return tmp_map;
}

void InttDeadMap::addDeadChannelIntt(const int layer,
                                     const int ladder_phi, const int ladder_z,
                                     const int strip_z, const int strip_phi)
{
  addDeadChannel(getInttKey(layer,
                            ladder_phi, ladder_z,
                            strip_z, strip_phi));
}

bool InttDeadMap::isDeadChannelIntt(const int layer,
                                    const int ladder_phi, const int ladder_z,
                                    const int strip_z, const int strip_phi) const
{
  if (isDeadChannel(getInttKey(layer,
                               ladder_phi, ladder_z,
                               strip_z, strip_phi)))
    return true;
  else if (isDeadChannel(getInttKey(layer,
                                    ladder_phi, ladder_z,
                                    strip_z, s_wildCardID)))
    return true;
  else if (isDeadChannel(getInttKey(layer,
                                    ladder_phi, ladder_z,
                                    s_wildCardID, s_wildCardID)))
    return true;
  else if (isDeadChannel(getInttKey(layer,
                                    ladder_phi, s_wildCardID,
                                    s_wildCardID, s_wildCardID)))
    return true;
  else if (isDeadChannel(getInttKey(layer,
                                    s_wildCardID, s_wildCardID,
                                    s_wildCardID, s_wildCardID)))
    return true;
  else
    return false;
}

int InttDeadMap::isValid() const
{
  return size() > 0;
}

void InttDeadMap::identify(std::ostream& os) const
{
  os << "InttDeadMap base class" << std::endl;
}

PHG4CellDefs::keytype InttDeadMap::getInttKey(int layer,
                                              int ladder_phi, int ladder_z,
                                              int strip_z, int strip_phi)
{
  static const int layer_bit = 8;
  static const int ladder_phi_bit = 16;
  static const int ladder_z_bit = 8;
  static const int strip_z_bit = 16;
  static const int strip_phi_bit = 16;

  bool wildcard = false;

  if (layer == s_wildCardID) wildcard = true;
  if (wildcard) layer = (1 << layer_bit) - 1;

  if (ladder_phi == s_wildCardID) wildcard = true;
  if (wildcard) ladder_phi = (1 << ladder_phi_bit) - 1;

  if (ladder_z == s_wildCardID) wildcard = true;
  if (wildcard) ladder_z = (1 << ladder_z_bit) - 1;

  if (strip_z == s_wildCardID) wildcard = true;
  if (wildcard) strip_z = (1 << strip_z_bit) - 1;

  if (strip_phi == s_wildCardID) wildcard = true;
  if (wildcard) strip_phi = (1 << strip_phi_bit) - 1;

  //  bit sum check
  assert(layer_bit + ladder_phi_bit + ladder_z_bit + strip_z_bit + strip_phi_bit == numeric_limits<PHG4CellDefs::keytype>::digits);

  //max range check
  assert(layer < (1 << layer_bit));
  assert(ladder_phi < (1 << ladder_phi_bit));
  assert(ladder_z < (1 << ladder_z_bit));
  assert(strip_z < (1 << strip_z_bit));
  assert(strip_phi < (1 << strip_phi_bit));

  //min range check
  assert(layer >= 0);
  assert(ladder_phi >= 0);
  assert(ladder_z >= 0);
  assert(strip_z >= 0);
  assert(strip_phi >= 0);

  PHG4CellDefs::keytype key = 0;

  key += layer;

  key <<= ladder_phi_bit;
  key += ladder_phi;

  key <<= ladder_z_bit;
  key += ladder_z;

  key <<= strip_z_bit;
  key += strip_z;

  key <<= strip_phi_bit;
  key += strip_phi;

  return key;
}