23.4. Adding the Material

Table of Contents

Choosing Materials for the Gamma-Ray Detector

In the gamma ray detector example the material choice is what turns a simple piece of geometry into a working detector. In this chapter you focus on selecting and defining the materials that will be used in the world, the detector housing, and especially the scintillator crystal.

You do not redesign the geometry here. Instead, you connect that geometry to realistic materials using Geant4 classes and, where possible, the NIST material database.

Role of Materials in the Example

Every physical volume in your detector must be associated with a material. The material affects how gamma rays interact, how much energy they deposit, and how many secondary particles are produced. For a scintillation detector the most important properties at this stage are:

Atomic composition and density, which determine interaction probabilities such as photoelectric absorption and Compton scattering.
Whether it represents air, vacuum, or a dense crystal, which affects attenuation and detector efficiency.

Optical properties and scintillation yields are handled later, when you add optical physics. Here you prepare the base material definitions on which those properties will be added.

In the gamma detector example you typically define at least three materials: the world medium, the scintillator crystal, and optionally a housing or encapsulation material.

Using the NIST Material Database

Geant4 provides a convenient database of many common elements and materials through the G4NistManager class. Whenever possible, you should use this database instead of redefining standard materials by hand. It reduces mistakes and keeps your simulation consistent with recommended values.

You usually access the manager once in your DetectorConstruction::Construct() method:

cpp

auto nist = G4NistManager::Instance();

To retrieve a predefined material, you call:

cpp

G4Material* worldMat = nist->FindOrBuildMaterial("G4_AIR");
G4Material* crystalMat = nist->FindOrBuildMaterial("G4_CESIUM_IODIDE");

If the material exists in the NIST list, FindOrBuildMaterial will create it the first time it is requested and return the pointer. Subsequent calls simply return the already created material. Many standard detector materials exist with names starting with G4_, such as G4_WATER, G4_Al, G4_Pb, and several scintillators.

Always use G4NistManager::Instance()->FindOrBuildMaterial("G4_NAME") for standard materials. Do not redefine common materials manually unless you have a specific reason and understand the consequences.

In the example, you will typically get air for the world and a scintillator material such as CsI or NaI from this database.

World Material

The world material should represent the environment that surrounds your detector. For a laboratory detector simulation, this is often normal air. If you want to approximate vacuum you can either use a predefined low density material or define your own near vacuum material.

For this example a simple and realistic choice is standard air:

cpp

auto nist = G4NistManager::Instance();
G4Material* worldMat = nist->FindOrBuildMaterial("G4_AIR");

You then assign this material to the logical volume of the world, which you already created when building the world geometry:

cpp

auto worldLV = new G4LogicalVolume(worldSolid, worldMat, "WorldLV");

This choice is sufficient for a basic gamma detector study. If you later want to investigate effects of shielding gases or vacuum, you can change worldMat without modifying your geometry construction logic.

Choosing a Scintillator Material

The core of a gamma ray detector is the scintillator crystal that converts gamma energy into light. In this example you represent it as a single solid volume with a dedicated material. The choice of material affects the detection efficiency, energy resolution, and gamma energy response.

Some common inorganic scintillators included in Geant4’s NIST list are:

NIST name	Typical use
`G4_SODIUM_IODIDE`	NaI(Tl) gamma spectrometers
`G4_CESIUM_IODIDE`	CsI detectors, calorimeters
`G4_BGO`	Bismuth germanate, high Z, dense
`G4_LSO`	Lutetium oxyorthosilicate, PET-like

For a simple gamma ray detector, NaI(Tl) or CsI are good introductory choices because they have relatively high light yield and are historically well studied. Geant4’s NIST database provides base NaI and CsI compositions. The activator (Tl) is normally treated through optical properties and scintillation yields, which you will handle when you add optical physics. For now you only need the base material:

cpp

auto nist = G4NistManager::Instance();
// Example choice: CsI
G4Material* crystalMat = nist->FindOrBuildMaterial("G4_CESIUM_IODIDE");
// Alternatively, NaI
// G4Material* crystalMat = nist->FindOrBuildMaterial("G4_SODIUM_IODIDE");

You then assign this crystalMat to the logical volume of your scintillator, for example:

cpp

auto crystalLV = new G4LogicalVolume(crystalSolid, crystalMat, "CrystalLV");

By changing the NIST material name only, you can reuse the same geometry and quickly compare different scintillator types in later studies.

Keep the scintillator material consistent throughout your simulation. If you change the crystal material, ensure that all related settings such as production cuts, optical properties, and calibration assumptions are updated accordingly.

Defining Optional Housing and Other Materials

Real detectors often include a housing or encapsulation, such as aluminum or stainless steel, that surrounds the scintillator. For a beginner gamma detector simulation this housing can often be omitted or simplified, but if your geometry includes such a shell you need to assign it a material as well.

Again, you can use G4NistManager:

cpp

auto nist = G4NistManager::Instance();
G4Material* housingMat = nist->FindOrBuildMaterial("G4_Al"); // Aluminum

After that, you attach the material to the corresponding logical volume, for example:

cpp

auto housingLV = new G4LogicalVolume(housingSolid, housingMat, "HousingLV");

Other optional materials you might introduce later include:

A light guide between the crystal and photodetector, often a plastic or glass material such as G4_PLEXIGLASS or G4_GLASS_PLATE.
A coupling medium like optical grease, which you can approximate using a light material such as G4_WATER or a custom low Z compound.
A photodetector window, which can use a glass material from the NIST database.

For the core beginner example, it is enough to pick a single housing material if the geometry includes it and to keep the rest as air.

Creating a Custom Scintillator Material (Optional)

If you do not find an exact material in the NIST list or want to define a specific mixture, you can create a custom G4Material. This requires defining the constituent elements and material density. While this is slightly more advanced, it can be helpful to understand the pattern here so you recognize such code when reading examples.

A typical custom definition uses the following pattern:

Retrieve or define elements:

cpp

auto nist = G4NistManager::Instance();
G4Element* elNa = nist->FindOrBuildElement("Na");
G4Element* elI  = nist->FindOrBuildElement("I");

Build the material:

cpp

G4double density = 3.67 * g/cm3;
G4Material* NaI = new G4Material("NaI", density, 2);
NaI->AddElement(elNa, 1);
NaI->AddElement(elI, 1);

The arguments in G4Material are the name, density, and number of components. In AddElement, the last parameter is the number of atoms per formula unit when you use integer counts like this.

For the introductory gamma detector, using the predefined G4_SODIUM_IODIDE or G4_CESIUM_IODIDE is strongly recommended, and you can ignore this manual construction until you need a truly custom material.

When creating custom materials, always use correct densities and stoichiometry. Small mistakes in density or composition can significantly change interaction probabilities and lead to misleading detector performance.

Connecting Materials to Geometry in the Example

At this point your gamma detector example has three key steps inside DetectorConstruction::Construct():

Create the world solid and logical volume and assign G4_AIR as the world material.
Create the crystal solid and logical volume and assign a scintillator material such as G4_CESIUM_IODIDE.
Optionally create a housing solid and logical volume and assign a material such as G4_Al.

The code structure typically looks like this in simplified form:

cpp

auto nist = G4NistManager::Instance();
// World
G4Material* worldMat = nist->FindOrBuildMaterial("G4_AIR");
auto worldLV = new G4LogicalVolume(worldSolid, worldMat, "WorldLV");
// Scintillator crystal
G4Material* crystalMat = nist->FindOrBuildMaterial("G4_CESIUM_IODIDE");
auto crystalLV = new G4LogicalVolume(crystalSolid, crystalMat, "CrystalLV");
// Optional housing
// G4Material* housingMat = nist->FindOrBuildMaterial("G4_Al");
// auto housingLV = new G4LogicalVolume(housingSolid, housingMat, "HousingLV");

Once these associations are in place and the physical placements are defined, the physics list can apply gamma interaction processes based on realistic material properties. This directly influences the energy deposition that you will record later and the shapes of the spectra you will analyze.

By separating the material choice from the geometry and by relying on the NIST database, you create a gamma detector example that is both simple to understand and easy to modify for future studies.