How to input Lu-177 spectrum information

Hi all,

I want to measure Lu-177 Dose in G4Box comprised of bins (10x10x10)

I made the model and made a isotropic source and I want to input Lu-177 Beta minus (e-) spectrum value to measure dose after simulation ( BeamOn 1000000)

But I don’t know how to input Spectrum value on Geant4

Please help me

This is my code below

#include “DetectorConstruction.hh”

#include “G4Material.hh”
#include “G4SystemOfUnits.hh”
#include “G4NistManager.hh”
#include “G4Box.hh”
#include “G4LogicalVolume.hh”
#include “G4PVPlacement.hh”
#include “G4PVReplica.hh”
#include “G4AutoDelete.hh”

#include “G4GeometryManager.hh”
#include “G4PhysicalVolumeStore.hh”
#include “G4LogicalVolumeStore.hh”
#include “G4SolidStore.hh”

#include “G4VisAttributes.hh”
#include “G4Colour.hh”

: G4VUserDetectorConstruction()


G4VPhysicalVolume DetectorConstruction::Construct()
// World
sol_World = new G4Box(“World”, 0.5m, 0.5m, 0.5m);
AIR = G4NistManager::Instance()->FindOrBuildMaterial(“G4_AIR”);
G4LogicalVolume* lv_World = new G4LogicalVolume(sol_World, AIR, “World”);
G4VPhysicalVolume* pv_World = new G4PVPlacement(0, G4ThreeVector(0.0, 0.0, 0.0), lv_World, “World”, 0, false, 0);

// User geometry
// Water box (x: 10 cm, y: 10 cm, z: 10 cm, at (0, 0, 0))
G4VSolid* sol_Geom = new G4Box("Geom", 5.0*cm, 5.0*cm, 5.0*cm);
G4Material* WATER = G4NistManager::Instance()->FindOrBuildMaterial("G4_WATER");
G4LogicalVolume* lv_Geom = new G4LogicalVolume(sol_Geom, WATER, "Geom");
G4VPhysicalVolume* pv_Geom = new G4PVPlacement(0, G4ThreeVector(0.0, 0.0, 0.0*cm), lv_Geom, "Geom", lv_World, false, 10);

//Replica of Water box (x-dir 0.1cm slices)
G4VSolid* sol_RepX = new G4Box("RepX",0.5*cm,5.0*cm,5.0*cm);
G4LogicalVolume* lv_RepX = new G4LogicalVolume(sol_RepX,WATER,"RepX");
G4VPhysicalVolume* pv_RepX = new G4PVReplica("RepX",lv_RepX,lv_Geom,kXAxis,10,1.0*cm);

//Replica of water box (z-dir 0.1cm slices)
G4VSolid* sol_RepZ = new G4Box("RepZ",0.5*cm,5.0*cm,0.5*cm);
G4LogicalVolume* lv_RepZ = new G4LogicalVolume(sol_RepZ,WATER,"RepZ");
G4VPhysicalVolume* pv_RepZ = new G4PVReplica("RepZ",lv_RepZ,lv_RepX,kZAxis,10,1.0*cm);

//Replica of water box (y-dir 0.1cm slices)
G4VSolid* sol_RepY = new G4Box("RepY",0.5*cm,0.5*cm,0.5*cm);
G4LogicalVolume* lv_RepY = new G4LogicalVolume(sol_RepY,WATER,"RepY");
G4VPhysicalVolume* pv_RepY = new G4PVReplica("RepY",lv_RepY,lv_RepZ,kYAxis,10,1.0*cm);

G4VisAttributes* va_Geom = new G4VisAttributes(G4Colour(0.0,0.0,1.0));
va_Geom -> SetForceWireframe(true);
lv_Geom -> SetVisAttributes(va_Geom);

G4VisAttributes* va_Rep = new G4VisAttributes(G4Colour(1.0,1.0,1.0,0.3));
va_Rep -> SetForceWireframe(true);
lv_RepX -> SetVisAttributes(va_Rep);
    lv_RepZ -> SetVisAttributes(va_Rep);
    lv_RepY -> SetVisAttributes(va_Rep);        
    return pv_World;


void DetectorConstruction::ConstructSDandField()


#include “G4ParticleGun.hh”
#include “G4SystemOfUnits.hh”
#include “G4RandomTools.hh”
#include “G4Electron.hh”

#include “PrimaryGeneratorAction.hh”

: G4VUserPrimaryGeneratorAction()
fPrimary = new G4ParticleGun();

delete fPrimary;

void PrimaryGeneratorAction::GeneratePrimaries(G4Event *anEvent)



Thanks for your help


why not to start directly with Lu-177 ?
Here a macro for examples/extended/radioactivedecay/rdecay01
Lu177.mac.txt (299 Bytes)