include/Geant4/G4TwistTrapParallelSide.hh

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 // G4TwistTrapParallelSide
0027 //
0028 // Class description:
0029 //
0030 // Class describing a twisted boundary surface for a trapezoid.
0031
0032 // Author: 27-Oct-2004 - O.Link (Oliver.Link@cern.ch)
0033 // --------------------------------------------------------------------
0034 #ifndef G4TWISTTRAPPARALLELSIDE_HH
0035 #define G4TWISTTRAPPARALLELSIDE_HH
0036
0037 #include "G4VTwistSurface.hh"
0038
0039 #include <vector>
0040
0041 class G4TwistTrapParallelSide : public G4VTwistSurface
0042 {
0043   public:
0044
0045     G4TwistTrapParallelSide(const G4String& name,
0046                               G4double  PhiTwist, // twist angle
0047                               G4double  pDz,      // half z lenght
0048                               G4double  pTheta, // direction between end planes
0049                               G4double  pPhi,   // by polar and azimutal angles
0050                               G4double  pDy1,     // half y length at -pDz
0051                               G4double  pDx1,     // half x length at -pDz,-pDy
0052                               G4double  pDx2,     // half x length at -pDz,+pDy
0053                               G4double  pDy2,     // half y length at +pDz
0054                               G4double  pDx3,     // half x length at +pDz,-pDy
0055                               G4double  pDx4,     // half x length at +pDz,+pDy
0056                               G4double  pAlph,    // tilt angle at +pDz
0057                               G4double  AngleSide // parity
0058                             );
0059
0060     ~G4TwistTrapParallelSide() override;
0061
0062     G4ThreeVector  GetNormal(const G4ThreeVector& xx,
0063                                    G4bool isGlobal = false) override ;
0064
0065     G4int DistanceToSurface(const G4ThreeVector& gp,
0066                             const G4ThreeVector& gv,
0067                                   G4ThreeVector  gxx[],
0068                                   G4double  distance[],
0069                                   G4int     areacode[],
0070                                   G4bool    isvalid[],
0071                             EValidate validate = kValidateWithTol) override;
0072
0073     G4int DistanceToSurface(const G4ThreeVector& gp,
0074                                   G4ThreeVector  gxx[],
0075                                   G4double       distance[],
0076                                   G4int          areacode[]) override;
0077
0078     G4TwistTrapParallelSide(__void__&);
0079       // Fake default constructor for usage restricted to direct object
0080       // persistency for clients requiring preallocation of memory for
0081       // persistifiable objects.
0082
0083   private:
0084
0085     G4int GetAreaCode(const G4ThreeVector& xx,
0086                             G4bool         withTol = true) override;
0087     void SetCorners() override;
0088     void SetBoundaries() override;
0089
0090     void GetPhiUAtX(const G4ThreeVector& p, G4double& phi, G4double& u);
0091     G4ThreeVector ProjectPoint(const G4ThreeVector& p,
0092                                      G4bool isglobal = false);
0093
0094     inline G4ThreeVector SurfacePoint(G4double phi, G4double u,
0095                                       G4bool isGlobal = false) override;
0096     inline G4double GetBoundaryMin(G4double phi) override;
0097     inline G4double GetBoundaryMax(G4double phi) override;
0098     inline G4double GetSurfaceArea() override;
0099     void GetFacets( G4int m, G4int n, G4double xyz[][3],
0100                     G4int faces[][4], G4int iside ) override;
0101
0102     inline G4ThreeVector NormAng(G4double phi, G4double u);
0103     inline G4double GetValueB(G4double phi) ;
0104     inline G4double Xcoef(G4double u);
0105       // To calculate the w(u) function
0106
0107   private:
0108
0109     G4double fTheta;
0110     G4double fPhi ;
0111
0112     G4double fDy1;
0113     G4double fDx1;
0114     G4double fDx2;
0115
0116     G4double fDy2;
0117     G4double fDx3;
0118     G4double fDx4;
0119
0120     G4double fDz;         // Half-length along the z axis
0121
0122     G4double fAlph;
0123     G4double fTAlph;      // std::tan(fAlph)
0124
0125     G4double fPhiTwist;   // twist angle ( dphi in surface equation)
0126
0127     G4double fAngleSide;
0128
0129     G4double fdeltaX;
0130     G4double fdeltaY;
0131
0132     G4double fDx4plus2;  // fDx4 + fDx2  == a2/2 + a1/2
0133     G4double fDx4minus2; // fDx4 - fDx2          -
0134     G4double fDx3plus1;  // fDx3 + fDx1  == d2/2 + d1/2
0135     G4double fDx3minus1; // fDx3 - fDx1          -
0136     G4double fDy2plus1;  // fDy2 + fDy1  == b2/2 + b1/2
0137     G4double fDy2minus1; // fDy2 - fDy1          -
0138     G4double fa1md1;     // 2 fDx2 - 2 fDx1  == a1 - d1
0139     G4double fa2md2;     // 2 fDx4 - 2 fDx3
0140 };
0141
0142 //========================================================
0143 // inline functions
0144 //========================================================
0145
0146
0147 inline
0148 G4double G4TwistTrapParallelSide::GetValueB(G4double phi)
0149 {
0150   return ( fDy2plus1 + fDy2minus1 * ( 2 * phi ) / fPhiTwist ) ;
0151 }
0152
0153 inline
0154 G4double G4TwistTrapParallelSide::Xcoef(G4double phi)
0155 {
0156   return GetValueB(phi)/2. ;
0157 }
0158
0159 inline G4ThreeVector
0160 G4TwistTrapParallelSide::
0161 SurfacePoint( G4double phi, G4double u, G4bool isGlobal )
0162 {
0163   // function to calculate a point on the surface, given by parameters phi,u
0164
0165   G4ThreeVector SurfPoint ( u*std::cos(phi) - Xcoef(phi)*std::sin(phi)
0166                           + fdeltaX*phi/fPhiTwist,
0167                             u*std::sin(phi) + Xcoef(phi)*std::cos(phi)
0168                           + fdeltaY*phi/fPhiTwist,
0169                             2*fDz*phi/fPhiTwist  );
0170   if (isGlobal) { return (fRot * SurfPoint + fTrans); }
0171   return SurfPoint;
0172 }
0173
0174
0175 inline
0176 G4double G4TwistTrapParallelSide::GetBoundaryMin(G4double phi)
0177 {
0178   return  -(fPhiTwist*(fDx2 + fDx4 - fDy2plus1*fTAlph)
0179           + 2*fDx4minus2*phi - 2*fDy2minus1*fTAlph*phi)/(2.*fPhiTwist) ;
0180 }
0181
0182 inline
0183 G4double G4TwistTrapParallelSide::GetBoundaryMax(G4double phi)
0184 {
0185   return (fDx2 + fDx4 + fDy2plus1*fTAlph)/ 2.
0186        + ((fDx4minus2 + fDy2minus1*fTAlph)*phi)/fPhiTwist ;
0187 }
0188
0189 inline
0190 G4double  G4TwistTrapParallelSide::GetSurfaceArea()
0191 {
0192   return 2*fDx4plus2*fDz ;
0193 }
0194
0195 inline
0196 G4ThreeVector G4TwistTrapParallelSide::NormAng( G4double phi, G4double u )
0197 {
0198   // function to calculate the norm at a given point on the surface
0199   // replace a1-d1
0200
0201   G4ThreeVector nvec(-2*fDz*std::sin(phi) ,
0202                       2*fDz*std::cos(phi) ,
0203                      -(fDy2minus1 + fPhiTwist*u + fdeltaY*std::cos(phi)
0204                      -fdeltaX*std::sin(phi))) ;
0205   return nvec.unit();
0206 }
0207
0208 #endif