B4DetectorConstruction Class Reference

#include <B4DetectorConstruction.hh>

Inheritance diagram for B4DetectorConstruction:

Public Member Functions
	B4DetectorConstruction ()
virtual	~B4DetectorConstruction ()
virtual G4VPhysicalVolume *	Construct ()
virtual void	ConstructSDandField ()
const G4VPhysicalVolume *	GetAbsorberPV () const
const G4VPhysicalVolume *	GetGapPV () const
	B4DetectorConstruction ()
virtual	~B4DetectorConstruction ()
virtual G4VPhysicalVolume *	Construct ()
virtual void	ConstructSDandField ()
const G4VPhysicalVolume *	GetAbsorberPV () const
const G4VPhysicalVolume *	GetGapPV () const
	B4DetectorConstruction ()
virtual	~B4DetectorConstruction ()
virtual G4VPhysicalVolume *	Construct ()
virtual void	ConstructSDandField ()
const G4VPhysicalVolume *	GetAbsorberPV () const
const G4VPhysicalVolume *	GetGapPV () const

Private Member Functions
void	DefineMaterials ()
G4VPhysicalVolume *	DefineVolumes ()
void	DefineMaterials ()
G4VPhysicalVolume *	DefineVolumes ()
void	DefineMaterials ()
G4VPhysicalVolume *	DefineVolumes ()

Private Attributes
G4VPhysicalVolume *	fAbsorberPV
G4VPhysicalVolume *	fGapPV
G4bool	fCheckOverlaps

Static Private Attributes
static G4ThreadLocal G4GlobalMagFieldMessenger *	fMagFieldMessenger = 0

Detailed Description

Detector construction class to define materials and geometry. The calorimeter is a box made of a given number of layers. A layer consists of an absorber plate and of a detection gap. The layer is replicated.

Four parameters define the geometry of the calorimeter :

• the thickness of an absorber plate,
• the thickness of a gap,
• the number of layers,
• the transverse size of the calorimeter (the input face is a square).

In addition a transverse uniform magnetic field is defined via G4GlobalMagFieldMessenger class.

Definition at line 54 of file B4DetectorConstruction.hh.

Constructor & Destructor Documentation

B4DetectorConstruction() [1/3]

B4DetectorConstruction::B4DetectorConstruction

(

)

Definition at line 64 of file B4DetectorConstruction.cxx.

  : G4VUserDetectorConstruction(),
    fAbsorberPV(0),
    fGapPV(0),
    fCheckOverlaps(true)
{}

~B4DetectorConstruction() [1/3]

B4DetectorConstruction::~B4DetectorConstruction ( )

virtual

Definition at line 73 of file B4DetectorConstruction.cxx.

{}

B4DetectorConstruction() [2/3]

B4DetectorConstruction::B4DetectorConstruction

(

)

~B4DetectorConstruction() [2/3]

virtual B4DetectorConstruction::~B4DetectorConstruction ( )

virtual

B4DetectorConstruction() [3/3]

B4DetectorConstruction::B4DetectorConstruction

(

)

~B4DetectorConstruction() [3/3]

virtual B4DetectorConstruction::~B4DetectorConstruction ( )

virtual

Member Function Documentation

Construct() [1/3]

G4VPhysicalVolume * B4DetectorConstruction::Construct ( )

virtual

Definition at line 77 of file B4DetectorConstruction.cxx.

{
  // Define materials
  DefineMaterials();
 
  // Define volumes
  return DefineVolumes();
}

ConstructSDandField() [1/3]

void B4DetectorConstruction::ConstructSDandField ( )

virtual

Definition at line 258 of file B4DetectorConstruction.cxx.

{
  // Added for VMC
  TG4GeometryManager::Instance()->ConstructSDandField();
  /*
    // Create global magnetic field messenger.
    // Uniform magnetic field is then created automatically if
    // the field value is not zero.
    G4ThreeVector fieldValue = G4ThreeVector();
    fMagFieldMessenger = new G4GlobalMagFieldMessenger(fieldValue);
    fMagFieldMessenger->SetVerboseLevel(1);
 
    // Register the field messenger for deleting
    G4AutoDelete::Register(fMagFieldMessenger);
  */
}

GetAbsorberPV() [1/3]

const G4VPhysicalVolume * B4DetectorConstruction::GetAbsorberPV ( ) const

inline

Definition at line 88 of file B4DetectorConstruction.hh.

                                                                            { 
  return fAbsorberPV; 
}

GetGapPV() [1/3]

const G4VPhysicalVolume * B4DetectorConstruction::GetGapPV ( ) const

inline

Definition at line 92 of file B4DetectorConstruction.hh.

                                                                        { 
  return fGapPV; 
}

DefineMaterials() [1/3]

void B4DetectorConstruction::DefineMaterials ( )

private

Definition at line 88 of file B4DetectorConstruction.cxx.

{
  // Lead material defined using NIST Manager
  G4NistManager* nistManager = G4NistManager::Instance();
  nistManager->FindOrBuildMaterial("G4_Pb");
 
  // Liquid argon material
  G4double a; // mass of a mole;
  G4double z; // z=mean number of protons;
  G4double density;
  new G4Material(
    "liquidArgon", z = 18., a = 39.95 * g / mole, density = 1.390 * g / cm3);
  // The argon by NIST Manager is a gas with a different density
 
  // Vacuum
  new G4Material("Galactic", z = 1., a = 1.01 * g / mole,
    density = universe_mean_density, kStateGas, 2.73 * kelvin, 3.e-18 * pascal);
 
  // Print materials
  G4cout << *(G4Material::GetMaterialTable()) << G4endl;
}

DefineVolumes() [1/3]

G4VPhysicalVolume * B4DetectorConstruction::DefineVolumes ( )

private

Definition at line 112 of file B4DetectorConstruction.cxx.

{
  // Geometry parameters
  G4int nofLayers = 10;
  G4double absoThickness = 10. * mm;
  G4double gapThickness = 5. * mm;
  G4double calorSizeXY = 10. * cm;
 
  G4double layerThickness = absoThickness + gapThickness;
  G4double calorThickness = nofLayers * layerThickness;
  G4double worldSizeXY = 1.2 * calorSizeXY;
  G4double worldSizeZ = 1.2 * calorThickness;
 
  // Get materials
  G4Material* defaultMaterial = G4Material::GetMaterial("Galactic");
  G4Material* absorberMaterial = G4Material::GetMaterial("G4_Pb");
  G4Material* gapMaterial = G4Material::GetMaterial("liquidArgon");
 
  if (!defaultMaterial || !absorberMaterial || !gapMaterial) {
    G4ExceptionDescription msg;
    msg << "Cannot retrieve materials already defined.";
    G4Exception("B4DetectorConstruction::DefineVolumes()", "MyCode0001",
      FatalException, msg);
  }
 
  //
  // World
  //
  G4VSolid* worldS = new G4Box("World",                // its name
    worldSizeXY / 2, worldSizeXY / 2, worldSizeZ / 2); // its size
 
  G4LogicalVolume* worldLV = new G4LogicalVolume(worldS, // its solid
    defaultMaterial,                                     // its material
    "WRLD");                                             // its name
 
  G4VPhysicalVolume* worldPV = new G4PVPlacement(0, // no rotation
    G4ThreeVector(),                                // at (0,0,0)
    worldLV,                                        // its logical volume
    "WRLD",                                         // its name
    0,                                              // its mother  volume
    false,                                          // no boolean operation
    0,                                              // copy number
    fCheckOverlaps);                                // checking overlaps
 
  //
  // Calorimeter
  //
  G4VSolid* calorimeterS = new G4Box("Calorimeter",        // its name
    calorSizeXY / 2, calorSizeXY / 2, calorThickness / 2); // its size
 
  G4LogicalVolume* calorLV = new G4LogicalVolume(calorimeterS, // its solid
    defaultMaterial,                                           // its material
    "CALO");                                                   // its name
 
  new G4PVPlacement(0, // no rotation
    G4ThreeVector(),   // at (0,0,0)
    calorLV,           // its logical volume
    "CALO",            // its name
    worldLV,           // its mother  volume
    false,             // no boolean operation
    0,                 // copy number
    fCheckOverlaps);   // checking overlaps
 
  //
  // Layer
  //
  G4VSolid* layerS = new G4Box("Layer",                    // its name
    calorSizeXY / 2, calorSizeXY / 2, layerThickness / 2); // its size
 
  G4LogicalVolume* layerLV = new G4LogicalVolume(layerS, // its solid
    defaultMaterial,                                     // its material
    "LAYE");                                             // its name
 
  new G4PVReplica("LAYE", // its name
    layerLV,              // its logical volume
    calorLV,              // its mother
    kZAxis,               // axis of replication
    nofLayers,            // number of replica
    layerThickness);      // witdth of replica
 
  //
  // Absorber
  //
  G4VSolid* absorberS = new G4Box("Abso",                 // its name
    calorSizeXY / 2, calorSizeXY / 2, absoThickness / 2); // its size
 
  G4LogicalVolume* absorberLV = new G4LogicalVolume(absorberS, // its solid
    absorberMaterial,                                          // its material
    "ABSO");                                                   // its name
 
  fAbsorberPV = new G4PVPlacement(0,          // no rotation
    G4ThreeVector(0., 0., -gapThickness / 2), // its position
    absorberLV,                               // its logical volume
    "ABSO",                                   // its name
    layerLV,                                  // its mother  volume
    false,                                    // no boolean operation
    0,                                        // copy number
    fCheckOverlaps);                          // checking overlaps
 
  //
  // Gap
  //
  G4VSolid* gapS = new G4Box("Gap",                      // its name
    calorSizeXY / 2, calorSizeXY / 2, gapThickness / 2); // its size
 
  G4LogicalVolume* gapLV = new G4LogicalVolume(gapS, // its solid
    gapMaterial,                                     // its material
    "GAPX");                                         // its name
 
  fGapPV = new G4PVPlacement(0,               // no rotation
    G4ThreeVector(0., 0., absoThickness / 2), // its position
    gapLV,                                    // its logical volume
    "GAPX",                                   // its name
    layerLV,                                  // its mother  volume
    false,                                    // no boolean operation
    0,                                        // copy number
    fCheckOverlaps);                          // checking overlaps
 
  //
  // print parameters
  //
  G4cout << "\n------------------------------------------------------------"
         << "\n---> The calorimeter is " << nofLayers << " layers of: [ "
         << absoThickness / mm << "mm of " << absorberMaterial->GetName()
         << " + " << gapThickness / mm << "mm of " << gapMaterial->GetName()
         << " ] "
         << "\n------------------------------------------------------------\n";
 
  //
  // Visualization attributes
  //
  worldLV->SetVisAttributes(G4VisAttributes::GetInvisible());
 
  G4VisAttributes* simpleBoxVisAtt =
    new G4VisAttributes(G4Colour(1.0, 1.0, 1.0));
  simpleBoxVisAtt->SetVisibility(true);
  calorLV->SetVisAttributes(simpleBoxVisAtt);
 
  //
  // Always return the physical World
  //
  return worldPV;
}

Construct() [2/3]

virtual G4VPhysicalVolume * B4DetectorConstruction::Construct ( )

virtual

ConstructSDandField() [2/3]

virtual void B4DetectorConstruction::ConstructSDandField ( )

virtual

GetAbsorberPV() [2/3]

const G4VPhysicalVolume * B4DetectorConstruction::GetAbsorberPV

(

)

const

GetGapPV() [2/3]

const G4VPhysicalVolume * B4DetectorConstruction::GetGapPV

(

)

const

DefineMaterials() [2/3]

void B4DetectorConstruction::DefineMaterials ( )

private

DefineVolumes() [2/3]

G4VPhysicalVolume * B4DetectorConstruction::DefineVolumes ( )

private

Construct() [3/3]

virtual G4VPhysicalVolume * B4DetectorConstruction::Construct ( )

virtual

ConstructSDandField() [3/3]

virtual void B4DetectorConstruction::ConstructSDandField ( )

virtual

GetAbsorberPV() [3/3]

const G4VPhysicalVolume * B4DetectorConstruction::GetAbsorberPV

(

)

const

GetGapPV() [3/3]

const G4VPhysicalVolume * B4DetectorConstruction::GetGapPV

(

)

const

DefineMaterials() [3/3]

void B4DetectorConstruction::DefineMaterials ( )

private

DefineVolumes() [3/3]

G4VPhysicalVolume * B4DetectorConstruction::DefineVolumes ( )

private

Member Data Documentation

fMagFieldMessenger

G4ThreadLocal G4GlobalMagFieldMessenger * B4DetectorConstruction::fMagFieldMessenger = 0

staticprivate

Definition at line 77 of file B4DetectorConstruction.hh.

fAbsorberPV

G4VPhysicalVolume * B4DetectorConstruction::fAbsorberPV

private

Definition at line 80 of file B4DetectorConstruction.hh.

fGapPV

G4VPhysicalVolume * B4DetectorConstruction::fGapPV

private

Definition at line 81 of file B4DetectorConstruction.hh.

fCheckOverlaps

G4bool B4DetectorConstruction::fCheckOverlaps

private

Definition at line 83 of file B4DetectorConstruction.hh.

---

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

• examples/E03/E03a/geant4/include/B4DetectorConstruction.hh
• examples/E03/E03b/geant4/include/B4DetectorConstruction.hh
• examples/E03/E03c/geant4/include/B4DetectorConstruction.hh
• examples/E03/E03a/geant4/src/B4DetectorConstruction.cxx
• examples/E03/E03b/geant4/src/B4DetectorConstruction.cxx
• examples/E03/E03c/geant4/src/B4DetectorConstruction.cxx