#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)

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< 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 ( edm::ParameterSet const & p )

Definition at line 32 of file PrintGeomSummary.cc.

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, edm::ParameterSet::getUntrackedParameter(), cuy::ii, nodeNames_, and solidShape_.

                                                            : 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";
 }

PrintGeomSummary::~PrintGeomSummary ( )

override

Definition at line 62 of file PrintGeomSummary.cc.

62 {}

Member Function Documentation

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

private

Definition at line 179 of file PrintGeomSummary.cc.

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

Referenced by dumpSummary().

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

void PrintGeomSummary::addSolid ( const DDLogicalPart & part )

private

Definition at line 95 of file PrintGeomSummary.cc.

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

Referenced by update().

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

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

private

Definition at line 136 of file PrintGeomSummary.cc.

References addName(), kount_, lvs_, printSummary(), sls_, AlCaHLTBitMon_QueryRunRegistry::string, and touch_.

Referenced by update().

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

void PrintGeomSummary::fillLV ( G4LogicalVolume * lv )

private

Definition at line 126 of file PrintGeomSummary.cc.

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

Referenced by update().

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

G4VPhysicalVolume * PrintGeomSummary::getTopPV ( )

private

Definition at line 174 of file PrintGeomSummary.cc.

Referenced by update().

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

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

private

Definition at line 197 of file PrintGeomSummary.cc.

References plotBeamSpotDB::first, gen::k, kount_, genParticles_cff::map, solidShape_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by dumpSummary().

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

void PrintGeomSummary::update ( const BeginOfJob * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 64 of file PrintGeomSummary.cc.

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

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

void PrintGeomSummary::update ( const BeginOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 104 of file PrintGeomSummary.cc.

References dumpSummary(), fillLV(), ecalTB2006H4_GenSimDigiReco_cfg::G4cout, getTopPV(), gen::k, lvs_, dataset::name, nodeNames_, sls_, AlCaHLTBitMon_QueryRunRegistry::string, theTopPV_, and touch_.

                                                     {
   theTopPV_ = getTopPV();
   if (theTopPV_) {
     lvs_.clear(); sls_.clear(); touch_.clear();
     fillLV(theTopPV_->GetLogicalVolume()); 
     std::string name = theTopPV_->GetName();
     dumpSummary(G4cout,name);
 
     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]);
     }
       }
     }
   }
 }