I instantiate a water material
G4NistManager* nist = G4NistManager::Instance();
G4Material* mymat = nist->FindOrBuildMaterial("G4_WATER")
and then, using
G4PenelopeRayleighModel::DumpFormFactorTable(mymat), I retrieve the Rayleigh scattering form factors which are calculated according to the independent atom model(IAM).
I then edit the form factors to add molecular interference (MI) effects and retrieve them back again to plot the MI form factors in comparison with the IAM form factors.
Common practice has it that the plots be normalised(divided) by sqrt(W), where W is the molecular weight of water but if I do that using 18.01528 g/mol then the results(first figure) are half of what they are reported in the literature(second figure)***SEE EDIT BELOW Geant4 implementation of inter-atomic interference effect in small-angle coherent X-ray scattering for materials of medical interest.
I notice units of “mol” missing from the y axis(W has units of g/mol and F has no units) in the literature so I am wondering why they disappear and if this is the reason why both Geant4 IAM and my MI form factors are half of those in the literature.
Results are not half. It is that if I use 9 g/mol instead of 18 g/mol that makes the results be equal.