ZdcNumberingScheme Class Reference

#include <ZdcNumberingScheme.h>

Public Member Functions
int	detectorLevel (const G4Step *)
void	detectorLevel (const G4Step *, int &, std::vector< int > &, std::vector< G4String > &)
unsigned int	getUnitID (const G4Step *aStep)
void	setVerbosity (const int)
	ZdcNumberingScheme (int)
	~ZdcNumberingScheme ()=default

Static Public Member Functions
static unsigned int	packZdcIndex (int subDet, int layer, int fiber, int channel, int z)
static void	unpackZdcIndex (const unsigned int &idx, int &subDet, int &layer, int &fiber, int &channel, int &z)

Private Attributes
int	verbosity

Detailed Description

Definition at line 12 of file ZdcNumberingScheme.h.

Constructor & Destructor Documentation

ZdcNumberingScheme()

ZdcNumberingScheme::ZdcNumberingScheme

(

int

iv

)

Definition at line 15 of file ZdcNumberingScheme.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv, and verbosity.

                                             {
  verbosity = iv;
  if (verbosity > 0)
    edm::LogVerbatim("ForwardSim") << "Creating ZDCNumberingScheme";
}

~ZdcNumberingScheme()

ZdcNumberingScheme::~ZdcNumberingScheme ( )

default

Member Function Documentation

detectorLevel() [1/2]

int ZdcNumberingScheme::detectorLevel

(

const G4Step *

aStep

)

Definition at line 144 of file ZdcNumberingScheme.cc.

References personalPlayback::level.

Referenced by getUnitID().

                                                         {
  //Find number of levels
  const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
  int level = 0;
  if (touch)
    level = ((touch->GetHistoryDepth()) + 1);
  return level;
}

detectorLevel() [2/2]

void ZdcNumberingScheme::detectorLevel	(	const G4Step *	aStep,
		int &	level,
		std::vector< int > &	copyno,
		std::vector< G4String > &	name
	)

Definition at line 153 of file ZdcNumberingScheme.cc.

References mps_fire::i, cuy::ii, personalPlayback::level, and Skims_PA_cff::name.

                                                                  {
  //Get name and copy numbers
  if (level > 0) {
    const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
    for (int ii = 0; ii < level; ii++) {
      int i = level - ii - 1;
      name.emplace_back(touch->GetVolume(i)->GetName());
      copyno.emplace_back(touch->GetReplicaNumber(i));
    }
  }
}

getUnitID()

unsigned int ZdcNumberingScheme::getUnitID

(

const G4Step *

aStep

)

Definition at line 23 of file ZdcNumberingScheme.cc.

References detectorLevel(), HcalZDCDetId::EM, FCDTask_cfi::fiber, HcalZDCDetId::HAD, nano_mu_digi_cff::layer, personalPlayback::level, HcalZDCDetId::LUM, Skims_PA_cff::name, packZdcIndex(), DetId::rawId(), hgcalPlots::section, HcalZDCDetId::Unknown, and ecaldqm::zside().

                                                              {
  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 = (zside != -1);

    HcalZDCDetId zdcId(section, true_for_positive_eta, channel);
    index = zdcId.rawId();
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("ForwardSim") << "DetectorId:\n" << 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: " << index;

#endif
  }

  return index;
}

packZdcIndex()

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

static

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

References FCDTask_cfi::fiber, heavyIonCSV_trainingSettings::idx, nano_mu_digi_cff::layer, hgcalPlots::section, unpackZdcIndex(), and z.

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: " << idx;
  int newsubdet, newlayer, newfiber, newchannel, newz;
  unpackZdcIndex(idx, newsubdet, newlayer, newfiber, newchannel, newz);
#endif

  return idx;
}

setVerbosity()

void ZdcNumberingScheme::setVerbosity

(

const int

iv

)

Definition at line 21 of file ZdcNumberingScheme.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv, and verbosity.

{ verbosity = iv; }

unpackZdcIndex()

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

static

Definition at line 130 of file ZdcNumberingScheme.cc.

References FCDTask_cfi::fiber, heavyIonCSV_trainingSettings::idx, nano_mu_digi_cff::layer, hgcalPlots::section, and z.

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
}

Member Data Documentation

verbosity

int ZdcNumberingScheme::verbosity

private

Definition at line 34 of file ZdcNumberingScheme.h.

Referenced by setVerbosity(), and ZdcNumberingScheme().