FP420NumberingScheme Class Reference

#include <FP420NumberingScheme.h>

Public Member Functions
virtual int	detectorLevel (const G4Step *) const
virtual void	detectorLevel (const G4Step *, int &, int *, G4String *) const
	FP420NumberingScheme ()
virtual unsigned int	getUnitID (const G4Step *aStep) const
virtual	~FP420NumberingScheme ()

Static Public Member Functions
static unsigned int	packFP420Index (int det, int zside, int station, int superplane)
static unsigned	packMYIndex (int rn0, int pn0, int sn0, int det, int zside, int sector, int zmodule)
static int	realzside (int rn0, int zsideinorder)
static int	unpackCopyIndex (int rn0, int zside)
static void	unpackFP420Index (const unsigned int &idx, int &det, int &zside, int &station, int &superplane)
static int	unpackLayerIndex (int rn0, int zside)
static void	unpackMYIndex (const int &idx, int rn0, int pn0, int sn0, int &det, int &zside, int &sector, int &zmodule)
static int	unpackOrientation (int rn0, int zside)

Detailed Description

Definition at line 15 of file FP420NumberingScheme.h.

Constructor & Destructor Documentation

FP420NumberingScheme()

FP420NumberingScheme::FP420NumberingScheme

(

)

Definition at line 17 of file FP420NumberingScheme.cc.

                                           {
  //  sn0=3, pn0=6, rn0=7;
}

~FP420NumberingScheme()

FP420NumberingScheme::~FP420NumberingScheme ( )

virtual

Definition at line 21 of file FP420NumberingScheme.cc.

                                            {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("FP420") << " Deleting FP420NumberingScheme";
#endif
}

Member Function Documentation

detectorLevel() [1/2]

int FP420NumberingScheme::detectorLevel ( const G4Step *  aStep ) const

virtual

Definition at line 27 of file FP420NumberingScheme.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 FP420NumberingScheme::detectorLevel ( const G4Step *  aStep, int &  level, int *  copyno, G4String *  name  ) const

virtual

Definition at line 36 of file FP420NumberingScheme.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[ii] = touch->GetVolume(i)->GetName();
      copyno[ii] = touch->GetReplicaNumber(i);
    }
  }
}

getUnitID()

unsigned int FP420NumberingScheme::getUnitID ( const G4Step *  aStep ) const

virtual

Definition at line 48 of file FP420NumberingScheme.cc.

References detectorLevel(), personalPlayback::level, Skims_PA_cff::name, packFP420Index(), relativeConstraints::station, and ecaldqm::zside().

Referenced by FP420SD::setDetUnitId().

                                                                      {
  unsigned intindex = 0;
  int level = detectorLevel(aStep);

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("FP420") << "FP420NumberingScheme number of levels= " << level;
#endif
  //  unsigned int intIndex = 0;
  if (level > 0) {
    int* copyno = new int[level];
    G4String* name = new G4String[level];
    detectorLevel(aStep, level, copyno, name);

    //    int det   = static_cast<int>(FP420);;

    int det = 1;
    int stationgen = 0;
    int zside = 0;
    int station = 0;
    int plane = 0;
    for (int ich = 0; ich < level; ich++) {
      /*
      // old set up configuration with equidistant stations 
      if(name[ich] == "FP420Ex") {
      stationgen   = copyno[ich];
      } else if(name[ich] == "SISTATION") {
      station   = stationgen;
      } else if(name[ich] == "SIPLANE") {
      plane   = copyno[ich];
      } else if(name[ich] == "SIDETL") {
      zside   = 1;
      } else if(name[ich] == "SIDETR") {
      zside   = 2;
      }
      */
      // new and old set up configurations are possible:
      if (name[ich] == "FP420E") {
        det = copyno[ich];
      } else if (name[ich] == "HPS240E") {
        det = copyno[ich] + 2;
      } else if (name[ich] == "FP420Ex1" || name[ich] == "HPS240Ex1") {
        stationgen = 1;
        //    } else if(name[ich] == "FP420Ex2") {
        //  stationgen   = 2;
      } else if (name[ich] == "FP420Ex3" || name[ich] == "HPS240Ex3") {
        stationgen = 2;  // was =3
      } else if (name[ich] == "SISTATION" || name[ich] == "HPS240SISTATION") {
        station = stationgen;
      } else if (name[ich] == "SIPLANE" || name[ich] == "HPS240SIPLANE") {
        plane = copyno[ich];
        // SIDETL (or R) can be ether X or Y type in next schemes of readout
        //        !!! (=...) zside
        //
        //      1                  2     <---copyno
        //   Front(=2) Empty(=4) Back(=6)     <--SIDETR OR SENSOR2
        //      1         2              <---copyno
        //   Front(=1) Back(=3) Empty(=5)     <--SIDETL OR SENSOR1
        //
      } else if (name[ich] == "SENSOR2" || name[ich] == "HPS240SENSOR2") {
        //    } else if(name[ich] == "SIDETR") {
        //  zside   = 4 * copyno[ich] - 2 ;//= 2   6 (copyno=1,2)
        zside = 3 * copyno[ich];  //= 3   6 (copyno=1,2)
      } else if (name[ich] == "SENSOR1" || name[ich] == "HPS240SENSOR1") {
        //     } else if(name[ich] == "SIDETL") {
        //  zside   = 2 * copyno[ich] - 1 ;//= 1   3 (copyno=1,2)
        zside = copyno[ich];  //= 1   2 (copyno=1,2)
      }
      //
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("FP420") << "FP420NumberingScheme  "
                                << "ich=" << ich << "copyno" << copyno[ich] << "name=" << name[ich];
#endif
    }
    // det = 1 for +FP420 ,  = 2 for -FP420  / (det-1) = 0,1
    // det = 3 for +HPS240 , = 4 for -HPS240 / (det-1) = 2,3
    // 0 is as default for every below:
    // Z index
    // station number 1 - 8   (in reality just 2 ones)
    // superplane(superlayer) number  1 - 16 (in reality just 5 ones)

    // intindex = myPacker.packEcalIndex (det, zside, station, plane);// see examples
    // intindex = myPacker.packCastorIndex (det, zside, station, plane);// see examples
    intindex = packFP420Index(det, zside, station, plane);
    //
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("FP420") << "FP420NumberingScheme det=" << det << " zside=" << zside << " station=" << station
                              << " plane=" << plane;
    for (int ich = 0; ich < level; ich++) {
      edm::LogVerbatim("FP420") << " name = " << name[ich] << " copy = " << copyno[ich];
      edm::LogVerbatim("FP420") << " packed index = intindex" << intindex;
    }
#endif

    delete[] copyno;
    delete[] name;
  }

  return intindex;
}

packFP420Index()

unsigned FP420NumberingScheme::packFP420Index ( int  det, int  zside, int  station, int  superplane  )

static

Definition at line 148 of file FP420NumberingScheme.cc.

References heavyIonCSV_trainingSettings::idx, relativeConstraints::station, and ecaldqm::zside().

Referenced by getUnitID().

                                                                                        {
  //
  unsigned int idx = ((det - 1) & 3) << 19;  //bit 19-20 (det-1):0- 3 = 4-->2**2 = 4  -> 4-1  ->((det-1)&3)  2 bit: 0-1
  // unsigned int idx = ((det-1)&1)<<20;//bit 20-20 (det-1):0- 1 = 2-->2**1 = 2  -> 2-1  ->((det-1)&1)  1 bit: 0
  idx += (zside & 7) << 7;  //bits  7- 9   zside:0- 7 = 8-->2**3 = 8  -> 8-1  ->  (zside&7)  3 bits:0-2
  //  idx += (zside&3)<<7;             //bits  7- 8   zside:0- 3 = 4-->2**2 = 4  -> 4-1  ->  (zside&3)  2 bits:0-1
  idx += (station & 7) << 4;  //bits  4- 6 station:0- 7 = 8-->2**3 = 8  -> 8-1  ->(station&7)  3 bits:0-2
  idx += (plane & 15);        //bits  0- 3   plane:0-15 =16-->2**4 =16  ->16-1  -> (plane&15)  4 bits:0-3

  //

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("FP420") << "FP420 packing: det " << det << " zside  " << zside << " station " << station
                            << " plane " << plane << " idx " << idx;
#endif
  //  int newdet, newzside, newstation,newplane;
  //  unpackFP420Index(idx, newdet, newzside, newstation,newplane);

  //

  return idx;
}

packMYIndex()

static unsigned FP420NumberingScheme::packMYIndex ( int  rn0, int  pn0, int  sn0, int  det, int  zside, int  sector, int  zmodule  )

inlinestatic

Definition at line 35 of file FP420NumberingScheme.h.

References hgcalTBTopologyTester_cfi::sector, and ecaldqm::zside().

Referenced by FP420Test::update().

                                                                                                      {
    int zScale = (rn0 - 1);           // rn0=3 - current --> for update  rn0=7
    int sScale = zScale * (pn0 - 1);  //pn0=6
    int dScale = sScale * (sn0 - 1);  //sn0=3

    unsigned int intindex = dScale * (det - 1) + sScale * (sector - 1) + zScale * (zmodule - 1) + zside;
    //
    return intindex;
  }

realzside()

static int FP420NumberingScheme::realzside ( int  rn0, int  zsideinorder  )

inlinestatic

Definition at line 97 of file FP420NumberingScheme.h.

References ecaldqm::zside().

Referenced by FP420Test::update().

                                                  {
    // zsideinorder:1 2 3 4 5 6
    //sensorsold    1 0 0 2 0 0
    //sensorsnew    1 0 5 2 4 0 ???
    //sensorsnew    1 3 0 2 0 6
    //zside         1 3 5 2 4 6  over layers 1 and 2
    int zside, zsidereal;
    if (zsideinorder < 0) {
      zside = 0;
    } else if (zsideinorder < 4) {
      zside = 2 * zsideinorder - 1;
    } else if (zsideinorder < 7) {
      zside = 2 * zsideinorder - 6;
    } else {
      zside = 0;
    }
    zsidereal = zside;
    //
    //old:
    if (rn0 == 3) {
      if (zside > 2)
        zsidereal = 0;
    }
    //new:
    if (rn0 == 7) {
      //    if(zside==3 || zside==6) zsidereal=0; // ????
      if (zside == 4 || zside == 5)
        zsidereal = 0;
    }
    //
    return zsidereal;
  }

unpackCopyIndex()

static int FP420NumberingScheme::unpackCopyIndex ( int  rn0, int  zside  )

inlinestatic

Definition at line 67 of file FP420NumberingScheme.h.

References a, b, createfilelist::int, and ecaldqm::zside().

Referenced by FP420Test::update().

                                                 {
    // 1,2,3
    int copyIndex = 0;
    if (zside > (rn0 - 1) || zside < 1) {
    } else {
      int layerIndex = 1, b;
      float a = (zside + 1) / 2.;
      b = (int)a;
      if (a - b != 0.)
        layerIndex = 2;
      if (layerIndex == 2)
        copyIndex = zside / 2;
      if (layerIndex == 1)
        copyIndex = (zside + 1) / 2;
    }

    return copyIndex;
  }

unpackFP420Index()

void FP420NumberingScheme::unpackFP420Index ( const unsigned int &  idx, int &  det, int &  zside, int &  station, int &  superplane  )

static

Definition at line 171 of file FP420NumberingScheme.cc.

References heavyIonCSV_trainingSettings::idx, relativeConstraints::station, and ecaldqm::zside().

Referenced by FP420Test::update().

                                                                                                                   {
  det = (idx >> 19) & 3;
  //det  = (idx>>20)&1;
  det += 1;
  zside = (idx >> 7) & 7;
  //  zside   = (idx>>7)&3;
  station = (idx >> 4) & 7;
  plane = idx & 15;
  //

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("FP420") << " FP420unpacking: idx=" << idx << " zside  " << zside << " station " << station
                            << " plane " << plane;
#endif
  //
}

unpackLayerIndex()

static int FP420NumberingScheme::unpackLayerIndex ( int  rn0, int  zside  )

inlinestatic

Definition at line 54 of file FP420NumberingScheme.h.

References a, b, createfilelist::int, and ecaldqm::zside().

Referenced by FP420Test::update().

                                                  {
    // 1,2
    int layerIndex = 1, b;
    float a = (zside + 1) / 2.;
    b = (int)a;
    if (a - b != 0.)
      layerIndex = 2;
    //
    if (zside > (rn0 - 1) || zside < 1)
      layerIndex = 0;
    return layerIndex;
  }

unpackMYIndex()

static void FP420NumberingScheme::unpackMYIndex ( const int &  idx, int  rn0, int  pn0, int  sn0, int &  det, int &  zside, int &  sector, int &  zmodule  )

inlinestatic

Definition at line 45 of file FP420NumberingScheme.h.

References heavyIonCSV_trainingSettings::idx, hgcalTBTopologyTester_cfi::sector, and ecaldqm::zside().

                                                                                                                        {
    int zScale = (rn0 - 1), sScale = (rn0 - 1) * (pn0 - 1), dScale = (rn0 - 1) * (pn0 - 1) * (sn0 - 1);

    det = (idx - 1) / dScale + 1;
    sector = (idx - 1 - dScale * (det - 1)) / sScale + 1;
    zmodule = (idx - 1 - dScale * (det - 1) - sScale * (sector - 1)) / zScale + 1;
    zside = idx - dScale * (det - 1) - sScale * (sector - 1) - zScale * (zmodule - 1);
  }

unpackOrientation()

static int FP420NumberingScheme::unpackOrientation ( int  rn0, int  zside  )

inlinestatic

Definition at line 86 of file FP420NumberingScheme.h.

References ecaldqm::zside().

Referenced by FP420Test::update().

                                                   {
    // Front: Orientation= 1; Back: Orientation= 2
    int Orientation = 2;
    if (zside > (rn0 - 1) || zside < 1)
      Orientation = 0;
    if (zside == 1 || zside == 2)
      Orientation = 1;
    //
    return Orientation;
  }