Optimal physics list classes for Sb-124 BeO photoneutron source

Physics Processes, Models and Cross Sections
1

I am currently creating an Sb-124 BeO photo neutron source. This is built upon the Neutron source example in geant4 version 10.06.
I am using QGSP_BIC_HP for inelastic hadron physics and for the GammaNuclearPhysics I use:
G4ProcessManager* pManager = G4Gamma::Gamma()->GetProcessManager();
G4PhotoNuclearProcess* process = new G4PhotoNuclearProcess();
G4LENDInelastic* lend = new G4LENDInelastic(G4Gamma::Gamma());
lend->SetMaxEnergy(10MeV);
process->RegisterMe(lend);
G4LENDCombinedCrossSection lendXS = new G4LENDCombinedCrossSection(G4Gamma::Gamma());
process->AddDataSet(lendXS);
pManager->AddDiscreteProcess(process);

While I expect to get the emerging neutrons to have a peak energy around 24 keV, I am getting a peak around 0.1 MeV. Am I missing something when defining the physics processes?