00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038 inline
00039 G4VSolid* G4ParameterisedNavigation::
00040 IdentifyAndPlaceSolid( G4int num,
00041 G4VPhysicalVolume *apparentPhys,
00042 G4VPVParameterisation *curParam )
00043 {
00044 G4VSolid *sampleSolid;
00045
00046 sampleSolid = curParam->ComputeSolid(num, apparentPhys);
00047 sampleSolid->ComputeDimensions(curParam, num, apparentPhys);
00048 curParam->ComputeTransformation(num, apparentPhys);
00049
00050 return sampleSolid;
00051 }
00052
00053
00054
00055
00056
00057 inline
00058 G4SmartVoxelNode*
00059 G4ParameterisedNavigation::ParamVoxelLocate( G4SmartVoxelHeader* pHead,
00060 const G4ThreeVector& localPoint )
00061 {
00062
00063
00064
00065 if ( pHead->GetParamAxis()==kUndefined )
00066 {
00067 fVoxelNode = G4VoxelNavigation::VoxelLocate(pHead,localPoint);
00068 }
00069 else
00070 {
00071 G4double targetHeaderMin, targetHeaderNodeWidth;
00072 G4int targetHeaderNoSlices, targetNodeNo;
00073 EAxis targetHeaderAxis;
00074
00075 targetHeaderAxis = pHead->GetAxis();
00076 targetHeaderNoSlices = pHead->GetNoSlices();
00077 targetHeaderMin = pHead->GetMinExtent();
00078 targetHeaderNodeWidth = (pHead->GetMaxExtent()-targetHeaderMin)
00079 / targetHeaderNoSlices;
00080 targetNodeNo = G4int ( (localPoint(targetHeaderAxis)-targetHeaderMin)
00081 / targetHeaderNodeWidth );
00082
00083
00084 if ( targetNodeNo<0 )
00085 {
00086 targetNodeNo = 0;
00087 }
00088 else if ( targetNodeNo>=targetHeaderNoSlices )
00089 {
00090 targetNodeNo = targetHeaderNoSlices-1;
00091 }
00092 fVoxelAxis = targetHeaderAxis;
00093 fVoxelNoSlices = targetHeaderNoSlices;
00094 fVoxelSliceWidth = targetHeaderNodeWidth;
00095 fVoxelNodeNo = targetNodeNo;
00096 fVoxelHeader = pHead;
00097 fVoxelNode = pHead->GetSlice(targetNodeNo)->GetNode();
00098 }
00099 return fVoxelNode;
00100 }
00101