47 #ifndef __G4TWISTTUBSHYPESIDE__
48 #define __G4TWISTTUBSHYPESIDE__
61 const G4int handedness,
101 G4bool isGlobal =
false) ;
105 G4bool isglobal =
false)
const ;
108 G4bool isGlobal =
false) ;
128 virtual void SetCorners();
130 virtual void SetCorners(
G4double EndInnerRadius[2],
135 virtual void SetBoundaries();
170 return std::sqrt(fR02 + tmpp.
z() * tmpp.
z() * fTan2Stereo);
177 G4double rho = std::sqrt(fR02 + z * z * fTan2Stereo) ;
179 G4ThreeVector SurfPoint (rho*std::cos(phi), rho*std::sin(phi), z) ;
181 if (isGlobal) {
return (
fRot * SurfPoint +
fTrans); }
191 return std::atan2( lowerlimit.
y(), lowerlimit.
x() ) ;
200 return std::atan2( upperlimit.
y(), upperlimit.
x() ) ;
virtual G4double GetSurfaceArea()
virtual void GetFacets(G4int m, G4int n, G4double xyz[][3], G4int faces[][4], G4int iside)
virtual G4double GetBoundaryMin(G4double phi)
virtual G4ThreeVector GetBoundaryAtPZ(G4int areacode, const G4ThreeVector &p) const
virtual G4int DistanceToSurface(const G4ThreeVector &gp, const G4ThreeVector &gv, G4ThreeVector gxx[], G4double distance[], G4int areacode[], G4bool isvalid[], EValidate validate=kValidateWithTol)
virtual G4double GetBoundaryMax(G4double phi)
HepRotation inverse() const
virtual EInside Inside(const G4ThreeVector &gp)
virtual G4ThreeVector GetNormal(const G4ThreeVector &xx, G4bool isGlobal=false)
static const G4int sAxis0
static const G4int sAxisMin
virtual G4double GetRhoAtPZ(const G4ThreeVector &p, G4bool isglobal=false) const
G4TwistTubsHypeSide(const G4String &name, const G4RotationMatrix &rot, const G4ThreeVector &tlate, const G4int handedness, const G4double kappa, const G4double tanstereo, const G4double r0, const EAxis axis0=kPhi, const EAxis axis1=kZAxis, G4double axis0min=-kInfinity, G4double axis1min=-kInfinity, G4double axis0max=kInfinity, G4double axis1max=kInfinity)
static const G4int sAxisMax
virtual ~G4TwistTubsHypeSide()
virtual G4ThreeVector SurfacePoint(G4double, G4double, G4bool isGlobal=false)