larg4::AuxDetReadoutGeometry Class Reference

#include "AuxDetReadoutGeometry.h"

Inheritance diagram for larg4::AuxDetReadoutGeometry:

Public Member Functions
	AuxDetReadoutGeometry (geo::AuxDetGeometryCore const *auxDetGeom, G4String name="AuxDetReadoutGeometry")
void	Construct () override

Private Member Functions
void	FindAndMakeAuxDet (std::vector< const G4VPhysicalVolume * > &path, unsigned int depth, G4Transform3D DepthToWorld)
void	FindAndMakeAuxDetSensitive (std::vector< const G4VPhysicalVolume * > &path, unsigned int depth, G4Transform3D DepthToWorld)

Private Attributes
geo::AuxDetGeometryCore const *	fAuxDetGeom
	Handle to the geometry. More...
uint32_t	fNumSensitiveVol
	number of sensitive volumes More...

Detailed Description

Definition at line 39 of file AuxDetReadoutGeometry.h.

Constructor & Destructor Documentation

larg4::AuxDetReadoutGeometry::AuxDetReadoutGeometry	(	geo::AuxDetGeometryCore const *	auxDetGeom,
		G4String	name = "AuxDetReadoutGeometry"
	)

Definition at line 25 of file AuxDetReadoutGeometry.cxx.

References fNumSensitiveVol.

    : G4VUserParallelWorld(name), fAuxDetGeom{auxDetGeom}, fNumSensitiveVol(0)
  {}

Member Function Documentation

void larg4::AuxDetReadoutGeometry::Construct ( )

override

The key method in this class; creates a parallel world view of those volumes relevant to the auxiliary detector readout. Required of any class that inherits from G4VUserParallelWorld

Definition at line 31 of file AuxDetReadoutGeometry.cxx.

References FindAndMakeAuxDet(), FindAndMakeAuxDetSensitive(), fNumSensitiveVol, and g4b::DetectorConstruction::GetWorld().

  {
    // We want to find all of the AuxDets that the Geometry service would find and make
    // each one a G4 sensitive detector.

    // Call initial case of a function that will rucursively run through the volume tree
    // down to max depth. Start at 0 depth with World, where the initial translation and
    // rotation should be 0 as well
    unsigned int MaxDepth = 8; // should be plenty
    std::vector<const G4VPhysicalVolume*> path(MaxDepth);
    path[0] = g4b::DetectorConstruction::GetWorld();
    G4Transform3D InitTransform(path[0]->GetObjectRotationValue(), path[0]->GetObjectTranslation());

    // first try to make sensitive volumes, if those are not in the gdml file (ie file was
    // made before the introduction of sensitive volumes) fall back to looking for the aux
    // dets
    FindAndMakeAuxDetSensitive(path, 0, InitTransform);
    if (fNumSensitiveVol < 1) FindAndMakeAuxDet(path, 0, InitTransform);
  }

void larg4::AuxDetReadoutGeometry::FindAndMakeAuxDet ( std::vector< const G4VPhysicalVolume * > &  path, unsigned int  depth, G4Transform3D  DepthToWorld  )

private

Definition at line 109 of file AuxDetReadoutGeometry.cxx.

References d, fAuxDetGeom, geo::AuxDetGeometryCore::FindAuxDetAtPosition(), fNumSensitiveVol, MF_LOG_DEBUG, and util::to_string().

Referenced by Construct().

  {
    G4LogicalVolume* LogicalVolumeAtDepth = path[depth]->GetLogicalVolume();

    std::string volName(path[depth]->GetName());
    if (volName.find("volAuxDet") != std::string::npos) {

      // find world coordinate of the AuxDet origin in cm
      G4Point3D local(0., 0., 0.);
      G4Point3D world = DepthToWorld * local; // G4 works in mm
      geo::Point_t const worldPos{
        world.x() / CLHEP::cm, world.y() / CLHEP::cm, world.z() / CLHEP::cm};

      auto const adNum = fAuxDetGeom->FindAuxDetAtPosition(worldPos);
      //  N.B. This name is expected by code in LArG4:
      std::string SDName = "AuxDetSD_AuxDet" + std::to_string(adNum) + "_0";
      AuxDetReadout* adReadout = new larg4::AuxDetReadout(SDName, adNum, 0);

      MF_LOG_DEBUG("AuxDetReadoutGeometry")
        << "found" << path[depth]->GetName() << ", number " << adNum << ":0";

      // Tell Geant4's sensitive-detector manager about the AuxDetReadout class
      (G4SDManager::GetSDMpointer())->AddNewDetector(adReadout);
      LogicalVolumeAtDepth->SetSensitiveDetector(adReadout);
      ++fNumSensitiveVol;
      return;
    }

    // Explore the next layer down -- unless it is a very deep geometry,
    // recursion should end before exception is thrown.
    unsigned int deeper = depth + 1;
    if (deeper >= path.size()) {
      throw cet::exception("AuxDetReadoutGeometry") << "exceeded maximum TGeoNode depth\n";
    }

    // Note that there will be nd different branches off of path[depth]
    G4int nd = LogicalVolumeAtDepth->GetNoDaughters();
    for (int d = 0; d < nd; ++d) {

      // get the physvol daughter in the logicalvol
      path[deeper] = LogicalVolumeAtDepth->GetDaughter(d);

      // keep track of the transform to world coordinates for PositionToAuxDet
      G4Transform3D DeeperToMother(path[deeper]->GetObjectRotationValue(),
                                   path[deeper]->GetObjectTranslation());
      G4Transform3D DeeperToWorld = DepthToWorld * DeeperToMother;
      FindAndMakeAuxDet(path, deeper, DeeperToWorld);
    }
  }

void larg4::AuxDetReadoutGeometry::FindAndMakeAuxDetSensitive ( std::vector< const G4VPhysicalVolume * > &  path, unsigned int  depth, G4Transform3D  DepthToWorld  )

private

Definition at line 52 of file AuxDetReadoutGeometry.cxx.

References d, fAuxDetGeom, geo::AuxDetGeometryCore::FindAuxDetSensitiveAtPosition(), fNumSensitiveVol, MF_LOG_DEBUG, and util::to_string().

Referenced by Construct().

  {
    G4LogicalVolume* LogicalVolumeAtDepth = path[depth]->GetLogicalVolume();

    std::string volName(path[depth]->GetName());
    if (volName.find("volAuxDet") != std::string::npos &&
        volName.find("Sensitive") != std::string::npos) {

      // find world coordinate of the AuxDet origin in cm
      G4Point3D local(0., 0., 0.);
      G4Point3D world = DepthToWorld * local; // G4 works in mm
      geo::Point_t const worldPos{
        world.x() / CLHEP::cm, world.y() / CLHEP::cm, world.z() / CLHEP::cm};

      size_t adNum = 0;
      size_t svNum = 0;
      fAuxDetGeom->FindAuxDetSensitiveAtPosition(worldPos, adNum, svNum);
      //  N.B. This name is expected by code in LArG4:
      std::string SDName = "AuxDetSD_AuxDet" + std::to_string(adNum) + "_" + std::to_string(svNum);
      AuxDetReadout* adReadout = new larg4::AuxDetReadout(SDName, adNum, svNum);

      MF_LOG_DEBUG("AuxDetReadoutGeometry")
        << "found" << path[depth]->GetName() << ", number " << adNum << ":" << svNum;

      // Tell Geant4's sensitive-detector manager about the AuxDetReadout class
      (G4SDManager::GetSDMpointer())->AddNewDetector(adReadout);
      LogicalVolumeAtDepth->SetSensitiveDetector(adReadout);
      ++fNumSensitiveVol;
      return;
    }

    // Explore the next layer down -- unless it is a very deep geometry,
    // recursion should end before exception is thrown.
    unsigned int deeper = depth + 1;
    if (deeper >= path.size()) {
      throw cet::exception("AuxDetReadoutGeometry") << "exceeded maximum TGeoNode depth\n";
    }

    // Note that there will be nd different branches off of path[depth]
    G4int nd = LogicalVolumeAtDepth->GetNoDaughters();
    for (int d = 0; d < nd; ++d) {

      // get the physvol daughter in the logicalvol
      path[deeper] = LogicalVolumeAtDepth->GetDaughter(d);

      // keep track of the transform to world coordinates for PositionToAuxDet
      G4Transform3D DeeperToMother(path[deeper]->GetObjectRotationValue(),
                                   path[deeper]->GetObjectTranslation());
      G4Transform3D DeeperToWorld = DepthToWorld * DeeperToMother;
      FindAndMakeAuxDetSensitive(path, deeper, DeeperToWorld);
    }
  }

Member Data Documentation

geo::AuxDetGeometryCore const* larg4::AuxDetReadoutGeometry::fAuxDetGeom

private

Handle to the geometry.

Definition at line 58 of file AuxDetReadoutGeometry.h.

Referenced by FindAndMakeAuxDet(), and FindAndMakeAuxDetSensitive().

uint32_t larg4::AuxDetReadoutGeometry::fNumSensitiveVol

private

number of sensitive volumes

Definition at line 59 of file AuxDetReadoutGeometry.h.

Referenced by AuxDetReadoutGeometry(), Construct(), FindAndMakeAuxDet(), and FindAndMakeAuxDetSensitive().

---

The documentation for this class was generated from the following files:

• larsim/v10_02_02/source/larsim/LegacyLArG4/AuxDetReadoutGeometry.h
• larsim/v10_02_02/source/larsim/LegacyLArG4/AuxDetReadoutGeometry.cxx