FWGeometry Class Reference

#include <FWGeometry.h>

Classes
struct	GeomDetInfo
struct	Range
class	VersionInfo

Public Types
enum	Detector {   Tracker = 1, Muon = 2, Ecal = 3, Hcal = 4,   Calo = 5, HGCalEE = 8, HGCalHSi = 9, HGCalHSc = 10,   HGCalTrigger = 11 }
typedef std::vector < FWGeometry::GeomDetInfo >	IdToInfo
typedef std::vector < FWGeometry::GeomDetInfo > ::const_iterator	IdToInfoItr
enum	SubDetector {   PixelBarrel = 1, PixelEndcap = 2, TIB = 3, TID = 4,   TOB = 5, TEC = 6, CSC = 7, DT = 8,   RPCBarrel = 9, RPCEndcap = 10, GEM = 11, ME0 = 12 }

Public Member Functions
void	applyGlobalTag (const std::string &gt)
void	clear (void)
bool	contains (unsigned int id) const
IdToInfoItr	find (unsigned int) const
	FWGeometry (void)
const float *	getCorners (unsigned int id) const
TEveGeoShape *	getEveShape (unsigned int id) const
TEveGeoShape *	getHGCScintillatorEveShape (unsigned int id) const
TEveGeoShape *	getHGCSiliconEveShape (unsigned int id) const
std::vector< unsigned int >	getMatchedIds (Detector det, SubDetector subdet) const
std::vector< unsigned int >	getMatchedIds (Detector det) const
const TGeoMatrix *	getMatrix (unsigned int id) const
const float *	getParameters (unsigned int id) const
int	getProducerVersion () const
TGeoShape *	getShape (unsigned int id) const
TGeoShape *	getShape (const GeomDetInfo &info) const
const float *	getShapePars (unsigned int id) const
const TrackerTopology *	getTrackerTopology () const
void	initMap (const FWRecoGeom::InfoMap &map)
bool	isEmpty () const
void	loadMap (const char *fileName)
void	localToGlobal (unsigned int id, const float *local, float *global, bool translatep=true) const
void	localToGlobal (unsigned int id, const float *local1, float *global1, const float *local2, float *global2, bool translatep=true) const
void	localToGlobal (const GeomDetInfo &info, const float *local, float *global, bool translatep=true) const
IdToInfoItr	mapEnd () const
bool	match_id (const GeomDetInfo &o, unsigned int mask) const
const VersionInfo &	versionInfo () const
	~FWGeometry (void)

Static Public Member Functions
static TFile *	findFile (const char *fileName)

Static Public Attributes
static const int	kDetOffset = 28
static const int	kSubdetOffset = 25

Private Attributes
std::string	m_fileName
IdToInfo	m_idToInfo
std::map< unsigned int, TGeoMatrix * >	m_idToMatrix
std::string	m_prodTag
int	m_producerVersion
std::unique_ptr< TrackerTopology >	m_trackerTopology
VersionInfo	m_versionInfo

Detailed Description

Definition at line 27 of file FWGeometry.h.

Member Typedef Documentation

typedef std::vector<FWGeometry::GeomDetInfo> FWGeometry::IdToInfo

Definition at line 136 of file FWGeometry.h.

typedef std::vector<FWGeometry::GeomDetInfo>::const_iterator FWGeometry::IdToInfoItr

Definition at line 137 of file FWGeometry.h.

Member Enumeration Documentation

enum FWGeometry::Detector

Enumerator
Tracker
Muon
Ecal
Hcal
Calo
HGCalEE
HGCalHSi
HGCalHSc
HGCalTrigger

Definition at line 32 of file FWGeometry.h.

                {
    Tracker = 1,
    Muon = 2,
    Ecal = 3,
    Hcal = 4,
    Calo = 5,
    HGCalEE = 8,
    HGCalHSi = 9,
    HGCalHSc = 10,
    HGCalTrigger = 11
  };

enum FWGeometry::SubDetector

Enumerator
PixelBarrel
PixelEndcap
TIB
TID
TOB
TEC
CSC
DT
RPCBarrel
RPCEndcap
GEM
ME0

Definition at line 43 of file FWGeometry.h.

                   {
    PixelBarrel = 1,
    PixelEndcap = 2,
    TIB = 3,
    TID = 4,
    TOB = 5,
    TEC = 6,
    CSC = 7,
    DT = 8,
    RPCBarrel = 9,
    RPCEndcap = 10,
    GEM = 11,
    ME0 = 12
  };

Constructor & Destructor Documentation

FWGeometry::FWGeometry

(

void

)

Definition at line 27 of file FWGeometry.cc.

: m_producerVersion(0) {}

FWGeometry::~FWGeometry

(

void

)

Definition at line 29 of file FWGeometry.cc.

{}

Member Function Documentation

void FWGeometry::applyGlobalTag

(

const std::string &

gt

)

Definition at line 59 of file FWGeometry.cc.

References fwLog, fwlog::kInfo, fwlog::kWarning, loadMap(), m_fileName, alignCSCRings::r, AlCaHLTBitMon_QueryRunRegistry::string, and ctpps_dqm_sourceclient-live_cfg::test.

Referenced by CmsShowMain::fileChangedSlot().

                                                          {
  const std::string fnRun2 = "cmsGeomRun2.root";
  const std::string fnRun3 = "cmsGeom2021.root";
  const std::string fnSLHC = "cmsGeom2026.root";

  TPMERegexp year_re("^[^_]+_[a-zA-Z]*20(\\d\\d)_");
  TPMERegexp run_re("^[^_]+_[a-zA-Z]*Run(\\d)_");

  TString test = globalTag.c_str();
  std::string cfn;
  if (year_re.Match(test)) {
    TString r = year_re[1];
    int year = atoi(r.Data());
    if (year < 18) {
      cfn = fnRun2;
    } else if (year < 21) {
      cfn = fnRun3;
    } else {
      cfn = fnSLHC;
    }
  } else if (run_re.Match(test)) {
    TString rn = run_re[1];
    if (rn == "1") {
      fwLog(fwlog::kWarning) << "Run1 geometry not included. Using Run2 geometry." << std::endl;
      cfn = fnRun2;
    } else if (rn == "2") {
      cfn = fnRun2;
    } else if (rn == "4") {
      cfn = fnSLHC;
    } else {
      fwLog(fwlog::kWarning) << "Detected Run" << rn << ". Using geometry scenario 2021.\n";
      cfn = fnRun3;
    }
  } else {
    fwLog(fwlog::kWarning) << "Could not guess geometry from global tag.  Using geometry scenario 2021.\n";
    cfn = fnRun3;
  }

  fwLog(fwlog::kInfo) << "Guessed geometry " << cfn << " from global tag " << globalTag << std::endl;
  if (cfn.compare(m_fileName)) {
    loadMap(cfn.c_str());
  }
}

void FWGeometry::clear ( void  )

inline

Definition at line 143 of file FWGeometry.h.

References m_idToInfo, and m_idToMatrix.

                   {
    m_idToInfo.clear();
    m_idToMatrix.clear();
  }

bool FWGeometry::contains ( unsigned int  id ) const

inline

Definition at line 139 of file FWGeometry.h.

References find(), and m_idToInfo.

Referenced by LumiList.LumiList::__contains__(), fireworks::addSiStripClusters(), FWTrackTrackingRecHitProxyBuilder::build(), FWME0SegmentProxyBuilder::build(), FWPhase2TrackerCluster1DDetProxyBuilder::build(), FWME0DigiProxyBuilder::build(), FWSiPixelClusterDetProxyBuilder::build(), FWCSCRecHitProxyBuilder::build(), FWSiStripDigiProxyBuilder::build(), FWCSCStripDigiProxyBuilder::build(), FWPhase2TrackerCluster1DProxyBuilder::build(), FWRPCDigiProxyBuilder::build(), FWCSCSegmentProxyBuilder::build(), FWCSCWireDigiProxyBuilder::build(), FWSiPixelDigiProxyBuilder::build(), FWGEMDigiProxyBuilder::build(), FWSiPixelClusterProxyBuilder::build(), FWTrajectorySeedProxyBuilder::build(), FWGEMPadDigiProxyBuilder::build(), FWME0RecHitProxyBuilder::buildViewType(), FWGEMSegmentProxyBuilder::buildViewType(), FWPSimHitProxyBuilder::buildViewType(), FWRPCRecHitProxyBuilder::buildViewType(), FWDTRecHitProxyBuilder::buildViewType(), FWGEMRecHitProxyBuilder::buildViewType(), FWDTSegmentProxyBuilder::buildViewType(), FWDTDigiProxyBuilder::buildViewType(), fireworks::pushNearbyPixelHits(), and fireworks::pushPixelHits().

{ return FWGeometry::find(id) != m_idToInfo.end(); }

FWGeometry::IdToInfoItr FWGeometry::find

(

unsigned int

id

)

const

Definition at line 514 of file FWGeometry.cc.

References SplitLinear::begin, dataset::end, cuda_std::lower_bound(), and m_idToInfo.

Referenced by FWCSCRecHitProxyBuilder::build(), FWSiStripClusterProxyBuilder::build(), FWME0RecHitProxyBuilder::buildViewType(), FWRPCRecHitProxyBuilder::buildViewType(), FWGEMRecHitProxyBuilder::buildViewType(), FWDTDigiProxyBuilder::buildViewType(), contains(), getCorners(), getEveShape(), getHGCScintillatorEveShape(), getHGCSiliconEveShape(), getMatrix(), getParameters(), getShape(), getShapePars(), localToGlobal(), and FWRPZViewGeometry::makeMuonGeometryRhoZ().

                                                            {
  FWGeometry::IdToInfoItr begin = m_idToInfo.begin();
  FWGeometry::IdToInfoItr end = m_idToInfo.end();
  return std::lower_bound(begin, end, id);
}

TFile * FWGeometry::findFile ( const char *  fileName )

static

Definition at line 38 of file FWGeometry.cc.

References MillePedeFileConverter_cfg::fileName, personalPlayback::fn, personalPlayback::fp, mps_fire::i, fed_dqm_sourceclient-live_cfg::path, mergeAndRegister::paths, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by loadMap(), and FWGeometryTableViewManager::setGeoManagerFromFile().

                                                {
  std::string searchPath = ".";

  if (gSystem->Getenv("CMSSW_SEARCH_PATH")) {
    TString paths = gSystem->Getenv("CMSSW_SEARCH_PATH");

    TObjArray* tokens = paths.Tokenize(":");
    for (int i = 0; i < tokens->GetEntries(); ++i) {
      TObjString* path = (TObjString*)tokens->At(i);
      searchPath += ":";
      searchPath += static_cast<const char*>(path->GetString());
      if (gSystem->Getenv("CMSSW_VERSION"))
        searchPath += "/Fireworks/Geometry/data/";
    }
  }

  TString fn = fileName;
  const char* fp = gSystem->FindFile(searchPath.c_str(), fn, kFileExists);
  return fp ? TFile::Open(fp) : nullptr;
}

const float * FWGeometry::getCorners

(

unsigned int

id

)

const

Definition at line 461 of file FWGeometry.cc.

References find(), fwLog, fwlog::kWarning, and m_idToInfo.

Referenced by FWPFCandidateWithHitsProxyBuilder::addHitsForCandidate(), FWCaloRecHitDigitSetProxyBuilder::build(), FWPRCaloTowerProxyBuilder::build(), FWHGCalMultiClusterProxyBuilder::build(), FWECaloParticleProxyBuilder::build(), FWHGCalTriggerCellProxyBuilder::build(), FWPCaloHitProxyBuilder::build(), FWHGCalTriggerClusterProxyBuilder::build(), FWCaloParticleProxyBuilder::build(), FWCaloClusterProxyBuilder::build(), FWTracksterHitsProxyBuilder::build(), FWPFEcalRecHitRPProxyBuilder::build(), FWPFEcalRecHitLegoProxyBuilder::build(), FWHGCRecHitProxyBuilder::build(), FWPhotonProxyBuilder::buildViewType(), FWECALCaloDataDetailViewBuilder::fillData(), FWECALDetailViewBuilder::fillEtaPhi(), FWHGTowerSliceSelector::findBinFromId(), FWHFTowerSliceSelector::findBinFromId(), fireworks::makeRhoPhiSuperCluster(), fireworks::makeRhoZSuperCluster(), FWCaloRecHitDigitSetProxyBuilder::scaleProduct(), FWPFCandidateWithHitsProxyBuilder::scaleProduct(), FW3DViewBase::showEcalBarrel(), and ValidateGeometry::validateCaloGeometry().

                                                         {
  // reco geometry points
  IdToInfoItr it = FWGeometry::find(id);
  if (it == m_idToInfo.end()) {
    fwLog(fwlog::kWarning) << "no reco geometry found for id " << id << std::endl;
    return nullptr;
  } else {
    return (*it).points;
  }
}

TEveGeoShape * FWGeometry::getEveShape

(

unsigned int

id

)

const

Definition at line 318 of file FWGeometry.cc.

References mps_check::array, find(), fwLog, getShape(), info(), fwlog::kWarning, m_idToInfo, FWGeometry::GeomDetInfo::matrix, and FWGeometry::GeomDetInfo::translation.

Referenced by FWConvTrackHitsDetailView::addModules(), FWTrackHitsDetailView::addModules(), FWME0SegmentProxyBuilder::build(), FWPhase2TrackerCluster1DDetProxyBuilder::build(), FWSiPixelClusterDetProxyBuilder::build(), FWTracksModulesProxyBuilder::build(), FWCSCSegmentProxyBuilder::build(), FWME0RecHitProxyBuilder::buildViewType(), FWGEMSegmentProxyBuilder::buildViewType(), FWRPCRecHitProxyBuilder::buildViewType(), FWGEMRecHitProxyBuilder::buildViewType(), FWDTSegmentProxyBuilder::buildViewType(), FWRPZViewGeometry::makeMuonGeometryRhoPhi(), FWRPZViewGeometry::makeMuonGeometryRhoZ(), FWRPZViewGeometry::showGEM(), FWRPZViewGeometry::showME0(), FW3DViewGeometry::showMuonBarrel(), FW3DViewGeometry::showMuonBarrelFull(), FW3DViewGeometry::showMuonEndcap(), FW3DViewGeometry::showPixelBarrel(), FWRPZViewGeometry::showPixelBarrel(), FW3DViewGeometry::showPixelEndcap(), FWRPZViewGeometry::showPixelEndcap(), FWRPZViewGeometry::showRpcEndcap(), FW3DViewGeometry::showTrackerBarrel(), FWRPZViewGeometry::showTrackerBarrel(), FW3DViewGeometry::showTrackerEndcap(), and FWRPZViewGeometry::showTrackerEndcap().

                                                           {
  IdToInfoItr it = FWGeometry::find(id);
  if (it == m_idToInfo.end()) {
    fwLog(fwlog::kWarning) << "no reco geoemtry found for id " << id << std::endl;
    return nullptr;
  } else {
    const GeomDetInfo& info = *it;
    double array[16] = {info.matrix[0],
                        info.matrix[3],
                        info.matrix[6],
                        0.,
                        info.matrix[1],
                        info.matrix[4],
                        info.matrix[7],
                        0.,
                        info.matrix[2],
                        info.matrix[5],
                        info.matrix[8],
                        0.,
                        info.translation[0],
                        info.translation[1],
                        info.translation[2],
                        1.};
    TEveGeoManagerHolder gmgr(TEveGeoShape::GetGeoMangeur());
    TEveGeoShape* shape = new TEveGeoShape(TString::Format("RecoGeom Id=%u", id));
    TGeoShape* geoShape = getShape(info);
    shape->SetShape(geoShape);
    // Set transformation matrix from a column-major array
    shape->SetTransMatrix(array);
    return shape;
  }
}

TEveGeoShape * FWGeometry::getHGCScintillatorEveShape

(

unsigned int

id

)

const

Definition at line 417 of file FWGeometry.cc.

References mps_check::array, find(), fwLog, info(), fwlog::kWarning, m_idToInfo, FWGeometry::GeomDetInfo::matrix, FWGeometry::GeomDetInfo::points, FWGeometry::GeomDetInfo::shape, FWGeometry::GeomDetInfo::translation, x, and y.

Referenced by FW3DViewGeometry::showHGCalHSc().

                                                                          {
  IdToInfoItr it = FWGeometry::find(id);
  if (it == m_idToInfo.end()) {
    fwLog(fwlog::kWarning) << "no reco geometry found for id " << id << std::endl;
    return nullptr;
  }

  GeomDetInfo info = *it;

  TEveGeoManagerHolder gmgr(TEveGeoShape::GetGeoMangeur());
  TEveGeoShape* shape = new TEveGeoShape(TString::Format("RecoGeom Id=%u", id));

  TGeoXtru* geoShape = new TGeoXtru(2);
  Double_t x[4] = {info.points[0], info.points[3], info.points[6], info.points[9]};
  Double_t y[4] = {info.points[1], info.points[4], info.points[7], info.points[10]};

  bool isNeg = info.shape[3] < 0;
  geoShape->DefinePolygon(4, x, y);
  geoShape->DefineSection(0, isNeg * info.shape[3]);
  geoShape->DefineSection(1, !isNeg * info.shape[3]);
  info.translation[2] = info.points[2];

  shape->SetShape(geoShape);
  double array[16] = {info.matrix[0],
                      info.matrix[3],
                      info.matrix[6],
                      0.,
                      info.matrix[1],
                      info.matrix[4],
                      info.matrix[7],
                      0.,
                      info.matrix[2],
                      info.matrix[5],
                      info.matrix[8],
                      0.,
                      info.translation[0],
                      info.translation[1],
                      info.translation[2],
                      1.};
  // Set transformation matrix from a column-major array
  shape->SetTransMatrix(array);
  return shape;
}

TEveGeoShape * FWGeometry::getHGCSiliconEveShape

(

unsigned int

id

)

const

Definition at line 351 of file FWGeometry.cc.

References mps_check::array, PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, validate-o2o-wbm::f, find(), fwLog, info(), fwlog::kWarning, m_idToInfo, FWGeometry::GeomDetInfo::matrix, FWGeometry::GeomDetInfo::points, mathSSE::sqrt(), FWGeometry::GeomDetInfo::translation, HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), x, and y.

Referenced by FW3DViewGeometry::showHGCalEE(), and FW3DViewGeometry::showHGCalHSi().

                                                                     {
#if 0
   const unsigned int type = (id>>26)&0x3;
   // select the middle cell of each waifer
   id &= ~0x3FF;
   id |= (type == 0) ? 0x16B : 0xE7;
#else
  float sideToSideWaferSize = 16.7441f;
  float dx = sideToSideWaferSize / 2;
  float sidey = dx / sqrt(3);
  float dy = 2 * sidey;

  int waferUint = (id >> 10) & 0xF;
  int waferVint = (id >> 15) & 0xF;
  float waferU = ((id >> 14) & 0x1) ? -sideToSideWaferSize * waferUint : sideToSideWaferSize * waferUint;
  float waferV = ((id >> 19) & 0x1) ? -sideToSideWaferSize * waferVint : sideToSideWaferSize * waferVint;

  float waferX = (-2 * waferU + waferV) / 2;
  float waferY = waferV * sqrt(3) / 2;
#endif
  IdToInfoItr it = FWGeometry::find(id);
  if (it == m_idToInfo.end()) {
    fwLog(fwlog::kWarning) << "no reco geometry found for id " << id << std::endl;
    return nullptr;
  }

  GeomDetInfo info = *it;

  TEveGeoManagerHolder gmgr(TEveGeoShape::GetGeoMangeur());
  TEveGeoShape* shape = new TEveGeoShape(TString::Format("RecoGeom Id=%u", id));

  float dz = fabs(info.points[14] - info.points[2]) * 0.5;

  info.translation[2] = (info.points[14] + info.points[2]) / 2.0f;
  info.translation[0] = waferX * ((0 < info.translation[2]) - (info.translation[2] < 0));
  info.translation[1] = waferY;

  TGeoXtru* geoShape = new TGeoXtru(2);
  Double_t x[6] = {-dx, -dx, 0.0, dx, dx, 0.0};
  Double_t y[6] = {-sidey, sidey, dy, sidey, -sidey, -dy};
  geoShape->DefinePolygon(6, x, y);
  geoShape->DefineSection(0, -dz);
  geoShape->DefineSection(1, dz);

  shape->SetShape(geoShape);
  double array[16] = {info.matrix[0],
                      info.matrix[3],
                      info.matrix[6],
                      0.,
                      info.matrix[1],
                      info.matrix[4],
                      info.matrix[7],
                      0.,
                      info.matrix[2],
                      info.matrix[5],
                      info.matrix[8],
                      0.,
                      info.translation[0],
                      info.translation[1],
                      info.translation[2],
                      1.};
  // Set transformation matrix from a column-major array
  shape->SetTransMatrix(array);
  return shape;
}

std::vector< unsigned int > FWGeometry::getMatchedIds	(	Detector	det,
		SubDetector	subdet
	)		const

Definition at line 254 of file FWGeometry.cc.

References m_idToInfo, and match_id().

Referenced by FW3DViewBase::showEcalBarrel(), FW3DViewGeometry::showHGCalEE(), FW3DViewGeometry::showHGCalHSc(), FW3DViewGeometry::showHGCalHSi(), FW3DViewGeometry::showPixelBarrel(), FWRPZViewGeometry::showPixelBarrel(), FW3DViewGeometry::showPixelEndcap(), FWRPZViewGeometry::showPixelEndcap(), FW3DViewGeometry::showTrackerBarrel(), FWRPZViewGeometry::showTrackerBarrel(), FW3DViewGeometry::showTrackerEndcap(), and FWRPZViewGeometry::showTrackerEndcap().

                                                                                        {
  std::vector<unsigned int> ids;
  unsigned int mask = (det << 4) | (subdet);
  for (IdToInfoItr it = m_idToInfo.begin(), itEnd = m_idToInfo.end(); it != itEnd; ++it) {
    if (FWGeometry::match_id(*it, mask))
      ids.push_back((*it).id);
  }

  return ids;
}

std::vector< unsigned int > FWGeometry::getMatchedIds

(

Detector

det

)

const

Definition at line 265 of file FWGeometry.cc.

References HGCalHSc, kDetOffset, submitPVResolutionJobs::key, and m_idToInfo.

                                                                    {
  std::vector<unsigned int> ids;
  std::set<unsigned int> cached_wafers;
  for (const auto& it : m_idToInfo) {
    if (((it.id >> kDetOffset) & 0xF) != det)
      continue;

    // select only the first cell of each wafer
    if (det != HGCalHSc) {
      auto key = ~0x3FF & it.id;
      if (cached_wafers.find(key) == cached_wafers.end()) {
        cached_wafers.insert(key);
        ids.push_back(it.id);
      }
      continue;
    }

    ids.push_back(it.id);
  }
  return ids;
}

const TGeoMatrix * FWGeometry::getMatrix

(

unsigned int

id

)

const

Definition at line 225 of file FWGeometry.cc.

References find(), fwLog, gpuClustering::id, info(), fwlog::kWarning, m_idToInfo, m_idToMatrix, FWGeometry::GeomDetInfo::matrix, makeMuonMisalignmentScenario::matrix, idealTransformation::rotation, and FWGeometry::GeomDetInfo::translation.

Referenced by FWTrackResidualDetailView::getSignedResidual(), CSCGeometryValidate::validateCSCChamberGeometry(), ValidateGeometry::validateCSChamberGeometry(), ValidateGeometry::validateCSCLayerGeometry(), DTGeometryValidate::validateDTChamberGeometry(), ValidateGeometry::validateDTChamberGeometry(), DTGeometryValidate::validateDTLayerGeometry(), ValidateGeometry::validateDTLayerGeometry(), GEMGeometryValidate::validateGEMChamberGeometry(), ME0GeometryValidate::validateME0ChamberGeometry(), ME0GeometryValidate::validateME0EtaPartitionGeometry2(), RPCGeometryValidate::validateRPCChamberGeometry(), ValidateGeometry::validateRPCGeometry(), and ValidateGeometry::validateTrackerGeometry().

                                                             {
  std::map<unsigned int, TGeoMatrix*>::iterator mit = m_idToMatrix.find(id);
  if (mit != m_idToMatrix.end())
    return mit->second;

  IdToInfoItr it = FWGeometry::find(id);
  if (it == m_idToInfo.end()) {
    fwLog(fwlog::kWarning) << "no reco geometry found for id " << id << std::endl;
    return nullptr;
  } else {
    const GeomDetInfo& info = *it;
    TGeoTranslation trans(info.translation[0], info.translation[1], info.translation[2]);
    TGeoRotation rotation;
    const Double_t matrix[9] = {info.matrix[0],
                                info.matrix[1],
                                info.matrix[2],
                                info.matrix[3],
                                info.matrix[4],
                                info.matrix[5],
                                info.matrix[6],
                                info.matrix[7],
                                info.matrix[8]};
    rotation.SetMatrix(matrix);

    m_idToMatrix[id] = new TGeoCombiTrans(trans, rotation);
    return m_idToMatrix[id];
  }
}

const float * FWGeometry::getParameters

(

unsigned int

id

)

const

Definition at line 472 of file FWGeometry.cc.

References find(), fwLog, fwlog::kWarning, and m_idToInfo.

Referenced by fireworks::addSiStripClusters(), FWHGCalMultiClusterProxyBuilder::build(), FWHGCalTriggerCellProxyBuilder::build(), FWHGCalTriggerClusterProxyBuilder::build(), FWSiStripClusterProxyBuilder::build(), FWCaloParticleProxyBuilder::build(), FWCSCStripDigiProxyBuilder::build(), FWSiStripDigiProxyBuilder::build(), FWRPCDigiProxyBuilder::build(), FWPhase2TrackerCluster1DProxyBuilder::build(), FWCaloClusterProxyBuilder::build(), FWSiPixelDigiProxyBuilder::build(), FWGEMDigiProxyBuilder::build(), FWSiPixelClusterProxyBuilder::build(), FWTracksterHitsProxyBuilder::build(), FWTrajectorySeedProxyBuilder::build(), FWTracksterLayersProxyBuilder::build(), FWHGCRecHitProxyBuilder::build(), FWGEMPadDigiProxyBuilder::build(), FWDTDigiProxyBuilder::buildViewType(), fireworks::pushNearbyPixelHits(), fireworks::pushPixelHits(), FW3DViewGeometry::showHGCalEE(), FW3DViewGeometry::showHGCalHSc(), FW3DViewGeometry::showHGCalHSi(), CSCGeometryValidate::validateCSCLayerGeometry(), ValidateGeometry::validateCSCLayerGeometry(), DTGeometryValidate::validateDTLayerGeometry(), ValidateGeometry::validateDTLayerGeometry(), GEMGeometryValidate::validateGEMEtaPartitionGeometry(), ME0GeometryValidate::validateME0EtaPartitionGeometry(), ValidateGeometry::validatePixelTopology(), ValidateGeometry::validateRPCGeometry(), RPCGeometryValidate::validateRPCStripsGeometry(), and ValidateGeometry::validateStripTopology().

                                                            {
  // reco geometry parameters
  IdToInfoItr it = FWGeometry::find(id);
  if (it == m_idToInfo.end()) {
    fwLog(fwlog::kWarning) << "no reco geometry found for id " << id << std::endl;
    return nullptr;
  } else {
    return (*it).parameters;
  }
}

int FWGeometry::getProducerVersion ( ) const

inline

Definition at line 152 of file FWGeometry.h.

References m_producerVersion.

{ return m_producerVersion; }

TGeoShape * FWGeometry::getShape

(

unsigned int

id

)

const

Definition at line 287 of file FWGeometry.cc.

References find(), fwLog, fwlog::kWarning, and m_idToInfo.

Referenced by FWSiStripClusterProxyBuilder::build(), and getEveShape().

                                                     {
  IdToInfoItr it = FWGeometry::find(id);
  if (it == m_idToInfo.end()) {
    fwLog(fwlog::kWarning) << "no reco geoemtry found for id " << id << std::endl;
    return nullptr;
  } else {
    return getShape(*it);
  }
}

TGeoShape * FWGeometry::getShape

(

const GeomDetInfo &

info

)

const

Definition at line 297 of file FWGeometry.cc.

References FWGeometry::GeomDetInfo::shape.

                                                             {
  TEveGeoManagerHolder gmgr(TEveGeoShape::GetGeoMangeur());
  TGeoShape* geoShape = nullptr;
  if (info.shape[0] == 1) {
    geoShape = new TGeoTrap(info.shape[3],  //dz
                            0,              //theta
                            0,              //phi
                            info.shape[4],  //dy1
                            info.shape[1],  //dx1
                            info.shape[2],  //dx2
                            0,              //alpha1
                            info.shape[4],  //dy2
                            info.shape[1],  //dx3
                            info.shape[2],  //dx4
                            0);             //alpha2
  } else
    geoShape = new TGeoBBox(info.shape[1], info.shape[2], info.shape[3]);

  return geoShape;
}

const float * FWGeometry::getShapePars

(

unsigned int

id

)

const

Definition at line 483 of file FWGeometry.cc.

References find(), fwLog, fwlog::kWarning, and m_idToInfo.

Referenced by FWHGCalMultiClusterProxyBuilder::build(), FWHGCalTriggerCellProxyBuilder::build(), FWHGCalTriggerClusterProxyBuilder::build(), FWCaloParticleProxyBuilder::build(), FWCSCStripDigiProxyBuilder::build(), FWPhase2TrackerCluster1DProxyBuilder::build(), FWCaloClusterProxyBuilder::build(), FWCSCWireDigiProxyBuilder::build(), FWTracksterHitsProxyBuilder::build(), FWHGCRecHitProxyBuilder::build(), CSCGeometryValidate::validateCSCChamberGeometry(), ValidateGeometry::validateCSChamberGeometry(), ValidateGeometry::validateCSCLayerGeometry(), DTGeometryValidate::validateDTChamberGeometry(), ValidateGeometry::validateDTChamberGeometry(), DTGeometryValidate::validateDTLayerGeometry(), ValidateGeometry::validateDTLayerGeometry(), GEMGeometryValidate::validateGEMChamberGeometry(), ME0GeometryValidate::validateME0ChamberGeometry(), ME0GeometryValidate::validateME0EtaPartitionGeometry2(), RPCGeometryValidate::validateRPCChamberGeometry(), ValidateGeometry::validateRPCGeometry(), and ValidateGeometry::validateTrackerGeometry().

                                                           {
  // reco geometry parameters
  IdToInfoItr it = FWGeometry::find(id);
  if (it == m_idToInfo.end()) {
    fwLog(fwlog::kWarning) << "no reco geometry found for id " << id << std::endl;
    return nullptr;
  } else {
    return (*it).shape;
  }
}

const TrackerTopology* FWGeometry::getTrackerTopology ( ) const

inline

Definition at line 156 of file FWGeometry.h.

References m_trackerTopology.

Referenced by FWTrackHitsDetailView::addModules(), fireworks::info(), FW3DViewGeometry::showPixelBarrel(), and FW3DViewGeometry::showPixelEndcap().

{ return m_trackerTopology.get(); }

void FWGeometry::initMap

(

const FWRecoGeom::InfoMap &

map

)

Definition at line 203 of file FWGeometry.cc.

References SplitLinear::begin, dataset::end, for(), mps_fire::i, dqmiolumiharvest::j, m_idToInfo, and makeMuonMisalignmentScenario::matrix.

Referenced by FWFFLooper::beginRun().

                                                     {
  FWRecoGeom::InfoMapItr begin = map.begin();
  FWRecoGeom::InfoMapItr end = map.end();
  unsigned int mapSize = map.size();
  if (m_idToInfo.size() != mapSize)
    m_idToInfo.resize(mapSize);
  unsigned int i = 0;
  for (FWRecoGeom::InfoMapItr it = begin; it != end; ++it, ++i) {
    m_idToInfo[i].id = it->id;
    for (unsigned int j = 0; j < 24; ++j)
      m_idToInfo[i].points[j] = it->points[j];
    for (unsigned int j = 0; j < 9; ++j)
      m_idToInfo[i].parameters[j] = it->topology[j];
    for (unsigned int j = 0; j < 5; ++j)
      m_idToInfo[i].shape[j] = it->shape[j];
    for (unsigned int j = 0; j < 3; ++j)
      m_idToInfo[i].translation[j] = it->translation[j];
    for (unsigned int j = 0; j < 9; ++j)
      m_idToInfo[i].matrix[j] = it->matrix[j];
  }
}

bool FWGeometry::isEmpty

(

void

)

const

Definition at line 31 of file FWGeometry.cc.

References m_idToInfo.

Referenced by plotting.Plot::clone(), and FWRPZViewGeometry::initStdGeoElements().

                               {
  // AMT this is a check if geomtery is not loaded
  // e.g. cmsShow starts with no data file and without given explicit argument ( --geometry-file option )

  return m_idToInfo.empty();
}

void FWGeometry::loadMap

(

const char *

fileName

)

Definition at line 103 of file FWGeometry.cc.

References FWGeometry::VersionInfo::cmsswVersion, FWGeometry::VersionInfo::extraDetectors, mergeVDriftHistosByStation::file, findFile(), StandaloneTrackerTopology::fromTrackerParametersXMLString(), fwLog, mps_fire::i, gpuClustering::id, dqmiolumiharvest::j, fwlog::kInfo, m_fileName, m_idToInfo, m_idToMatrix, m_producerVersion, m_trackerTopology, m_versionInfo, makeMuonMisalignmentScenario::matrix, AlCaHLTBitMon_ParallelJobs::p, fed_dqm_sourceclient-live_cfg::path, FWGeometry::VersionInfo::productionTag, AlCaHLTBitMon_QueryRunRegistry::string, and MainPageGenerator::tree.

Referenced by applyGlobalTag(), CSCGeometryValidate::CSCGeometryValidate(), DTGeometryValidate::DTGeometryValidate(), GEMGeometryValidate::GEMGeometryValidate(), CmsShowMainBase::loadGeometry(), ME0GeometryValidate::ME0GeometryValidate(), RPCGeometryValidate::RPCGeometryValidate(), and ValidateGeometry::ValidateGeometry().

                                              {
  TFile* file = findFile(iFileName);
  if (!file) {
    throw std::runtime_error("ERROR: failed to find geometry file. Initialization failed.");
    return;
  }
  m_fileName = iFileName;
  TTree* tree = static_cast<TTree*>(file->Get("idToGeo"));
  if (!tree) {
    throw std::runtime_error("ERROR: cannot find detector id map in the file. Initialization failed.");
    return;
  }

  unsigned int id;
  Float_t points[24];
  Float_t topology[9];
  Float_t shape[5];
  Float_t translation[3];
  Float_t matrix[9];
  bool loadPoints = tree->GetBranch("points") != nullptr;
  bool loadParameters = tree->GetBranch("topology") != nullptr;
  bool loadShape = tree->GetBranch("shape") != nullptr;
  bool loadTranslation = tree->GetBranch("translation") != nullptr;
  bool loadMatrix = tree->GetBranch("matrix") != nullptr;
  tree->SetBranchAddress("id", &id);
  if (loadPoints)
    tree->SetBranchAddress("points", &points);
  if (loadParameters)
    tree->SetBranchAddress("topology", &topology);
  if (loadShape)
    tree->SetBranchAddress("shape", &shape);
  if (loadTranslation)
    tree->SetBranchAddress("translation", &translation);
  if (loadMatrix)
    tree->SetBranchAddress("matrix", &matrix);

  // reset previous values
  m_idToInfo.clear();
  for (const auto& p : m_idToMatrix)
    delete p.second;
  m_trackerTopology.reset();

  unsigned int treeSize = tree->GetEntries();
  if (m_idToInfo.size() != treeSize)
    m_idToInfo.resize(treeSize);
  for (unsigned int i = 0; i < treeSize; ++i) {
    tree->GetEntry(i);

    m_idToInfo[i].id = id;
    if (loadPoints) {
      for (unsigned int j = 0; j < 24; ++j)
        m_idToInfo[i].points[j] = points[j];
    }
    if (loadParameters) {
      for (unsigned int j = 0; j < 9; ++j)
        m_idToInfo[i].parameters[j] = topology[j];
    }
    if (loadShape) {
      for (unsigned int j = 0; j < 5; ++j)
        m_idToInfo[i].shape[j] = shape[j];
    }
    if (loadTranslation) {
      for (unsigned int j = 0; j < 3; ++j)
        m_idToInfo[i].translation[j] = translation[j];
    }
    if (loadMatrix) {
      for (unsigned int j = 0; j < 9; ++j)
        m_idToInfo[i].matrix[j] = matrix[j];
    }
  }

  m_versionInfo.productionTag = static_cast<TNamed*>(file->Get("tag"));
  m_versionInfo.cmsswVersion = static_cast<TNamed*>(file->Get("CMSSW_VERSION"));
  m_versionInfo.extraDetectors = static_cast<TObjArray*>(file->Get("ExtraDetectors"));

  TString path = file->GetPath();
  if (path.EndsWith(":/"))
    path.Resize(path.Length() - 2);

  if (m_versionInfo.productionTag)
    fwLog(fwlog::kInfo) << Form(
        "Load %s %s from %s\n", tree->GetName(), m_versionInfo.productionTag->GetTitle(), path.Data());
  else
    fwLog(fwlog::kInfo) << Form("Load %s from %s\n", tree->GetName(), path.Data());

  TNamed* producerInfo = static_cast<TNamed*>(file->Get("PRODUCER_VERSION"));
  if (producerInfo) {
    m_producerVersion = atoi(producerInfo->GetTitle());
  }

  TNamed* ttopology = static_cast<TNamed*>(file->Get("TrackerTopology"));
  if (ttopology) {
    std::string xml = ttopology->GetTitle();
    m_trackerTopology =
        std::make_unique<TrackerTopology>(StandaloneTrackerTopology::fromTrackerParametersXMLString(xml));
  }

  file->Close();
}

void FWGeometry::localToGlobal	(	unsigned int	id,
		const float *	local,
		float *	global,
		bool	translatep = true
	)		const

Definition at line 494 of file FWGeometry.cc.

References find(), fwLog, fwlog::kWarning, and m_idToInfo.

Referenced by fireworks::addSiStripClusters(), FWTrackTrackingRecHitProxyBuilder::build(), FWME0SegmentProxyBuilder::build(), FWME0DigiProxyBuilder::build(), FWCSCRecHitProxyBuilder::build(), FWSiStripClusterProxyBuilder::build(), FWSiStripDigiProxyBuilder::build(), FWCSCStripDigiProxyBuilder::build(), FWRPCDigiProxyBuilder::build(), FWPhase2TrackerCluster1DProxyBuilder::build(), FWCSCSegmentProxyBuilder::build(), FWSiPixelDigiProxyBuilder::build(), FWCSCWireDigiProxyBuilder::build(), FWGEMDigiProxyBuilder::build(), FWSiPixelClusterProxyBuilder::build(), FWTrackingParticleProxyBuilderFullFramework::build(), FWTrajectorySeedProxyBuilder::build(), FWGEMPadDigiProxyBuilder::build(), FWME0RecHitProxyBuilder::buildViewType(), FWGEMSegmentProxyBuilder::buildViewType(), FWPSimHitProxyBuilder::buildViewType(), FWRPCRecHitProxyBuilder::buildViewType(), FWDTRecHitProxyBuilder::buildViewType(), FWGEMRecHitProxyBuilder::buildViewType(), FWDTSegmentProxyBuilder::buildViewType(), FWDTDigiProxyBuilder::buildViewType(), FWRPZViewGeometry::estimateProjectionSizeDT(), localToGlobal(), and fireworks::pushPixelCluster().

                                                                                                        {
  IdToInfoItr it = FWGeometry::find(id);
  if (it == m_idToInfo.end()) {
    fwLog(fwlog::kWarning) << "no reco geometry found for id " << id << std::endl;
  } else {
    localToGlobal(*it, local, global, translatep);
  }
}

void FWGeometry::localToGlobal	(	unsigned int	id,
		const float *	local1,
		float *	global1,
		const float *	local2,
		float *	global2,
		bool	translatep = true
	)		const

Definition at line 503 of file FWGeometry.cc.

References find(), fwLog, fwlog::kWarning, localToGlobal(), and m_idToInfo.

                                                                                                                      {
  IdToInfoItr it = FWGeometry::find(id);
  if (it == m_idToInfo.end()) {
    fwLog(fwlog::kWarning) << "no reco geometry found for id " << id << std::endl;
  } else {
    localToGlobal(*it, local1, global1, translatep);
    localToGlobal(*it, local2, global2, translatep);
  }
}

void FWGeometry::localToGlobal	(	const GeomDetInfo &	info,
		const float *	local,
		float *	global,
		bool	translatep = true
	)		const

Definition at line 520 of file FWGeometry.cc.

References mps_fire::i, FWGeometry::GeomDetInfo::matrix, and FWGeometry::GeomDetInfo::translation.

                                                                                                                {
  for (int i = 0; i < 3; ++i) {
    global[i] = translatep ? info.translation[i] : 0;
    global[i] += local[0] * info.matrix[3 * i] + local[1] * info.matrix[3 * i + 1] + local[2] * info.matrix[3 * i + 2];
  }
}

IdToInfoItr FWGeometry::mapEnd ( ) const

inline

Definition at line 141 of file FWGeometry.h.

References m_idToInfo.

Referenced by FWRPZViewGeometry::makeMuonGeometryRhoZ().

{ return m_idToInfo.end(); }

bool FWGeometry::match_id ( const GeomDetInfo &  o, unsigned int  mask  ) const

inline

Definition at line 131 of file FWGeometry.h.

References FWGeometry::GeomDetInfo::id, kDetOffset, and kSubdetOffset.

Referenced by getMatchedIds().

                                                               {
    unsigned int id = o.id;
    return (((((id >> kDetOffset) & 0xF) << 4) | ((id >> kSubdetOffset) & 0x7)) == mask);
  }

const VersionInfo& FWGeometry::versionInfo ( ) const

inline

Definition at line 150 of file FWGeometry.h.

References m_versionInfo.

Referenced by FWRPZView::populateController(), FW3DViewGeometry::showMuonEndcap(), and FWRPZViewGeometry::showRpcEndcap().

{ return m_versionInfo; }

Member Data Documentation

const int FWGeometry::kDetOffset = 28

static

Definition at line 29 of file FWGeometry.h.

Referenced by getMatchedIds(), and match_id().

const int FWGeometry::kSubdetOffset = 25

static

Definition at line 30 of file FWGeometry.h.

Referenced by match_id().

std::string FWGeometry::m_fileName

private

Definition at line 173 of file FWGeometry.h.

Referenced by applyGlobalTag(), and loadMap().

IdToInfo FWGeometry::m_idToInfo

private

Definition at line 165 of file FWGeometry.h.

Referenced by clear(), contains(), find(), getCorners(), getEveShape(), getHGCScintillatorEveShape(), getHGCSiliconEveShape(), getMatchedIds(), getMatrix(), getParameters(), getShape(), getShapePars(), initMap(), isEmpty(), loadMap(), localToGlobal(), and mapEnd().

std::map<unsigned int, TGeoMatrix*> FWGeometry::m_idToMatrix

mutableprivate

Definition at line 163 of file FWGeometry.h.

Referenced by clear(), getMatrix(), and loadMap().

std::string FWGeometry::m_prodTag

private

Definition at line 167 of file FWGeometry.h.

int FWGeometry::m_producerVersion

private

Definition at line 171 of file FWGeometry.h.

Referenced by getProducerVersion(), and loadMap().

std::unique_ptr<TrackerTopology> FWGeometry::m_trackerTopology

private

Definition at line 175 of file FWGeometry.h.

Referenced by getTrackerTopology(), and loadMap().

VersionInfo FWGeometry::m_versionInfo

private

Definition at line 169 of file FWGeometry.h.

Referenced by loadMap(), and versionInfo().