Tool for building the MuonWall detector. More...

#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 ()=default
virtual void	BuildGeometry () override final
	Override DetectorGeometryBase::BuildGeometry method.
StatusCode	initialize () override
	Athena method.
virtual void	Build () override
	purely virtual methods being implemented here
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
Gaudi::Property< double >	m_zLength {this, "ZLength", 0.}
Gaudi::Property< double >	m_yLength {this, "YLength", 0.}
Gaudi::Property< double >	m_xLength {this, "XLength", 0.}
Gaudi::Property< bool >	m_backWall {this, "backWall", true}
Gaudi::Property< bool >	m_sideWall {this, "sideWall", false}

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

◆ ~MuonWallTool()

MuonWallTool::~MuonWallTool ( )

default

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 39 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

Gaudi::Property<bool> MuonWallTool::m_backWall {this, "backWall", true}

private

Definition at line 32 of file MuonWallTool.h.

32{this, "backWall", true};

◆ 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.

71{this, "DetectorName", "", "Detector name (same as the Tool name if not set"};

◆ 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.

72{false};

◆ 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.

66{this, "GeometryNotifierSvc", "G4GeometryNotifierSvc", "Detector name (same as the Tool name if not set"};

◆ 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.

76{this, "OffsetX" , 0.0, "Offset in the X-direction"};

◆ 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.

77{this, "OffsetY" , 0.0, "Offset in the Y-direction"};

◆ 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.

78{this, "OffsetZ" , 0.0, "Offset in the Z-direction"};

◆ 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.

73{this, "RotateX" , 0.0, "Rotation around the X-axis"};

◆ 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.

74{this, "RotateY" , 0.0, "Rotation around the Y-axis"};

◆ 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.

75{this, "RotateZ" , 0.0, "Rotation around the Z-axis"};

◆ m_sideWall

Gaudi::Property<bool> MuonWallTool::m_sideWall {this, "sideWall", false}

private

Definition at line 33 of file MuonWallTool.h.

33{this, "sideWall", false};

◆ 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.

64{this, "SubDetectors", {}, "Tool handle array of all subdetector tools"};

◆ m_theParent

IDetectorGeometryTool* DetectorGeometryBase::m_theParent {}

protectedinherited

Definition at line 68 of file DetectorGeometryBase.h.

68{};

◆ m_xLength

Gaudi::Property<double> MuonWallTool::m_xLength {this, "XLength", 0.}

private

Definition at line 31 of file MuonWallTool.h.

31{this, "XLength", 0.};

◆ m_yLength

Gaudi::Property<double> MuonWallTool::m_yLength {this, "YLength", 0.}

private

Definition at line 30 of file MuonWallTool.h.

30{this, "YLength", 0.};

◆ m_zLength

Gaudi::Property<double> MuonWallTool::m_zLength {this, "ZLength", 0.}

private

Definition at line 29 of file MuonWallTool.h.

29{this, "ZLength", 0.};

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