Message: About the MedLinac Advanced Example!!!!! Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question About the MedLinac Advanced Example!!!!! 

Keywords: MedLinac, Phantom, Region of Interest
Forum: Documentation and Examples
Date: 12 May, 2005
From: Francisco Garcia <Francisco Garcia>

To whom can clear out this question!!!

This is cocerning the Advanced example MedLinac. In the MedLinacPhantomROGeometry.cc file there is next code:
.
.
.
.
 // variables for x division ...
  G4double halfXVoxelDimensionX = PhantomDimensionX/NumberOfVoxelsAlongX;
  G4double halfXVoxelDimensionZ = PhantomDimensionZ;
  G4double voxelXThickness = 2*halfXVoxelDimensionX;

 // variables for y division ...
  G4double halfZVoxelDimensionX = PhantomDimensionX;
  G4double halfZVoxelDimensionZ = PhantomDimensionZ/NumberOfVoxelsAlongZ;
  G4double voxelZThickness = 2*halfZVoxelDimensionZ;
.
.
.
But not for the y division!!!

Then with the RO definition for X, I am agree

 // X division first...
  G4Box *ROPhantomXDivision = new G4Box("ROPhantomXDivision",
					halfXVoxelDimensionX,
					halfXVoxelDimensionZ,
					halfXVoxelDimensionZ);

  G4LogicalVolume *ROPhantomXDivisionLog = new G4LogicalVolume(ROPhantomXDivision,H2O,"ROPhantomXDivisionLog",
                                        	       0,0,0);

  G4VPhysicalVolume *ROPhantomXDivisionPhys = new G4PVReplica("ROPhantomXDivisionPhys",ROPhantomXDivisionLog,
                 ROPhantomPhys, kXAxis,NumberOfVoxelsAlongX,                                                    
                                            voxelXThickness);

With the RO definition for Z, I am agree

  // ...then Z division

  G4Box *ROPhantomZDivision = new G4Box("ROPhantomZDivision",
					halfZVoxelDimensionX,
					halfZVoxelDimensionX,
					halfZVoxelDimensionZ);

  G4LogicalVolume *ROPhantomZDivisionLog = new G4LogicalVolume(ROPhantomZDivision,H2O,"ROPhantomZDivisionLog",
						       0,0,0);

  G4VPhysicalVolume *ROPhantomZDivisionPhys = new G4PVReplica("ROPhantomZDivisionPhys",ROPhantomZDivisionLog,
		ROPhantomXDivisionPhys,kZAxis,NumberOfVoxelsAlongZ,
						  voxelZThickness);

But for the RO definition for Y, I don't understand next thin, why it is used halfZVoxelDimensionX for all coordinates of the ROPhantomYDivision and then replicate this for NumberOfVoxelsAlongZ times???

Some figures shown that halfZVoxelDimensionX = 15 cm and NumberOfVoxelsAlongZ = 300 so this will means you are creating voxels out of phantom itself in that direction.

  // ...then Y  division

  G4Box *ROPhantomYDivision = new G4Box("ROPhantomYDivision",
					halfZVoxelDimensionX,
					halfZVoxelDimensionX,
					halfZVoxelDimensionX);

  G4LogicalVolume *ROPhantomYDivisionLog = new G4LogicalVolume(ROPhantomYDivision,H2O,"ROPhantomYDivisionLog",
						       0,0,0);

  ROPhantomYDivisionPhys = new G4PVReplica("ROPhantomYDivisionPhys",
				     ROPhantomYDivisionLog,
				     ROPhantomZDivisionPhys,kYAxis,
			      NumberOfVoxelsAlongZ,voxelZThickness);

For this run I am using Geant4.7.0.patch01 with gcc 3.2 200.. (Red Hat 8.0 3.2-7).

Best regards,

Francisco

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

 Add Message Add Message
to: "About the MedLinac Advanced Example!!!!!"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]