CMS 3D CMS Logo
List of all members | Classes | Public Types | Public Member Functions | Static Public Member Functions | Static Public Attributes | Private Attributes
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) const
 
std::vector< unsigned int > getMatchedIds (Detector det, SubDetector subdet) const
 
const TGeoMatrix * getMatrix (unsigned int id) const
 
const float * getParameters (unsigned int id) const
 
int getProducerVersion () const
 
TGeoShape * getShape (const GeomDetInfo &info) const
 
TGeoShape * getShape (unsigned int id) 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 (const GeomDetInfo &info, const float *local, float *global, bool translatep=true) const
 
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
 
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 136 of file FWGeometry.h.

◆ IdToInfoItr

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

Definition at line 137 of file FWGeometry.h.

Member Enumeration Documentation

◆ Detector

enum FWGeometry::Detector
Enumerator
Tracker 
Muon 
Ecal 
Hcal 
Calo 
HGCalEE 
HGCalHSi 
HGCalHSc 
HGCalTrigger 

Definition at line 32 of file FWGeometry.h.

32  {
33  Tracker = 1,
34  Muon = 2,
35  Ecal = 3,
36  Hcal = 4,
37  Calo = 5,
38  HGCalEE = 8,
39  HGCalHSi = 9,
40  HGCalHSc = 10,
41  HGCalTrigger = 11
42  };

◆ SubDetector

enum FWGeometry::SubDetector
Enumerator
PixelBarrel 
PixelEndcap 
TIB 
TID 
TOB 
TEC 
CSC 
DT 
RPCBarrel 
RPCEndcap 
GEM 
ME0 

Definition at line 43 of file FWGeometry.h.

43  {
44  PixelBarrel = 1,
45  PixelEndcap = 2,
46  TIB = 3,
47  TID = 4,
48  TOB = 5,
49  TEC = 6,
50  CSC = 7,
51  DT = 8,
52  RPCBarrel = 9,
53  RPCEndcap = 10,
54  GEM = 11,
55  ME0 = 12
56  };

Constructor & Destructor Documentation

◆ FWGeometry()

FWGeometry::FWGeometry ( void  )

Definition at line 24 of file FWGeometry.cc.

24 : m_producerVersion(0) {}

◆ ~FWGeometry()

FWGeometry::~FWGeometry ( void  )

Definition at line 26 of file FWGeometry.cc.

26 {}

Member Function Documentation

◆ applyGlobalTag()

void FWGeometry::applyGlobalTag ( const std::string &  gt)

Definition at line 56 of file FWGeometry.cc.

56  {
57  const std::string fnRun2 = "cmsGeomRun2.root";
58  const std::string fnRun3 = "cmsGeom2021.root";
59  const std::string fnSLHC = "cmsGeom2026.root";
60 
61  TPMERegexp year_re("^[^_]+_[a-zA-Z]*20(\\d\\d)_");
62  TPMERegexp run_re("^[^_]+_[a-zA-Z]*Run(\\d)_");
63 
64  TString test = globalTag.c_str();
65  std::string cfn;
66  if (year_re.Match(test)) {
67  TString r = year_re[1];
68  int year = atoi(r.Data());
69  if (year < 18) {
70  cfn = fnRun2;
71  } else if (year < 21) {
72  cfn = fnRun3;
73  } else {
74  cfn = fnSLHC;
75  }
76  } else if (run_re.Match(test)) {
77  TString rn = run_re[1];
78  if (rn == "1") {
79  fwLog(fwlog::kWarning) << "Run1 geometry not included. Using Run2 geometry." << std::endl;
80  cfn = fnRun2;
81  } else if (rn == "2") {
82  cfn = fnRun2;
83  } else if (rn == "4") {
84  cfn = fnSLHC;
85  } else {
86  fwLog(fwlog::kWarning) << "Detected Run" << rn << ". Using geometry scenario 2021.\n";
87  cfn = fnRun3;
88  }
89  } else {
90  fwLog(fwlog::kWarning) << "Could not guess geometry from global tag. Using geometry scenario 2021.\n";
91  cfn = fnRun3;
92  }
93 
94  fwLog(fwlog::kInfo) << "Guessed geometry " << cfn << " from global tag " << globalTag << std::endl;
95  if (cfn.compare(m_fileName)) {
96  loadMap(cfn.c_str());
97  }
98 }

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

Referenced by CmsShowMain::fileChangedSlot().

◆ clear()

void FWGeometry::clear ( void  )
inline

Definition at line 143 of file FWGeometry.h.

143  {
144  m_idToInfo.clear();
145  m_idToMatrix.clear();
146  }

References m_idToInfo, and m_idToMatrix.

◆ contains()

bool FWGeometry::contains ( unsigned int  id) const
inline

Definition at line 139 of file FWGeometry.h.

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

References find(), and m_idToInfo.

Referenced by LumiList.LumiList::__contains__().

◆ find()

FWGeometry::IdToInfoItr FWGeometry::find ( unsigned int  id) const

Definition at line 517 of file FWGeometry.cc.

517  {
518  FWGeometry::IdToInfoItr begin = m_idToInfo.begin();
519  FWGeometry::IdToInfoItr end = m_idToInfo.end();
520  return std::lower_bound(begin, end, id);
521 }

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

Referenced by FWSiStripClusterProxyBuilder::build(), contains(), getCorners(), getEveShape(), getHGCScintillatorEveShape(), getHGCSiliconEveShape(), getMatrix(), getParameters(), getShape(), getShapePars(), localToGlobal(), and FWRPZViewGeometry::makeMuonGeometryRhoZ().

◆ findFile()

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

Definition at line 35 of file FWGeometry.cc.

35  {
36  std::string searchPath = ".";
37 
38  if (gSystem->Getenv("CMSSW_SEARCH_PATH")) {
39  TString paths = gSystem->Getenv("CMSSW_SEARCH_PATH");
40 
41  TObjArray* tokens = paths.Tokenize(":");
42  for (int i = 0; i < tokens->GetEntries(); ++i) {
43  TObjString* path = (TObjString*)tokens->At(i);
44  searchPath += ":";
45  searchPath += static_cast<const char*>(path->GetString());
46  if (gSystem->Getenv("CMSSW_VERSION"))
47  searchPath += "/Fireworks/Geometry/data/";
48  }
49  }
50 
51  TString fn = fileName;
52  const char* fp = gSystem->FindFile(searchPath.c_str(), fn, kFileExists);
53  return fp ? TFile::Open(fp) : nullptr;
54 }

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().

◆ getCorners()

const float * FWGeometry::getCorners ( unsigned int  id) const

Definition at line 464 of file FWGeometry.cc.

464  {
465  // reco geometry points
466  IdToInfoItr it = FWGeometry::find(id);
467  if (it == m_idToInfo.end()) {
468  fwLog(fwlog::kWarning) << "no reco geometry found for id " << id << std::endl;
469  return nullptr;
470  } else {
471  return (*it).points;
472  }
473 }

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(), FWPFEcalRecHitRPProxyBuilder::build(), FWPFEcalRecHitLegoProxyBuilder::build(), FWHGCRecHitProxyBuilder::build(), FWECALCaloDataDetailViewBuilder::fillData(), FWECALDetailViewBuilder::fillEtaPhi(), FWHGTowerSliceSelector::findBinFromId(), FWHFTowerSliceSelector::findBinFromId(), fireworks::makeRhoPhiSuperCluster(), fireworks::makeRhoZSuperCluster(), FWCaloRecHitDigitSetProxyBuilder::scaleProduct(), FWPFCandidateWithHitsProxyBuilder::scaleProduct(), FW3DViewBase::showEcalBarrel(), and ValidateGeometry::validateCaloGeometry().

◆ getEveShape()

TEveGeoShape * FWGeometry::getEveShape ( unsigned int  id) const

Definition at line 321 of file FWGeometry.cc.

321  {
322  IdToInfoItr it = FWGeometry::find(id);
323  if (it == m_idToInfo.end()) {
324  fwLog(fwlog::kWarning) << "no reco geoemtry found for id " << id << std::endl;
325  return nullptr;
326  } else {
327  const GeomDetInfo& info = *it;
328  double array[16] = {info.matrix[0],
329  info.matrix[3],
330  info.matrix[6],
331  0.,
332  info.matrix[1],
333  info.matrix[4],
334  info.matrix[7],
335  0.,
336  info.matrix[2],
337  info.matrix[5],
338  info.matrix[8],
339  0.,
340  info.translation[0],
341  info.translation[1],
342  info.translation[2],
343  1.};
344  TEveGeoManagerHolder gmgr(TEveGeoShape::GetGeoMangeur());
345  TEveGeoShape* shape = new TEveGeoShape(TString::Format("RecoGeom Id=%u", id));
346  TGeoShape* geoShape = getShape(info);
347  shape->SetShape(geoShape);
348  // Set transformation matrix from a column-major array
349  shape->SetTransMatrix(array);
350  return shape;
351  }
352 }

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

Referenced by FWTrackHitsDetailView::addModules(), FWConvTrackHitsDetailView::addModules(), FWME0SegmentProxyBuilder::build(), FWCSCSegmentProxyBuilder::build(), 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().

◆ getHGCScintillatorEveShape()

TEveGeoShape * FWGeometry::getHGCScintillatorEveShape ( unsigned int  id) const

Definition at line 420 of file FWGeometry.cc.

420  {
421  IdToInfoItr it = FWGeometry::find(id);
422  if (it == m_idToInfo.end()) {
423  fwLog(fwlog::kWarning) << "no reco geometry found for id " << id << std::endl;
424  return nullptr;
425  }
426 
427  GeomDetInfo info = *it;
428 
429  TEveGeoManagerHolder gmgr(TEveGeoShape::GetGeoMangeur());
430  TEveGeoShape* shape = new TEveGeoShape(TString::Format("RecoGeom Id=%u", id));
431 
432  TGeoXtru* geoShape = new TGeoXtru(2);
433  Double_t x[4] = {info.points[0], info.points[3], info.points[6], info.points[9]};
434  Double_t y[4] = {info.points[1], info.points[4], info.points[7], info.points[10]};
435 
436  bool isNeg = info.shape[3] < 0;
437  geoShape->DefinePolygon(4, x, y);
438  geoShape->DefineSection(0, isNeg * info.shape[3]);
439  geoShape->DefineSection(1, !isNeg * info.shape[3]);
440  info.translation[2] = info.points[2];
441 
442  shape->SetShape(geoShape);
443  double array[16] = {info.matrix[0],
444  info.matrix[3],
445  info.matrix[6],
446  0.,
447  info.matrix[1],
448  info.matrix[4],
449  info.matrix[7],
450  0.,
451  info.matrix[2],
452  info.matrix[5],
453  info.matrix[8],
454  0.,
455  info.translation[0],
456  info.translation[1],
457  info.translation[2],
458  1.};
459  // Set transformation matrix from a column-major array
460  shape->SetTransMatrix(array);
461  return shape;
462 }

References mps_check::array, find(), fwLog, info(), fwlog::kWarning, m_idToInfo, x, and y.

Referenced by FW3DViewGeometry::showHGCalHSc().

◆ getHGCSiliconEveShape()

TEveGeoShape * FWGeometry::getHGCSiliconEveShape ( unsigned int  id) const

Definition at line 354 of file FWGeometry.cc.

354  {
355 #if 0
356  const unsigned int type = (id>>26)&0x3;
357  // select the middle cell of each waifer
358  id &= ~0x3FF;
359  id |= (type == 0) ? 0x16B : 0xE7;
360 #else
361  float sideToSideWaferSize = 16.7441f;
362  float dx = sideToSideWaferSize / 2;
363  float sidey = dx / sqrt(3);
364  float dy = 2 * sidey;
365 
366  int waferUint = (id >> 10) & 0xF;
367  int waferVint = (id >> 15) & 0xF;
368  float waferU = ((id >> 14) & 0x1) ? -sideToSideWaferSize * waferUint : sideToSideWaferSize * waferUint;
369  float waferV = ((id >> 19) & 0x1) ? -sideToSideWaferSize * waferVint : sideToSideWaferSize * waferVint;
370 
371  float waferX = (-2 * waferU + waferV) / 2;
372  float waferY = waferV * sqrt(3) / 2;
373 #endif
374  IdToInfoItr it = FWGeometry::find(id);
375  if (it == m_idToInfo.end()) {
376  fwLog(fwlog::kWarning) << "no reco geometry found for id " << id << std::endl;
377  return nullptr;
378  }
379 
380  GeomDetInfo info = *it;
381 
382  TEveGeoManagerHolder gmgr(TEveGeoShape::GetGeoMangeur());
383  TEveGeoShape* shape = new TEveGeoShape(TString::Format("RecoGeom Id=%u", id));
384 
385  float dz = fabs(info.points[14] - info.points[2]) * 0.5;
386 
387  info.translation[2] = (info.points[14] + info.points[2]) / 2.0f;
388  info.translation[0] = waferX * ((0 < info.translation[2]) - (info.translation[2] < 0));
389  info.translation[1] = waferY;
390 
391  TGeoXtru* geoShape = new TGeoXtru(2);
392  Double_t x[6] = {-dx, -dx, 0.0, dx, dx, 0.0};
393  Double_t y[6] = {-sidey, sidey, dy, sidey, -sidey, -dy};
394  geoShape->DefinePolygon(6, x, y);
395  geoShape->DefineSection(0, -dz);
396  geoShape->DefineSection(1, dz);
397 
398  shape->SetShape(geoShape);
399  double array[16] = {info.matrix[0],
400  info.matrix[3],
401  info.matrix[6],
402  0.,
403  info.matrix[1],
404  info.matrix[4],
405  info.matrix[7],
406  0.,
407  info.matrix[2],
408  info.matrix[5],
409  info.matrix[8],
410  0.,
411  info.translation[0],
412  info.translation[1],
413  info.translation[2],
414  1.};
415  // Set transformation matrix from a column-major array
416  shape->SetTransMatrix(array);
417  return shape;
418 }

References mps_check::array, PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, f, find(), fwLog, info(), fwlog::kWarning, m_idToInfo, mathSSE::sqrt(), x, testProducerWithPsetDescEmpty_cfi::x1, and y.

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

◆ getMatchedIds() [1/2]

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

Definition at line 262 of file FWGeometry.cc.

262  {
263  std::vector<unsigned int> ids;
264 
265  for (const auto& it : m_idToInfo) {
266  if (((it.id >> kDetOffset) & 0xF) != det)
267  continue;
268  // select only the 1st layer
269  // if(((it->id>>17)&0x1F) != 12) continue;
270 
271  // select only the first cell of each wafer
272  if (det != HGCalHSc) {
273  bool flag(false);
274 
275  for (const auto& id_it : ids) {
276  flag = (~0x3FF & it.id) == (~0x3FF & id_it);
277  if (flag)
278  break;
279  }
280 
281  if (flag)
282  continue;
283  }
284 
285  ids.push_back(it.id);
286  }
287  return ids;
288 }

References RemoveAddSevLevel::flag, HGCalHSc, kDetOffset, and m_idToInfo.

◆ getMatchedIds() [2/2]

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

Definition at line 251 of file FWGeometry.cc.

251  {
252  std::vector<unsigned int> ids;
253  unsigned int mask = (det << 4) | (subdet);
254  for (IdToInfoItr it = m_idToInfo.begin(), itEnd = m_idToInfo.end(); it != itEnd; ++it) {
255  if (FWGeometry::match_id(*it, mask))
256  ids.push_back((*it).id);
257  }
258 
259  return ids;
260 }

References m_idToInfo, and match_id().

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

◆ getMatrix()

const TGeoMatrix * FWGeometry::getMatrix ( unsigned int  id) const

Definition at line 222 of file FWGeometry.cc.

222  {
223  std::map<unsigned int, TGeoMatrix*>::iterator mit = m_idToMatrix.find(id);
224  if (mit != m_idToMatrix.end())
225  return mit->second;
226 
227  IdToInfoItr it = FWGeometry::find(id);
228  if (it == m_idToInfo.end()) {
229  fwLog(fwlog::kWarning) << "no reco geometry found for id " << id << std::endl;
230  return nullptr;
231  } else {
232  const GeomDetInfo& info = *it;
233  TGeoTranslation trans(info.translation[0], info.translation[1], info.translation[2]);
234  TGeoRotation rotation;
235  const Double_t matrix[9] = {info.matrix[0],
236  info.matrix[1],
237  info.matrix[2],
238  info.matrix[3],
239  info.matrix[4],
240  info.matrix[5],
241  info.matrix[6],
242  info.matrix[7],
243  info.matrix[8]};
244  rotation.SetMatrix(matrix);
245 
246  m_idToMatrix[id] = new TGeoCombiTrans(trans, rotation);
247  return m_idToMatrix[id];
248  }
249 }

References find(), fwLog, triggerObjects_cff::id, info(), 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().

◆ getParameters()

const float * FWGeometry::getParameters ( unsigned int  id) const

Definition at line 475 of file FWGeometry.cc.

475  {
476  // reco geometry parameters
477  IdToInfoItr it = FWGeometry::find(id);
478  if (it == m_idToInfo.end()) {
479  fwLog(fwlog::kWarning) << "no reco geometry found for id " << id << std::endl;
480  return nullptr;
481  } else {
482  return (*it).parameters;
483  }
484 }

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

Referenced by FWHGCalMultiClusterProxyBuilder::build(), FWHGCalTriggerCellProxyBuilder::build(), FWSiStripClusterProxyBuilder::build(), FWHGCalTriggerClusterProxyBuilder::build(), FWCaloParticleProxyBuilder::build(), FWCaloClusterProxyBuilder::build(), FWHGCRecHitProxyBuilder::build(), CSCGeometryValidate::validateCSCLayerGeometry(), ValidateGeometry::validateCSCLayerGeometry(), DTGeometryValidate::validateDTLayerGeometry(), ValidateGeometry::validateDTLayerGeometry(), GEMGeometryValidate::validateGEMEtaPartitionGeometry(), ME0GeometryValidate::validateME0EtaPartitionGeometry(), ValidateGeometry::validatePixelTopology(), ValidateGeometry::validateRPCGeometry(), RPCGeometryValidate::validateRPCStripsGeometry(), and ValidateGeometry::validateStripTopology().

◆ getProducerVersion()

int FWGeometry::getProducerVersion ( ) const
inline

Definition at line 152 of file FWGeometry.h.

152 { return m_producerVersion; }

References m_producerVersion.

◆ getShape() [1/2]

TGeoShape * FWGeometry::getShape ( const GeomDetInfo &  info) const

Definition at line 300 of file FWGeometry.cc.

300  {
301  TEveGeoManagerHolder gmgr(TEveGeoShape::GetGeoMangeur());
302  TGeoShape* geoShape = nullptr;
303  if (info.shape[0] == 1) {
304  geoShape = new TGeoTrap(info.shape[3], //dz
305  0, //theta
306  0, //phi
307  info.shape[4], //dy1
308  info.shape[1], //dx1
309  info.shape[2], //dx2
310  0, //alpha1
311  info.shape[4], //dy2
312  info.shape[1], //dx3
313  info.shape[2], //dx4
314  0); //alpha2
315  } else
316  geoShape = new TGeoBBox(info.shape[1], info.shape[2], info.shape[3]);
317 
318  return geoShape;
319 }

References info().

◆ getShape() [2/2]

TGeoShape * FWGeometry::getShape ( unsigned int  id) const

Definition at line 290 of file FWGeometry.cc.

290  {
291  IdToInfoItr it = FWGeometry::find(id);
292  if (it == m_idToInfo.end()) {
293  fwLog(fwlog::kWarning) << "no reco geoemtry found for id " << id << std::endl;
294  return nullptr;
295  } else {
296  return getShape(*it);
297  }
298 }

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

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

◆ getShapePars()

const float * FWGeometry::getShapePars ( unsigned int  id) const

Definition at line 486 of file FWGeometry.cc.

486  {
487  // reco geometry parameters
488  IdToInfoItr it = FWGeometry::find(id);
489  if (it == m_idToInfo.end()) {
490  fwLog(fwlog::kWarning) << "no reco geometry found for id " << id << std::endl;
491  return nullptr;
492  } else {
493  return (*it).shape;
494  }
495 }

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

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

◆ getTrackerTopology()

const TrackerTopology* FWGeometry::getTrackerTopology ( ) const
inline

Definition at line 156 of file FWGeometry.h.

156 { return m_trackerTopology.get(); }

References m_trackerTopology.

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

◆ initMap()

void FWGeometry::initMap ( const FWRecoGeom::InfoMap &  map)

Definition at line 200 of file FWGeometry.cc.

200  {
201  FWRecoGeom::InfoMapItr begin = map.begin();
202  FWRecoGeom::InfoMapItr end = map.end();
203  unsigned int mapSize = map.size();
204  if (m_idToInfo.size() != mapSize)
205  m_idToInfo.resize(mapSize);
206  unsigned int i = 0;
207  for (FWRecoGeom::InfoMapItr it = begin; it != end; ++it, ++i) {
208  m_idToInfo[i].id = it->id;
209  for (unsigned int j = 0; j < 24; ++j)
210  m_idToInfo[i].points[j] = it->points[j];
211  for (unsigned int j = 0; j < 9; ++j)
212  m_idToInfo[i].parameters[j] = it->topology[j];
213  for (unsigned int j = 0; j < 5; ++j)
214  m_idToInfo[i].shape[j] = it->shape[j];
215  for (unsigned int j = 0; j < 3; ++j)
216  m_idToInfo[i].translation[j] = it->translation[j];
217  for (unsigned int j = 0; j < 9; ++j)
218  m_idToInfo[i].matrix[j] = it->matrix[j];
219  }
220 }

References begin, end, cms::cuda::for(), mps_fire::i, dqmiolumiharvest::j, m_idToInfo, genParticles_cff::map, makeMuonMisalignmentScenario::matrix, and HLT_2018_cff::points.

Referenced by FWFFLooper::beginRun(), and FWFFService::postBeginRun().

◆ isEmpty()

bool FWGeometry::isEmpty ( void  ) const

Definition at line 28 of file FWGeometry.cc.

28  {
29  // AMT this is a check if geomtery is not loaded
30  // e.g. cmsShow starts with no data file and without given explicit argument ( --geometry-file option )
31 
32  return m_idToInfo.empty();
33 }

References m_idToInfo.

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

◆ loadMap()

void FWGeometry::loadMap ( const char *  fileName)

Definition at line 100 of file FWGeometry.cc.

100  {
101  TFile* file = findFile(iFileName);
102  if (!file) {
103  throw std::runtime_error("ERROR: failed to find geometry file. Initialization failed.");
104  return;
105  }
106  m_fileName = iFileName;
107  TTree* tree = static_cast<TTree*>(file->Get("idToGeo"));
108  if (!tree) {
109  throw std::runtime_error("ERROR: cannot find detector id map in the file. Initialization failed.");
110  return;
111  }
112 
113  unsigned int id;
114  Float_t points[24];
115  Float_t topology[9];
116  Float_t shape[5];
117  Float_t translation[3];
118  Float_t matrix[9];
119  bool loadPoints = tree->GetBranch("points") != nullptr;
120  bool loadParameters = tree->GetBranch("topology") != nullptr;
121  bool loadShape = tree->GetBranch("shape") != nullptr;
122  bool loadTranslation = tree->GetBranch("translation") != nullptr;
123  bool loadMatrix = tree->GetBranch("matrix") != nullptr;
124  tree->SetBranchAddress("id", &id);
125  if (loadPoints)
126  tree->SetBranchAddress("points", &points);
127  if (loadParameters)
128  tree->SetBranchAddress("topology", &topology);
129  if (loadShape)
130  tree->SetBranchAddress("shape", &shape);
131  if (loadTranslation)
132  tree->SetBranchAddress("translation", &translation);
133  if (loadMatrix)
134  tree->SetBranchAddress("matrix", &matrix);
135 
136  // reset previous values
137  m_idToInfo.clear();
138  for (const auto& p : m_idToMatrix)
139  delete p.second;
140  m_trackerTopology.reset();
141 
142  unsigned int treeSize = tree->GetEntries();
143  if (m_idToInfo.size() != treeSize)
144  m_idToInfo.resize(treeSize);
145  for (unsigned int i = 0; i < treeSize; ++i) {
146  tree->GetEntry(i);
147 
148  m_idToInfo[i].id = id;
149  if (loadPoints) {
150  for (unsigned int j = 0; j < 24; ++j)
151  m_idToInfo[i].points[j] = points[j];
152  }
153  if (loadParameters) {
154  for (unsigned int j = 0; j < 9; ++j)
155  m_idToInfo[i].parameters[j] = topology[j];
156  }
157  if (loadShape) {
158  for (unsigned int j = 0; j < 5; ++j)
159  m_idToInfo[i].shape[j] = shape[j];
160  }
161  if (loadTranslation) {
162  for (unsigned int j = 0; j < 3; ++j)
163  m_idToInfo[i].translation[j] = translation[j];
164  }
165  if (loadMatrix) {
166  for (unsigned int j = 0; j < 9; ++j)
167  m_idToInfo[i].matrix[j] = matrix[j];
168  }
169  }
170 
171  m_versionInfo.productionTag = static_cast<TNamed*>(file->Get("tag"));
172  m_versionInfo.cmsswVersion = static_cast<TNamed*>(file->Get("CMSSW_VERSION"));
173  m_versionInfo.extraDetectors = static_cast<TObjArray*>(file->Get("ExtraDetectors"));
174 
175  TString path = file->GetPath();
176  if (path.EndsWith(":/"))
177  path.Resize(path.Length() - 2);
178 
179  if (m_versionInfo.productionTag)
180  fwLog(fwlog::kInfo) << Form(
181  "Load %s %s from %s\n", tree->GetName(), m_versionInfo.productionTag->GetTitle(), path.Data());
182  else
183  fwLog(fwlog::kInfo) << Form("Load %s from %s\n", tree->GetName(), path.Data());
184 
185  TNamed* producerInfo = static_cast<TNamed*>(file->Get("PRODUCER_VERSION"));
186  if (producerInfo) {
187  m_producerVersion = atoi(producerInfo->GetTitle());
188  }
189 
190  TNamed* ttopology = static_cast<TNamed*>(file->Get("TrackerTopology"));
191  if (ttopology) {
192  std::string xml = ttopology->GetTitle();
193  m_trackerTopology = std::unique_ptr<TrackerTopology>(
194  new TrackerTopology(StandaloneTrackerTopology::fromTrackerParametersXMLString(xml)));
195  }
196 
197  file->Close();
198 }

References FWGeometry::VersionInfo::cmsswVersion, FWGeometry::VersionInfo::extraDetectors, FrontierConditions_GlobalTag_cff::file, findFile(), StandaloneTrackerTopology::fromTrackerParametersXMLString(), fwLog, mps_fire::i, triggerObjects_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, HLT_2018_cff::points, FWGeometry::VersionInfo::productionTag, AlCaHLTBitMon_QueryRunRegistry::string, and ecaldqm::topology().

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

◆ localToGlobal() [1/3]

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

Definition at line 523 of file FWGeometry.cc.

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

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

◆ localToGlobal() [2/3]

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

Definition at line 497 of file FWGeometry.cc.

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

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

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

◆ localToGlobal() [3/3]

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

Definition at line 506 of file FWGeometry.cc.

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

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

◆ mapEnd()

IdToInfoItr FWGeometry::mapEnd ( ) const
inline

Definition at line 141 of file FWGeometry.h.

141 { return m_idToInfo.end(); }

References m_idToInfo.

Referenced by FWRPZViewGeometry::makeMuonGeometryRhoZ().

◆ match_id()

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

Definition at line 131 of file FWGeometry.h.

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

References kDetOffset, kSubdetOffset, and EcalTangentSkim_cfg::o.

Referenced by getMatchedIds().

◆ versionInfo()

const VersionInfo& FWGeometry::versionInfo ( ) const
inline

Definition at line 150 of file FWGeometry.h.

150 { return m_versionInfo; }

References m_versionInfo.

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

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 173 of file FWGeometry.h.

Referenced by applyGlobalTag(), and loadMap().

◆ m_idToInfo

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().

◆ m_idToMatrix

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

Definition at line 163 of file FWGeometry.h.

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

◆ m_prodTag

std::string FWGeometry::m_prodTag
private

Definition at line 167 of file FWGeometry.h.

◆ m_producerVersion

int FWGeometry::m_producerVersion
private

Definition at line 171 of file FWGeometry.h.

Referenced by getProducerVersion(), and loadMap().

◆ m_trackerTopology

std::unique_ptr<TrackerTopology> FWGeometry::m_trackerTopology
private

Definition at line 175 of file FWGeometry.h.

Referenced by getTrackerTopology(), and loadMap().

◆ m_versionInfo

VersionInfo FWGeometry::m_versionInfo
private

Definition at line 169 of file FWGeometry.h.

Referenced by loadMap(), and versionInfo().

HLT_2018_cff.points
points
Definition: HLT_2018_cff.py:20125
DDAxes::y
fwLog
#define fwLog(_level_)
Definition: fwLog.h:45
mps_fire.i
i
Definition: mps_fire.py:355
FWGeometry::match_id
bool match_id(const GeomDetInfo &o, unsigned int mask) const
Definition: FWGeometry.h:131
makeMuonMisalignmentScenario.matrix
list matrix
Definition: makeMuonMisalignmentScenario.py:141
FWGeometry::loadMap
void loadMap(const char *fileName)
Definition: FWGeometry.cc:100
FWGeometry::TID
Definition: FWGeometry.h:47
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
FWGeometry::VersionInfo::productionTag
TNamed * productionTag
Definition: FWGeometry.h:68
TrackerTopology
Definition: TrackerTopology.h:16
FWGeometry::RPCBarrel
Definition: FWGeometry.h:52
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
tree
Definition: tree.py:1
FWGeometry::localToGlobal
void localToGlobal(unsigned int id, const float *local, float *global, bool translatep=true) const
Definition: FWGeometry.cc:497
FWGeometry::IdToInfoItr
std::vector< FWGeometry::GeomDetInfo >::const_iterator IdToInfoItr
Definition: FWGeometry.h:137
info
static const TGPicture * info(bool iBackgroundIsBlack)
Definition: FWCollectionSummaryWidget.cc:152
personalPlayback.fp
fp
Definition: personalPlayback.py:523
FWGeometry::DT
Definition: FWGeometry.h:51
mps_check.array
array
Definition: mps_check.py:216
FWGeometry::RPCEndcap
Definition: FWGeometry.h:53
DDAxes::x
MillePedeFileConverter_cfg.fileName
fileName
Definition: MillePedeFileConverter_cfg.py:32
FWGeometry::m_versionInfo
VersionInfo m_versionInfo
Definition: FWGeometry.h:169
Muon
Definition: Muon.py:1
end
#define end
Definition: vmac.h:39
EcalTangentSkim_cfg.o
o
Definition: EcalTangentSkim_cfg.py:36
parameters
parameters
Definition: BeamSpot_PayloadInspector.cc:14
FWGeometry::kDetOffset
static const int kDetOffset
Definition: FWGeometry.h:29
StandaloneTrackerTopology::fromTrackerParametersXMLString
TrackerTopology fromTrackerParametersXMLString(const std::string &xmlContent)
Definition: StandaloneTrackerTopology.cc:180
testProducerWithPsetDescEmpty_cfi.x1
x1
Definition: testProducerWithPsetDescEmpty_cfi.py:33
test
Definition: SmallWORMDict.h:13
ecaldqm::topology
const CaloTopology * topology(nullptr)
FWGeometry::Tracker
Definition: FWGeometry.h:33
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
cuda_std::lower_bound
__host__ constexpr __device__ RandomIt lower_bound(RandomIt first, RandomIt last, const T &value, Compare comp={})
Definition: cudastdAlgorithm.h:27
FWGeometry::Hcal
Definition: FWGeometry.h:36
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
FWGeometry::TOB
Definition: FWGeometry.h:48
idealTransformation.rotation
dictionary rotation
Definition: idealTransformation.py:1
FWGeometry::VersionInfo::extraDetectors
TObjArray * extraDetectors
Definition: FWGeometry.h:70
FWGeometry::CSC
Definition: FWGeometry.h:50
FrontierConditions_GlobalTag_cff.file
file
Definition: FrontierConditions_GlobalTag_cff.py:13
FWGeometry::m_idToInfo
IdToInfo m_idToInfo
Definition: FWGeometry.h:165
FWGeometry::Ecal
Definition: FWGeometry.h:35
FWGeometry::m_producerVersion
int m_producerVersion
Definition: FWGeometry.h:171
fwlog::kWarning
Definition: fwLog.h:35
PVValHelper::dy
Definition: PVValidationHelpers.h:49
FWGeometry::m_trackerTopology
std::unique_ptr< TrackerTopology > m_trackerTopology
Definition: FWGeometry.h:175
FWGeometry::TIB
Definition: FWGeometry.h:46
FWRecoGeom::InfoMapItr
std::vector< FWRecoGeom::Info >::const_iterator InfoMapItr
Definition: FWRecoGeom.h:24
FWGeometry::HGCalHSi
Definition: FWGeometry.h:39
cms::cuda::for
for(int i=first, nt=offsets[nh];i< nt;i+=gridDim.x *blockDim.x)
Definition: HistoContainer.h:27
FWGeometry::HGCalHSc
Definition: FWGeometry.h:40
FWGeometry::TEC
Definition: FWGeometry.h:49
FWGeometry::m_idToMatrix
std::map< unsigned int, TGeoMatrix * > m_idToMatrix
Definition: FWGeometry.h:163
FWGeometry::PixelBarrel
Definition: FWGeometry.h:44
alignCSCRings.r
r
Definition: alignCSCRings.py:93
FWGeometry::ME0
Definition: FWGeometry.h:55
fwlog::kInfo
Definition: fwLog.h:35
FWGeometry::VersionInfo::cmsswVersion
TNamed * cmsswVersion
Definition: FWGeometry.h:69
FWGeometry::HGCalTrigger
Definition: FWGeometry.h:41
type
type
Definition: HCALResponse.h:21
FWGeometry::findFile
static TFile * findFile(const char *fileName)
Definition: FWGeometry.cc:35
download_sqlite_cfg.globalTag
globalTag
Definition: download_sqlite_cfg.py:4
PVValHelper::dz
Definition: PVValidationHelpers.h:50
triggerObjects_cff.id
id
Definition: triggerObjects_cff.py:31
FWGeometry::Calo
Definition: FWGeometry.h:37
personalPlayback.fn
fn
Definition: personalPlayback.py:515
Skims_PA_cff.paths
paths
Definition: Skims_PA_cff.py:18
FWGeometry::GEM
Definition: FWGeometry.h:54
FWGeometry::getShape
TGeoShape * getShape(unsigned int id) const
Definition: FWGeometry.cc:290
FWGeometry::kSubdetOffset
static const int kSubdetOffset
Definition: FWGeometry.h:30
castor_dqm_sourceclient_file_cfg.path
path
Definition: castor_dqm_sourceclient_file_cfg.py:37
DTRecHitClients_cfi.local
local
Definition: DTRecHitClients_cfi.py:10
FWGeometry::m_fileName
std::string m_fileName
Definition: FWGeometry.h:173
genParticles_cff.map
map
Definition: genParticles_cff.py:11
FWGeometry::PixelEndcap
Definition: FWGeometry.h:45
dqmiolumiharvest.j
j
Definition: dqmiolumiharvest.py:66
begin
#define begin
Definition: vmac.h:32
PVValHelper::dx
Definition: PVValidationHelpers.h:48
FWGeometry::HGCalEE
Definition: FWGeometry.h:38
RemoveAddSevLevel.flag
flag
Definition: RemoveAddSevLevel.py:116
FWGeometry::find
IdToInfoItr find(unsigned int) const
Definition: FWGeometry.cc:517
Generated for CMSSW Reference Manual by doxygen 1.8.16