larg4::AuxDetReadoutGeometry Class Reference

#include <AuxDetReadoutGeometry.h>

Inheritance diagram for larg4::AuxDetReadoutGeometry:

Public Member Functions
	AuxDetReadoutGeometry (const G4String name="AuxDetReadoutGeometry")
	Constructor and destructor. More...
virtual	~AuxDetReadoutGeometry ()
virtual void	Construct ()

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
art::ServiceHandle < geo::Geometry const >	fGeo
	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

(

const G4String

name = "AuxDetReadoutGeometry"

)

Constructor and destructor.

Definition at line 24 of file AuxDetReadoutGeometry.cxx.

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

larg4::AuxDetReadoutGeometry::~AuxDetReadoutGeometry ( )

virtual

Definition at line 28 of file AuxDetReadoutGeometry.cxx.

  {}

Member Function Documentation

void larg4::AuxDetReadoutGeometry::Construct ( )

virtual

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 32 of file AuxDetReadoutGeometry.cxx.

  {

    // 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
    this->FindAndMakeAuxDetSensitive(path, 0, InitTransform);
    if(fNumSensitiveVol < 1)
      this->FindAndMakeAuxDet(path, 0, InitTransform);


    return;
  }// end Construct

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

private

Definition at line 116 of file AuxDetReadoutGeometry.cxx.

  {

    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
      double worldPos[3] = { world.x() / CLHEP::cm, world.y() / CLHEP::cm, world.z() / CLHEP::cm };

      unsigned int adNum;
      fGeo->PositionToAuxDet(worldPos, adNum);
      //  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;
      this->FindAndMakeAuxDet(path, deeper, DeeperToWorld);
    }

  }

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

private

Definition at line 59 of file AuxDetReadoutGeometry.cxx.

  {

    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
      double worldPos[3] = { world.x() / CLHEP::cm, world.y() / CLHEP::cm, world.z() / CLHEP::cm };

      size_t adNum = 0;
      size_t svNum = 0;
      fGeo->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;
      this->FindAndMakeAuxDetSensitive(path, deeper, DeeperToWorld);
    }

  }

Member Data Documentation

art::ServiceHandle<geo::Geometry const> larg4::AuxDetReadoutGeometry::fGeo

private

Handle to the geometry.

Definition at line 61 of file AuxDetReadoutGeometry.h.

uint32_t larg4::AuxDetReadoutGeometry::fNumSensitiveVol

private

number of sensitive volumes

Definition at line 62 of file AuxDetReadoutGeometry.h.

---

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

• AuxDetReadoutGeometry.h
• AuxDetReadoutGeometry.cxx