Detector construction. More...

#include <Doxymodules_parameterisations.h>

Inheritance diagram for Par03DetectorConstruction:

Public Member Functions
	Par03DetectorConstruction ()

virtual	~Par03DetectorConstruction ()

virtual G4VPhysicalVolume *	Construct () final

virtual void	ConstructSDandField () final

void	SetRadius (G4double aRadius)

G4double	GetRadius () const

void	SetLength (G4double aLength)

G4double	GetLength () const

void	SetMaterial (const G4String &aMaterial)

G4String	GetMaterial () const

void	SetNbOfLayers (G4int aNumber)

G4int	GetNbOfLayers () const

void	SetNbOfRhoCells (G4int aNumber)

G4int	GetNbOfRhoCells () const

void	SetNbOfPhiCells (G4int aNumber)

G4int	GetNbOfPhiCells () const

void	Print () const

Private Attributes
Par03DetectorMessenger *	fDetectorMessenger
	Messenger that allows to modify geometry.

G4LogicalVolume *	fLogicCell = nullptr
	Logical volume of replicated cell.

G4double	fWorldSize = 10 * m
	World size (in each X, Y, Z dimension)

G4double	fDetectorRadius = 10 * cm
	Radius of the cylindrical detector.

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

G4Material *	fDetectorMaterial = nullptr
	Material of the detector.

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

G4int	fNbOfRhoCells = 10
	Number of cells along radius.

G4int	fNbOfPhiCells = 10
	Number of cells in azimuthal angle.

Detailed Description

Detector construction.

Creates a cylindrical detector, with cylinder axis along Z-axis. It is placed in the world volume so that its bases are located at z=0 and z=Length. Dimensions of the detector (Radius and Length) and material can be set using the UI commands. Readout geometry of the detector is created, and can be set by UI commands. Cells are created along z-axis, azimuthal angle, and radius (cylindrical segmentation). Sensitive detector Par03SensitiveDetector is attached to the cell volume. Region for the detector is created as an envelope of the fast simulation.

Definition at line 60 of file Doxymodules_parameterisations.h.

Constructor & Destructor Documentation

◆ Par03DetectorConstruction()

Par03DetectorConstruction::Par03DetectorConstruction ( )

Definition at line 50 of file Par03DetectorConstruction.cc.

  : G4VUserDetectorConstruction()
{
  fDetectorMessenger = new Par03DetectorMessenger(this);
 
  G4NistManager* nistManager = G4NistManager::Instance();
  fDetectorMaterial          = nistManager->FindOrBuildMaterial("G4_Fe");
}

◆ ~Par03DetectorConstruction()

Par03DetectorConstruction::~Par03DetectorConstruction ( )

virtualdefault

Member Function Documentation

◆ Construct()

G4VPhysicalVolume * Par03DetectorConstruction::Construct ( )

finalvirtual

Definition at line 65 of file Par03DetectorConstruction.cc.

{
  //--------- Material definition ---------
  G4NistManager* nistManager = G4NistManager::Instance();
  G4Material* air            = nistManager->FindOrBuildMaterial("G4_AIR");
 
  //--------- Derived dimensions ---------
  G4double full2Pi        = 2. * CLHEP::pi * rad;
  G4double layerThickness = fDetectorLength / fNbOfLayers;
  G4double cellPhi        = full2Pi / fNbOfPhiCells;
  G4double cellDR         = fDetectorRadius / fNbOfRhoCells;
 
  //--------- World ---------
  auto fSolidWorld  = new G4Box("World",                  // name
                               fWorldSize / 2.,          // half-width in X
                               fWorldSize / 2.,          // half-width in Y
                               fWorldSize / 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
                                        999,              // copy number
                                        true);            // copy number
 
  //--------- Detector envelope ---------
  auto fSolidDetector = new G4Tubs("Detector",               // name
                                   0,                        // inner radius
                                   fDetectorRadius,          // outer radius
                                   fDetectorLength / 2.,     // half-width in Z
                                   0,                        // start angle
                                   full2Pi);                 // delta angle
  auto fLogicDetector = new G4LogicalVolume(fSolidDetector,  // solid
                                            fDetectorMaterial,  // material
                                            "Detector");        // name
  new G4PVPlacement(
    0,  // no rotation
    G4ThreeVector(0, 0,
                  fDetectorLength / 2),  // detector face starts at (0,0,0)
    fLogicDetector,                      // logical volume
    "Detector",                          // name
    fLogicWorld,                         // mother volume
    false,                               // not used
    99,                                  // copy number
    true);                               // check overlaps
 
  // Region for fast simulation
  auto detectorRegion = new G4Region("DetectorRegion");
  detectorRegion->AddRootLogicalVolume(fLogicDetector);
 
  //--------- Readout geometry ---------
  // Layers (along z)
  auto fSolidLayer = new G4Tubs("Layer",               // name
                                0,                     // inner radius
                                fDetectorRadius,       // outer radius
                                layerThickness / 2.,   // half-width in Z
                                0,                     // start angle
                                full2Pi);              // delta angle
  auto fLogicLayer = new G4LogicalVolume(fSolidLayer,  // solid
                                         air,          // material
                                         "Layer");     // name
  if(fNbOfLayers > 1)
    new G4PVReplica("Layer",          // name
                    fLogicLayer,      // logical volume
                    fLogicDetector,   // mother volume
                    kZAxis,           // axis of replication
                    fNbOfLayers,      // number of replicas
                    layerThickness);  // width of single replica
  else
    new G4PVPlacement(0,                // no rotation
                      G4ThreeVector(),  // place at centre of mother volume
                      fLogicLayer,      // logical volume
                      "Layer",          // name
                      fLogicDetector,   // mother volume
                      false,            // not used
                      0,                // copy number
                      true);            // check overlaps
 
  // Layer segment (division in phi)
  auto fSolidRow = new G4Tubs("Row",                // name
                              0,                    // inner radius
                              fDetectorRadius,      // outer radius
                              layerThickness / 2.,  // half-width in Z
                              0,                    // start angle
                              cellPhi);             // delta angle
 
  auto fLogicRow = new G4LogicalVolume(fSolidRow,   // solid
                                       air,         // material
                                       "Segment");  // name
  if(fNbOfPhiCells > 1)
    new G4PVReplica("Segment",      // name
                    fLogicRow,      // logical volume
                    fLogicLayer,    // mother volume
                    kPhi,           // axis of replication
                    fNbOfPhiCells,  // number of replicas
                    cellPhi);       // width of single replica
  else
    new G4PVPlacement(0,                // no rotation
                      G4ThreeVector(),  // place at centre of mother volume
                      fLogicRow,        // logical volume
                      "Row",            // name
                      fLogicLayer,      // mother volume
                      false,            // not used
                      0,                // copy number
                      true);            // check overlaps
 
  // Final cells (segment slices in radius)
  // No volume can be placed inside a radial replication
  auto fSolidCell = new G4Tubs("Cell",               // name
                               0,                    // inner radius
                               cellDR,               // outer radius
                               layerThickness / 2.,  // half-width in Z
                               0,                    // start angle
                               cellPhi);             // delta angle
 
  fLogicCell = new G4LogicalVolume(fSolidCell,         // solid
                                   fDetectorMaterial,  // material
                                   "Cell");            // name
  if(fNbOfRhoCells > 1)
    new G4PVReplica("Cell",         // name
                    fLogicCell,     // logical volume
                    fLogicRow,      // mother volume
                    kRho,           // axis of replication
                    fNbOfRhoCells,  // number of replicas
                    cellDR);        // width of single replica
  else
    new G4PVPlacement(0,                // no rotation
                      G4ThreeVector(),  // place at centre of mother volume
                      fLogicCell,       // logical volume
                      "Cell",           // name
                      fLogicRow,        // mother volume
                      false,            // not used
                      0,                // copy number
                      true);            // check overlaps
 
  //--------- Visualisation settings ---------
  fLogicWorld->SetVisAttributes(G4VisAttributes::GetInvisible());
  fLogicLayer->SetVisAttributes(G4VisAttributes::GetInvisible());
  fLogicRow->SetVisAttributes(G4VisAttributes::GetInvisible());
  G4VisAttributes attribs;
  attribs.SetColour(G4Colour(0, 0, 1, 0.3));
  attribs.SetForceSolid(true);
  fLogicCell->SetVisAttributes(attribs);
 
  Print();
  return fPhysicWorld;
}

◆ ConstructSDandField()

void Par03DetectorConstruction::ConstructSDandField ( )

finalvirtual

Definition at line 219 of file Par03DetectorConstruction.cc.

{
  Par03SensitiveDetector* caloSD = new Par03SensitiveDetector(
    "sensitiveDetector", fNbOfLayers, fNbOfPhiCells, fNbOfRhoCells);
  G4SDManager::GetSDMpointer()->AddNewDetector(caloSD);
  SetSensitiveDetector(fLogicCell, caloSD);
 
  auto detectorRegion =
    G4RegionStore::GetInstance()->GetRegion("DetectorRegion");
  new Par03EMShowerModel("model", detectorRegion);
}

◆ SetRadius()

void Par03DetectorConstruction::SetRadius ( G4double aRadius )

Definition at line 261 of file Par03DetectorConstruction.cc.

{
  // check if fits within world volume
  if(aRadius >= fWorldSize / 2.)
    G4Exception("Par03DetectorConstruction::SetRadius()", "InvalidSetup",
                FatalException,
                ("Detector radius cannot be larger than the world size (" +
                 G4String(G4BestUnit(fWorldSize / 2., "Length")) + ")")
                  .c_str());
  fDetectorRadius = aRadius;
}

◆ GetRadius()

G4double Par03DetectorConstruction::GetRadius ( ) const

inline

Definition at line 64 of file Par03DetectorConstruction.hh.

64{ return fDetectorRadius; };

◆ SetLength()

void Par03DetectorConstruction::SetLength ( G4double aLength )

Definition at line 274 of file Par03DetectorConstruction.cc.

{
  // check if fits within world volume
  if(aLength >= fWorldSize / 2.)
    G4Exception("Par03DetectorConstruction::SetLength()", "InvalidSetup",
                FatalException,
                ("Detector length cannot be larger than the world size (" +
                 G4String(G4BestUnit(fWorldSize / 2., "Length")) + ")")
                  .c_str());
  fDetectorLength = aLength;
}

◆ GetLength()

G4double Par03DetectorConstruction::GetLength ( ) const

inline

Definition at line 68 of file Par03DetectorConstruction.hh.

68{ return fDetectorLength; };

◆ SetMaterial()

void Par03DetectorConstruction::SetMaterial ( const G4String & aMaterial )

Definition at line 247 of file Par03DetectorConstruction.cc.

{
  // search material by its name
  G4Material* material = G4NistManager::Instance()->FindOrBuildMaterial(aName);
  if(material)
    fDetectorMaterial = material;
  else
    G4Exception("Par03DetectorConstruction::SetMaterial()", "InvalidSetup",
                FatalException, ("Unknown material name: " + aName).c_str());
  G4RunManager::GetRunManager()->PhysicsHasBeenModified();
}

◆ GetMaterial()

G4String Par03DetectorConstruction::GetMaterial ( ) const

inline

Definition at line 72 of file Par03DetectorConstruction.hh.

72{ return fDetectorMaterial->GetName(); };

◆ SetNbOfLayers()

void Par03DetectorConstruction::SetNbOfLayers ( G4int aNumber )

inline

Definition at line 75 of file Par03DetectorConstruction.hh.

75{ fNbOfLayers = aNumber; };

◆ GetNbOfLayers()

G4int Par03DetectorConstruction::GetNbOfLayers ( ) const

inline

Definition at line 77 of file Par03DetectorConstruction.hh.

77{ return fNbOfLayers; };

◆ SetNbOfRhoCells()

void Par03DetectorConstruction::SetNbOfRhoCells ( G4int aNumber )

inline

Definition at line 79 of file Par03DetectorConstruction.hh.

79{ fNbOfRhoCells = aNumber; };

◆ GetNbOfRhoCells()

G4int Par03DetectorConstruction::GetNbOfRhoCells ( ) const

inline

Definition at line 81 of file Par03DetectorConstruction.hh.

81{ return fNbOfRhoCells; };

◆ SetNbOfPhiCells()

void Par03DetectorConstruction::SetNbOfPhiCells ( G4int aNumber )

inline

Definition at line 83 of file Par03DetectorConstruction.hh.

83{ fNbOfPhiCells = aNumber; };

◆ GetNbOfPhiCells()

G4int Par03DetectorConstruction::GetNbOfPhiCells ( ) const

inline

Definition at line 85 of file Par03DetectorConstruction.hh.

85{ return fNbOfPhiCells; };

◆ Print()

void Par03DetectorConstruction::Print ( ) const

Definition at line 233 of file Par03DetectorConstruction.cc.

{
  G4cout << "\n------------------------------------------------------"
         << "\n--- Detector material:\t" << fDetectorMaterial->GetName()
         << "\n--- Detector length:\t" << G4BestUnit(fDetectorLength, "Length")
         << "\n--- Detector radius:\t" << G4BestUnit(fDetectorRadius, "Length")
         << "\n--- Number of layers:\t" << fNbOfLayers
         << "\n--- Number of R-cells:\t" << fNbOfRhoCells
         << "\n--- Number of phi-cells:\t" << fNbOfPhiCells << G4endl;
  G4cout << "-----------------------------------------------------" << G4endl;
}

Member Data Documentation

◆ fDetectorMessenger

Par03DetectorMessenger* Par03DetectorConstruction::fDetectorMessenger

private

Messenger that allows to modify geometry.

Definition at line 92 of file Par03DetectorConstruction.hh.

◆ fLogicCell

G4LogicalVolume* Par03DetectorConstruction::fLogicCell = nullptr

private

Logical volume of replicated cell.

Definition at line 94 of file Par03DetectorConstruction.hh.

◆ fWorldSize

G4double Par03DetectorConstruction::fWorldSize = 10 * m

private

World size (in each X, Y, Z dimension)

Definition at line 96 of file Par03DetectorConstruction.hh.

◆ fDetectorRadius

G4double Par03DetectorConstruction::fDetectorRadius = 10 * cm

private

Radius of the cylindrical detector.

Definition at line 98 of file Par03DetectorConstruction.hh.

◆ fDetectorLength

G4double Par03DetectorConstruction::fDetectorLength = 30 * cm

private

Length of the cylindrical detector (along z axis)

Definition at line 100 of file Par03DetectorConstruction.hh.

◆ fDetectorMaterial

G4Material* Par03DetectorConstruction::fDetectorMaterial = nullptr

private

Material of the detector.

Definition at line 102 of file Par03DetectorConstruction.hh.

◆ fNbOfLayers

G4int Par03DetectorConstruction::fNbOfLayers = 10

private

Number of layers = slices along z axis.

Definition at line 104 of file Par03DetectorConstruction.hh.

◆ fNbOfRhoCells

G4int Par03DetectorConstruction::fNbOfRhoCells = 10

private

Number of cells along radius.

Definition at line 106 of file Par03DetectorConstruction.hh.

◆ fNbOfPhiCells

G4int Par03DetectorConstruction::fNbOfPhiCells = 10

private

Number of cells in azimuthal angle.

Definition at line 108 of file Par03DetectorConstruction.hh.

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

.doxygen/Doxymodules_parameterisations.h
extended/parameterisations/Par03/include/Par03DetectorConstruction.hh
extended/parameterisations/Par03/src/Par03DetectorConstruction.cc

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ Par03DetectorConstruction()

◆ ~Par03DetectorConstruction()

Member Function Documentation

◆ Construct()

◆ ConstructSDandField()

◆ SetRadius()

◆ GetRadius()

◆ SetLength()

◆ GetLength()

◆ SetMaterial()

◆ GetMaterial()

◆ SetNbOfLayers()

◆ GetNbOfLayers()

◆ SetNbOfRhoCells()

◆ GetNbOfRhoCells()

◆ SetNbOfPhiCells()

◆ GetNbOfPhiCells()

◆ Print()

Member Data Documentation

◆ fDetectorMessenger

◆ fLogicCell

◆ fWorldSize

◆ fDetectorRadius

◆ fDetectorLength

◆ fDetectorMaterial

◆ fNbOfLayers

◆ fNbOfRhoCells

◆ fNbOfPhiCells