ZdcNumberingScheme Namespace Reference

Functions
unsigned int	getUnitID (const G4Step *aStep)
unsigned int	packZdcIndex (int subDet, int layer, int fiber, int channel, int z)
void	unpackZdcIndex (const unsigned int &idx, int &subDet, int &layer, int &fiber, int &channel, int &z)

Function Documentation

getUnitID()

unsigned int ZdcNumberingScheme::getUnitID

(

const G4Step *

aStep

)

Definition at line 42 of file ZdcNumberingScheme.cc.

References TauDecayModes::dec, HcalZDCDetId::EM, HcalZDCDetId::HAD, personalPlayback::level, HcalZDCDetId::LUM, Skims_PA_cff::name, packZdcIndex(), DetId::rawId(), hgcalPlots::section, HcalZDCDetId::Unknown, and ecaldqm::zside().

Referenced by ZdcSD::setDetUnitId().

                                              {
    uint32_t index = 0;
    int level = detectorLevel(aStep);

    if (level > 0) {
      std::vector<int> copyno;
      std::vector<G4String> name;

      detectorLevel(aStep, level, copyno, name);

      int zside = 0;
      int channel = 0;
      int fiber = 0;
      int layer = 0;
      HcalZDCDetId::Section section = HcalZDCDetId::Unknown;

      for (int ich = 0; ich < level; ich++) {
        if (name[ich] == "ZDC") {
          if (copyno[ich] == 1)
            zside = 1;
          if (copyno[ich] == 2)
            zside = -1;
        } else if (name[ich] == "ZDC_EMLayer") {
          section = HcalZDCDetId::EM;
#ifdef EDM_ML_DEBUG
          layer = copyno[ich];
#endif
        } else if (name[ich] == "ZDC_EMFiber") {
          fiber = copyno[ich];
          if (fiber < 20)
            channel = 1;
          else if (fiber < 39)
            channel = 2;
          else if (fiber < 58)
            channel = 3;
          else if (fiber < 77)
            channel = 4;
          else
            channel = 5;
        } else if (name[ich] == "ZDC_LumLayer") {
          section = HcalZDCDetId::LUM;
          layer = copyno[ich];
          channel = layer;
        } else if (name[ich] == "ZDC_HadLayer") {
          section = HcalZDCDetId::HAD;
          layer = copyno[ich];
          if (layer < 6)
            channel = 1;
          else if (layer < 12)
            channel = 2;
          else if (layer < 18)
            channel = 3;
          else
            channel = 4;
        }
#ifdef EDM_ML_DEBUG
        else if (name[ich] == "ZDC_LumGas") {
          fiber = 1;
        } else if (name[ich] == "ZDC_HadFiber") {
          fiber = copyno[ich];
        }
#endif
      }

#ifdef EDM_ML_DEBUG
      unsigned intindex = 0;
      intindex = packZdcIndex(section, layer, fiber, channel, zside);
#endif

      bool true_for_positive_eta = true;
      if (zside == -1)
        true_for_positive_eta = false;

      HcalZDCDetId zdcId(section, true_for_positive_eta, channel);
      index = zdcId.rawId();

#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("ForwardSim") << "DetectorId: " << zdcId;

      edm::LogVerbatim("ForwardSim") << "ZdcNumberingScheme:"
                                     << "  getUnitID - # of levels = " << level;
      for (int ich = 0; ich < level; ich++)
        edm::LogVerbatim("ForwardSim") << "  " << ich << ": copyno " << copyno[ich] << " name=" << name[ich]
                                       << "  section " << section << " zside " << zside << " layer " << layer
                                       << " fiber " << fiber << " channel " << channel << "packedIndex =" << intindex
                                       << " detId raw: " << std::hex << index << std::dec;

#endif
    }

    return index;
  }

packZdcIndex()

unsigned ZdcNumberingScheme::packZdcIndex	(	int	subDet,
		int	layer,
		int	fiber,
		int	channel,
		int	z
	)

pack the Unit ID for Zdc
z = 1,2 = -z,+z; subDet = 1,2,3 = EM,Lum,HAD; fiber = 1-96 (EM,HAD), 1 (Lum); channel = 1-5 (EM), layer# (Lum), 1-3 (HAD)

Definition at line 135 of file ZdcNumberingScheme.cc.

References TauDecayModes::dec, heavyIonCSV_trainingSettings::idx, hgcalPlots::section, and unpackZdcIndex().

Referenced by getUnitID().

                                                                               {
    unsigned int idx = ((z - 1) & 1) << 20;  //bit 20
    idx += (channel & 7) << 17;              //bits 17-19
    idx += (fiber & 255) << 9;               //bits 9-16
    idx += (layer & 127) << 2;               //bits 2-8
    idx += (section & 3);                    //bits 0-1

#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("ForwardSim") << "ZDC packing: section " << section << " layer  " << layer << " fiber " << fiber
                                   << " channel " << channel << " zside " << z << "idx: " << std::hex << idx
                                   << std::dec;
    int newsubdet, newlayer, newfiber, newchannel, newz;
    unpackZdcIndex(idx, newsubdet, newlayer, newfiber, newchannel, newz);
#endif

    return idx;
  }

unpackZdcIndex()

void ZdcNumberingScheme::unpackZdcIndex	(	const unsigned int &	idx,
		int &	subDet,
		int &	layer,
		int &	fiber,
		int &	channel,
		int &	z
	)

Definition at line 153 of file ZdcNumberingScheme.cc.

References heavyIonCSV_trainingSettings::idx, and hgcalPlots::section.

Referenced by packZdcIndex(), and ZdcTestAnalysis::update().

                                                                                                           {
    z = 1 + ((idx >> 20) & 1);
    channel = (idx >> 17) & 7;
    fiber = (idx >> 9) & 255;
    layer = (idx >> 2) & 127;
    section = idx & 3;

#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("ForwardSim") << "ZDC unpacking: idx:" << idx << " -> section " << section << " layer " << layer
                                   << " fiber " << fiber << " channel " << channel << " zside " << z;
#endif
  }