Problem 1303 - the new density is not correct in DicomDetectorConstruction.cc
Summary: the new density is not correct in DicomDetectorConstruction.cc
Status: RESOLVED INVALID
Alias: None
Product: Examples/Extended
Classification: Unclassified
Component: medical/DICOM (show other problems)
Version: 9.5
Hardware: All All
: P5 trivial
Assignee: Pedro.Arce
URL:
Depends on:
Blocks:
 
Reported: 2012-04-13 22:53 CEST by Fada Guan
Modified: 2012-10-21 10:33 CEST (History)
0 users

See Also:

Attachments
Add an attachment (proposed patch, testcase, etc.)

Note You need to log in before you can comment on or make changes to this problem.
Description Fada Guan 2012-04-13 22:53:57 CEST 
on line 372  in DicomDetectorConstruction.cc

densityBin = fDensityDiffs[mateID] * (G4int(density/fDensityDiffs[mateID])+0.5);

this new density calculated is not correct.

For example, the density of lungInhale should be lower than lungExhale.
However, the above conversion can result in a result that a new lungInhale material has the same density as a new lungExhale material.

part of the output information are:

Material : LungInhale__0.4499999881

Material : LungExhale__0.4499999881
Comment 1 Pedro.Arce 2012-08-23 13:55:26 CEST 
I am not the author of this code, but I do not see any problem in it. My suspicion is that in your DICOM file you have two voxels assigned to different materials but the densities are the same (or close). Can you please check it?

Pedro
Comment 2 Pedro.Arce 2012-10-21 10:33:29 CEST 
I am closing this problem, as I have received no answer to my request of more information.