Computed tomography
From UCL HEP PBT Wiki
m |
m |
||
Line 58: | Line 58: | ||
</pre> | </pre> | ||
- | ; | + | ; Run macro run.mac |
- | + | ||
- | + | ||
- | + | ||
- | + | ||
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
Line 68: | Line 64: | ||
</pre> | </pre> | ||
- | * | + | == <span style="color:#000080"> How to analyze data </span> == |
+ | |||
+ | This is a [http://geant4.web.cern.ch/geant4/UserDocumentation/Doxygen/examples_doc/html/ExampleDICOM.html link] to the official DICOM explanation notes. We recommend to read it before proceeding with this tutorial. | ||
+ | |||
+ | === <span style="color:#000080"> Text files </span> === | ||
+ | |||
+ | The macro run.mac produces a text file '''dicom.out''' which includes the dose deposition in several voxels. | ||
+ | |||
+ | === <span style="color:#000080"> Run with different settings </span> === | ||
+ | |||
+ | ; Chose among several voxelization algorithms | ||
+ | |||
+ | * <span style="color:#ff0000"> G4RegularNavigation </span>: | ||
+ | |||
+ | '''G4RegularNavigation''' class is the default class for this tutorial. This algorithm skips frontiers between voxels when they have the same material i.e. "replacing group of voxels with a smaller number of larger voxels". | ||
+ | |||
+ | * <span style="color:#ff0000"> G4NestedParameterization </span>: | ||
+ | |||
+ | To run with this voxelization algorithm you need to set the variable '''DICOM_NESTED_PARAM''' | ||
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
Line 74: | Line 88: | ||
</pre> | </pre> | ||
- | You can check if | + | You can check if it is set correctly by typing in the command line |
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
Line 80: | Line 94: | ||
</pre> | </pre> | ||
- | Then, run the code: | + | Then, compile and run the code: |
<pre style="color: #800000; background-color: #dcdcdc"> | <pre style="color: #800000; background-color: #dcdcdc"> | ||
+ | [username@plus1 DICOM_build]$ make | ||
+ | |||
[username@plus1 DICOM_build]$ ./DICOM run.mac | [username@plus1 DICOM_build]$ ./DICOM run.mac | ||
</pre> | </pre> | ||
- | * | + | * <span style="color:#ff0000"> G4SmartVoxel / G4VoxelNavigation </span>: |
Using smart voxels required a huge amount of memory. At /home/username/DICOMFolder/DICOM_source/src/DicomRegularDetectorConstruction.cc set | Using smart voxels required a huge amount of memory. At /home/username/DICOMFolder/DICOM_source/src/DicomRegularDetectorConstruction.cc set | ||
Line 102: | Line 118: | ||
</pre> | </pre> | ||
- | + | ; Change the type, energy and position of the incident particle | |
- | + | At /home/username/DICOMFolder/DICOM_source/src/DicomPrimaryGeneratorAction.cc change the following lines: | |
- | + | <pre style="color: #800000; background-color: #dcdcdc"> | |
+ | G4ParticleDefinition* particle = particleTable->FindParticle(particleName="e-"); | ||
- | + | fParticleGun->SetParticleEnergy(100.*MeV); | |
+ | |||
+ | fParticleGun->SetParticlePosition(G4ThreeVector(0.,0.,0.)); | ||
+ | </pre> | ||
+ | |||
+ | You can chose among several particles like "gamma", "e+", "alpha" and "He3". Then, compile and run the code | ||
+ | |||
+ | <pre style="color: #800000; background-color: #dcdcdc"> | ||
+ | [username@plus1 DICOM_build]$ make | ||
+ | |||
+ | [username@plus1 DICOM_build]$ ./DICOM run.mac | ||
+ | </pre> | ||
=== <span style="color:#000080"> Visualisation </span> === | === <span style="color:#000080"> Visualisation </span> === |