TrackingTools/TrackAssociator/src/DetIdInfo.cc - master - githubmirror/cmssw

TrackingTools/TrackAssociator/src/DetIdInfo.cc
179 строк6 KB
Cms Build
Clang-Format
17 июн 2019, 12:15
17 июн 2019, 12:15514e187
#include "DetIdInfo.h"

#include "DataFormats/MuonDetId/interface/MuonSubdetId.h"
#include "DataFormats/MuonDetId/interface/DTChamberId.h"
#include "DataFormats/MuonDetId/interface/CSCDetId.h"
#include "DataFormats/MuonDetId/interface/RPCDetId.h"
#include "DataFormats/MuonDetId/interface/GEMDetId.h"
#include "DataFormats/MuonDetId/interface/ME0DetId.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"

#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"

#include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"
#include "DataFormats/SiStripDetId/interface/StripSubdetector.h"

#include "DataFormats/CaloTowers/interface/CaloTowerDetId.h"

#include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
#include "DataFormats/HcalDetId/interface/HcalDetId.h"

#include <sstream>

std::string DetIdInfo::info(const DetId& id, const TrackerTopology* tTopo) {
  std::ostringstream oss;

  oss << "DetId: " << id.rawId() << "\n";

  switch (id.det()) {
    case DetId::Tracker:
      switch (id.subdetId()) {
        case StripSubdetector::TIB: {
          oss << "TIB ";
        } break;
        case StripSubdetector::TOB: {
          oss << "TOB ";
        } break;
        case StripSubdetector::TEC: {
          oss << "TEC ";
        } break;
        case StripSubdetector::TID: {
          oss << "TID ";
        } break;
        case (int)PixelSubdetector::PixelBarrel: {
          oss << "PixBarrel ";
        } break;
        case (int)PixelSubdetector::PixelEndcap: {
          oss << "PixEndcap ";
        } break;
      }
      if (tTopo != nullptr)
        oss << tTopo->layer(id);
      break;

    case DetId::Muon:
      switch (id.subdetId()) {
        case MuonSubdetId::DT: {
          DTChamberId detId(id.rawId());
          oss << "DT chamber (wheel, station, sector): " << detId.wheel() << ", " << detId.station() << ", "
              << detId.sector();
        } break;
        case MuonSubdetId::CSC: {
          CSCDetId detId(id.rawId());
          oss << "CSC chamber (endcap, station, ring, chamber, layer): " << detId.endcap() << ", " << detId.station()
              << ", " << detId.ring() << ", " << detId.chamber() << ", " << detId.layer();
        } break;
        case MuonSubdetId::RPC: {
          RPCDetId detId(id.rawId());
          oss << "RPC chamber ";
          switch (detId.region()) {
            case 0:
              oss << "/ barrel / (wheel, station, sector, layer, subsector, roll): " << detId.ring() << ", "
                  << detId.station() << ", " << detId.sector() << ", " << detId.layer() << ", " << detId.subsector()
                  << ", " << detId.roll();
              break;
            case 1:
              oss << "/ forward endcap / (wheel, station, sector, layer, subsector, roll): " << detId.ring() << ", "
                  << detId.station() << ", " << detId.sector() << ", " << detId.layer() << ", " << detId.subsector()
                  << ", " << detId.roll();
              break;
            case -1:
              oss << "/ backward endcap / (wheel, station, sector, layer, subsector, roll): " << detId.ring() << ", "
                  << detId.station() << ", " << detId.sector() << ", " << detId.layer() << ", " << detId.subsector()
                  << ", " << detId.roll();
              break;
          }
        } break;
        case MuonSubdetId::GEM: {
          GEMDetId detId(id.rawId());
          oss << "GEM chamber (endcap, station, ring, chamber, layer): " << detId.region() << ", " << detId.station()
              << ", " << detId.ring() << ", " << detId.chamber() << ", " << detId.layer();
        } break;
        case MuonSubdetId::ME0: {
          ME0DetId detId(id.rawId());
          oss << "ME0 chamber (endcap, station, ring, chamber, layer): " << detId.region() << ", " << detId.station()
              << ", " << detId.chamber() << ", " << detId.layer();
        } break;
      }
      break;

    case DetId::Calo: {
      CaloTowerDetId detId(id.rawId());
      oss << "CaloTower (ieta, iphi): " << detId.ieta() << ", " << detId.iphi();
    } break;

    case DetId::Ecal:
      switch (id.subdetId()) {
        case EcalBarrel: {
          EBDetId detId(id);
          oss << "EcalBarrel (ieta, iphi, tower_ieta, tower_iphi): " << detId.ieta() << ", " << detId.iphi() << ", "
              << detId.tower_ieta() << ", " << detId.tower_iphi();
        } break;
        case EcalEndcap: {
          EEDetId detId(id);
          oss << "EcalEndcap (ix, iy, SuperCrystal, crystal, quadrant): " << detId.ix() << ", " << detId.iy() << ", "
              << detId.isc() << ", " << detId.ic() << ", " << detId.iquadrant();
        } break;
        case EcalPreshower:
          oss << "EcalPreshower";
          break;
        case EcalTriggerTower:
          oss << "EcalTriggerTower";
          break;
        case EcalLaserPnDiode:
          oss << "EcalLaserPnDiode";
          break;
      }
      break;

    case DetId::Hcal: {
      HcalDetId detId(id);
      switch (detId.subdet()) {
        case HcalEmpty:
          oss << "HcalEmpty ";
          break;
        case HcalBarrel:
          oss << "HcalBarrel ";
          break;
        case HcalEndcap:
          oss << "HcalEndcap ";
          break;
        case HcalOuter:
          oss << "HcalOuter ";
          break;
        case HcalForward:
          oss << "HcalForward ";
          break;
        case HcalTriggerTower:
          oss << "HcalTriggerTower ";
          break;
        case HcalOther:
          oss << "HcalOther ";
          break;
      }
      oss << "(ieta, iphi, depth):" << detId.ieta() << ", " << detId.iphi() << ", " << detId.depth();
    } break;
    default:;
  }
  return oss.str();
}

std::string DetIdInfo::info(const std::set<DetId>& idSet, const TrackerTopology* tTopo) {
  std::string text;
  for (std::set<DetId>::const_iterator id = idSet.begin(); id != idSet.end(); id++) {
    text += info(*id, tTopo);
    text += "\n";
  }
  return text;
}

std::string DetIdInfo::info(const std::vector<DetId>& idSet, const TrackerTopology* tTopo) {
  std::string text;
  for (std::vector<DetId>::const_iterator id = idSet.begin(); id != idSet.end(); id++) {
    text += info(*id, tTopo);
    text += "\n";
  }
  return text;
}