FastSimulation/Tracking/src/TrackingLayer.cc - master - githubmirror/cmssw

FastSimulation/Tracking/src/TrackingLayer.cc
223 строки10 KB
Sunanda
Code check
10 окт 2022, 08:38
10 окт 2022, 08:38d87e7da
#include "FastSimulation/Tracking/interface/TrackingLayer.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "DataFormats/SiStripDetId/interface/StripSubdetector.h"

const TrackingLayer::eqfct TrackingLayer::_eqfct;

TrackingLayer::TrackingLayer() : _subDet(Det::UNKNOWN), _side(Side::BARREL), _layerNumber(0) {}

TrackingLayer TrackingLayer::createFromDetId(const DetId& detId, const TrackerTopology& trackerTopology) {
  TrackingLayer trackingLayer;
  uint32_t subdet = detId.subdetId();
  //BPix
  if (subdet == PixelSubdetector::PixelBarrel) {
    trackingLayer._subDet = TrackingLayer::Det::PXB;
    trackingLayer._side = TrackingLayer::Side::BARREL;
    trackingLayer._layerNumber = trackerTopology.pxbLayer(detId);
  }
  //FPix
  else if (subdet == PixelSubdetector::PixelEndcap) {
    trackingLayer._subDet = TrackingLayer::Det::PXD;
    if (trackerTopology.pxfSide(detId) == 1) {
      trackingLayer._side = TrackingLayer::Side::NEG_ENDCAP;
    } else if (trackerTopology.pxfSide(detId) == 2) {
      trackingLayer._side = TrackingLayer::Side::POS_ENDCAP;
    } else {
      throw cms::Exception("FastSimulation/Tracking")
          << "Tracker hit for seeding in FPix seems neither on positive nor on negative disk side: "
          << trackerTopology.print(detId).c_str();
    }
    trackingLayer._layerNumber = trackerTopology.pxfDisk(detId);
  }
  //TIB
  else if (subdet == StripSubdetector::TIB) {
    trackingLayer._subDet = TrackingLayer::Det::TIB;
    trackingLayer._side = TrackingLayer::Side::BARREL;
    trackingLayer._layerNumber = trackerTopology.tibLayer(detId);
  }
  //TID
  else if (subdet == StripSubdetector::TID) {
    trackingLayer._subDet = TrackingLayer::Det::TID;
    if (trackerTopology.tidSide(detId) == 1) {
      trackingLayer._side = TrackingLayer::Side::NEG_ENDCAP;
    } else if (trackerTopology.tidSide(detId) == 2) {
      trackingLayer._side = TrackingLayer::Side::POS_ENDCAP;
    } else {
      throw cms::Exception("FastSimulation/Tracking")
          << "Tracker hit for seeding in TID seems neither on positive nor on negative disk side: "
          << trackerTopology.print(detId).c_str();
    }
    trackingLayer._layerNumber = trackerTopology.tidWheel(detId);
  }
  //TOB
  else if (subdet == StripSubdetector::TOB) {
    trackingLayer._subDet = TrackingLayer::Det::TOB;
    trackingLayer._side = TrackingLayer::Side::BARREL;
    trackingLayer._layerNumber = trackerTopology.tobLayer(detId);
  }
  //TEC
  else if (subdet == StripSubdetector::TEC) {
    trackingLayer._subDet = TrackingLayer::Det::TEC;
    if (trackerTopology.tecSide(detId) == 1) {
      trackingLayer._side = TrackingLayer::Side::NEG_ENDCAP;
    } else if (trackerTopology.tecSide(detId) == 2) {
      trackingLayer._side = TrackingLayer::Side::POS_ENDCAP;
    } else {
      throw cms::Exception("FastSimulation/Tracking")
          << "Tracker hit for seeding in TEC seems neither on positive nor on negative disk side: "
          << trackerTopology.print(detId).c_str();
    }
    trackingLayer._layerNumber = trackerTopology.tecWheel(detId);
  } else {
    throw cms::Exception("FastSimulation/Tracking")
        << "Cannot determine seeding layer from DetId:" << trackerTopology.print(detId).c_str() << std::endl;
  }
  edm::LogVerbatim("TrackingLayer") << "LayerSpec::createFromDetId: " << trackerTopology.print(detId).c_str()
                                    << ", parsed=" << trackingLayer._layerNumber;
  return trackingLayer;
}

TrackingLayer TrackingLayer::createFromString(std::string layerSpecification) {
  TrackingLayer trackingLayer;
  if (layerSpecification.substr(0, 4) == "BPix") {
    trackingLayer._subDet = TrackingLayer::Det::PXB;
    trackingLayer._side = TrackingLayer::Side::BARREL;
    trackingLayer._layerNumber = std::atoi(layerSpecification.substr(4, 1).c_str());
  } else if (layerSpecification.substr(0, 4) == "FPix") {
    trackingLayer._subDet = TrackingLayer::Det::PXD;
    trackingLayer._layerNumber = std::atoi(layerSpecification.substr(4, 1).c_str());
    if (layerSpecification.substr(layerSpecification.size() - 3) == "pos") {
      trackingLayer._side = TrackingLayer::Side::POS_ENDCAP;
    } else if (layerSpecification.substr(layerSpecification.size() - 3) == "neg") {
      trackingLayer._side = TrackingLayer::Side::NEG_ENDCAP;
    } else {
      throw cms::Exception("FastSimulation/Tracking/python")
          << "FPix seeding layer configuration '" << layerSpecification.c_str()
          << "' does not specify the side correctly!";
    }

  } else if (layerSpecification.substr(0, 3) == "TIB") {
    trackingLayer._subDet = TrackingLayer::Det::TIB;
    trackingLayer._side = TrackingLayer::Side::BARREL;
    trackingLayer._layerNumber = std::atoi(layerSpecification.substr(3, 1).c_str());
  } else if (layerSpecification.substr(0, 4) == "MTIB") {
    trackingLayer._subDet = TrackingLayer::Det::TIB;
    trackingLayer._side = TrackingLayer::Side::BARREL;
    trackingLayer._layerNumber = std::atoi(layerSpecification.substr(4, 1).c_str());
  }

  else if (layerSpecification.substr(0, 3) == "TID") {
    trackingLayer._subDet = TrackingLayer::Det::TID;
    trackingLayer._layerNumber = std::atoi(layerSpecification.substr(3, 1).c_str());
    if (layerSpecification.substr(layerSpecification.size() - 3) == "pos") {
      trackingLayer._side = TrackingLayer::Side::POS_ENDCAP;
    } else if (layerSpecification.substr(layerSpecification.size() - 3) == "neg") {
      trackingLayer._side = TrackingLayer::Side::NEG_ENDCAP;
    } else {
      throw cms::Exception("FastSimulation/Tracking/python")
          << "TID seeding layer configuration '" << layerSpecification.c_str()
          << "' does not specify the side correctly!";
    }
  } else if (layerSpecification.substr(0, 4) == "MTID") {
    trackingLayer._subDet = TrackingLayer::Det::TID;
    trackingLayer._layerNumber = std::atoi(layerSpecification.substr(4, 1).c_str());
    if (layerSpecification.substr(layerSpecification.size() - 3) == "pos") {
      trackingLayer._side = TrackingLayer::Side::POS_ENDCAP;
    } else if (layerSpecification.substr(layerSpecification.size() - 3) == "neg") {
      trackingLayer._side = TrackingLayer::Side::NEG_ENDCAP;
    } else {
      throw cms::Exception("FastSimulation/Tracking/python")
          << "MTID seeding layer configuration '" << layerSpecification.c_str()
          << "' does not specify the side correctly!";
    }
  } else if (layerSpecification.substr(0, 3) == "TOB") {
    trackingLayer._subDet = TrackingLayer::Det::TOB;
    trackingLayer._side = TrackingLayer::Side::BARREL;
    trackingLayer._layerNumber = std::atoi(layerSpecification.substr(3, 1).c_str());
  } else if (layerSpecification.substr(0, 4) == "MTOB") {
    trackingLayer._subDet = TrackingLayer::Det::TOB;
    trackingLayer._side = TrackingLayer::Side::BARREL;
    trackingLayer._layerNumber = std::atoi(layerSpecification.substr(4, 1).c_str());
  } else if (layerSpecification.substr(0, 3) == "TEC") {
    trackingLayer._subDet = TrackingLayer::Det::TEC;
    trackingLayer._layerNumber = std::atoi(layerSpecification.substr(3, 1).c_str());
    if (layerSpecification.substr(layerSpecification.size() - 3) == "pos") {
      trackingLayer._side = TrackingLayer::Side::POS_ENDCAP;
    } else if (layerSpecification.substr(layerSpecification.size() - 3) == "neg") {
      trackingLayer._side = TrackingLayer::Side::NEG_ENDCAP;
    } else {
      throw cms::Exception("FastSimulation/Tracking/python")
          << "TEC seeding layer configuration '" << layerSpecification.c_str()
          << "' does not specify the side correctly!";
    }

  } else if (layerSpecification.substr(0, 4) == "MTEC") {
    trackingLayer._subDet = TrackingLayer::Det::TEC;
    trackingLayer._layerNumber = std::atoi(layerSpecification.substr(4, 1).c_str());
    if (layerSpecification.substr(layerSpecification.size() - 3) == "pos") {
      trackingLayer._side = TrackingLayer::Side::POS_ENDCAP;
    } else if (layerSpecification.substr(layerSpecification.size() - 3) == "neg") {
      trackingLayer._side = TrackingLayer::Side::NEG_ENDCAP;
    } else {
      throw cms::Exception("FastSimulation/Tracking/python")
          << "MTEC seeding layer configuration '" << layerSpecification.c_str()
          << "' does not specify the side correctly!";
    }
  } else {
    throw cms::Exception("FastSimulation/Tracking/python")
        << "Bad data naming in seeding layer configuration."
        << "no case sensitive name of ['BPix','FPix','TIB','MTIB','TID','MTID','TOB','TEC','MTEC'] matches '"
        << layerSpecification.c_str() << "'";
  }
  edm::LogVerbatim("TrackingLayer") << "LayerSpec::createFromString: " << layerSpecification.c_str()
                                    << ", parsed=" << trackingLayer._layerNumber;

  return trackingLayer;
}

std::string TrackingLayer::toString() const {
  std::stringstream ss;
  switch (_subDet) {
    case TrackingLayer::Det::UNKNOWN:
      ss << "unknown";
      break;
    case TrackingLayer::Det::PXB:
      ss << " BPix";
      break;
    case TrackingLayer::Det::PXD:
      ss << " FPix";
      break;
    case TrackingLayer::Det::TIB:
      ss << " TIB";
      break;
    case TrackingLayer::Det::TID:
      ss << " TID";
      break;
    case TrackingLayer::Det::TOB:
      ss << " TOB";
      break;
    case TrackingLayer::Det::TEC:
      ss << " TEC";
      break;
  }
  ss << _layerNumber;
  switch (_side) {
    case TrackingLayer::Side::BARREL:
      break;
    case Side::NEG_ENDCAP:
      ss << "_neg";
      break;
    case TrackingLayer::Side::POS_ENDCAP:
      ss << "_pos";
      break;
  }

  return ss.str();
}

std::string TrackingLayer::toIdString() const {
  std::stringstream ss;
  ss << getSubDetNumber() << ":" << getLayerNumber() << ":" << getSideNumber();
  return ss.str();
}