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

None 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)

1 Feedback: Re: About the MedLinac Advanced Example!!!!!   (Michela Piergentili - 02 Jun, 2005)
 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 ]