#include "MuonWall/MuonWallTool.h"

Inheritance diagram for MuonWallTool:

Collaboration diagram for MuonWallTool:

Public Member Functions
	MuonWallTool (const std::string &type, const std::string &name, const IInterface *parent)

	~MuonWallTool ()

virtual void	BuildGeometry () override final
	Override DetectorGeometryBase::BuildGeometry method. More...

StatusCode	initialize () override
	Athena method. More...

virtual void	Build () override
	purely virtual methods being implemented here More...

virtual void	PositionInParent () override

virtual void	SetRotationAndOffset ()

virtual void	BuildSubDetectors () override

virtual void	SetEnvelope () override

virtual void	ResetEnvelope () override

virtual void	SetAsWorld () override

std::string	GetDetectorName () const override

void	SetDetectorName (const std::string &) override

void	SetParent (IDetectorGeometryTool *) override

Envelope &	GetEnvelope () override final

G4VPhysicalVolume *	GetWorldVolume () override final

Protected Attributes
ToolHandleArray< IDetectorGeometryTool >	m_subDetTools {this, "SubDetectors", {}, "Tool handle array of all subdetector tools"}

ServiceHandle< IG4GeometryNotifierSvc >	m_notifierSvc {this, "GeometryNotifierSvc", "G4GeometryNotifierSvc", "Detector name (same as the Tool name if not set"}

IDetectorGeometryTool *	m_theParent {}

Envelope	m_envelope

Gaudi::Property< std::string >	m_detectorName {this, "DetectorName", "", "Detector name (same as the Tool name if not set"}

bool	m_isWorld {false}

Gaudi::Property< double >	m_rotateX {this, "RotateX" , 0.0, "Rotation around the X-axis"}

Gaudi::Property< double >	m_rotateY {this, "RotateY" , 0.0, "Rotation around the Y-axis"}

Gaudi::Property< double >	m_rotateZ {this, "RotateZ" , 0.0, "Rotation around the Z-axis"}

Gaudi::Property< double >	m_offsetX {this, "OffsetX" , 0.0, "Offset in the X-direction"}

Gaudi::Property< double >	m_offsetY {this, "OffsetY" , 0.0, "Offset in the Y-direction"}

Gaudi::Property< double >	m_offsetZ {this, "OffsetZ" , 0.0, "Offset in the Z-direction"}

Private Attributes
double	m_zLength

double	m_yLength

double	m_xLength

bool	m_backWall

bool	m_sideWall

Detailed Description

Tool for building the MuonWall detector.

Definition at line 19 of file MuonWallTool.h.

Constructor & Destructor Documentation

◆ MuonWallTool()

MuonWallTool::MuonWallTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Definition at line 34 of file MuonWallTool.cxx.

   : DetectorGeometryBase(type,name,parent),
     m_zLength(0.),
     m_yLength(0.),
     m_xLength(0.),
     m_backWall(true),
     m_sideWall(false)
 {
   ATH_MSG_DEBUG( "MuonWallTool constructor for " << name );
   declareProperty("ZLength", m_zLength, "");
   declareProperty("YLength", m_yLength, "");
   declareProperty("XLength", m_xLength, "");
   declareProperty("backWall", m_backWall, "");
   declareProperty("sideWall", m_sideWall, "");
 }

◆ ~MuonWallTool()

MuonWallTool::~MuonWallTool ( )

Definition at line 50 of file MuonWallTool.cxx.

51 {

52 }

Member Function Documentation

◆ Build()

void DetectorGeometryBase::Build ( )

overridevirtualinherited

purely virtual methods being implemented here

Definition at line 53 of file DetectorGeometryBase.cxx.

 {
   ATH_MSG_VERBOSE( name() << "::Build() (Base class method): Starting. Number of registered volumes "<<G4LogicalVolumeStore::GetInstance()->size() );
   
   SetEnvelope();
   ATH_MSG_VERBOSE( name() << "::Build() - Envelope set. Number of registered volumes "<<G4LogicalVolumeStore::GetInstance()->size() );
  
   m_notifierSvc->SetCurrentDetectorName(m_detectorName.value());
   BuildGeometry();
   ATH_MSG_VERBOSE( name() << "::Build() - Geometry built. Number of registered volumes "<<G4LogicalVolumeStore::GetInstance()->size() );
  
   SetRotationAndOffset();
   ATH_MSG_VERBOSE( name() << "::Build() - Volume moved around. Number of registered volumes "<<G4LogicalVolumeStore::GetInstance()->size() );
  
   PositionInParent();
   ATH_MSG_VERBOSE( name() << "::Build() - Connected with parent. Number of registered volumes "<<G4LogicalVolumeStore::GetInstance()->size() );
  
   BuildSubDetectors();
   ATH_MSG_VERBOSE( name() << "::Build() (Base class method): Finished. Number of registered volumes "<<G4LogicalVolumeStore::GetInstance()->size() );
   return;
 }

◆ BuildGeometry()

void MuonWallTool::BuildGeometry ( )

finaloverridevirtual

Override DetectorGeometryBase::BuildGeometry method.

Reimplemented from DetectorGeometryBase.

Definition at line 54 of file MuonWallTool.cxx.

                                  {
  
   ATH_MSG_DEBUG( "Building Geometry back muon wall: "
                  << (m_backWall ? "true" : "false")
                  << "  side muon wall: "
                  << (m_sideWall ? "true" : "false") );
  
   // MuonWall description :
  
   // Scintillator : 400*200*20, calculating half-size
   const double zScintillator = m_zLength / 2.;
   const double yScintillator = m_yLength / 2.;
   const double xScintillator = m_xLength / 2.;
   const double dzmuonwall = 750. * CLHEP::mm;
   const double dymuonwall = 425. * CLHEP::mm;
   const double dxmuonwall = xScintillator;
  
   G4Box *wall = new G4Box("MuonWall", dxmuonwall, dymuonwall, dzmuonwall);
  
   // Get the materials
   StoredMaterialManager* materialManager = nullptr;
   if (StatusCode::SUCCESS != detStore()->retrieve(materialManager, std::string("MATERIALS"))) {
     ATH_MSG_ERROR( "Could not find Material Manager MATERIALS" );
     return;
   }
  
   const GeoMaterial *geoAir = materialManager->getMaterial("tile::Air");
   const GeoMaterial *geoScintillator = materialManager->getMaterial("tile::Scintillator");
   Geo2G4MaterialFactory theMaterialFactory;
   G4Material *airMaterial = theMaterialFactory.Build(geoAir);
   G4Material *scintillatorMaterial = theMaterialFactory.Build(geoScintillator);
  
   G4LogicalVolume *wallV = new G4LogicalVolume(wall, airMaterial, "MuonWall");
  
   // ------- Create lead layers and place them inside the mother box --------
   double zLayer, yLayer, xLayer = 0.;
  
   if (m_backWall) {
  
     const int nrOfLayers = 12;
  
     for (int j = 0; j < nrOfLayers; j++) {
  
       // Scintillator
       G4Box *scintillatorLayer = new G4Box("ScintillatorLayer", xScintillator, yScintillator, zScintillator);
  
       G4LogicalVolume *scintillatorLayerV = new G4LogicalVolume(scintillatorLayer, scintillatorMaterial,
           "MuScintillatorLayer");
  
       // scintillatorLayerV->GetLogicalVolume()->SetSensitiveDetector(muonwallSD);
  
       if (j < 6) {
         yLayer = yScintillator + 25. * CLHEP::mm;
         zLayer = -450. * CLHEP::mm + (2 * j + 1) * zScintillator;
       } else {
         yLayer = -yScintillator + 25. * CLHEP::mm;
         zLayer = -450. * CLHEP::mm + (2 * (j - 6) + 1) * zScintillator;
       }
  
       G4PVPlacement* scintillatorLayerVPhys __attribute__((unused)) =
           new G4PVPlacement(0,
                             G4ThreeVector(xLayer,yLayer,zLayer),
                             scintillatorLayerV,
                             "MuScintillatorLayer",
                             wallV,
                             false,
                             j+1);
       }
   }
  
   if (m_sideWall) {
  
     // Scintillator : 500*100*20 instead of 400*200*20
     const double yScintillator1 = yScintillator + 50. * CLHEP::mm;
     const double zScintillator1 = zScintillator - 50. * CLHEP::mm;
  
     const int nScintLayers = 3;
  
     for (int j = 0; j < nScintLayers; j++) {
  
       G4Box* scintillatorLayer = new G4Box("SideScintiLayer", xScintillator, yScintillator1, zScintillator1);
  
       G4LogicalVolume* scintillatorLayerV = new G4LogicalVolume(scintillatorLayer, scintillatorMaterial,
           "SideMuScintiLayer");
  
       yLayer = -175. * CLHEP::mm;
       zLayer = -750. * CLHEP::mm + (2 * j + 1) * zScintillator1;
  
       G4PVPlacement* scintillatorLayerVPhys __attribute__((unused)) =
           new G4PVPlacement(0,
                             G4ThreeVector(xLayer,yLayer,zLayer),
                             scintillatorLayerV,
                             "SideMuScintiLayer",
                             wallV,
                             false,
                             j+1);
       }
   }
  
   // FINAL STEP
   m_envelope.theEnvelope = wallV;
 }

◆ BuildSubDetectors()

void DetectorGeometryBase::BuildSubDetectors ( )

overridevirtualinherited

Definition at line 155 of file DetectorGeometryBase.cxx.

 {
   ATH_MSG_VERBOSE( name() << "::BuildSubDetectors() (Base class method): Starting");
   for (auto& subDetTool: m_subDetTools)
     {
       ATH_MSG_VERBOSE(name() << "::BuildSubDetectors() (Base class method):  Positioning "<<subDetTool->GetDetectorName()<<" within "<<m_detectorName.value());
       subDetTool->SetParent(this);
       subDetTool->Build();
     }
   ATH_MSG_VERBOSE( name() << "::BuildSubDetectors() (Base class method): Finished");
 }

◆ GetDetectorName()

std::string DetectorGeometryBase::GetDetectorName ( ) const

overrideinherited

Definition at line 180 of file DetectorGeometryBase.cxx.

 {
   return m_detectorName.value();
 }

◆ GetEnvelope()

Envelope & DetectorGeometryBase::GetEnvelope ( )

finaloverrideinherited

Definition at line 195 of file DetectorGeometryBase.cxx.

 {
   return m_envelope;
 }

◆ GetWorldVolume()

G4VPhysicalVolume * DetectorGeometryBase::GetWorldVolume ( )

finaloverrideinherited

Definition at line 200 of file DetectorGeometryBase.cxx.

 {
   if (m_isWorld)
     return m_envelope.thePositionedVolume;
   else
     {
       ATH_MSG_ERROR("trying to get World from a DetectorTool which World is not!");
       return 0;
     }
 }

◆ initialize()

StatusCode DetectorGeometryBase::initialize ( )

overrideinherited

Athena method.

called at initialization time, being customized here

Definition at line 19 of file DetectorGeometryBase.cxx.

 {
   ATH_MSG_VERBOSE( name() << "::initialize(): starting." );
   if(m_detectorName.empty())
     {
       m_detectorName = this->name();
       // re-initialize m_detectorName in order to take the real detector name rather than the path to it
       size_t ipos=m_detectorName.value().find_last_of('.');
       size_t length=m_detectorName.value().size();
       if (ipos<length)
         {
           ATH_MSG_VERBOSE( "m_detectorName: " << m_detectorName.value() << " needs to be reset.");
           m_detectorName=m_detectorName.value().substr(ipos+1,length-ipos-1);
           ATH_MSG_VERBOSE( "m_detectorName default value reset to " << m_detectorName.value());
         }
     }
   ATH_MSG_DEBUG( name() << "::initialize() (Base class method): Detector name = " << m_detectorName.value() );
   CHECK(m_notifierSvc.retrieve());
  
   //  This fires initialize() for each of those tools
   if (m_subDetTools.size())
     {
       ATH_MSG_DEBUG( name() << "::initialize(): Initializing list of " << m_subDetTools.size() << " detectors." );
       CHECK( m_subDetTools.retrieve() );
     }
   else
     {
       ATH_MSG_DEBUG( name() << "::initialize(): no sub-detectors to initialize." );
     }
  
   ATH_MSG_VERBOSE( name() << "::initialize(): finished." );
   return StatusCode::SUCCESS;
 }

◆ PositionInParent()

void DetectorGeometryBase::PositionInParent ( )

overridevirtualinherited

Reimplemented in GeoDetectorTool.

Definition at line 120 of file DetectorGeometryBase.cxx.

 {
   ATH_MSG_VERBOSE( name() << "::PositionInParent() (Base class method)");
   if (m_isWorld)
     {
       // check that the detector is built
       if (m_envelope.IsBuilt())
         {
           G4VPhysicalVolume* physWorld= new G4PVPlacement(0,G4ThreeVector(),
                                                           m_envelope.theEnvelope,m_envelope.theEnvelope->GetName(),0,false,0,false);
           m_envelope.thePositionedVolume=physWorld;
         }
     }
   else
     {
       // check that there is a parent
       if (!m_theParent)
         {
           ATH_MSG_ERROR("Parent not set for "<<m_detectorName.value()<<"!!!!!!!!!!");
         }
       else
         {
           if (m_theParent->GetEnvelope().IsBuilt())
             {
               // G4VPhysicalVolume *physVol = new G4PVPlacement(0,G4ThreeVector(),
               //                                                m_envelope.theEnvelope,m_envelope.theEnvelope->GetName(),m_theParent->GetEnvelope().theEnvelope,false,0,false);
               // TODO: implement a possible rotation/displacement - something like this based on the old code?
               G4VPhysicalVolume *physVol = new G4PVPlacement(m_envelope.theRotation,m_envelope.thePosition,
                                                              m_envelope.theEnvelope,m_envelope.theEnvelope->GetName(),m_theParent->GetEnvelope().theEnvelope,false,0);
               m_envelope.thePositionedVolume=physVol;
             }
         }
     }
 }

◆ ResetEnvelope()

void DetectorGeometryBase::ResetEnvelope ( )

overridevirtualinherited

Definition at line 171 of file DetectorGeometryBase.cxx.

172 {

173 }

◆ SetAsWorld()

void DetectorGeometryBase::SetAsWorld ( )

overridevirtualinherited

Definition at line 185 of file DetectorGeometryBase.cxx.

 {
   m_isWorld=true;
 }

◆ SetDetectorName()

void DetectorGeometryBase::SetDetectorName ( const std::string & s )

overrideinherited

Definition at line 175 of file DetectorGeometryBase.cxx.

 {
   m_detectorName=s;
 }

◆ SetEnvelope()

void DetectorGeometryBase::SetEnvelope ( )

overridevirtualinherited

Definition at line 167 of file DetectorGeometryBase.cxx.

168 {

169 }

◆ SetParent()

void DetectorGeometryBase::SetParent ( IDetectorGeometryTool * p )

overrideinherited

Definition at line 190 of file DetectorGeometryBase.cxx.

 {
   m_theParent=p;
 }

◆ SetRotationAndOffset()

void DetectorGeometryBase::SetRotationAndOffset ( )

virtualinherited

Definition at line 80 of file DetectorGeometryBase.cxx.

 {
   ATH_MSG_VERBOSE( name() << "::SetRotationAndOffset() (Base class method)");
   // Firstly do the rotation
   if (!m_envelope.theRotation)
     {
       // m_envelope.theRotation is null, so create an identity
       // rotation first.
       // FIXME probably a neater way to do this part.
       m_envelope.theRotation=new G4RotationMatrix;
       // add the extra rotations.
       m_envelope.theRotation->rotateX(m_rotateX);
       m_envelope.theRotation->rotateY(m_rotateY);
       m_envelope.theRotation->rotateZ(m_rotateZ);
       if (m_envelope.thePositionedVolume)
         {
           // Override the rotation for m_envelope.thePositionedVolume.
           m_envelope.thePositionedVolume->SetRotation(m_envelope.theRotation);
         }
     }
   else
     {
       // m_envelope.theRotation already exists, so just add the
       // extra rotations.
       m_envelope.theRotation->rotateX(m_rotateX);
       m_envelope.theRotation->rotateY(m_rotateY);
       m_envelope.theRotation->rotateZ(m_rotateZ);
     }
   // Secondly add the additional position offset to the existing
   // m_envelope.thePosition vector.
   m_envelope.thePosition+=G4ThreeVector(m_offsetX,m_offsetY,m_offsetZ);
   if (m_envelope.thePositionedVolume)
     {
       // Override the translation for m_envelope.thePositionedVolume.
       m_envelope.thePositionedVolume->SetTranslation(m_envelope.thePosition);
     }
  
   ATH_MSG_VERBOSE( name() << "::SetRotationAndOffset() (Base class method): Finished" );
   return;
 }

Member Data Documentation

◆ m_backWall

bool MuonWallTool::m_backWall

private

Definition at line 32 of file MuonWallTool.h.

◆ m_detectorName

Gaudi::Property<std::string> DetectorGeometryBase::m_detectorName {this, "DetectorName", "", "Detector name (same as the Tool name if not set"}

protectedinherited

Definition at line 71 of file DetectorGeometryBase.h.

◆ m_envelope

Envelope DetectorGeometryBase::m_envelope

protectedinherited

Definition at line 70 of file DetectorGeometryBase.h.

◆ m_isWorld

bool DetectorGeometryBase::m_isWorld {false}

protectedinherited

Definition at line 72 of file DetectorGeometryBase.h.

◆ m_notifierSvc

ServiceHandle<IG4GeometryNotifierSvc> DetectorGeometryBase::m_notifierSvc {this, "GeometryNotifierSvc", "G4GeometryNotifierSvc", "Detector name (same as the Tool name if not set"}

protectedinherited

Definition at line 66 of file DetectorGeometryBase.h.

◆ m_offsetX

Gaudi::Property<double> DetectorGeometryBase::m_offsetX {this, "OffsetX" , 0.0, "Offset in the X-direction"}

protectedinherited

Definition at line 76 of file DetectorGeometryBase.h.

◆ m_offsetY

Gaudi::Property<double> DetectorGeometryBase::m_offsetY {this, "OffsetY" , 0.0, "Offset in the Y-direction"}

protectedinherited

Definition at line 77 of file DetectorGeometryBase.h.

◆ m_offsetZ

Gaudi::Property<double> DetectorGeometryBase::m_offsetZ {this, "OffsetZ" , 0.0, "Offset in the Z-direction"}

protectedinherited

Definition at line 78 of file DetectorGeometryBase.h.

◆ m_rotateX

Gaudi::Property<double> DetectorGeometryBase::m_rotateX {this, "RotateX" , 0.0, "Rotation around the X-axis"}

protectedinherited

Definition at line 73 of file DetectorGeometryBase.h.

◆ m_rotateY

Gaudi::Property<double> DetectorGeometryBase::m_rotateY {this, "RotateY" , 0.0, "Rotation around the Y-axis"}

protectedinherited

Definition at line 74 of file DetectorGeometryBase.h.

◆ m_rotateZ

Gaudi::Property<double> DetectorGeometryBase::m_rotateZ {this, "RotateZ" , 0.0, "Rotation around the Z-axis"}

protectedinherited

Definition at line 75 of file DetectorGeometryBase.h.

◆ m_sideWall

bool MuonWallTool::m_sideWall

private

Definition at line 33 of file MuonWallTool.h.

◆ m_subDetTools

ToolHandleArray<IDetectorGeometryTool> DetectorGeometryBase::m_subDetTools {this, "SubDetectors", {}, "Tool handle array of all subdetector tools"}

protectedinherited

Definition at line 64 of file DetectorGeometryBase.h.

◆ m_theParent

IDetectorGeometryTool* DetectorGeometryBase::m_theParent {}

protectedinherited

Definition at line 68 of file DetectorGeometryBase.h.

◆ m_xLength

double MuonWallTool::m_xLength

private

Definition at line 31 of file MuonWallTool.h.

◆ m_yLength

double MuonWallTool::m_yLength

private

Definition at line 30 of file MuonWallTool.h.

◆ m_zLength

double MuonWallTool::m_zLength

private

Definition at line 29 of file MuonWallTool.h.

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

Public Member Functions

Protected Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ MuonWallTool()

◆ ~MuonWallTool()

Member Function Documentation

◆ Build()

◆ BuildGeometry()

◆ BuildSubDetectors()

◆ GetDetectorName()

◆ GetEnvelope()

◆ GetWorldVolume()

◆ initialize()

◆ PositionInParent()

◆ ResetEnvelope()

◆ SetAsWorld()

◆ SetDetectorName()

◆ SetEnvelope()

◆ SetParent()

◆ SetRotationAndOffset()

Member Data Documentation

◆ m_backWall

◆ m_detectorName

◆ m_envelope

◆ m_isWorld

◆ m_notifierSvc

◆ m_offsetX

◆ m_offsetY

◆ m_offsetZ

◆ m_rotateX

◆ m_rotateY

◆ m_rotateZ

◆ m_sideWall

◆ m_subDetTools

◆ m_theParent

◆ m_xLength

◆ m_yLength

◆ m_zLength