PrintG4Touch Class Reference

Inheritance diagram for PrintG4Touch:

Public Member Functions
	PrintG4Touch (edm::ParameterSet const &p)
	~PrintG4Touch () override=default
Public Member Functions inherited from SimWatcher
virtual void	beginRun (edm::EventSetup const &)
bool	isMT () const
const SimWatcher &	operator= (const SimWatcher &)=delete
virtual void	registerConsumes (edm::ConsumesCollector)
	SimWatcher ()
	SimWatcher (const SimWatcher &)=delete
virtual	~SimWatcher ()
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	dumpSummary (std::ostream &out=G4cout)
void	dumpTouch (G4VPhysicalVolume *pv, unsigned int leafDepth, std::ostream &out=G4cout)
G4LogicalVolume *	getTopLV ()
G4VPhysicalVolume *	getTopPV ()
void	getTouch (G4VPhysicalVolume *pv, unsigned int leafDepth, unsigned int copym, std::vector< std::string > &touches)
void	update (const BeginOfRun *run) override
	This routine will be called when the appropriate signal arrives. More...

Private Attributes
bool	dd4hep_
G4NavigationHistory	fHistory_
G4VPhysicalVolume *	theTopPV_
bool	verbosity_

Additional Inherited Members
Protected Member Functions inherited from SimWatcher
void	setMT (bool val)

Detailed Description

Definition at line 35 of file PrintG4Touch.cc.

Constructor & Destructor Documentation

PrintG4Touch()

PrintG4Touch::PrintG4Touch

(

edm::ParameterSet const &

p

)

Definition at line 56 of file PrintG4Touch.cc.

References dd4hep_, ecalTB2006H4_GenSimDigiReco_cfg::G4cout, AlCaHLTBitMon_ParallelJobs::p, and verbosity_.

                                                   {
  dd4hep_ = p.getUntrackedParameter<bool>("dd4hep", false);
  verbosity_ = p.getUntrackedParameter<bool>("verbosity", false);
  G4cout << "PrintG4Touch:: initialised for dd4hep " << dd4hep_ << " with verbosity levels:" << verbosity_ << G4endl;
}

~PrintG4Touch()

PrintG4Touch::~PrintG4Touch ( )

overridedefault

Member Function Documentation

add1touchable()

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

private

Definition at line 112 of file PrintG4Touch.cc.

References cuy::ii.

Referenced by countNoTouchables().

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

countNoTouchables()

int PrintG4Touch::countNoTouchables ( )

private

Definition at line 105 of file PrintG4Touch.cc.

References add1touchable(), and getTopLV().

Referenced by dumpSummary().

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

dumpSummary()

void PrintG4Touch::dumpSummary ( std::ostream &  out = G4cout )

private

Definition at line 79 of file PrintG4Touch.cc.

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

Referenced by update().

                                              {
  //---------- 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: "
      << DD4hep2DDDName::namePV(static_cast<std::string>(theTopPV_->GetName()), dd4hep_) << 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;
  const G4RegionStore *regs = G4RegionStore::GetInstance();
  out << " Number of G4Region's: " << regs->size() << G4endl;
}

dumpTouch()

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

private

Definition at line 118 of file PrintG4Touch.cc.

References angle_units::operators::convertRadToDeg(), dd4hep_, fHistory_, DD4hep2DDDName::nameSolid(), MillePedeFileConverter_cfg::out, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by update().

                                                                                           {
  if (leafDepth == 0)
    fHistory_.SetFirstEntry(pv);
  else
    fHistory_.NewLevel(pv, kNormal, pv->GetCopyNo());

  G4LogicalVolume *lv = pv->GetLogicalVolume();

  G4ThreeVector globalpoint = fHistory_.GetTopTransform().Inverse().TransformPoint(G4ThreeVector(0, 0, 0));
  std::string mother = (pv->GetMotherLogical())
                           ? (DD4hep2DDDName::nameSolid(
                                 static_cast<std::string>(pv->GetMotherLogical()->GetSolid()->GetName()), dd4hep_))
                           : "World";
  std::string lvname = DD4hep2DDDName::nameSolid(static_cast<std::string>(lv->GetName()), dd4hep_);
  out << leafDepth << "### VOLUME = " << lvname << " 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();
}

getTopLV()

G4LogicalVolume * PrintG4Touch::getTopLV ( )

private

Definition at line 185 of file PrintG4Touch.cc.

References theTopPV_.

Referenced by countNoTouchables().

{ return theTopPV_->GetLogicalVolume(); }

getTopPV()

G4VPhysicalVolume * PrintG4Touch::getTopPV ( )

private

Definition at line 181 of file PrintG4Touch.cc.

Referenced by update().

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

getTouch()

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

private

Definition at line 147 of file PrintG4Touch.cc.

References filterCSVwithJSON::copy, dd4hep_, fHistory_, Skims_PA_cff::name, DD4hep2DDDName::nameSolid(), AlCaHLTBitMon_QueryRunRegistry::string, and to_string().

Referenced by update().

                                                           {
  if (leafDepth == 0)
    fHistory_.SetFirstEntry(pv);
  else
    fHistory_.NewLevel(pv, kNormal, pv->GetCopyNo());

  std::string mother = (pv->GetMotherLogical())
                           ? (DD4hep2DDDName::nameSolid(
                                 static_cast<std::string>(pv->GetMotherLogical()->GetSolid()->GetName()), dd4hep_))
                           : "World";

  G4LogicalVolume *lv = pv->GetLogicalVolume();
  std::string lvname = DD4hep2DDDName::nameSolid(static_cast<std::string>(lv->GetSolid()->GetName()), dd4hep_);
  unsigned int copy = static_cast<unsigned int>(pv->GetCopyNo());

  std::string type = static_cast<std::string>(lv->GetSolid()->GetEntityType());

  std::string name = lvname + " " + std::to_string(copy) + " " + mother + " " + std::to_string(copym) + " " + type;
  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();
}

update()

void PrintG4Touch::update ( const BeginOfRun *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun *>.

Definition at line 62 of file PrintG4Touch.cc.

References dumpSummary(), dumpTouch(), ecalTB2006H4_GenSimDigiReco_cfg::G4cout, getTopPV(), getTouch(), jetUpdater_cfi::sort, theTopPV_, and verbosity_.

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

                                               {
  //Now take action
  theTopPV_ = getTopPV();

  dumpSummary(G4cout);

  std::vector<std::string> touches;
  getTouch(theTopPV_, 0, 1, touches);
  std::sort(touches.begin(), touches.end());
  for (const auto &touch : touches)
    G4cout << touch << G4endl;

  //---------- Dump LV and PV information
  if (verbosity_)
    dumpTouch(theTopPV_, 0, G4cout);
}

Member Data Documentation

dd4hep_

bool PrintG4Touch::dd4hep_

private

Definition at line 51 of file PrintG4Touch.cc.

Referenced by dumpSummary(), dumpTouch(), getTouch(), and PrintG4Touch().

fHistory_

G4NavigationHistory PrintG4Touch::fHistory_

private

Definition at line 53 of file PrintG4Touch.cc.

Referenced by dumpTouch(), and getTouch().

theTopPV_

G4VPhysicalVolume* PrintG4Touch::theTopPV_

private

Definition at line 52 of file PrintG4Touch.cc.

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

verbosity_

bool PrintG4Touch::verbosity_

private

Definition at line 51 of file PrintG4Touch.cc.

Referenced by PrintG4Touch(), and update().