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

/////////////////////// DetectorConstruction.cc////////////////////////////////
#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”

DetectorConstruction::DetectorConstruction()
: G4VUserDetectorConstruction()
{
}

DetectorConstruction::~DetectorConstruction()
{
}

G4VPhysicalVolume DetectorConstruction::Construct()
{
// World
G4VSolid
sol_World = new G4Box(“World”, 0.5m, 0.5m, 0.5m);
G4Material
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()
{

}

///////////////////PrimaryGeneratorAction.cc///////////////////
#include “G4ParticleGun.hh”
#include “G4SystemOfUnits.hh”
#include “G4RandomTools.hh”
#include “G4Electron.hh”

#include “PrimaryGeneratorAction.hh”

PrimaryGeneratorAction::PrimaryGeneratorAction()
: G4VUserPrimaryGeneratorAction()
{
fPrimary = new G4ParticleGun();
}

PrimaryGeneratorAction::~PrimaryGeneratorAction()
{
delete fPrimary;
}

void PrimaryGeneratorAction::GeneratePrimaries(G4Event *anEvent)
{
fPrimary->SetParticleDefinition(G4Electron::Definition());
fPrimary->SetParticleEnergy(*MeV);
fPrimary->SetParticlePosition(G4ThreeVector());
fPrimary->SetParticleMomentumDirection(G4RandomDirection());

fPrimary->GeneratePrimaryVertex(anEvent);

}

Thanks for your help

Junsu

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