Detector construction. More...

#include <Doxymodules_parameterisations.h>

Inheritance diagram for Par04DetectorConstruction:

Public Member Functions
	Par04DetectorConstruction ()

virtual	~Par04DetectorConstruction ()

virtual G4VPhysicalVolume *	Construct () final

virtual void	ConstructSDandField () final

void	SetInnerRadius (G4double aInnerRadius)
	Set inner radius of the cylindrical detector.

G4double	GetInnerRadius () const
	Get inner radius of the cylindrical detector.

void	SetLength (G4double aLength)
	Set length radius of the cylindrical detector.

G4double	GetLength () const
	Get length of the cylindrical detector (along z-axis)

void	SetNbOfLayers (G4int aNumber)
	Set number of layers.

G4int	GetNbOfLayers () const
	Get number of layers.

void	SetAbsorberMaterial (const std::size_t aLayer, const G4String &aMaterial)
	Set material of the layer (from NIST materials)

G4String	GetAbsorberMaterial (const std::size_t aLayer) const
	Get name of the material of the layer.

void	SetAbsorberThickness (const std::size_t aLayer, const G4double aThickness)
	Set thickness of the layer.

G4double	GetAbsorberThickness (const std::size_t aLayer) const
	Get thickness of the layer.

void	SetAbsorberSensitivity (const std::size_t aLayer, const G4bool aSensitivity)
	Set sensitivity of the layer.

G4bool	GetAbsorberSensitivity (const std::size_t aLayer) const
	Get sensitivity of the layer.

void	SetMeshNbOfCells (G4ThreeVector aNb)
	Set number of Mesh cells in cylindrical coordinates (r, phi, z)

void	SetMeshNbOfCells (std::size_t aIndex, G4double aNb)
	Set number of Mesh cells in cylindrical coordinates along one of the axis.

G4ThreeVector	GetMeshNbOfCells () const
	Get number of Mesh cells in cylindrical coordinates (r, phi, z)

void	SetMeshSizeOfCells (G4ThreeVector aNb)
	Set size of Mesh cells in cylindrical coordinates (r, phi, z)

void	SetMeshSizeOfCells (std::size_t aIndex, G4double aNb)
	Set size of Mesh cells in cylindrical coordinates along one of the axis.

G4ThreeVector	GetMeshSizeOfCells () const
	Get size of Mesh cells in cylindrical coordinates (r, phi, z)

void	Print () const
	Print detector information.

Private Attributes
Par04DetectorMessenger *	fDetectorMessenger = nullptr
	Messenger that allows to modify geometry.

G4double	fDetectorInnerRadius = 80 * cm
	Inner radius of the cylindrical detector.

G4double	fDetectorLength = 24 * cm
	Length of the cylindrical detector (along z axis)

std::vector< G4LogicalVolume * >	fLayerLogical
	Logical volume(s) of the sensitive absorbers.

std::array< G4Material *, 2 >	fAbsorberMaterial = { nullptr, nullptr }
	Material(s) of the layers.

std::array< G4double, 2 >	fAbsorberThickness = { 1 * cm, 0 }
	Thickness(es) of the layers.

std::array< G4bool, 2 >	fAbsorberSensitivity = { true, 0 }
	Sensitivity of the layers.

G4int	fNbOfLayers = 24
	Number of layers = slices along z axis.

G4ThreeVector	fMeshNbOfCells = { 40, 50, 48 }
	Mesh number of cells (Nr, Nphi, Nz)

G4ThreeVector	fMeshSizeOfCells = { 5 * mm, 2 * CLHEP::pi / 50 * CLHEP::rad, 5 * mm }
	Mesh size of cells (dr, dphi, dz).

Detailed Description

Detector construction.

Creates a cylindrical detector, with cylinder axis along Z-axis. It is placed in the centre of the world volume. Dimensions of the detector (inner radius and its length) as well as composition (number of radial layers, number of absorbers, its thicknesses and materials) can be set using the UI commands. There may be up to two absorbers used to build the layers.

TODO Extend to allow use of more than two different absorbers.

Each absorber may have differnt thickness (along radial direction), material, and may be either sensitive to particle passage (will register deposited energy) or not (passive). Readout geometry of the detector is dynamic, and its size can be set by UI commands. Cells are created along z-axis, azimuthal angle, and radius (cylindrical segmentation). The z axis is parallel to the direction of the particle entering the detector volume. The mesh also starts at the entrance to the detector volume. In order to define this enrance position and direction, a fast simulation model is used. For the detector volume, if particle has entered, sets the particle direction and position in the event information. Those vectors define the readout mesh for the event.

TODO In order to speed up the simulation, fast simulation that checks for the entrance properites should be defined in a very thin region instead of a whole region of the detector.

Sensitive detector Par04SensitiveDetector is attached to the detector volume Region for the detector is created as an envelope of the fast simulation.

Definition at line 79 of file Doxymodules_parameterisations.h.

Constructor & Destructor Documentation

◆ Par04DetectorConstruction()

Par04DetectorConstruction::Par04DetectorConstruction ( )

Definition at line 60 of file Par04DetectorConstruction.cc.

  : G4VUserDetectorConstruction()
{
  fDetectorMessenger         = new Par04DetectorMessenger(this);
  G4NistManager* nistManager = G4NistManager::Instance();
  fAbsorberMaterial[0]       = nistManager->FindOrBuildMaterial("G4_PbWO4");
}

◆ ~Par04DetectorConstruction()

Par04DetectorConstruction::~Par04DetectorConstruction ( )

virtualdefault

Member Function Documentation

◆ Construct()

G4VPhysicalVolume * Par04DetectorConstruction::Construct ( )

finalvirtual

useful variable

Definition at line 74 of file Par04DetectorConstruction.cc.

{
  //--------- Material definition ---------
  G4NistManager* nistManager = G4NistManager::Instance();
  G4Material* air            = nistManager->FindOrBuildMaterial("G4_AIR");
 
  //--------- Derived dimensions ---------
  G4double full2Pi = 2. * CLHEP::pi * rad;
  G4double layerThickness =
    std::accumulate(fAbsorberThickness.begin(), fAbsorberThickness.end(), 0.);
  G4double detectorOuterRadius = fDetectorInnerRadius + fNbOfLayers * layerThickness;
  G4double worldSizeXY         = detectorOuterRadius * 4.;
  G4double worldSizeZ          = fDetectorLength * 2;
  // check number of materials: (1 = homogeneous calo, 2 = sampling calo)
  G4int nbOfMaterials = 0;
  for(const auto material : fAbsorberMaterial)
  {
    if(material != nullptr)
      nbOfMaterials++;
  }
 
  //--------- World ---------
  auto fSolidWorld  = new G4Box("World",                  // name
                                worldSizeXY / 2.,         // half-width in X
                                worldSizeXY / 2.,         // half-width in Y
                                worldSizeZ / 2.);         // half-width in Z
  auto fLogicWorld  = new G4LogicalVolume(fSolidWorld,    // solid
                                          air,            // material
                                          "World");       // name
  auto fPhysicWorld = new G4PVPlacement(0,                // no rotation
                                        G4ThreeVector(),  // at (0,0,0)
                                        fLogicWorld,      // logical volume
                                        "World",          // name
                                        0,                // mother volume
                                        false,            // not used
                                        99999,            // copy number
                                        false);           // check overlaps
  fLogicWorld->SetVisAttributes(G4VisAttributes::GetInvisible());
 
  //--------- Detector envelope ---------
  auto fSolidDetector = new G4Tubs("Detector",               // name
                                   fDetectorInnerRadius,     // inner radius
                                   detectorOuterRadius,      // outer radius
                                   fDetectorLength / 2.,     // half-width in Z
                                   0,                        // start angle
                                   full2Pi);                 // delta angle
  auto fLogicDetector = new G4LogicalVolume(fSolidDetector,  // solid
                                            air,             // material
                                            "Detector");     // name
  new G4PVPlacement(0,                                       // no rotation
                    G4ThreeVector(0, 0, 0),                  // detector centre at (0,0,0)
                    fLogicDetector,                          // logical volume
                    "Detector",                              // name
                    fLogicWorld,                             // mother volume
                    false,                                   // not used
                    999,                                    // copy number
                    false);                                  // check overlaps
 
  // Region for fast simulation
  auto detectorRegion = new G4Region("DetectorRegion");
  detectorRegion->AddRootLogicalVolume(fLogicDetector);
 
  //--------- Detector layers: material ---------
  std::array<G4VisAttributes, 2> attribs;
  attribs[0].SetColour(G4Colour(0, 0, 1, 0.1));
  attribs[0].SetForceSolid(true);
  attribs[1].SetColour(G4Colour(1, 0, 0, 0.1));
  attribs[1].SetForceSolid(true);
  /// useful variable
  G4double innerRadius = fDetectorInnerRadius;
  for(G4int iLayer = 0; iLayer < fNbOfLayers; iLayer++)
  {
    for(G4int iMaterial = 0; iMaterial < nbOfMaterials; iMaterial++)
    {
      auto fSolidLayer = new G4Tubs("Layer",      // name
                                    innerRadius,  // inner radius
                                    innerRadius + fAbsorberThickness[iMaterial],  // outer radius
                                    fDetectorLength / 2.,  // half-width in Z
                                    0,                     // start angle
                                    full2Pi);              // delta angle
      G4LogicalVolume* logical = new G4LogicalVolume(fSolidLayer,  // solid
                                                     fAbsorberMaterial[iMaterial],  // material
                                                     "Layer");                      // name
      new G4PVPlacement(0,                                                          // no rotation
                        G4ThreeVector(),                     // place at centre of mother volume
                        logical,                             // logical volume
                        "Layer",                             // name
                        fLogicDetector,                      // mother volume
                        false,                               // not used
                        iLayer * nbOfMaterials + iMaterial,  // copy number
                        false);                              // check overlaps
      logical->SetVisAttributes(attribs[iMaterial]);
      innerRadius += fAbsorberThickness[iMaterial];
      if(fAbsorberSensitivity[iMaterial])
      {
        fLayerLogical.push_back(logical);
      }
    }
  }
 
  Print();
  return fPhysicWorld;
}

◆ ConstructSDandField()

void Par04DetectorConstruction::ConstructSDandField ( )

finalvirtual

Definition at line 180 of file Par04DetectorConstruction.cc.

{
  Par04SensitiveDetector* caloSD =
    new Par04SensitiveDetector("sensitiveDetector", fMeshNbOfCells, fMeshSizeOfCells);
  G4SDManager::GetSDMpointer()->AddNewDetector(caloSD);
  for(const auto logical : fLayerLogical)
  {
    SetSensitiveDetector(logical, caloSD);
  }
 
  auto detectorRegion = G4RegionStore::GetInstance()->GetRegion("DetectorRegion");
  // Par04DefineMeshModel needs to be first model to call
  new Par04DefineMeshModel("defineMesh", detectorRegion);
#ifdef USE_INFERENCE
  new Par04MLFastSimModel("inferenceModel", detectorRegion);
#endif
}

◆ SetInnerRadius()

void Par04DetectorConstruction::SetInnerRadius ( G4double aInnerRadius )

Set inner radius of the cylindrical detector.

Definition at line 264 of file Par04DetectorConstruction.cc.

264{ fDetectorInnerRadius = aRadius; }

◆ GetInnerRadius()

G4double Par04DetectorConstruction::GetInnerRadius ( ) const

inline

Get inner radius of the cylindrical detector.

Definition at line 84 of file Par04DetectorConstruction.hh.

84{ return fDetectorInnerRadius; };

◆ SetLength()

void Par04DetectorConstruction::SetLength ( G4double aLength )

Set length radius of the cylindrical detector.

Definition at line 268 of file Par04DetectorConstruction.cc.

268{ fDetectorLength = aLength; }

◆ GetLength()

G4double Par04DetectorConstruction::GetLength ( ) const

inline

Get length of the cylindrical detector (along z-axis)

Definition at line 88 of file Par04DetectorConstruction.hh.

88{ return fDetectorLength; };

◆ SetNbOfLayers()

void Par04DetectorConstruction::SetNbOfLayers ( G4int aNumber )

inline

Set number of layers.

Definition at line 90 of file Par04DetectorConstruction.hh.

90{ fNbOfLayers = aNumber; };

◆ GetNbOfLayers()

G4int Par04DetectorConstruction::GetNbOfLayers ( ) const

inline

Get number of layers.

Definition at line 92 of file Par04DetectorConstruction.hh.

92{ return fNbOfLayers; };

◆ SetAbsorberMaterial()

void Par04DetectorConstruction::SetAbsorberMaterial	(	const std::size_t	aLayer,
		const G4String &	aMaterial
	)

Set material of the layer (from NIST materials)

Definition at line 218 of file Par04DetectorConstruction.cc.

{
  // search material by its name
  G4Material* material = G4NistManager::Instance()->FindOrBuildMaterial(aName);
  if(material)
    fAbsorberMaterial[aLayer] = material;
  else
    G4Exception("Par04DetectorConstruction::SetAbsorberMaterial()", "InvalidSetup", FatalException,
                ("Unknown material name: " + aName).c_str());
  G4RunManager::GetRunManager()->PhysicsHasBeenModified();
}

◆ GetAbsorberMaterial()

G4String Par04DetectorConstruction::GetAbsorberMaterial ( const std::size_t aLayer ) const

inline

Get name of the material of the layer.

Definition at line 97 of file Par04DetectorConstruction.hh.

  {
    return fAbsorberMaterial[aLayer]->GetName();
  };

◆ SetAbsorberThickness()

void Par04DetectorConstruction::SetAbsorberThickness	(	const std::size_t	aLayer,
		const G4double	aThickness
	)

Set thickness of the layer.

Definition at line 232 of file Par04DetectorConstruction.cc.

{
  if(aLayer < fAbsorberThickness.size())
    fAbsorberThickness[aLayer] = aThickness;
  else
    G4Exception("Par04DetectorConstruction::SetAbsorberThickness()", "InvalidSetup", FatalException,
                ("Requested layer " + std::to_string(aLayer) +
                 " is larger than number of available layers (" +
                 std::to_string(fAbsorberThickness.size()) + ").")
                  .c_str());
  G4RunManager::GetRunManager()->PhysicsHasBeenModified();
}

◆ GetAbsorberThickness()

G4double Par04DetectorConstruction::GetAbsorberThickness ( const std::size_t aLayer ) const

inline

Get thickness of the layer.

Definition at line 104 of file Par04DetectorConstruction.hh.

  {
    return fAbsorberThickness[aLayer];
  };

◆ SetAbsorberSensitivity()

void Par04DetectorConstruction::SetAbsorberSensitivity	(	const std::size_t	aLayer,
		const G4bool	aSensitivity
	)

Set sensitivity of the layer.

Definition at line 248 of file Par04DetectorConstruction.cc.

{
  if(aLayer < fAbsorberSensitivity.size())
    fAbsorberSensitivity[aLayer] = aSensitivity;
  else
    G4Exception(
      "Par04DetectorConstruction::SetAbsorberSensitivity()", "InvalidSetup", FatalException,
      ("Requested layer " + std::to_string(aLayer) + " is larger than number of available layers ("+
       std::to_string(fAbsorberSensitivity.size()) + ").")
        .c_str());
  G4RunManager::GetRunManager()->PhysicsHasBeenModified();
}

◆ GetAbsorberSensitivity()

G4bool Par04DetectorConstruction::GetAbsorberSensitivity ( const std::size_t aLayer ) const

inline

Get sensitivity of the layer.

Definition at line 111 of file Par04DetectorConstruction.hh.

  {
    return fAbsorberSensitivity[aLayer];
  };

◆ SetMeshNbOfCells() [1/2]

void Par04DetectorConstruction::SetMeshNbOfCells ( G4ThreeVector aNb )

inline

Set number of Mesh cells in cylindrical coordinates (r, phi, z)

Definition at line 117 of file Par04DetectorConstruction.hh.

117{ fMeshNbOfCells = aNb; };

◆ SetMeshNbOfCells() [2/2]

void Par04DetectorConstruction::SetMeshNbOfCells	(	std::size_t	aIndex,
		G4double	aNb
	)

inline

Set number of Mesh cells in cylindrical coordinates along one of the axis.

Parameters

[in] aIndex index of cylindrical axis (0,1,2) = (r, phi, z)

Definition at line 120 of file Par04DetectorConstruction.hh.

120{ fMeshNbOfCells[aIndex] = aNb; };

◆ GetMeshNbOfCells()

G4ThreeVector Par04DetectorConstruction::GetMeshNbOfCells ( ) const

inline

Get number of Mesh cells in cylindrical coordinates (r, phi, z)

Definition at line 122 of file Par04DetectorConstruction.hh.

122{ return fMeshNbOfCells; };

◆ SetMeshSizeOfCells() [1/2]

void Par04DetectorConstruction::SetMeshSizeOfCells ( G4ThreeVector aNb )

inline

Set size of Mesh cells in cylindrical coordinates (r, phi, z)

Definition at line 124 of file Par04DetectorConstruction.hh.

124{ fMeshSizeOfCells = aNb; };

◆ SetMeshSizeOfCells() [2/2]

void Par04DetectorConstruction::SetMeshSizeOfCells	(	std::size_t	aIndex,
		G4double	aNb
	)

inline

Set size of Mesh cells in cylindrical coordinates along one of the axis.

Parameters

[in] aIndex index of cylindrical axis (0,1,2) = (r, phi, z)

Definition at line 127 of file Par04DetectorConstruction.hh.

128 { fMeshSizeOfCells[aIndex] = aNb; };

◆ GetMeshSizeOfCells()

G4ThreeVector Par04DetectorConstruction::GetMeshSizeOfCells ( ) const

inline

Get size of Mesh cells in cylindrical coordinates (r, phi, z)

Definition at line 130 of file Par04DetectorConstruction.hh.

130{ return fMeshSizeOfCells; };

◆ Print()

void Par04DetectorConstruction::Print ( ) const

Print detector information.

Definition at line 200 of file Par04DetectorConstruction.cc.

{
  G4cout << "\n------------------------------------------------------"
         << "\n--- Detector length:\t" << G4BestUnit(fDetectorLength, "Length")
         << "\n--- Detector inner radius:\t" << G4BestUnit(fDetectorInnerRadius, "Length")
         << "\n--- Number of layers:\t" << fNbOfLayers << G4endl << "\n--- 1st layer: \t"
         << G4BestUnit(fAbsorberThickness[0], "Length") << " of "
         << (fAbsorberSensitivity[0] ? "active " : "passive ") << fAbsorberMaterial[0]->GetName()
         << G4endl;
  if(fAbsorberMaterial[1] != nullptr)
    G4cout << "--- 2nd layer: \t" << G4BestUnit(fAbsorberThickness[1], "Length") << " of "
           << (fAbsorberSensitivity[1] ? "active " : "passive ") << fAbsorberMaterial[1]->GetName()
           << G4endl;
  G4cout << "-----------------------------------------------------" << G4endl;
}

Member Data Documentation

◆ fDetectorMessenger

Par04DetectorMessenger* Par04DetectorConstruction::fDetectorMessenger = nullptr

private

Messenger that allows to modify geometry.

Definition at line 137 of file Par04DetectorConstruction.hh.

◆ fDetectorInnerRadius

G4double Par04DetectorConstruction::fDetectorInnerRadius = 80 * cm

private

Inner radius of the cylindrical detector.

Definition at line 139 of file Par04DetectorConstruction.hh.

◆ fDetectorLength

G4double Par04DetectorConstruction::fDetectorLength = 24 * cm

private

Length of the cylindrical detector (along z axis)

Definition at line 141 of file Par04DetectorConstruction.hh.

◆ fLayerLogical

std::vector<G4LogicalVolume*> Par04DetectorConstruction::fLayerLogical

private

Logical volume(s) of the sensitive absorbers.

Definition at line 143 of file Par04DetectorConstruction.hh.

◆ fAbsorberMaterial

std::array<G4Material*, 2> Par04DetectorConstruction::fAbsorberMaterial = { nullptr, nullptr }

private

Material(s) of the layers.

Definition at line 145 of file Par04DetectorConstruction.hh.

145{ nullptr, nullptr };

◆ fAbsorberThickness

std::array<G4double, 2> Par04DetectorConstruction::fAbsorberThickness = { 1 * cm, 0 }

private

Thickness(es) of the layers.

Definition at line 147 of file Par04DetectorConstruction.hh.

147{ 1 * cm, 0 };

◆ fAbsorberSensitivity

std::array<G4bool, 2> Par04DetectorConstruction::fAbsorberSensitivity = { true, 0 }

private

Sensitivity of the layers.

Definition at line 149 of file Par04DetectorConstruction.hh.

149{ true, 0 };

◆ fNbOfLayers

G4int Par04DetectorConstruction::fNbOfLayers = 24

private

Number of layers = slices along z axis.

Definition at line 151 of file Par04DetectorConstruction.hh.

◆ fMeshNbOfCells

G4ThreeVector Par04DetectorConstruction::fMeshNbOfCells = { 40, 50, 48 }

private

Mesh number of cells (Nr, Nphi, Nz)

Definition at line 153 of file Par04DetectorConstruction.hh.

153{ 40, 50, 48 };

◆ fMeshSizeOfCells

G4ThreeVector Par04DetectorConstruction::fMeshSizeOfCells = { 5 * mm, 2 * CLHEP::pi / 50 * CLHEP::rad, 5 * mm }

private

Mesh size of cells (dr, dphi, dz).

Definition at line 155 of file Par04DetectorConstruction.hh.

155{ 5 * mm, 2 * CLHEP::pi / 50 * CLHEP::rad, 5 * mm };

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

.doxygen/Doxymodules_parameterisations.h
extended/parameterisations/Par04/include/Par04DetectorConstruction.hh
extended/parameterisations/Par04/src/Par04DetectorConstruction.cc

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Par04DetectorConstruction()

◆ ~Par04DetectorConstruction()

Member Function Documentation

◆ Construct()

◆ ConstructSDandField()

◆ SetInnerRadius()

◆ GetInnerRadius()

◆ SetLength()

◆ GetLength()

◆ SetNbOfLayers()

◆ GetNbOfLayers()

◆ SetAbsorberMaterial()

◆ GetAbsorberMaterial()

◆ SetAbsorberThickness()

◆ GetAbsorberThickness()

◆ SetAbsorberSensitivity()

◆ GetAbsorberSensitivity()

◆ SetMeshNbOfCells() [1/2]

◆ SetMeshNbOfCells() [2/2]

◆ GetMeshNbOfCells()

◆ SetMeshSizeOfCells() [1/2]

◆ SetMeshSizeOfCells() [2/2]

◆ GetMeshSizeOfCells()

◆ Print()

Member Data Documentation

◆ fDetectorMessenger

◆ fDetectorInnerRadius

◆ fDetectorLength

◆ fLayerLogical

◆ fAbsorberMaterial

◆ fAbsorberThickness

◆ fAbsorberSensitivity

◆ fNbOfLayers

◆ fMeshNbOfCells

◆ fMeshSizeOfCells