HcalTBNumberingScheme Class Reference

#include <SimG4CMS/HcalTestBeam/interface/HcalTBNumberingScheme.h>

Public Member Functions
	HcalTBNumberingScheme ()
virtual	~HcalTBNumberingScheme ()

Static Public Member Functions
static uint32_t	getUnitID (const uint32_t id, const int mode)
static std::vector< uint32_t >	getUnitIDs (const int type, const int mode)

Detailed Description

Description: Numbering scheme for hadron calorimeter in test beam

Usage: Sets up tower ID's of all towers in the 2004 Hcal test beam setup

Definition at line 25 of file HcalTBNumberingScheme.h.

Constructor & Destructor Documentation

HcalTBNumberingScheme::HcalTBNumberingScheme ( )

inline

Definition at line 27 of file HcalTBNumberingScheme.h.

{}

virtual HcalTBNumberingScheme::~HcalTBNumberingScheme ( )

inlinevirtual

Definition at line 28 of file HcalTBNumberingScheme.h.

{}

Member Function Documentation

uint32_t HcalTBNumberingScheme::getUnitID ( const uint32_t  id, const int  mode  )

static

Definition at line 27 of file HcalTBNumberingScheme.cc.

References TauDecayModes::dec, watchdog::group, HcalBarrel, HcalTestNumbering::packHcalIndex(), HcalTestNumbering::unpackHcalIndex(), and ecaldqm::zside().

Referenced by HcalTB04Analysis::fillBuffer().

                                                                              {
  int subdet, zside, group, ieta, iphi, lay;
  HcalTestNumbering::unpackHcalIndex(idHit, subdet, zside, group, ieta, iphi, lay);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HcalTBSim") << "HcalTBNumberingScheme: i/p ID 0x" << std::hex << idHit << std::dec << " det "
                                << zside << " group " << group << " layer " << lay << " eta " << ieta << " phi "
                                << iphi;
#endif
  uint32_t idunit;
  if (subdet == static_cast<int>(HcalBarrel)) {
    if (lay <= 17)
      group = 1;
    else
      group = 2;
  }
  if (mode > 0) {
    if (subdet == static_cast<int>(HcalBarrel) && iphi > 4) {
      if (lay <= 17) {
        // HB2 (unmasked and masked)
        if (ieta > 4 && ieta < 10) {
          idunit = HcalTestNumbering::packHcalIndex(0, 0, 1, 0, iphi, lay);
        } else {
          idunit = HcalTestNumbering::packHcalIndex(0, 0, 1, ieta, iphi, group);
        }
      } else {
        // HO behind HB2
        idunit = HcalTestNumbering::packHcalIndex(0, 0, 1, ieta, iphi, 18);
      }
    } else {
      // HB1, HE, HO behind HB1
      idunit = HcalTestNumbering::packHcalIndex(0, 0, 1, ieta, iphi, group);
    }
  } else {
    idunit = HcalTestNumbering::packHcalIndex(0, 0, 1, ieta, iphi, group);
  }

  HcalTestNumbering::unpackHcalIndex(idunit, subdet, zside, group, ieta, iphi, lay);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HcalTBSim") << "HcalTBNumberingScheme: idHit 0x" << std::hex << idHit << " idunit 0x" << idunit
                                << std::dec << "\n"
                                << "HcalTBNumberingScheme: o/p ID 0x" << std::hex << idunit << std::dec << " det "
                                << zside << " group " << group << " layer " << lay << " eta " << ieta << " phi "
                                << iphi;
#endif
  return idunit;
}

std::vector< uint32_t > HcalTBNumberingScheme::getUnitIDs ( const int  type, const int  mode  )

static

Definition at line 74 of file HcalTBNumberingScheme.cc.

References gpuClustering::id, HcalTestNumbering::packHcalIndex(), and createJobs::tmp.

Referenced by HcalTB04Analysis::init().

                                                                                    {
  std::vector<uint32_t> tmp;
  int iphi, ieta, lay;
  uint32_t id;

  if (type != 1) {
    // Include HB and HO id's
    if (mode > 0) {
      // HB1 and masked part of HB2
      for (ieta = 1; ieta < 17; ieta++) {
        for (iphi = 1; iphi < 9; iphi++) {
          id = HcalTestNumbering::packHcalIndex(0, 0, 1, ieta, iphi, 1);
          tmp.push_back(id);
        }
      }
      // HO behind HB1
      for (ieta = 1; ieta < 16; ieta++) {
        for (iphi = 2; iphi < 5; iphi++) {
          id = HcalTestNumbering::packHcalIndex(0, 0, 1, ieta, iphi, 2);
          tmp.push_back(id);
        }
      }
      // HB2
      for (lay = 1; lay < 18; lay++) {
        for (iphi = 5; iphi < 9; iphi++) {
          id = HcalTestNumbering::packHcalIndex(0, 0, 1, 0, iphi, lay);
          tmp.push_back(id);
        }
      }
      // HO behind HB2
      lay = 18;
      for (ieta = 1; ieta < 16; ieta++) {
        for (iphi = 5; iphi < 8; iphi++) {
          id = HcalTestNumbering::packHcalIndex(0, 0, 1, ieta, iphi, lay);
          tmp.push_back(id);
        }
      }
    } else {
      // HB1 & HB2
      for (ieta = 1; ieta < 17; ieta++) {
        for (iphi = 1; iphi < 9; iphi++) {
          id = HcalTestNumbering::packHcalIndex(0, 0, 1, ieta, iphi, 1);
          tmp.push_back(id);
        }
      }
      // HO behind HB
      for (ieta = 1; ieta < 16; ieta++) {
        for (iphi = 2; iphi < 8; iphi++) {
          id = HcalTestNumbering::packHcalIndex(0, 0, 1, ieta, iphi, 2);
          tmp.push_back(id);
        }
      }
    }
  }

  if (type > 0) {
    // Include HE id's
    for (ieta = 15; ieta < 17; ieta++) {
      for (iphi = 3; iphi < 7; iphi++) {
        id = HcalTestNumbering::packHcalIndex(0, 0, 1, ieta, iphi, 3);
        tmp.push_back(id);
      }
    }
    for (iphi = 3; iphi < 7; iphi++) {
      id = HcalTestNumbering::packHcalIndex(0, 0, 1, 17, iphi, 1);
      tmp.push_back(id);
    }
    for (ieta = 18; ieta < 21; ieta++) {
      for (iphi = 3; iphi < 7; iphi++) {
        for (int idep = 1; idep < 3; idep++) {
          id = HcalTestNumbering::packHcalIndex(0, 0, 1, ieta, iphi, idep);
          tmp.push_back(id);
        }
      }
    }
    for (ieta = 21; ieta < 26; ieta++) {
      for (iphi = 2; iphi < 4; iphi++) {
        for (int idep = 1; idep < 3; idep++) {
          id = HcalTestNumbering::packHcalIndex(0, 0, 1, ieta, iphi, idep);
          tmp.push_back(id);
        }
      }
    }
  }

  return tmp;
}