PixelEndcapName.cc

Go to the documentation of this file.

#include "DataFormats/TrackerCommon/interface/PixelEndcapName.h"

#include <sstream>
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"

using namespace std;

namespace {
  //const bool phase1 = false;
  const bool pilot_blade = false;
}  // namespace

// Decodes the pixel name from the cmssw DetID, uses tracker topology
PixelEndcapName::PixelEndcapName(const DetId& id, const TrackerTopology* tt, bool phase)
    : PixelModuleName(false), thePart(mO), theDisk(0), theBlade(0), thePannel(0), thePlaquette(0), phase1(phase) {
  //PXFDetId cmssw_numbering(id);
  int side = tt->pxfSide(id);
  int tmpBlade = tt->pxfBlade(id);
  theDisk = tt->pxfDisk(id);
  thePannel = tt->pxfPanel(id);
  bool outer = false;  // outer means with respect to the LHC ring (x - axis) Bm(p)I/O

  if (phase1) {  // phase 1

    int ring = 0;  // ring number , according to radius, 1-lower, 2- higher

    if (tmpBlade >= 7 && tmpBlade <= 17) {  // ring 1, O
      outer = true;
      theBlade = tmpBlade - 6;  //7...17-->1...11
      ring = 1;                 // lower radius

    } else if (tmpBlade <= 6) {  // ring 1, I
      theBlade = 7 - tmpBlade;   //1...6-->6...1
      ring = 1;

    } else if (tmpBlade >= 18 && tmpBlade <= 22) {  // ring 1, I
      theBlade = 29 - tmpBlade;                     //18...22-->11...7
      ring = 1;

    } else if (tmpBlade >= 32 && tmpBlade <= 48) {  // ring 2, O
      outer = true;
      theBlade = tmpBlade - 31;  //32...48-->28...12
      ring = 2;                  // higher radius

    } else if (tmpBlade >= 23 && tmpBlade <= 31) {  // ring 2, I
      theBlade = 32 - tmpBlade;                     //23...31-->9...1
      ring = 2;

    } else if (tmpBlade >= 49 && tmpBlade <= 56) {  // ring 2, I
      theBlade = 66 - tmpBlade;                     //49...56-->28...21
      ring = 2;

    }  //  end

    thePlaquette = ring;

  } else {  // phase0

    // hack for the pilot blade
    if (pilot_blade && theDisk == 3) {  // do only for disk 3
      //cout<<tmpBlade<<" "<<theDisk<<endl;
      if (tmpBlade >= 1 && tmpBlade <= 4) {
        // convert the sequential counting of pilot blades to std. cmssw convention
        if (tmpBlade < 3)
          tmpBlade += 3;
        else
          tmpBlade += 13;
      } else {
        edm::LogError("Bad PilotBlade blade number ") << tmpBlade;
      }
    }

    if (tmpBlade >= 7 && tmpBlade <= 18) {
      outer = true;
      theBlade = tmpBlade - 6;
    } else if (tmpBlade <= 6) {
      theBlade = 7 - tmpBlade;
    } else if (tmpBlade >= 19) {
      theBlade = 31 - tmpBlade;
    }

    thePlaquette = tt->pxfModule(id);  // ring number is coded as plaqutte (plaqs are unused)

  }  // end phase1

  if (side == 1 && outer)
    thePart = mO;
  else if (side == 1 && !outer)
    thePart = mI;
  else if (side == 2 && outer)
    thePart = pO;
  else if (side == 2 && !outer)
    thePart = pI;
}

// Decodes the pixel name from the cmssw DetID, uses old pixel name classes
PixelEndcapName::PixelEndcapName(const DetId& id, bool phase)
    : PixelModuleName(false), thePart(mO), theDisk(0), theBlade(0), thePannel(0), thePlaquette(0), phase1(phase) {
  PXFDetId cmssw_numbering(id);
  int side = cmssw_numbering.side();
  int tmpBlade = cmssw_numbering.blade();
  theDisk = cmssw_numbering.disk();
  thePannel = cmssw_numbering.panel();
  bool outer = false;  // outer means with respect to the LHC ring (x - axis)

  if (phase1) {  // phase1

    int ring = 0;  // ring number , according to radius, 1-lower, 2- higher

    if (tmpBlade >= 7 && tmpBlade <= 17) {  // ring 1, O
      outer = true;
      theBlade = tmpBlade - 6;  //7...17-->1...11
      ring = 1;                 // lower radius

    } else if (tmpBlade <= 6) {  // ring 1, I
      theBlade = 7 - tmpBlade;   //1...6-->6...1
      ring = 1;

    } else if (tmpBlade >= 18 && tmpBlade <= 22) {  // ring 1, I
      theBlade = 29 - tmpBlade;                     //18...22-->11...7
      ring = 1;

    } else if (tmpBlade >= 32 && tmpBlade <= 48) {  // ring 2, O
      outer = true;
      theBlade = tmpBlade - 31;  //32...48-->28...12
      ring = 2;                  // higher radius

    } else if (tmpBlade >= 23 && tmpBlade <= 31) {  // ring 2, I
      theBlade = 32 - tmpBlade;                     //23...31-->9...1
      ring = 2;

    } else if (tmpBlade >= 49 && tmpBlade <= 56) {  // ring 2, I
      theBlade = 66 - tmpBlade;                     //49...56-->28...21
      ring = 2;

    }  //  end

    thePlaquette = ring;

  } else {  // phase 0

    // hack for the pilot blade
    if (pilot_blade && theDisk == 3) {  // do only for disk 3
      if (tmpBlade >= 1 && tmpBlade <= 4) {
        // convert the sequential counting of pilot blades to std. cmssw convention
        if (tmpBlade < 3)
          tmpBlade += 3;
        else
          tmpBlade += 13;
      } else {
        edm::LogError("Bad PilotBlade blade number ") << tmpBlade;
      }
    }

    if (tmpBlade >= 7 && tmpBlade <= 18) {
      outer = true;  // outer means with respect to the LHC ring (x - axis)
      theBlade = tmpBlade - 6;
    } else if (tmpBlade <= 6) {
      theBlade = 7 - tmpBlade;
    } else if (tmpBlade >= 19) {
      theBlade = 31 - tmpBlade;
    }

    thePlaquette = cmssw_numbering.module();

  }  // end phase1

  if (side == 1 && outer)
    thePart = mO;
  else if (side == 1 && !outer)
    thePart = mI;
  else if (side == 2 && outer)
    thePart = pO;
  else if (side == 2 && !outer)
    thePart = pI;
}

// constructor from name string
PixelEndcapName::PixelEndcapName(std::string name, bool phase)
    : PixelModuleName(false), thePart(mO), theDisk(0), theBlade(0), thePannel(0), thePlaquette(0), phase1(phase) {
  // parse the name string
  // first, check to make sure this is an FPix name, should start with "FPix_"
  // also check to make sure the needed parts are present
  if ((name.substr(0, 5) != "FPix_") || (name.find("_B") == string::npos) || (name.find("_D") == string::npos) ||
      (name.find("_BLD") == string::npos) || (name.find("_PNL") == string::npos) ||
      ((phase1 && name.find("_RNG") == string::npos)) || ((!phase1 && name.find("_PLQ") == string::npos))) {
    edm::LogError("BadNameString|SiPixel")
        << "Bad name string in PixelEndcapName::PixelEndcapName(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"));

  // get the half cylinder
  string hcString = name.substr(name.find("_B") + 2, name.find("_D") - name.find("_B") - 2);
  if (hcString == "mO")
    thePart = mO;
  else if (hcString == "mI")
    thePart = mI;
  else if (hcString == "pO")
    thePart = pO;
  else if (hcString == "pI")
    thePart = pI;
  else {
    edm::LogError("BadNameString|SiPixel")
        << "Unable to determine half cylinder in PixelEndcapName::PixelEndcapName(std::string): " << name;
  }

  // get the disk
  string diskString = name.substr(name.find("_D") + 2, name.find("_BLD") - name.find("_D") - 2);
  if (diskString == "1")
    theDisk = 1;
  else if (diskString == "2")
    theDisk = 2;
  else if (diskString == "3")
    theDisk = 3;
  else {
    edm::LogError("BadNameString|SiPixel")
        << "Unable to determine disk number in PixelEndcapName::PixelEndcapName(std::string): " << name;
  }

  // get the blade
  string bladeString = name.substr(name.find("_BLD") + 4, name.find("_PNL") - name.find("_BLD") - 4);
  // since atoi() doesn't report errors, do it the long way
  if (bladeString == "1")
    theBlade = 1;
  else if (bladeString == "2")
    theBlade = 2;
  else if (bladeString == "3")
    theBlade = 3;
  else if (bladeString == "4")
    theBlade = 4;
  else if (bladeString == "5")
    theBlade = 5;
  else if (bladeString == "6")
    theBlade = 6;
  else if (bladeString == "7")
    theBlade = 7;
  else if (bladeString == "8")
    theBlade = 8;
  else if (bladeString == "9")
    theBlade = 9;
  else if (bladeString == "10")
    theBlade = 10;
  else if (bladeString == "11")
    theBlade = 11;
  else if (bladeString == "12")
    theBlade = 12;
  else if (bladeString == "13")
    theBlade = 13;
  else if (bladeString == "14")
    theBlade = 14;
  else if (bladeString == "15")
    theBlade = 15;
  else if (bladeString == "16")
    theBlade = 16;
  else if (bladeString == "17")
    theBlade = 17;
  else {
    edm::LogError("BadNameString|SiPixel")
        << "Unable to determine blade number in PixelEndcapName::PixelEndcapName(std::string): " << name;
  }

  // find the panel
  string panelString;
  if (phase1)
    panelString = name.substr(name.find("_PNL") + 4, name.find("_RNG") - name.find("_PNL") - 4);
  else
    panelString = name.substr(name.find("_PNL") + 4, name.find("_PLQ") - name.find("_PNL") - 4);

  if (panelString == "1")
    thePannel = 1;
  else if (panelString == "2")
    thePannel = 2;
  else {
    edm::LogError("BadNameString|SiPixel")
        << "Unable to determine panel number in PixelEndcapName::PixelEndcapName(std::string): " << name;
  }

  // find the plaquette, for phase 1 this is the rung number
  if (phase1) {  // phase1

    string ringString = name.substr(name.find("_RNG") + 4, name.size() - name.find("_RNG") - 4);
    if (ringString == "1")
      thePlaquette = 1;  // code ring in the plaquette
    else if (ringString == "2")
      thePlaquette = 2;
    else {
      edm::LogError("BadNameString|SiPixel")
          << "Unable to determine ring number in PixelEndcapName::PixelEndcapName(std::string): " << name;
    }

  } else {  // phase 0

    // find the plaquette
    string plaquetteString = name.substr(name.find("_PLQ") + 4, name.size() - name.find("_PLQ") - 4);
    if (plaquetteString == "1")
      thePlaquette = 1;
    else if (plaquetteString == "2")
      thePlaquette = 2;
    else if (plaquetteString == "3")
      thePlaquette = 3;
    else if (plaquetteString == "4")
      thePlaquette = 4;
    else {
      edm::LogError("BadNameString|SiPixel")
          << "Unable to determine plaquette number in PixelEndcapName::PixelEndcapName(std::string): " << name;
    }

  }  // end phase1

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

PixelModuleName::ModuleType PixelEndcapName::moduleType() const {
  ModuleType type = v1x2;

  if (phase1) {  // phase1

    type = v2x8;

  } else {  // phase 0

    if (pannelName() == 1) {
      if (plaquetteName() == 1) {
        type = v1x2;
      } else if (plaquetteName() == 2) {
        type = v2x3;
      } else if (plaquetteName() == 3) {
        type = v2x4;
      } else if (plaquetteName() == 4) {
        type = v1x5;
      }
    } else {
      if (plaquetteName() == 1) {
        type = v2x3;
      } else if (plaquetteName() == 2) {
        type = v2x4;
      } else if (plaquetteName() == 3) {
        type = v2x5;
      }
    }

    // hack for the pilot blade
    if (pilot_blade && theDisk == 3) {
      type = v2x8;
    }  // do only for disk 3

  }  // end phase1

  return type;
}

bool PixelEndcapName::operator==(const PixelModuleName& other) const {
  return other.isBarrel() ? false : (dynamic_cast<const PixelEndcapName&>(other) == *this);
}

string PixelEndcapName::name() const {
  std::ostringstream stm;

  if (phase1)
    stm << "FPix_B" << thePart << "_D" << theDisk << "_BLD" << theBlade << "_PNL" << thePannel << "_RNG"
        << thePlaquette;
  else
    stm << "FPix_B" << thePart << "_D" << theDisk << "_BLD" << theBlade << "_PNL" << thePannel << "_PLQ"
        << thePlaquette;

  return stm.str();
}

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

// return the DetId
DetId PixelEndcapName::getDetId(const TrackerTopology* tt) {
  uint32_t side = 0;
  uint32_t disk = 0;
  uint32_t blade = 0;
  uint32_t panel = 0;
  uint32_t module = 0;

  // figure out the side
  HalfCylinder hc = halfCylinder();
  if (hc == mO || hc == mI)
    side = 1;
  else if (hc == pO || hc == pI)
    side = 2;

  // get disk/blade/panel/module numbers from PixelEndcapName object
  disk = static_cast<uint32_t>(diskName());
  uint32_t tmpBlade = static_cast<uint32_t>(bladeName());
  panel = static_cast<uint32_t>(pannelName());

  // convert blade numbering to cmssw convention
  bool outer = false;
  outer = (hc == mO) || (hc == pO);

  if (phase1) {  // phase1

    module = 1;
    int ring = static_cast<uint32_t>(ringName());  // this is ring for phase1
    if (ring == 1) {                               // ring 1, lower radius
      if (outer) {
        if (tmpBlade >= 1 && tmpBlade <= 11)
          blade = tmpBlade + 6;
      } else {  //inner
        if (tmpBlade <= 6)
          blade = 7 - tmpBlade;
        else
          blade = 29 - tmpBlade;
      }

    } else if (ring == 2) {  //ring 2, upper radius
      if (outer) {
        if (tmpBlade >= 1 && tmpBlade <= 17)
          blade = tmpBlade + 31;
      } else {  //inner
        if (tmpBlade <= 9)
          blade = 32 - tmpBlade;
        else
          blade = 66 - tmpBlade;
      }
    }

  } else {  // phase 0

    if (outer) {
      blade = tmpBlade + 6;
    } else {  // inner
      if (tmpBlade <= 6)
        blade = 7 - tmpBlade;
      else if (tmpBlade <= 12)
        blade = 31 - tmpBlade;
    }

    // hack for the pilot blade
    if (pilot_blade && theDisk == 3) {  // do only for disk 3
      //cout<<tmpBlade<<" "<<blade<<endl;
      if (blade <= 5)
        blade -= 3;
      else
        blade -= 13;
      //cout<<tmpBlade<<" "<<blade<<endl;
    }

    module = static_cast<uint32_t>(plaquetteName());
  }  // end phase1

  // create and return the DetId
  DetId id = tt->pxfDetId(side, disk, blade, panel, module);

  return id;

}  // PXFDetId PixelEndcapName::getDetId()

// return the DetId
PXFDetId PixelEndcapName::getDetId() {
  uint32_t side = 0;
  uint32_t disk = 0;
  uint32_t blade = 0;
  uint32_t panel = 0;
  uint32_t module = 0;

  // figure out the side
  HalfCylinder hc = halfCylinder();
  if (hc == mO || hc == mI)
    side = 1;
  else if (hc == pO || hc == pI)
    side = 2;

  // get disk/blade/panel/module numbers from PixelEndcapName object
  disk = static_cast<uint32_t>(diskName());
  uint32_t tmpBlade = static_cast<uint32_t>(bladeName());
  panel = static_cast<uint32_t>(pannelName());

  // convert blade numbering to cmssw convention
  bool outer = false;
  outer = (hc == mO) || (hc == pO);

  if (phase1) {  // phase1

    module = 1;                                    // for phase 1 always 1,
    int ring = static_cast<uint32_t>(ringName());  // this is ring for phase1
    if (ring == 1) {                               // ring 1, lower radius
      if (outer) {
        if (tmpBlade >= 1 && tmpBlade <= 11)
          blade = tmpBlade + 6;
      } else {  //inner
        if (tmpBlade <= 6)
          blade = 7 - tmpBlade;
        else
          blade = 29 - tmpBlade;
      }

    } else if (ring == 2) {  //ring 2, upper radius
      if (outer) {
        if (tmpBlade >= 1 && tmpBlade <= 17)
          blade = tmpBlade + 31;
      } else {  //inner
        if (tmpBlade <= 9)
          blade = 32 - tmpBlade;
        else
          blade = 66 - tmpBlade;
      }
    }

  } else {  // phase 0

    if (outer) {
      blade = tmpBlade + 6;
    } else {  // inner
      if (tmpBlade <= 6)
        blade = 7 - tmpBlade;
      else if (tmpBlade <= 12)
        blade = 31 - tmpBlade;
    }

    // hack for the pilot blade
    if (pilot_blade && theDisk == 3) {  // do only for disk 3
      //cout<<tmpBlade<<" "<<blade<<endl;
      if (blade <= 5)
        blade -= 3;
      else
        blade -= 13;
      //cout<<tmpBlade<<" "<<blade<<endl;
    }

    module = static_cast<uint32_t>(plaquetteName());
  }  // end phase1

  // create and return the DetId
  return PXFDetId(side, disk, blade, panel, module);

}  // PXFDetId PixelEndcapName::getDetId()