2// ********************************************************************
3// * License and Disclaimer *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
26// class G4ParameterisedNavigation Inline implementation
28// --------------------------------------------------------------------
30// ********************************************************************
31// IdentifyAndPlaceSolid
32// ********************************************************************
34inline G4VSolid* G4ParameterisedNavigation::
35IdentifyAndPlaceSolid( G4int num,
36 G4VPhysicalVolume *apparentPhys, // PhysV or PhysT
37 G4VPVParameterisation *curParam )
39 G4VSolid *sampleSolid;
41 sampleSolid = curParam->ComputeSolid(num, apparentPhys);
42 sampleSolid->ComputeDimensions(curParam, num, apparentPhys);
43 curParam->ComputeTransformation(num, apparentPhys);
48// ********************************************************************
50// ********************************************************************
52inline G4SmartVoxelNode*
53G4ParameterisedNavigation::ParamVoxelLocate( G4SmartVoxelHeader* pHead,
54 const G4ThreeVector& localPoint )
56 // If no parameterisation axis is specified, adopt default
57 // location strategy as for placements
59 if ( pHead->GetParamAxis()==kUndefined )
61 fVoxelNode = G4VoxelNavigation::VoxelLocate(pHead,localPoint);
65 G4double targetHeaderMin, targetHeaderNodeWidth;
66 G4int targetHeaderNoSlices, targetNodeNo;
67 EAxis targetHeaderAxis;
69 targetHeaderAxis = pHead->GetAxis();
70 targetHeaderNoSlices = G4int(pHead->GetNoSlices());
71 targetHeaderMin = pHead->GetMinExtent();
72 targetHeaderNodeWidth = (pHead->GetMaxExtent()-targetHeaderMin)
73 / targetHeaderNoSlices;
74 targetNodeNo = G4int ( (localPoint(targetHeaderAxis)-targetHeaderMin)
75 / targetHeaderNodeWidth );
76 // Rounding protection
82 else if ( targetNodeNo>=targetHeaderNoSlices )
84 targetNodeNo = targetHeaderNoSlices-1;
86 fVoxelAxis = targetHeaderAxis;
87 fVoxelNoSlices = targetHeaderNoSlices;
88 fVoxelSliceWidth = targetHeaderNodeWidth;
89 fVoxelNodeNo = targetNodeNo;
91 fVoxelNode = pHead->GetSlice(targetNodeNo)->GetNode();