PixelBarrelName.cc

Go to the documentation of this file.

#include "DataFormats/SiPixelDetId/interface/PixelBarrelName.h"

#include <sstream>
#include <iostream>

#include "DataFormats/SiPixelDetId/interface/PXBDetId.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

using namespace std;


PixelBarrelName::PixelBarrelName(const DetId & id) 
  : PixelModuleName(true)
{
 
//  uint32_t rawId = id.rawId(); 
  PXBDetId cmssw_numbering(id);


  theLayer = cmssw_numbering.layer();

  int oldModule = cmssw_numbering.module() -4; if (oldModule<=0) oldModule--;

  int oldLadder = cmssw_numbering.ladder();
  if (theLayer == 1) {
    if (oldLadder <= 5) oldLadder = 6-oldLadder;
    else if (oldLadder >= 6 && oldLadder <= 15 ) oldLadder = 5-oldLadder;
    else if (oldLadder >= 16) oldLadder = 26-oldLadder;
  } 
  else if (theLayer == 2) {
    if (oldLadder <= 8) oldLadder = 9-oldLadder;
    else if (oldLadder >= 9 && oldLadder <= 24) oldLadder = 8-oldLadder;
    else if (oldLadder >= 25) oldLadder = 41-oldLadder; 
  } 
  else if (theLayer == 3) {
    if (oldLadder <= 11) oldLadder = 12-oldLadder;
    else if (oldLadder >= 12 && oldLadder <= 33) oldLadder = 11-oldLadder;
    else if (oldLadder >= 34) oldLadder = 56-oldLadder;
  } 

  //
  // part
  //
  if      (oldModule < 0 && oldLadder < 0) thePart = mO; 
  else if (oldModule > 0 && oldLadder < 0) thePart = pO;
  else if (oldModule < 0 && oldLadder > 0) thePart = mI;
  else if (oldModule > 0 && oldLadder > 0) thePart = pI;
  

  //
  // ladder
  //
  theLadder = abs(oldLadder);

  //
  // module
  //
  theModule = abs(oldModule);
 
}

// constructor from name string
PixelBarrelName::PixelBarrelName(std::string name) 
  : PixelModuleName(true), thePart(mO), theLayer(0),
    theModule(0), theLadder(0) {

  // parse the name string
  // first, check to make sure this is an BPix name, should start with "BPix_"
  // also check to make sure the needed parts are present
  if ( (name.substr(0, 5) != "BPix_") ||
       (name.find("_B") == string::npos) || 
       (name.find("_LYR") == string::npos) ||
       (name.find("_LDR") == string::npos) || 
       (name.find("_MOD") == string::npos) ) {
    edm::LogError ("BadNameString|SiPixel") 
      << "Bad name string in PixelBarrelName::PixelBarrelName(std::string): "
      << name;
    return;
  }

  // strip off ROC part if it's there
  if (name.find("_ROC") != string::npos)
    name = name.substr(0, name.find("_ROC"));

  // find shell
  string shellString = name.substr(name.find("_B")+2, name.find("_SEC")-name.find("_B")-2);
  if (shellString == "mO") thePart = mO;
  else if (shellString == "mI") thePart = mI;
  else if (shellString == "pO") thePart = pO;
  else if (shellString == "pI") thePart = pI;
  else {
    edm::LogError ("BadNameString|SiPixel")
      << "Unable to determine shell in PixelBarrelName::PixelBarrelName(std::string): "
      << name;
  }

  // find the layer
  string layerString = name.substr(name.find("_LYR")+4, name.find("_LDR")-name.find("_LYR")-4);
  if (layerString == "1") theLayer = 1;
  else if (layerString == "2") theLayer = 2;
  else if (layerString == "3") theLayer = 3;
  else {
    edm::LogError ("BadNameString|SiPixel")
      << "Unable to determine layer in PixelBarrelName::PixelBarrelName(std::string): "
      << name;
  }

  // find the ladder
  string ladderString = name.substr(name.find("_LDR")+4, name.find("_MOD")-name.find("_LDR")-4);
  if (ladderString == "1H") theLadder = 1;
  else if (ladderString == "10H" && theLayer == 1) theLadder = 10;
  else if (ladderString == "16H" && theLayer == 2) theLadder = 16;
  else if (ladderString == "22H" && theLayer == 3) theLadder = 22;
  else if (ladderString.substr(ladderString.size()-1, 1) == "F") {
    int ladderNum = atoi(ladderString.substr(0, ladderString.size() -1).c_str());
    if (theLayer == 1 && ladderNum > 1 && ladderNum < 10) theLadder = ladderNum;
    else if (theLayer == 2 && ladderNum > 1 && ladderNum < 16) theLadder = ladderNum;
    else if (theLayer == 3 && ladderNum > 1 && ladderNum < 22) theLadder = ladderNum;
    else {
      edm::LogError ("BadNameString|SiPixel")
    << "Unable to determine ladder in PixelBarrelName::PixelBarrelName(std::string): "
    << name;
    }
  } // full ladders
  else {
    edm::LogError ("BadNameString|SiPixel")
      << "Unable to determine ladder in PixelBarrelName::PixelBarrelName(std::string): "
      << name;
  }

  // find the module
  string moduleString = name.substr(name.find("_MOD")+4, name.size()-name.find("_MOD")-4);
  if (moduleString == "1") theModule = 1;
  else if (moduleString == "2") theModule = 2;
  else if (moduleString == "3") theModule = 3;
  else if (moduleString == "4") theModule = 4;
  else {
    edm::LogError ("BadNameString|SiPixel")
      << "Unable to determine module in PixelBarrelName::PixelBarrelName(std::string): "
      << name;
  }

} // PixelBarrelName::PixelBarrelName(std::string name)

int PixelBarrelName::sectorName() const
{
  int sector = 0;
  if (theLayer==1) {
    switch (theLadder) {
    case 1 : case 2: {sector = 1; break;}
    case 3 :         {sector = 2; break;}
    case 4 :         {sector = 3; break;}
    case 5 :         {sector = 4; break;}
    case 6 :         {sector = 5; break;}
    case 7 :         {sector = 6; break;}
    case 8 :         {sector = 7; break;}
    case 9 : case 10:{sector = 8; break;}
    default: ;
    };
  } else if (theLayer==2) {
    switch (theLadder) {
    case  1 : case  2: {sector = 1; break;}
    case  3 : case  4: {sector = 2; break;}
    case  5 : case  6: {sector = 3; break;}
    case  7 : case  8: {sector = 4; break;}
    case  9 : case 10: {sector = 5; break;}
    case 11 : case 12: {sector = 6; break;}
    case 13 : case 14: {sector = 7; break;}
    case 15 : case 16: {sector = 8; break;}
    default: ;
    };
  } else if (theLayer==3) {
    switch (theLadder) {
    case  1 : case  2: case  3: {sector = 1; break;}
    case  4 : case  5: case  6: {sector = 2; break;}
    case  7 : case  8: case  9: {sector = 3; break;}
    case 10 : case 11:          {sector = 4; break;}
    case 12 : case 13:          {sector = 5; break;}
    case 14 : case 15: case 16: {sector = 6; break;}
    case 17 : case 18: case 19: {sector = 7; break;}
    case 20 : case 21: case 22: {sector = 8; break;}
    default: ;
    };
  }
  return sector;

}

bool PixelBarrelName::isHalfModule() const 
{
  bool halfModule = false;
  if (theLadder == 1) halfModule = true;
  if (theLayer == 1 && theLadder == 10) halfModule = true;
  if (theLayer == 2 && theLadder == 16) halfModule = true;
  if (theLayer == 3 && theLadder == 22) halfModule = true;
  return halfModule;
}

PixelModuleName::ModuleType  PixelBarrelName::moduleType() const 
{
  return isHalfModule() ? PixelBarrelName::v1x8 : PixelBarrelName::v2x8;
}

bool PixelBarrelName::operator==(const PixelModuleName &o) const
{
  if ( o.isBarrel() ) { 
    const PixelBarrelName *other = dynamic_cast<const PixelBarrelName*>(&o);
    return (    other
             && thePart   == other->thePart 
             && theLayer  == other->theLayer 
             && theModule == other->theModule 
             && theLadder == other->theLadder);
  } else return false; 
}

string PixelBarrelName::name() const 
{
   std::ostringstream stm;
   
   stm<<"BPix_B"<<thePart<<"_SEC"<<sectorName()<<"_LYR"<<theLayer<<"_LDR"<<theLadder;
   if ( isHalfModule() ) stm <<"H"; else stm <<"F";
   stm << "_MOD" << theModule;

   return stm.str();
}

// return the DetId
PXBDetId PixelBarrelName::getDetId() {
  
  uint32_t layer = 0;
  uint32_t ladder = 0;
  uint32_t module = 0;

  layer = layerName();
  uint32_t tmpLadder = ladderName();
  uint32_t tmpModule = moduleName();

  // translate the ladder number from the naming convention to the cmssw convention
  bool outer = false;
  Shell shell = thePart;
  outer = (shell == mO) || (shell == pO);
  if (outer) {
    if (layer == 1)
      ladder = tmpLadder + 5;
    else if (layer == 2)
      ladder = tmpLadder + 8;
    else if (layer == 3)
      ladder = tmpLadder + 11;
  } // outer
  else { // inner
    if (layer == 1) {
      if (tmpLadder <= 5) ladder = 6 - tmpLadder;
      else if (tmpLadder <= 10) ladder = 26 - tmpLadder;
    } // layer 1
    else if (layer == 2) {
      if (tmpLadder <= 8) ladder = 9 - tmpLadder;
      else if (tmpLadder <= 16) ladder = 41 - tmpLadder;
    } // layer 2
    else if (layer == 3) {
      if (tmpLadder <= 11) ladder = 12 - tmpLadder;
      else if (tmpLadder <= 22) ladder = 56 - tmpLadder;
    } // layer 3
  } // inner

  // translate the module number from naming convention to cmssw convention
  // numbering starts at positive z
  if (shell == pO || shell == pI)
    module = tmpModule + 4;
  else // negative z side
    module = 5 - tmpModule;

  return PXBDetId(layer, ladder, module);

} // PXBDetId PixelBarrelName::getDetId()

std::ostream & operator<<( std::ostream& out, const PixelBarrelName::Shell& t)
{
  switch (t) {
    case(PixelBarrelName::pI) : {out << "pI"; break;}
    case(PixelBarrelName::pO) : {out << "pO"; break;}
    case(PixelBarrelName::mI) : {out << "mI"; break;}
    case(PixelBarrelName::mO) : {out << "mO"; break;}
    default: out << "unknown";
  };
  return out;
}