#include <PrintGeomInfoAction.h>

Inheritance diagram for PrintGeomInfoAction:

Public Member Functions
	PrintGeomInfoAction (edm::ParameterSet const &p)

	~PrintGeomInfoAction () override

Public Member Functions inherited from SimWatcher
const SimWatcher &	operator= (const SimWatcher &)=delete

	SimWatcher ()

	SimWatcher (const SimWatcher &)=delete

virtual	~SimWatcher ()

Public Member Functions inherited from Observer< const BeginOfJob * >
	Observer ()

void	slotForUpdate (const BeginOfJob * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const BeginOfRun * >
	Observer ()

void	slotForUpdate (const BeginOfRun * iT)

virtual	~Observer ()

Private Member Functions
void	add1touchable (G4LogicalVolume *lv, int &nTouch)

int	countNoTouchables ()

void	dumpG4LVLeaf (G4LogicalVolume *lv, unsigned int leafDepth, unsigned int count, std::ostream &out=std::cout)

void	dumpG4LVList (std::ostream &out=std::cout)

void	dumpG4LVTree (std::ostream &out=std::cout)

void	dumpHierarchyLeafPVLV (G4LogicalVolume *lv, unsigned int leafDepth, std::ostream &out=std::cout)

void	dumpHierarchyTreePVLV (std::ostream &out=std::cout)

void	dumpInFile ()

void	dumpLV (G4LogicalVolume *lv, unsigned int leafDepth, std::ostream &out=std::cout)

void	dumpMaterialList (std::ostream &out=std::cout)

void	dumpPV (G4VPhysicalVolume *pv, unsigned int leafDepth, std::ostream &out=std::cout)

void	dumpSolid (G4VSolid *sol, unsigned int leafDepth, std::ostream &out=std::cout)

void	dumpSummary (std::ostream &out=std::cout)

void	dumpTouch (G4VPhysicalVolume *pv, unsigned int leafDepth, std::ostream &out=std::cout)

G4LogicalVolume *	getTopLV ()

G4VPhysicalVolume *	getTopPV ()

void	getTouch (G4VPhysicalVolume *pv, unsigned int leafDepth, unsigned int copym, std::vector< std::string > &touches)

std::string	spacesFromLeafDepth (unsigned int leafDepth)

void	update (const BeginOfJob *job) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfRun *run) override
	This routine will be called when the appropriate signal arrives. More...

Private Attributes
bool	dd4hep_

bool	dumpAtts_

bool	dumpLV_

bool	dumpLVList_

bool	dumpLVTree_

bool	dumpMaterial_

bool	dumpPV_

bool	dumpReplica_

bool	dumpRotation_

bool	dumpSense_

bool	dumpSolid_

bool	dumpSummary_

bool	dumpTouch_

G4NavigationHistory	fHistory_

bool	fileDetail_

std::string	fileLV_

std::string	fileMat_

std::string	filePV_

std::string	fileSolid_

std::string	fileTouch_

std::string	name_

std::vector< std::string >	names_

int	nchar_

G4VPhysicalVolume *	theTopPV_

Additional Inherited Members

Detailed Description

Definition at line 23 of file PrintGeomInfoAction.h.

Constructor & Destructor Documentation

◆ PrintGeomInfoAction()

PrintGeomInfoAction::PrintGeomInfoAction ( edm::ParameterSet const & p )

Definition at line 42 of file PrintGeomInfoAction.cc.

                                                                  {
   dumpSummary_ = p.getUntrackedParameter<bool>("DumpSummary", true);
   dumpLVTree_ = p.getUntrackedParameter<bool>("DumpLVTree", true);
   dumpLVList_ = p.getUntrackedParameter<bool>("DumpLVList", false);
   dumpMaterial_ = p.getUntrackedParameter<bool>("DumpMaterial", false);
   dumpLV_ = p.getUntrackedParameter<bool>("DumpLV", false);
   dumpSolid_ = p.getUntrackedParameter<bool>("DumpSolid", false);
   dumpAtts_ = p.getUntrackedParameter<bool>("DumpAttributes", false);
   dumpPV_ = p.getUntrackedParameter<bool>("DumpPV", false);
   dumpRotation_ = p.getUntrackedParameter<bool>("DumpRotation", false);
   dumpReplica_ = p.getUntrackedParameter<bool>("DumpReplica", false);
   dumpTouch_ = p.getUntrackedParameter<bool>("DumpTouch", false);
   dumpSense_ = p.getUntrackedParameter<bool>("DumpSense", false);
   dd4hep_ = p.getUntrackedParameter<bool>("DD4Hep", false);
   name_ = p.getUntrackedParameter<std::string>("Name", "*");
   nchar_ = name_.find('*');
   name_.assign(name_, 0, nchar_);
   names_ = p.getUntrackedParameter<std::vector<std::string> >("Names");
   fileMat_ = p.getUntrackedParameter<std::string>("MaterialFileName", "");
   fileSolid_ = p.getUntrackedParameter<std::string>("SolidFileName", "");
   fileLV_ = p.getUntrackedParameter<std::string>("LVFileName", "");
   filePV_ = p.getUntrackedParameter<std::string>("PVFileName", "");
   fileTouch_ = p.getUntrackedParameter<std::string>("TouchFileName", "");
   fileDetail_ = p.getUntrackedParameter<bool>("FileDetail", false);
   G4cout << "PrintGeomInfoAction:: initialised for dd4hep " << dd4hep_ << " with verbosity levels:"
          << " Summary   " << dumpSummary_ << " LVTree   " << dumpLVTree_ << " LVList    " << dumpLVList_ << " Material "
          << dumpMaterial_ << "\n                                                        "
          << " LV        " << dumpLV_ << " Solid    " << dumpSolid_ << " Attribs   " << dumpAtts_
          << "\n                                                        "
          << " PV        " << dumpPV_ << " Rotation " << dumpRotation_ << " Replica   " << dumpReplica_
          << "\n                                                        "
          << " Touchable " << dumpTouch_ << " for names (0-" << nchar_ << ") = " << name_
          << "\n                                                        "
          << " Sensitive " << dumpSense_ << " Files " << fileMat_ << ":" << fileSolid_ << ":" << fileLV_ << ":"
          << filePV_ << ":" << fileTouch_ << " FileDetail " << fileDetail_
          << "\n                                                        for " << names_.size() << " names:";
   for (unsigned int i = 0; i < names_.size(); i++)
     G4cout << " " << names_[i];
   G4cout << G4endl;
 }

References dd4hep_, dumpAtts_, dumpLV_, dumpLVList_, dumpLVTree_, dumpMaterial_, dumpPV_, dumpReplica_, dumpRotation_, dumpSense_, dumpSolid_, dumpSummary_, dumpTouch_, fileDetail_, fileLV_, fileMat_, filePV_, fileSolid_, fileTouch_, ecalTB2006H4_GenSimDigiReco_cfg::G4cout, mps_fire::i, name_, names_, nchar_, AlCaHLTBitMon_ParallelJobs::p, and AlCaHLTBitMon_QueryRunRegistry::string.

◆ ~PrintGeomInfoAction()

PrintGeomInfoAction::~PrintGeomInfoAction ( )

override

Definition at line 83 of file PrintGeomInfoAction.cc.

83 {}

Member Function Documentation

◆ add1touchable()

void PrintGeomInfoAction::add1touchable	(	G4LogicalVolume *	lv,
		int &	nTouch
	)

private

Definition at line 246 of file PrintGeomInfoAction.cc.

                                                                         {
   int siz = lv->GetNoDaughters();
   for (int ii = 0; ii < siz; ii++)
     add1touchable(lv->GetDaughter(ii)->GetLogicalVolume(), ++nTouch);
 }

References cuy::ii.

Referenced by countNoTouchables().

◆ countNoTouchables()

int PrintGeomInfoAction::countNoTouchables ( )

private

Definition at line 239 of file PrintGeomInfoAction.cc.

                                            {
   int nTouch = 0;
   G4LogicalVolume *lv = getTopLV();
   add1touchable(lv, nTouch);
   return nTouch;
 }

References add1touchable(), and getTopLV().

Referenced by dumpSummary().

◆ dumpG4LVLeaf()

void PrintGeomInfoAction::dumpG4LVLeaf	(	G4LogicalVolume *	lv,
		unsigned int	leafDepth,
		unsigned int	count,
		std::ostream &	out = `std::cout`
	)

private

Definition at line 217 of file PrintGeomInfoAction.cc.

                                                         {
   for (unsigned int ii = 0; ii < leafDepth; ii++)
     out << "  ";
   out << " LV:(" << leafDepth << ") " << lv->GetName() << " (" << count << ")" << G4endl;
   //--- If a volume is placed n types as daughter of this LV, it should only be counted once
   std::map<G4LogicalVolume *, unsigned int> lvCount;
   std::map<G4LogicalVolume *, unsigned int>::const_iterator cite;
   int siz = lv->GetNoDaughters();
   for (int ii = 0; ii < siz; ii++) {
     cite = lvCount.find(lv->GetDaughter(ii)->GetLogicalVolume());
     if (cite != lvCount.end())
       lvCount[cite->first] = (cite->second) + 1;
     else
       lvCount.insert(std::pair<G4LogicalVolume *, unsigned int>(lv->GetDaughter(ii)->GetLogicalVolume(), 1));
   }
   for (cite = lvCount.begin(); cite != lvCount.end(); cite++)
     dumpG4LVLeaf((cite->first), leafDepth + 1, (cite->second), out);
 }

References submitPVResolutionJobs::count, cuy::ii, and MillePedeFileConverter_cfg::out.

Referenced by dumpG4LVTree().

◆ dumpG4LVList()

void PrintGeomInfoAction::dumpG4LVList ( std::ostream & out = std::cout )

private

Definition at line 194 of file PrintGeomInfoAction.cc.

                                                       {
   out << " @@@@@@@@@@@@@@@@ DUMPING G4LogicalVolume's List  " << G4endl;
   const G4LogicalVolumeStore *lvs = G4LogicalVolumeStore::GetInstance();
   std::vector<G4LogicalVolume *>::const_iterator lvcite;
   for (lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++)
     out << "LV:" << (*lvcite)->GetName() << "\tMaterial: " << (*lvcite)->GetMaterial()->GetName() << G4endl;
 }

References MillePedeFileConverter_cfg::out.

Referenced by update().

◆ dumpG4LVTree()

void PrintGeomInfoAction::dumpG4LVTree ( std::ostream & out = std::cout )

private

Definition at line 202 of file PrintGeomInfoAction.cc.

                                                       {
   out << " @@@@@@@@@@@@@@@@ DUMPING G4LogicalVolume's Tree  " << G4endl;
   G4LogicalVolume *lv = getTopLV();
   dumpG4LVLeaf(lv, 0, 1, out);
 }

References dumpG4LVLeaf(), getTopLV(), and MillePedeFileConverter_cfg::out.

Referenced by update().

◆ dumpHierarchyLeafPVLV()

void PrintGeomInfoAction::dumpHierarchyLeafPVLV	(	G4LogicalVolume *	lv,
		unsigned int	leafDepth,
		std::ostream &	out = `std::cout`
	)

private

Definition at line 270 of file PrintGeomInfoAction.cc.

                                                                                                             {
   //----- Dump this LV
   dumpLV(lv, leafDepth, out);
  
   //----- Get LV daughters from list of PV daughters
   mmlvpv lvpvDaughters;
   std::set<G4LogicalVolume *> lvDaughters;
   int NoDaughters = lv->GetNoDaughters();
   while ((NoDaughters--) > 0) {
     G4VPhysicalVolume *pvD = lv->GetDaughter(NoDaughters);
     lvpvDaughters.insert(mmlvpv::value_type(pvD->GetLogicalVolume(), pvD));
     lvDaughters.insert(pvD->GetLogicalVolume());
   }
  
   std::set<G4LogicalVolume *>::const_iterator scite;
   mmlvpv::const_iterator mmcite;
  
   //----- Dump daughters PV and LV
   for (scite = lvDaughters.begin(); scite != lvDaughters.end(); scite++) {
     std::pair<mmlvpv::iterator, mmlvpv::iterator> mmER = lvpvDaughters.equal_range(*scite);
     //----- Dump daughters PV of this LV
     for (mmcite = mmER.first; mmcite != mmER.second; mmcite++)
       dumpPV((*mmcite).second, leafDepth + 1, out);
     //----- Dump daughters LV
     dumpHierarchyLeafPVLV(*scite, leafDepth + 1, out);
   }
 }

References dumpLV(), dumpPV(), and MillePedeFileConverter_cfg::out.

Referenced by dumpHierarchyTreePVLV().

◆ dumpHierarchyTreePVLV()

void PrintGeomInfoAction::dumpHierarchyTreePVLV ( std::ostream & out = std::cout )

private

Definition at line 252 of file PrintGeomInfoAction.cc.

                                                                {
   //dumps in the following order:
   //    1) a LV with details
   //    2) list of PVs daughters of this LV with details
   //    3) list of LVs daughters of this LV and for each go to 1)
  
   //----- Get top PV
   G4LogicalVolume *topLV = getTopLV();
  
   //----- Dump this leaf (it will recursively dump all the tree)
   dumpHierarchyLeafPVLV(topLV, 0, out);
   dumpPV(theTopPV_, 0, out);
  
   //----- Dump the touchables (it will recursively dump all the tree)
   if (dumpTouch_)
     dumpTouch(theTopPV_, 0, out);
 }

References dumpHierarchyLeafPVLV(), dumpPV(), dumpTouch(), dumpTouch_, getTopLV(), MillePedeFileConverter_cfg::out, and theTopPV_.

Referenced by update().

◆ dumpInFile()

void PrintGeomInfoAction::dumpInFile ( )

private

Definition at line 430 of file PrintGeomInfoAction.cc.

                                      {
   //---------- Dump number objects of each class in a file
   if (theTopPV_ != nullptr) {
     if (!fileMat_.empty()) {
       const G4MaterialTable *matTab = G4Material::GetMaterialTable();
       std::ofstream fout(fileMat_.c_str());
       for (std::vector<G4Material *>::const_iterator matite = matTab->begin(); matite != matTab->end(); matite++) {
         if (!fileDetail_)
           fout << (*matite)->GetName() << "\n";
         else
           fout << (*matite)->GetName() << " " << (*matite)->GetRadlen() << " " << (*matite)->GetNuclearInterLength()
                << "\n";
       }
       fout.close();
     }
     const G4LogicalVolumeStore *lvs = G4LogicalVolumeStore::GetInstance();
     if (!fileSolid_.empty()) {
       std::ofstream fout(fileSolid_.c_str());
       for (std::vector<G4LogicalVolume *>::const_iterator lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++)
         if (!fileDetail_)
           fout << (*lvcite)->GetSolid()->GetName() << "\n";
         else
           fout << (*lvcite)->GetSolid()->GetName() << "  " << (*lvcite)->GetSolid()->GetCubicVolume() << "\n";
       fout.close();
     }
     if (!fileLV_.empty()) {
       std::ofstream fout(fileLV_.c_str());
       for (std::vector<G4LogicalVolume *>::const_iterator lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++)
         if (!fileDetail_)
           fout << (*lvcite)->GetName() << "\n";
         else
           fout << (*lvcite)->GetName() << "  " << (*lvcite)->GetMass(false, false) << "\n";
       fout.close();
     }
     if (!filePV_.empty()) {
       const G4PhysicalVolumeStore *pvs = G4PhysicalVolumeStore::GetInstance();
       std::ofstream fout(filePV_.c_str());
       for (std::vector<G4VPhysicalVolume *>::const_iterator pvcite = pvs->begin(); pvcite != pvs->end(); pvcite++) {
         if (fileDetail_) {
           if (dd4hep_)
             fout << (*pvcite)->GetName() << " " << (*pvcite)->GetTranslation().x() << " "
                  << (*pvcite)->GetTranslation().y() << " " << (*pvcite)->GetTranslation().z() << "\n";
           else
             fout << (*pvcite)->GetName() << "_" << (*pvcite)->GetCopyNo() << " " << (*pvcite)->GetTranslation().x()
                  << " " << (*pvcite)->GetTranslation().y() << " " << (*pvcite)->GetTranslation().z() << "\n";
         } else {
           if (dd4hep_)
             fout << (*pvcite)->GetName() << "\n";
           else
             fout << (*pvcite)->GetName() << "_" << (*pvcite)->GetCopyNo() << "\n";
         }
       }
       fout.close();
     }
     if (!fileTouch_.empty()) {
       std::ofstream fout(fileTouch_.c_str());
       std::vector<std::string> touches;
       getTouch(theTopPV_, 0, 1, touches);
       std::sort(touches.begin(), touches.end());
       for (const auto &touch : touches)
         fout << touch << "\n";
       fout.close();
     }
   }
 }

References dd4hep_, fileDetail_, fileLV_, fileMat_, filePV_, fileSolid_, fileTouch_, groupFilesInBlocks::fout, getTouch(), FSQDQM_cfi::pvs, jetUpdater_cfi::sort, and theTopPV_.

Referenced by update().

◆ dumpLV()

void PrintGeomInfoAction::dumpLV	(	G4LogicalVolume *	lv,
		unsigned int	leafDepth,
		std::ostream &	out = `std::cout`
	)

private

Definition at line 298 of file PrintGeomInfoAction.cc.

                                                                                              {
   std::string spaces = spacesFromLeafDepth(leafDepth);
  
   //----- dump name
   if (dumpLV_) {
     out << leafDepth << spaces << "$$$ VOLUME = " << lv->GetName() << "  Solid: " << lv->GetSolid()->GetName()
         << "  MATERIAL: " << lv->GetMaterial()->GetName() << G4endl;
     if (dumpSolid_)
       dumpSolid(lv->GetSolid(), leafDepth, out);  //----- dump solid
  
     //----- dump LV info
     //--- material
     if (dumpAtts_) {
       //--- Visualisation attributes
       const G4VisAttributes *fVA = lv->GetVisAttributes();
       if (fVA != nullptr) {
         out << spaces << "  VISUALISATION ATTRIBUTES: " << G4endl;
         out << spaces << "    IsVisible " << fVA->IsVisible() << G4endl;
         out << spaces << "    IsDaughtersInvisible " << fVA->IsDaughtersInvisible() << G4endl;
         out << spaces << "    Colour " << fVA->GetColour() << G4endl;
         out << spaces << "    LineStyle " << fVA->GetLineStyle() << G4endl;
         out << spaces << "    LineWidth " << fVA->GetLineWidth() << G4endl;
         out << spaces << "    IsForceDrawingStyle " << fVA->IsForceDrawingStyle() << G4endl;
         out << spaces << "    ForcedDrawingStyle " << fVA->GetForcedDrawingStyle() << G4endl;
       }
  
       //--- User Limits
       G4UserLimits *fUL = lv->GetUserLimits();
       G4Track dummy;
       if (fUL != nullptr) {
         out << spaces << "    MaxAllowedStep " << fUL->GetMaxAllowedStep(dummy) << G4endl;
         out << spaces << "    UserMaxTrackLength " << fUL->GetUserMaxTrackLength(dummy) << G4endl;
         out << spaces << "    UserMaxTime " << fUL->GetUserMaxTime(dummy) << G4endl;
         out << spaces << "    UserMinEkine " << fUL->GetUserMinEkine(dummy) << G4endl;
         out << spaces << "    UserMinRange " << fUL->GetUserMinRange(dummy) << G4endl;
       }
  
       //--- other LV info
       if (lv->GetSensitiveDetector())
         out << spaces << "  IS SENSITIVE DETECTOR " << G4endl;
       if (lv->GetFieldManager())
         out << spaces << "  FIELD ON " << G4endl;
  
       // Pointer (possibly NULL) to optimisation info objects.
       out << spaces << "        Quality for optimisation, average number of voxels to be spent per content "
           << lv->GetSmartless() << G4endl;
  
       // Pointer (possibly NULL) to G4FastSimulationManager object.
       if (lv->GetFastSimulationManager())
         out << spaces << "     Logical Volume is an envelope for a FastSimulationManager " << G4endl;
       out << spaces << "     Weight used in the event biasing technique = " << lv->GetBiasWeight() << G4endl;
     }
   }
 }

References dumpAtts_, dumpLV_, dumpSolid(), dumpSolid_, MillePedeFileConverter_cfg::out, spacesFromLeafDepth(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by dumpHierarchyLeafPVLV().

◆ dumpMaterialList()

void PrintGeomInfoAction::dumpMaterialList ( std::ostream & out = std::cout )

private

Definition at line 208 of file PrintGeomInfoAction.cc.

                                                           {
   out << " @@@@@@@@@@@@@@@@ DUMPING G4Material List ";
   const G4MaterialTable *matTab = G4Material::GetMaterialTable();
   out << " with " << matTab->size() << " materials " << G4endl;
   std::vector<G4Material *>::const_iterator matite;
   for (matite = matTab->begin(); matite != matTab->end(); matite++)
     out << "Material: " << (*matite) << G4endl;
 }

References MillePedeFileConverter_cfg::out.

Referenced by update().

◆ dumpPV()

void PrintGeomInfoAction::dumpPV	(	G4VPhysicalVolume *	pv,
		unsigned int	leafDepth,
		std::ostream &	out = `std::cout`
	)

private

Definition at line 353 of file PrintGeomInfoAction.cc.

                                                                                                {
   std::string spaces = spacesFromLeafDepth(leafDepth);
  
   //----- PV info
   if (dumpPV_) {
     std::string mother = "World";
     if (pv->GetMotherLogical())
       mother = pv->GetMotherLogical()->GetName();
     out << leafDepth << spaces << "### VOLUME = " << pv->GetName() << " Copy No " << pv->GetCopyNo() << " in " << mother
         << " at " << pv->GetTranslation();
   }
   if (!pv->IsReplicated()) {
     if (dumpPV_) {
       if (pv->GetRotation() == nullptr)
         out << " with no rotation" << G4endl;
       else if (!dumpRotation_)
         out << " with rotation" << G4endl;  //just rotation name
       else
         out << " with rotation " << *(pv->GetRotation()) << G4endl;
     }
   } else {
     if (dumpReplica_) {
       out << spaces << "    It is replica: " << G4endl;
       EAxis axis;
       int nReplicas;
       double width;
       double offset;
       bool consuming;
       pv->GetReplicationData(axis, nReplicas, width, offset, consuming);
       out << spaces << "     axis " << axis << G4endl << spaces << "     nReplicas " << nReplicas << G4endl;
       if (pv->GetParameterisation() != nullptr)
         out << spaces << "    It is parameterisation " << G4endl;
       else
         out << spaces << "     width " << width << G4endl << spaces << "     offset " << offset << G4endl << spaces
             << "     consuming" << consuming << G4endl;
       if (pv->GetParameterisation() != nullptr)
         out << spaces << "    It is parameterisation " << G4endl;
     }
   }
 }

References dumpPV_, dumpReplica_, dumpRotation_, hltrates_dqm_sourceclient-live_cfg::offset, MillePedeFileConverter_cfg::out, MetAnalyzer::pv(), spacesFromLeafDepth(), AlCaHLTBitMon_QueryRunRegistry::string, and ApeEstimator_cff::width.

Referenced by dumpHierarchyLeafPVLV(), and dumpHierarchyTreePVLV().

◆ dumpSolid()

void PrintGeomInfoAction::dumpSolid	(	G4VSolid *	sol,
		unsigned int	leafDepth,
		std::ostream &	out = `std::cout`
	)

private

Definition at line 394 of file PrintGeomInfoAction.cc.

                                                                                           {
   std::string spaces = spacesFromLeafDepth(leafDepth);
   out << spaces << *(sol) << G4endl;
 }

References MillePedeFileConverter_cfg::out, spacesFromLeafDepth(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by dumpLV().

◆ dumpSummary()

void PrintGeomInfoAction::dumpSummary ( std::ostream & out = std::cout )

private

Definition at line 171 of file PrintGeomInfoAction.cc.

                                                      {
   //---------- Dump number of objects of each class
   out << " @@@@@@@@@@@@@@@@@@ Dumping G4 geometry objects Summary " << G4endl;
   if (theTopPV_ == nullptr) {
     out << " No volume created " << G4endl;
     return;
   }
   out << " @@@ Geometry built inside world volume: " << theTopPV_->GetName() << G4endl;
   // Get number of solids (< # LV if several LV share a solid)
   const G4LogicalVolumeStore *lvs = G4LogicalVolumeStore::GetInstance();
   std::vector<G4LogicalVolume *>::const_iterator lvcite;
   std::set<G4VSolid *> theSolids;
   for (lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++)
     theSolids.insert((*lvcite)->GetSolid());
   out << " Number of G4VSolid's: " << theSolids.size() << G4endl;
   out << " Number of G4LogicalVolume's: " << lvs->size() << G4endl;
   const G4PhysicalVolumeStore *pvs = G4PhysicalVolumeStore::GetInstance();
   out << " Number of G4VPhysicalVolume's: " << pvs->size() << G4endl;
   out << " Number of Touchable's: " << countNoTouchables() << G4endl;
   const G4MaterialTable *matTab = G4Material::GetMaterialTable();
   out << " Number of G4Material's: " << matTab->size() << G4endl;
 }

References countNoTouchables(), MillePedeFileConverter_cfg::out, FSQDQM_cfi::pvs, and theTopPV_.

Referenced by update().

◆ dumpTouch()

void PrintGeomInfoAction::dumpTouch	(	G4VPhysicalVolume *	pv,
		unsigned int	leafDepth,
		std::ostream &	out = `std::cout`
	)

private

Definition at line 399 of file PrintGeomInfoAction.cc.

                                                                                                   {
   std::string spaces = spacesFromLeafDepth(leafDepth);
   if (leafDepth == 0)
     fHistory_.SetFirstEntry(pv);
   else
     fHistory_.NewLevel(pv, kNormal, pv->GetCopyNo());
  
   G4ThreeVector globalpoint = fHistory_.GetTopTransform().Inverse().TransformPoint(G4ThreeVector(0, 0, 0));
   G4LogicalVolume *lv = pv->GetLogicalVolume();
  
   std::string mother = "World";
   if (pv->GetMotherLogical())
     mother = pv->GetMotherLogical()->GetName();
   std::string lvname = lv->GetName();
   lvname.assign(lvname, 0, nchar_);
   if (lvname == name_)
     out << leafDepth << spaces << "### VOLUME = " << lv->GetName() << " Copy No " << pv->GetCopyNo() << " in " << mother
         << " global position of centre " << globalpoint << " (r = " << globalpoint.perp()
         << ", phi = " << convertRadToDeg(globalpoint.phi()) << ")" << G4endl;
  
   int NoDaughters = lv->GetNoDaughters();
   while ((NoDaughters--) > 0) {
     G4VPhysicalVolume *pvD = lv->GetDaughter(NoDaughters);
     if (!pvD->IsReplicated())
       dumpTouch(pvD, leafDepth + 1, out);
   }
  
   if (leafDepth > 0)
     fHistory_.BackLevel();
 }

References angle_units::operators::convertRadToDeg(), fHistory_, name_, nchar_, MillePedeFileConverter_cfg::out, MetAnalyzer::pv(), spacesFromLeafDepth(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by dumpHierarchyTreePVLV().

◆ getTopLV()

G4LogicalVolume * PrintGeomInfoAction::getTopLV ( )

private

Definition at line 540 of file PrintGeomInfoAction.cc.

540 { return theTopPV_->GetLogicalVolume(); }

References theTopPV_.

Referenced by countNoTouchables(), dumpG4LVTree(), and dumpHierarchyTreePVLV().

◆ getTopPV()

G4VPhysicalVolume * PrintGeomInfoAction::getTopPV ( )

private

Definition at line 536 of file PrintGeomInfoAction.cc.

                                                  {
   return G4TransportationManager::GetTransportationManager()->GetNavigatorForTracking()->GetWorldVolume();
 }

Referenced by update().

◆ getTouch()

void PrintGeomInfoAction::getTouch	(	G4VPhysicalVolume *	pv,
		unsigned int	leafDepth,
		unsigned int	copym,
		std::vector< std::string > &	touches
	)

private

Definition at line 496 of file PrintGeomInfoAction.cc.

                                                                   {
   if (leafDepth == 0)
     fHistory_.SetFirstEntry(pv);
   else
     fHistory_.NewLevel(pv, kNormal, pv->GetCopyNo());
  
   std::string mother = "World";
   if (pv->GetMotherLogical())
     mother = static_cast<std::string>(dd4hep::dd::noNamespace(pv->GetMotherLogical()->GetName()));
  
   G4LogicalVolume *lv = pv->GetLogicalVolume();
   std::string lvname = static_cast<std::string>(dd4hep::dd::noNamespace(lv->GetName()));
   unsigned int copy = static_cast<unsigned int>(pv->GetCopyNo());
  
   std::string name = lvname + ":" + std::to_string(copy) + "_" + mother + ":" + std::to_string(copym);
   touches.emplace_back(name);
  
   int NoDaughters = lv->GetNoDaughters();
   while ((NoDaughters--) > 0) {
     G4VPhysicalVolume *pvD = lv->GetDaughter(NoDaughters);
     if (!pvD->IsReplicated())
       getTouch(pvD, leafDepth + 1, copy, touches);
   }
  
   if (leafDepth > 0)
     fHistory_.BackLevel();
 }

References filterCSVwithJSON::copy, fHistory_, Skims_PA_cff::name, MetAnalyzer::pv(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by dumpInFile().

◆ spacesFromLeafDepth()

std::string PrintGeomInfoAction::spacesFromLeafDepth ( unsigned int leafDepth )

private

Definition at line 527 of file PrintGeomInfoAction.cc.

                                                                          {
   std::string spaces;
   unsigned int ii;
   for (ii = 0; ii < leafDepth; ii++) {
     spaces += "  ";
   }
   return spaces;
 }

References cuy::ii, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by dumpLV(), dumpPV(), dumpSolid(), dumpTouch(), and update().

◆ update() [1/2]

void PrintGeomInfoAction::update ( const BeginOfJob * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 85 of file PrintGeomInfoAction.cc.

                                                       {
   if (dumpSense_) {
     if (dd4hep_) {
       edm::ESTransientHandle<cms::DDCompactView> pDD;
       (*job)()->get<IdealGeometryRecord>().get(pDD);
  
       G4cout << "PrintGeomInfoAction::Get Printout of Sensitive Volumes "
              << "for " << names_.size() << " Readout Units" << G4endl;
       for (unsigned int i = 0; i < names_.size(); i++) {
         std::string sd = names_[i];
         const cms::DDFilter filter("ReadOutName", sd);
         cms::DDFilteredView fv(*pDD, filter);
         G4cout << "PrintGeomInfoAction:: Get Filtered view for ReadOutName = " << sd << G4endl;
         G4cout << "Lengths are in mm, angles in degrees" << G4endl;
  
         std::string spaces = spacesFromLeafDepth(1);
  
         while (fv.firstChild()) {
           auto tran = fv.translation() / dd4hep::mm;
           std::vector<int> copy = fv.copyNos();
           auto lvname = fv.name();
           unsigned int leafDepth = copy.size();
           G4cout << leafDepth << spaces << "### VOLUME = " << lvname << " Copy No";
           for (unsigned int k = 0; k < leafDepth; ++k)
             G4cout << " " << copy[k];
           G4cout << " Centre at " << tran << " (r = " << tran.Rho() << ", phi = " << convertRadToDeg(tran.phi()) << ")"
                  << G4endl;
         }
       }
     } else {
       edm::ESTransientHandle<DDCompactView> pDD;
       (*job)()->get<IdealGeometryRecord>().get(pDD);
  
       G4cout << "PrintGeomInfoAction::Get Printout of Sensitive Volumes "
              << "for " << names_.size() << " Readout Units" << G4endl;
       for (unsigned int i = 0; i < names_.size(); i++) {
         std::string attribute = "ReadOutName";
         std::string sd = names_[i];
         DDSpecificsMatchesValueFilter filter{DDValue(attribute, sd, 0)};
         DDFilteredView fv(*pDD, filter);
         G4cout << "PrintGeomInfoAction:: Get Filtered view for " << attribute << " = " << sd << G4endl;
         G4cout << "Lengths are in mm, angles in degrees" << G4endl;
         bool dodet = fv.firstChild();
  
         std::string spaces = spacesFromLeafDepth(1);
  
         while (dodet) {
           const DDLogicalPart &log = fv.logicalPart();
           std::string lvname = log.name().name();
           DDTranslation tran = fv.translation();
           std::vector<int> copy = fv.copyNumbers();
  
           unsigned int leafDepth = copy.size();
           G4cout << leafDepth << spaces << "### VOLUME = " << lvname << " Copy No";
           for (int k = leafDepth - 1; k >= 0; k--)
             G4cout << " " << copy[k];
           G4cout << " Centre at " << tran << " (r = " << tran.Rho() << ", phi = " << convertRadToDeg(tran.phi()) << ")"
                  << G4endl;
           dodet = fv.next();
         }
       }
     }
   }
 }

References angle_units::operators::convertRadToDeg(), filterCSVwithJSON::copy, cms::DDFilteredView::copyNos(), DDFilteredView::copyNumbers(), dd4hep_, dumpSense_, ALCARECOTkAlBeamHalo_cff::filter, DDFilteredView::firstChild(), cms::DDFilteredView::firstChild(), ecalTB2006H4_GenSimDigiReco_cfg::G4cout, mps_fire::i, dqmdumpme::k, dqm-mbProfile::log, DDFilteredView::logicalPart(), cms::DDFilteredView::name(), names_, DDFilteredView::next(), sd, spacesFromLeafDepth(), AlCaHLTBitMon_QueryRunRegistry::string, DDFilteredView::translation(), and cms::DDFilteredView::translation().

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

◆ update() [2/2]

void PrintGeomInfoAction::update ( const BeginOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 150 of file PrintGeomInfoAction.cc.

                                                       {
   theTopPV_ = getTopPV();
  
   if (dumpSummary_)
     dumpSummary(G4cout);
   if (dumpLVTree_)
     dumpG4LVTree(G4cout);
  
   //---------- Dump list of objects of each class with detail of parameters
   if (dumpMaterial_)
     dumpMaterialList(G4cout);
   if (dumpLVList_)
     dumpG4LVList(G4cout);
  
   //---------- Dump LV and PV information
   if (dumpLV_ || dumpPV_ || dumpTouch_)
     dumpHierarchyTreePVLV(G4cout);
  
   dumpInFile();
 }