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, Forward = 6, 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

IdToInfo

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

Definition at line 137 of file FWGeometry.h.

IdToInfoItr

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

Definition at line 138 of file FWGeometry.h.

Member Enumeration Documentation

Detector

enum FWGeometry::Detector

Enumerator
Tracker
Muon
Ecal
Hcal
Calo
Forward
HGCalEE
HGCalHSi
HGCalHSc
HGCalTrigger

Definition at line 32 of file FWGeometry.h.

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

SubDetector

enum FWGeometry::SubDetector

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

Definition at line 44 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::FWGeometry

(

void

)

Definition at line 27 of file FWGeometry.cc.

: m_producerVersion(0) {}

~FWGeometry()

FWGeometry::~FWGeometry

(

void

)

Definition at line 29 of file FWGeometry.cc.

{}

Member Function Documentation

applyGlobalTag()

void FWGeometry::applyGlobalTag

(

const std::string &

gt

)

Definition at line 59 of file FWGeometry.cc.

References fwLog, download_sqlite_cfg::globalTag, fwlog::kInfo, fwlog::kWarning, loadMap(), m_fileName, alignCSCRings::r, AlCaHLTBitMon_QueryRunRegistry::string, and HLT_2024v13_cff::year.

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());
  }
}

clear()

void FWGeometry::clear ( void  )

inline

Definition at line 144 of file FWGeometry.h.

References m_idToInfo, and m_idToMatrix.

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

contains()

bool FWGeometry::contains ( unsigned int  id ) const

inline

Definition at line 140 of file FWGeometry.h.

References find(), and m_idToInfo.

Referenced by LumiList.LumiList::__contains__().

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

find()

FWGeometry::IdToInfoItr FWGeometry::find

(

unsigned int

id

)

const

Definition at line 492 of file FWGeometry.cc.

References mps_fire::end, pfDeepBoostedJetPreprocessParams_cfi::lower_bound, and m_idToInfo.

Referenced by FWSiStripClusterProxyBuilder::build(), 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);
}

findFile()

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, castor_dqm_sourceclient_file_cfg::path, Skims_PA_cff::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;
}

getCorners()

const float * FWGeometry::getCorners

(

unsigned int

id

)

const

Definition at line 439 of file FWGeometry.cc.

References find(), fwLog, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, 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(), FWECALCaloDataDetailViewBuilder::fillData(), FWECALDetailViewBuilder::fillEtaPhi(), FWHFTowerSliceSelector::findBinFromId(), FWHGTowerSliceSelector::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;
  }
}

getEveShape()

TEveGeoShape * FWGeometry::getEveShape

(

unsigned int

id

)

const

Definition at line 317 of file FWGeometry.cc.

References mps_check::array, find(), fwLog, getShape(), info(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, fwlog::kWarning, m_idToInfo, and l1trig_cff::shape.

Referenced by FWConvTrackHitsDetailView::addModules(), FWTrackHitsDetailView::addModules(), FWME0SegmentProxyBuilder::build(), FWCSCSegmentProxyBuilder::build(), FWDTSegmentProxyBuilder::buildViewType(), FWRPZViewGeometry::makeMuonGeometryRhoPhi(), FWRPZViewGeometry::makeMuonGeometryRhoZ(), FWRPZViewGeometry::showGEM(), FWRPZViewGeometry::showME0(), FWRPZViewGeometry::showMtdBarrel(), FW3DViewGeometry::showMtdBarrel(), FW3DViewGeometry::showMtdEndcap(), 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;
  }
}

getHGCScintillatorEveShape()

TEveGeoShape * FWGeometry::getHGCScintillatorEveShape

(

unsigned int

id

)

const

Definition at line 395 of file FWGeometry.cc.

References mps_check::array, find(), fwLog, info(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, fwlog::kWarning, m_idToInfo, l1trig_cff::shape, 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;
}

getHGCSiliconEveShape()

TEveGeoShape * FWGeometry::getHGCSiliconEveShape

(

unsigned int

id

)

const

Definition at line 350 of file FWGeometry.cc.

References mps_check::array, find(), fwLog, mps_fire::i, info(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, fwlog::kWarning, m_idToInfo, l1trig_cff::shape, x, and y.

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

                                                                     {
  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[6];
  Double_t y[6];
  for (unsigned int i = 0; i < 6; ++i) {
    x[i] = info.points[i * 3];
    y[i] = info.points[3 * i + 1];
  }
  geoShape->DefinePolygon(6, x, y);
  geoShape->DefineSection(0, info.points[2] - 0.0150);  // First plane at the Z position of the wafer, minus 150um
  geoShape->DefineSection(1, info.points[2] + 0.0150);  // Second plane at the Z position of the wafer, minus 150um

  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.,
                      0.,  // translation x
                      0.,  // translation y
                      0.,  // translation z
                      1.};
  // Set transformation matrix from a column-major array
  shape->SetTransMatrix(array);
  return shape;
}

getMatchedIds() [1/2]

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

Definition at line 254 of file FWGeometry.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, m_idToInfo, ALPAKA_ACCELERATOR_NAMESPACE::pixelClustering::pixelStatus::mask, and match_id().

Referenced by FW3DViewGeometry::showHGCalEE(), FW3DViewGeometry::showHGCalHSc(), FW3DViewGeometry::showHGCalHSi(), FWRPZViewGeometry::showMtdBarrel(), FW3DViewGeometry::showMtdBarrel(), FW3DViewGeometry::showMtdEndcap(), 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;
}

getMatchedIds() [2/2]

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

(

Detector

det

)

const

Definition at line 265 of file FWGeometry.cc.

References HGCalHSc, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, kDetOffset, submitPVResolutionJobs::key, and m_idToInfo.

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

    // select only the fake DetIds that have all the (u,v) bits set at 1.  This
    // is used to draw the HGCal Geometry that is wafer-based for the silicon
    // part. The Scintillators are treated on a tile-basis.
    if (det == HGCalHSc) {
      ids.push_back(it.id);
    } else {
      auto key = 0x3FF;  // 10 bits mask of 1s.
      if ((it.id | key) == it.id) {
        ids.push_back(it.id);
      }
    }
  }
  return ids;
}

getMatrix()

const TGeoMatrix * FWGeometry::getMatrix

(

unsigned int

id

)

const

Definition at line 225 of file FWGeometry.cc.

References find(), fwLog, l1ctLayer2EG_cff::id, info(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, fwlog::kWarning, m_idToInfo, m_idToMatrix, makeMuonMisalignmentScenario::matrix, and idealTransformation::rotation.

Referenced by 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];
  }
}

getParameters()

const float * FWGeometry::getParameters

(

unsigned int

id

)

const

Definition at line 450 of file FWGeometry.cc.

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

Referenced by FWHGCalMultiClusterProxyBuilder::build(), FWHGCalTriggerCellProxyBuilder::build(), FWHGCalTriggerClusterProxyBuilder::build(), FWSiStripClusterProxyBuilder::build(), FWCaloParticleProxyBuilder::build(), FWCaloClusterProxyBuilder::build(), FWTracksterHitsProxyBuilder::build(), FWTracksterLayersProxyBuilder::build(), FWHGCRecHitProxyBuilder::build(), 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;
  }
}

getProducerVersion()

int FWGeometry::getProducerVersion ( ) const

inline

Definition at line 153 of file FWGeometry.h.

References m_producerVersion.

{ return m_producerVersion; }

getShape() [1/2]

TGeoShape * FWGeometry::getShape

(

unsigned int

id

)

const

Definition at line 286 of file FWGeometry.cc.

References find(), fwLog, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, 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);
  }
}

getShape() [2/2]

TGeoShape * FWGeometry::getShape

(

const GeomDetInfo &

info

)

const

Definition at line 296 of file FWGeometry.cc.

References info().

                                                             {
  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;
}

getShapePars()

const float * FWGeometry::getShapePars

(

unsigned int

id

)

const

Definition at line 461 of file FWGeometry.cc.

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

Referenced by FWHGCalMultiClusterProxyBuilder::build(), FWHGCalTriggerCellProxyBuilder::build(), FWHGCalTriggerClusterProxyBuilder::build(), FWCaloParticleProxyBuilder::build(), FWCaloClusterProxyBuilder::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;
  }
}

getTrackerTopology()

const TrackerTopology* FWGeometry::getTrackerTopology ( ) const

inline

Definition at line 157 of file FWGeometry.h.

References m_trackerTopology.

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

{ return m_trackerTopology.get(); }

initMap()

void FWGeometry::initMap

(

const FWRecoGeom::InfoMap &

map

)

Definition at line 203 of file FWGeometry.cc.

References mps_fire::end, cms::cuda::for(), mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, dqmiolumiharvest::j, m_idToInfo, genParticles_cff::map, makeMuonMisalignmentScenario::matrix, and l1trig_cff::shape.

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];
  }
}

isEmpty()

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();
}

loadMap()

void FWGeometry::loadMap

(

const char *

fileName

)

Definition at line 103 of file FWGeometry.cc.

References FWGeometry::VersionInfo::cmsswVersion, FWGeometry::VersionInfo::extraDetectors, geometryDiff::file, findFile(), StandaloneTrackerTopology::fromTrackerParametersXMLString(), fwLog, mps_fire::i, l1ctLayer2EG_cff::id, dqmiolumiharvest::j, fwlog::kInfo, m_fileName, m_idToInfo, m_idToMatrix, m_producerVersion, m_trackerTopology, m_versionInfo, makeMuonMisalignmentScenario::matrix, AlCaHLTBitMon_ParallelJobs::p, castor_dqm_sourceclient_file_cfg::path, FWGeometry::VersionInfo::productionTag, l1trig_cff::shape, AlCaHLTBitMon_QueryRunRegistry::string, and HLT_2024v13_cff::topology.

Referenced by applyGlobalTag(), CmsShowMainBase::loadGeometry(), 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();
}

localToGlobal() [1/3]

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

Definition at line 472 of file FWGeometry.cc.

References find(), fwLog, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, fwlog::kWarning, DTRecHitClients_cfi::local, and m_idToInfo.

Referenced by FWSiStripClusterProxyBuilder::build(), FWRPZViewGeometry::estimateProjectionSizeDT(), and localToGlobal().

                                                                                                        {
  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);
  }
}

localToGlobal() [2/3]

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

Definition at line 481 of file FWGeometry.cc.

References find(), fwLog, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, 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);
  }
}

localToGlobal() [3/3]

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

Definition at line 498 of file FWGeometry.cc.

References mps_fire::i, info(), and DTRecHitClients_cfi::local.

                                                                                                                {
  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];
  }
}

mapEnd()

IdToInfoItr FWGeometry::mapEnd ( ) const

inline

Definition at line 142 of file FWGeometry.h.

References m_idToInfo.

Referenced by FWRPZViewGeometry::makeMuonGeometryRhoZ().

{ return m_idToInfo.end(); }

match_id()

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

inline

Definition at line 132 of file FWGeometry.h.

References kDetOffset, kSubdetOffset, ALPAKA_ACCELERATOR_NAMESPACE::pixelClustering::pixelStatus::mask, and EcalTangentSkim_cfg::o.

Referenced by getMatchedIds().

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

versionInfo()

const VersionInfo& FWGeometry::versionInfo ( ) const

inline

Definition at line 151 of file FWGeometry.h.

References m_versionInfo.

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

{ return m_versionInfo; }

Member Data Documentation

kDetOffset

const int FWGeometry::kDetOffset = 28

static

Definition at line 29 of file FWGeometry.h.

Referenced by getMatchedIds(), and match_id().

kSubdetOffset

const int FWGeometry::kSubdetOffset = 25

static

Definition at line 30 of file FWGeometry.h.

Referenced by match_id().

m_fileName

std::string FWGeometry::m_fileName

private

Definition at line 174 of file FWGeometry.h.

Referenced by applyGlobalTag(), and loadMap().

m_idToInfo

IdToInfo FWGeometry::m_idToInfo

private

Definition at line 166 of file FWGeometry.h.

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

m_idToMatrix

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

mutableprivate

Definition at line 164 of file FWGeometry.h.

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

m_prodTag

std::string FWGeometry::m_prodTag

private

Definition at line 168 of file FWGeometry.h.

m_producerVersion

int FWGeometry::m_producerVersion

private

Definition at line 172 of file FWGeometry.h.

Referenced by getProducerVersion(), and loadMap().

m_trackerTopology

std::unique_ptr<TrackerTopology> FWGeometry::m_trackerTopology

private

Definition at line 176 of file FWGeometry.h.

Referenced by getTrackerTopology(), and loadMap().

m_versionInfo

VersionInfo FWGeometry::m_versionInfo

private

Definition at line 170 of file FWGeometry.h.

Referenced by loadMap(), and versionInfo().