PrintGeomSummary Class Reference

#include <PrintGeomSummary.h>

Inheritance diagram for PrintGeomSummary:

Public Member Functions
	PrintGeomSummary (edm::ParameterSet const &p)
	~PrintGeomSummary () override
Public Member Functions inherited from SimWatcher
	SimWatcher ()
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	addName (std::string name)
void	addSolid (const DDLogicalPart &part)
void	dumpSummary (std::ostream &out, std::string name)
void	fillLV (G4LogicalVolume *lv)
void	fillPV (G4VPhysicalVolume *pv)
G4VPhysicalVolume *	getTopPV ()
void	printSummary (std::ostream &out)
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
std::map< DDSolidShape, std::pair< int, int > >	kount_
std::vector< G4LogicalVolume * >	lvs_
std::vector< std::string >	nodeNames_
std::vector< G4VPhysicalVolume * >	pvs_
std::vector< G4VSolid * >	sls_
std::map< std::string, DDSolidShape >	solidMap_
std::map< DDSolidShape, std::string >	solidShape_
G4VPhysicalVolume *	theTopPV_
std::vector< G4LogicalVolume * >	touch_

Additional Inherited Members

Detailed Description

Definition at line 24 of file PrintGeomSummary.h.

Constructor & Destructor Documentation

PrintGeomSummary()

PrintGeomSummary::PrintGeomSummary

(

edm::ParameterSet const &

p

)

Definition at line 32 of file PrintGeomSummary.cc.

                                                           : theTopPV_(nullptr) {
  std::vector<std::string> defNames;
  nodeNames_ = p.getUntrackedParameter<std::vector<std::string>>("NodeNames", defNames);
  G4cout << "PrintGeomSummary:: initialised for " << nodeNames_.size() << " nodes:" << G4endl;
  for (unsigned int ii = 0; ii < nodeNames_.size(); ii++)
    G4cout << "Node[" << ii << "] : " << nodeNames_[ii] << G4endl;
 
  solidShape_[DDSolidShape::ddbox] = "Box";
  solidShape_[DDSolidShape::ddtubs] = "Tube";
  solidShape_[DDSolidShape::ddtrap] = "Trapezoid";
  solidShape_[DDSolidShape::ddcons] = "Cone";
  solidShape_[DDSolidShape::ddpolycone_rz] = "Polycone_rz";
  solidShape_[DDSolidShape::ddpolyhedra_rz] = "Polyhedra_rz";
  solidShape_[DDSolidShape::ddpolycone_rrz] = "Polycone_rrz";
  solidShape_[DDSolidShape::ddpolyhedra_rrz] = "Polyhedra_rrz";
  solidShape_[DDSolidShape::ddtorus] = "Torus";
  solidShape_[DDSolidShape::ddunion] = "UnionSolid";
  solidShape_[DDSolidShape::ddsubtraction] = "SubtractionSolid";
  solidShape_[DDSolidShape::ddintersection] = "IntersectionSolid";
  solidShape_[DDSolidShape::ddshapeless] = "ShapelessSolid";
  solidShape_[DDSolidShape::ddpseudotrap] = "PseudoTrapezoid";
  solidShape_[DDSolidShape::ddtrunctubs] = "TruncatedTube";
  solidShape_[DDSolidShape::ddsphere] = "Sphere";
  solidShape_[DDSolidShape::ddellipticaltube] = "EllipticalTube";
  solidShape_[DDSolidShape::ddcuttubs] = "CutTubs";
  solidShape_[DDSolidShape::ddextrudedpolygon] = "ExtrudedPolygon";
  solidShape_[DDSolidShape::dd_not_init] = "Unknown";
}

References dd_not_init, ddbox, ddcons, ddcuttubs, ddellipticaltube, ddextrudedpolygon, ddintersection, ddpolycone_rrz, ddpolycone_rz, ddpolyhedra_rrz, ddpolyhedra_rz, ddpseudotrap, ddshapeless, ddsphere, ddsubtraction, ddtorus, ddtrap, ddtrunctubs, ddtubs, ddunion, ecalTB2006H4_GenSimDigiReco_cfg::G4cout, cuy::ii, nodeNames_, AlCaHLTBitMon_ParallelJobs::p, and solidShape_.

~PrintGeomSummary()

PrintGeomSummary::~PrintGeomSummary ( )

override

Definition at line 61 of file PrintGeomSummary.cc.

{}

Member Function Documentation

addName()

void PrintGeomSummary::addName ( std::string  name )

private

Definition at line 194 of file PrintGeomSummary.cc.

                                             {
  bool refl(false);
  if (name.find("_refl") < name.size()) {
    refl = true;
    name = name.substr(0, (name.find("_refl")));
  }
  std::map<std::string, DDSolidShape>::const_iterator jt = solidMap_.find(name);
  DDSolidShape shape = (jt == solidMap_.end()) ? DDSolidShape::dd_not_init : jt->second;
  std::map<DDSolidShape, std::pair<int, int>>::iterator itr = kount_.find(shape);
  if (itr == kount_.end()) {
    kount_[shape] = (refl) ? std::pair<int, int>(0, 1) : std::pair<int, int>(1, 0);
  } else {
    kount_[shape] = (refl) ? std::pair<int, int>(((itr->second).first), ++((itr->second).second))
                           : std::pair<int, int>(++((itr->second).first), ((itr->second).second));
  }
}

References dd_not_init, kount_, genParticles_cff::map, Skims_PA_cff::name, truncPyr::refl(), edm::second(), and solidMap_.

Referenced by dumpSummary().

addSolid()

void PrintGeomSummary::addSolid ( const DDLogicalPart &  part )

private

Definition at line 90 of file PrintGeomSummary.cc.

                                                         {
  const DDSolid& solid = part.solid();
  std::map<DDSolidShape, std::string>::iterator it = solidShape_.find(solid.shape());
  std::string name = solid.name().name();
  if (it == solidShape_.end())
    solidMap_[name] = DDSolidShape::dd_not_init;
  else
    solidMap_[name] = it->first;
  //G4cout << "Solid " << name << " is of shape " << solidMap_[name] << G4endl;
}

References dd_not_init, Skims_PA_cff::name, DDName::name(), DDBase< N, C >::name(), DDSolid::shape(), solidMap_, solidShape_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by update().

dumpSummary()

void PrintGeomSummary::dumpSummary ( std::ostream &  out, std::string  name  )

private

Definition at line 158 of file PrintGeomSummary.cc.

                                                                  {
  //---------- Dump number of objects of each class
  out << G4endl << G4endl << "@@@@@@@@@@@@@@@@@@ Dumping Summary For Node " << name << G4endl;
  out << " Number of G4VSolid's: " << sls_.size() << G4endl;
  out << " Number of G4LogicalVolume's: " << lvs_.size() << G4endl;
  out << " Number of Touchable's: " << touch_.size() << G4endl;
  //First the solids
  out << G4endl << "Occurence of each type of shape among Solids" << G4endl;
  kount_.clear();
  for (std::vector<G4VSolid*>::iterator it = sls_.begin(); it != sls_.end(); ++it) {
    std::string name = (*it)->GetName();
    addName(name);
  }
  printSummary(out);
  //Then the logical volumes
  out << G4endl << "Occurence of each type of shape among Logical Volumes" << G4endl;
  kount_.clear();
  for (std::vector<G4LogicalVolume*>::iterator it = lvs_.begin(); it != lvs_.end(); ++it) {
    std::string name = ((*it)->GetSolid())->GetName();
    addName(name);
  }
  printSummary(out);
  //Finally the touchables
  out << G4endl << "Occurence of each type of shape among Touchables" << G4endl;
  kount_.clear();
  for (std::vector<G4LogicalVolume*>::iterator it = touch_.begin(); it != touch_.end(); ++it) {
    std::string name = ((*it)->GetSolid())->GetName();
    addName(name);
  }
  printSummary(out);
}

References addName(), kount_, lvs_, Skims_PA_cff::name, MillePedeFileConverter_cfg::out, printSummary(), sls_, AlCaHLTBitMon_QueryRunRegistry::string, and touch_.

Referenced by update().

fillLV()

void PrintGeomSummary::fillLV ( G4LogicalVolume *  lv )

private

Definition at line 140 of file PrintGeomSummary.cc.

                                                 {
  if (std::find(lvs_.begin(), lvs_.end(), lv) == lvs_.end())
    lvs_.push_back(lv);
  G4VSolid* sl = lv->GetSolid();
  if (std::find(sls_.begin(), sls_.end(), sl) == sls_.end())
    sls_.push_back(sl);
  touch_.push_back(lv);
  for (int ii = 0; ii < (int)(lv->GetNoDaughters()); ii++)
    fillLV(lv->GetDaughter(ii)->GetLogicalVolume());
}

References spr::find(), cuy::ii, createfilelist::int, lvs_, sls_, and touch_.

Referenced by update().

fillPV()

void PrintGeomSummary::fillPV ( G4VPhysicalVolume *  pv )

private

Definition at line 151 of file PrintGeomSummary.cc.

                                                   {
  if (std::find(pvs_.begin(), pvs_.end(), pv) == pvs_.end())
    pvs_.push_back(pv);
  for (int ii = 0; ii < (int)(pv->GetLogicalVolume()->GetNoDaughters()); ii++)
    fillPV(pv->GetLogicalVolume()->GetDaughter(ii));
}

References spr::find(), cuy::ii, createfilelist::int, MetAnalyzer::pv(), and pvs_.

Referenced by update().

getTopPV()

G4VPhysicalVolume * PrintGeomSummary::getTopPV ( )

private

Definition at line 190 of file PrintGeomSummary.cc.

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

Referenced by update().

printSummary()

void PrintGeomSummary::printSummary ( std::ostream &  out )

private

Definition at line 211 of file PrintGeomSummary.cc.

                                                   {
  int k(0);
  for (std::map<DDSolidShape, std::pair<int, int>>::iterator itr = kount_.begin(); itr != kount_.end(); ++itr, ++k) {
    std::string shape = solidShape_[itr->first];
    out << "Shape [" << k << "]  " << shape << " # " << (itr->second).first << " : " << (itr->second).second << G4endl;
  }
}

References dqmdumpme::first, dqmdumpme::k, kount_, genParticles_cff::map, MillePedeFileConverter_cfg::out, solidShape_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by dumpSummary().

update() [1/2]

void PrintGeomSummary::update ( const BeginOfJob *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 63 of file PrintGeomSummary.cc.

                                                   {
  edm::ESTransientHandle<DDCompactView> pDD;
  (*job)()->get<IdealGeometryRecord>().get(pDD);
  const DDCompactView* cpv = &(*pDD);
 
  const auto& gra = cpv->graph();
 
  using Graph = DDCompactView::Graph;
  using adjl_iterator = Graph::const_adj_iterator;
 
  Graph::index_type i = 0;
  solidMap_.clear();
  for (adjl_iterator git = gra.begin(); git != gra.end(); ++git) {
    const DDLogicalPart& ddLP = gra.nodeData(git);
    addSolid(ddLP);
    ++i;
    if (!git->empty()) {
      // ask for children of ddLP
      for (Graph::edge_list::const_iterator cit = git->begin(); cit != git->end(); ++cit) {
        const DDLogicalPart& ddcurLP = gra.nodeData(cit->first);
        addSolid(ddcurLP);
      }
    }
  }
  G4cout << "Finds " << solidMap_.size() << " different solids in the tree" << G4endl;
}

References addSolid(), ecalTB2006H4_GenSimDigiReco_cfg::G4cout, DDCompactView::graph(), mps_fire::i, and solidMap_.

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

update() [2/2]

void PrintGeomSummary::update ( const BeginOfRun *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 101 of file PrintGeomSummary.cc.

                                                   {
  theTopPV_ = getTopPV();
  if (theTopPV_) {
    lvs_.clear();
    sls_.clear();
    touch_.clear();
    fillLV(theTopPV_->GetLogicalVolume());
    std::string name = theTopPV_->GetName();
    dumpSummary(G4cout, name);
 
    pvs_.clear();
    fillPV(theTopPV_);
    G4cout << " Number of G4VPhysicalVolume's for " << name << ": " << pvs_.size() << G4endl;
 
    for (unsigned int k = 0; k < nodeNames_.size(); ++k) {
      const G4LogicalVolumeStore* lvs = G4LogicalVolumeStore::GetInstance();
      std::vector<G4LogicalVolume*>::const_iterator lvcite;
      for (lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++) {
        if ((*lvcite)->GetName() == (G4String)(nodeNames_[k])) {
          lvs_.clear();
          sls_.clear();
          touch_.clear();
          fillLV(*lvcite);
          dumpSummary(G4cout, nodeNames_[k]);
        }
      }
      const G4PhysicalVolumeStore* pvs = G4PhysicalVolumeStore::GetInstance();
      std::vector<G4VPhysicalVolume*>::const_iterator pvcite;
      for (pvcite = pvs->begin(); pvcite != pvs->end(); pvcite++) {
        if ((*pvcite)->GetName() == (G4String)(nodeNames_[k])) {
          pvs_.clear();
          fillPV(*pvcite);
          G4cout << " Number of G4VPhysicalVolume's for " << nodeNames_[k] << ": " << pvs_.size() << G4endl;
        }
      }
    }
  }
}

References dumpSummary(), fillLV(), fillPV(), ecalTB2006H4_GenSimDigiReco_cfg::G4cout, getTopPV(), dqmdumpme::k, lvs_, Skims_PA_cff::name, nodeNames_, FSQDQM_cfi::pvs, pvs_, sls_, AlCaHLTBitMon_QueryRunRegistry::string, theTopPV_, and touch_.

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

Member Data Documentation

kount_

std::map<DDSolidShape, std::pair<int, int> > PrintGeomSummary::kount_

private

Definition at line 48 of file PrintGeomSummary.h.

Referenced by addName(), dumpSummary(), and printSummary().

lvs_

std::vector<G4LogicalVolume *> PrintGeomSummary::lvs_

private

Definition at line 45 of file PrintGeomSummary.h.

Referenced by dumpSummary(), fillLV(), and update().

nodeNames_

std::vector<std::string> PrintGeomSummary::nodeNames_

private

Definition at line 41 of file PrintGeomSummary.h.

Referenced by PrintGeomSummary(), and update().

pvs_

std::vector<G4VPhysicalVolume *> PrintGeomSummary::pvs_

private

Definition at line 47 of file PrintGeomSummary.h.

Referenced by fillPV(), and update().

sls_

std::vector<G4VSolid *> PrintGeomSummary::sls_

private

Definition at line 46 of file PrintGeomSummary.h.

Referenced by dumpSummary(), fillLV(), and update().

solidMap_

std::map<std::string, DDSolidShape> PrintGeomSummary::solidMap_

private

Definition at line 43 of file PrintGeomSummary.h.

Referenced by addName(), addSolid(), and update().

solidShape_

std::map<DDSolidShape, std::string> PrintGeomSummary::solidShape_

private

Definition at line 42 of file PrintGeomSummary.h.

Referenced by addSolid(), PrintGeomSummary(), and printSummary().

theTopPV_

G4VPhysicalVolume* PrintGeomSummary::theTopPV_

private

Definition at line 44 of file PrintGeomSummary.h.

Referenced by update().

touch_

std::vector<G4LogicalVolume *> PrintGeomSummary::touch_

private

Definition at line 45 of file PrintGeomSummary.h.

Referenced by dumpSummary(), fillLV(), and update().