|Message: Re: DICOM: adding magnetic field produces NaN||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
Finding another issue: I decided to throw away all the NaNs before they get written into the EvtMap in the simulation where i have theGeometry->SetMagField(G4ThreeVector(0.,0.,0.)); and then compare the result with a simulation where i do not set any magnetic field at all (that runs without producing NaNs), so that i see if the difference is small enough that i can ignore the NaNs.
You find plots attached of the dose in a water phantom throwing away the NaNs, and the difference when subtracting the dose simulated without any magnetic field at all. The irritating part is that areas of positive values exist, while the dose where something was thrown away should always be smaller than where nothing is thrown away. - And this is obviously systematic in the buil-up region, not noise. The penumbras differ, too. The resolution you see in the plots is 0.5mmx0.5mmx3mm (512x512x70) voxels, and the radiation field 2cmx2cm, number of particles 10^9. Although it may be argued whether or not such a difference is clinically relevant, i thought that switching on the magnetic field transportation should not lead to so big effects..? Need to test what happens in more complex CT geometries before trusting the magnetic field results? Is this known? Is this effect amplified when B>0T?
The difference in integral dose is less than 0.5%, so that throwing away the NaNs seems to be basically possible.
All the best,
Attachment: http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/08/03/09.41-67441-dicomNaNTo0_xy.png http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/08/03/09.41-39258-dicomNaNTo0_xz.png http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/08/03/09.41-45356-dicomNaNTo0_yz.png http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/08/03/09.41-60019-om0MinusNoMagAtAll_xy.png http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/08/03/09.41-78036-om0MinusNoMagAtAll_xz.png http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/08/03/09.41-36915-om0MinusNoMagAtAll_yz.png