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 26 of file HcalTBNumberingScheme.h.

Constructor & Destructor Documentation

HcalTBNumberingScheme::HcalTBNumberingScheme ( )

inline

Definition at line 29 of file HcalTBNumberingScheme.h.

{}

virtual HcalTBNumberingScheme::~HcalTBNumberingScheme ( )

inlinevirtual

Definition at line 30 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 HcalBarrel, LogDebug, HcalTestNumbering::packHcalIndex(), and HcalTestNumbering::unpackHcalIndex().

Referenced by HcalTB04Analysis::fillBuffer().

                                      {

  int   subdet, zside, group, ieta, iphi, lay;
  HcalTestNumbering::unpackHcalIndex(idHit, subdet, zside, group,
                     ieta, iphi, lay);
  LogDebug("HcalTBSim") << "HcalTBNumberingScheme: i/p ID 0x" << std::hex
            << idHit << std::dec << " det " << zside << " group "
            << group << " layer " << lay << " eta " << ieta 
            << " phi " << iphi;

  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);
  LogDebug("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;
  
  return idunit;
}

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

static

Definition at line 76 of file HcalTBNumberingScheme.cc.

References ExpressReco_HICollisions_FallBack::id, HcalTestNumbering::packHcalIndex(), and 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;
}