Problem 1210

Summary: no dose gets written to dicom.out file
Product: Examples/Extended Reporter: Derek Dolney <Derek.Dolney>
Component: medical/DICOMAssignee: Pedro.Arce
Status: RESOLVED FIXED    
Severity: major    
Priority: P3    
Version: 9.4   
Hardware: All   
OS: All   

Description Derek Dolney 2011-04-28 23:52:00 CEST
Since 4.9.4, this example no longer works. Building and running with the defaults, using run.mac, should write out dose to the file dicom.out. Indeed, it works with 4.9.3.p02.

May be related to bug 1209.
Comment 1 Derek Dolney 2011-04-28 23:55:32 CEST
May be related to Problem 1209.
Comment 2 Pedro.Arce 2011-04-29 12:00:55 CEST
There is a bug int the DICOM example: the name of the SD starts with a capital in one place and not in the other. Could you check if it happens the same in your case? It should read:

./src/DicomRunAction.cc:  theSDName.push_back(G4String("phantomSD"));
./src/DicomDetectorConstruction.cc:  G4String concreteSDname = "phantomSD";

Pedro
Comment 3 Derek Dolney 2011-04-29 20:25:46 CEST
This fixes the problem. Thank you.
Comment 4 Pedro.Arce 2011-05-03 12:56:56 CEST
I do not see any problem. Without skipping I get a longer list of voxels with score:

=============================================================
 Number of event processed : 100
=============================================================
 opened file dicom.out for dose output
     Voxel |   DoseDeposit  |
    7992     2.07542e-21 Gy
    8120     3.30981e-22 Gy
    8121     4.42141e-21 Gy
    8122     2.19112e-22 Gy
    8127     5.64535e-20 Gy
    8128     7.92449e-20 Gy
    8253     3.52617e-22 Gy
    8254     3.03714e-21 Gy
    8255     4.89962e-20 Gy
    8256     6.87428e-20 Gy
    24506     2.21257e-21 Gy
    24511     8.49559e-20 Gy
    24512     9.07961e-20 Gy
    24634     1.31314e-21 Gy
    24635     3.12405e-21 Gy
    24636     4.54652e-21 Gy
    24637     2.98413e-21 Gy
    24638     5.59209e-21 Gy
    24639     3.49854e-20 Gy
    24640     6.50744e-20 Gy
    24764     7.11161e-23 Gy
    24765     2.28296e-21 Gy
    24766     4.65128e-22 Gy
    40641     3.64866e-22 Gy
    40768     3.07089e-22 Gy
    40769     1.32696e-21 Gy
    40895     7.54732e-20 Gy
    40896     7.56551e-20 Gy
    41022     3.02737e-21 Gy
   41023     3.09145e-20 Gy
    41024     7.61299e-20 Gy
    41147     5.1506e-22 Gy
    41148     3.60304e-21 Gy
=============================================
Comment 5 Derek Dolney 2011-05-04 00:30:07 CEST
Oh, sorry. I am posting comments to Problem 1209, but Pedro, you are commenting here. I would close this bug, because with the SD name correction, the DICOM example then works out-of-the-box, so to speak, since voxel boundaries are skipped by default. Problem 1209 is not specific to the DICOM example: dose is not scored in any G4PhantomParameterisation when using SetSkipEqualMaterials(false).
Comment 6 Pedro.Arce 2011-06-06 11:32:26 CEST
 fix will be included in the next patch and releases