When I tried to use the model G4LivermoreComptonModel or G4LivermorePolarizedComptonModel with the process G4ComptonScattering, it didn't work. The command which was the problem is : pos = dataMap.find(Z); of the function G4VCrossSectionHandler::FindValue called by G4double cs = crossSectionHandler->FindValue(G4int(Z), GammaEnergy) in the function G4LivermorePolarizedComptonModel::ComputeCrossSectionPerAtom. This function is called by the command InitialiseElementSelectors(particle,cuts); in the function G4LivermorePolarizedComptonModel::Initialise. If this command(InitialiseElementsSectors) is moved at the end of this function(Initialise) it works. Is it normal ?
Dear Helene, I can confirm that the problem is known and is fixed for the next Geant4 release, due this June. You solved it the correct way. Best regards, S. Incerti