Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/geometry/solids/specific/include/G4TwistTrapParallelSide.hh

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

  1 //
  2 // ********************************************************************
  3 // * License and Disclaimer                                           *
  4 // *                                                                  *
  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.                             *
 10 // *                                                                  *
 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.         *
 17 // *                                                                  *
 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 // ********************************************************************
 25 //
 26 // G4TwistTrapParallelSide
 27 //
 28 // Class description:
 29 //
 30 // Class describing a twisted boundary surface for a trapezoid.
 31 
 32 // Author: 27-Oct-2004 - O.Link (Oliver.Link@cern.ch)
 33 // --------------------------------------------------------------------
 34 #ifndef G4TWISTTRAPPARALLELSIDE_HH
 35 #define G4TWISTTRAPPARALLELSIDE_HH
 36 
 37 #include "G4VTwistSurface.hh"
 38 
 39 #include <vector>
 40 
 41 class G4TwistTrapParallelSide : public G4VTwistSurface
 42 {
 43   public:
 44    
 45     G4TwistTrapParallelSide(const G4String& name,
 46                               G4double  PhiTwist, // twist angle
 47                               G4double  pDz,      // half z lenght
 48                               G4double  pTheta, // direction between end planes
 49                               G4double  pPhi,   // by polar and azimutal angles
 50                               G4double  pDy1,     // half y length at -pDz
 51                               G4double  pDx1,     // half x length at -pDz,-pDy
 52                               G4double  pDx2,     // half x length at -pDz,+pDy
 53                               G4double  pDy2,     // half y length at +pDz
 54                               G4double  pDx3,     // half x length at +pDz,-pDy
 55                               G4double  pDx4,     // half x length at +pDz,+pDy
 56                               G4double  pAlph,    // tilt angle at +pDz
 57                               G4double  AngleSide // parity
 58                             );
 59   
 60     ~G4TwistTrapParallelSide() override;
 61    
 62     G4ThreeVector  GetNormal(const G4ThreeVector& xx,
 63                                    G4bool isGlobal = false) override ;   
 64    
 65     G4int DistanceToSurface(const G4ThreeVector& gp,
 66                             const G4ThreeVector& gv,
 67                                   G4ThreeVector  gxx[],
 68                                   G4double  distance[],
 69                                   G4int     areacode[],
 70                                   G4bool    isvalid[],
 71                             EValidate validate = kValidateWithTol) override;
 72                                                   
 73     G4int DistanceToSurface(const G4ThreeVector& gp,
 74                                   G4ThreeVector  gxx[],
 75                                   G4double       distance[],
 76                                   G4int          areacode[]) override;
 77 
 78     G4TwistTrapParallelSide(__void__&);
 79       // Fake default constructor for usage restricted to direct object
 80       // persistency for clients requiring preallocation of memory for
 81       // persistifiable objects.
 82 
 83   private:
 84 
 85     G4int GetAreaCode(const G4ThreeVector& xx, 
 86                             G4bool         withTol = true) override;
 87     void SetCorners() override;
 88     void SetBoundaries() override;
 89 
 90     void GetPhiUAtX(const G4ThreeVector& p, G4double& phi, G4double& u);
 91     G4ThreeVector ProjectPoint(const G4ThreeVector& p,
 92                                      G4bool isglobal = false);
 93 
 94     inline G4ThreeVector SurfacePoint(G4double phi, G4double u,
 95                                       G4bool isGlobal = false) override;
 96     inline G4double GetBoundaryMin(G4double phi) override;
 97     inline G4double GetBoundaryMax(G4double phi) override;
 98     inline G4double GetSurfaceArea() override;
 99     void GetFacets( G4int m, G4int n, G4double xyz[][3],
100                     G4int faces[][4], G4int iside ) override;
101 
102     inline G4ThreeVector NormAng(G4double phi, G4double u);
103     inline G4double GetValueB(G4double phi) ;
104     inline G4double Xcoef(G4double u);
105       // To calculate the w(u) function
106 
107   private:
108 
109     G4double fTheta;   
110     G4double fPhi ;
111 
112     G4double fDy1;   
113     G4double fDx1;     
114     G4double fDx2;     
115 
116     G4double fDy2;   
117     G4double fDx3;     
118     G4double fDx4;     
119 
120     G4double fDz;         // Half-length along the z axis
121 
122     G4double fAlph;
123     G4double fTAlph;      // std::tan(fAlph)
124     
125     G4double fPhiTwist;   // twist angle ( dphi in surface equation)
126 
127     G4double fAngleSide;
128 
129     G4double fdeltaX;
130     G4double fdeltaY;
131 
132     G4double fDx4plus2;  // fDx4 + fDx2  == a2/2 + a1/2
133     G4double fDx4minus2; // fDx4 - fDx2          -
134     G4double fDx3plus1;  // fDx3 + fDx1  == d2/2 + d1/2
135     G4double fDx3minus1; // fDx3 - fDx1          -
136     G4double fDy2plus1;  // fDy2 + fDy1  == b2/2 + b1/2
137     G4double fDy2minus1; // fDy2 - fDy1          -
138     G4double fa1md1;     // 2 fDx2 - 2 fDx1  == a1 - d1
139     G4double fa2md2;     // 2 fDx4 - 2 fDx3 
140 };   
141 
142 //========================================================
143 // inline functions
144 //========================================================
145 
146 
147 inline 
148 G4double G4TwistTrapParallelSide::GetValueB(G4double phi) 
149 {
150   return ( fDy2plus1 + fDy2minus1 * ( 2 * phi ) / fPhiTwist ) ;
151 }
152 
153 inline
154 G4double G4TwistTrapParallelSide::Xcoef(G4double phi)
155 {
156   return GetValueB(phi)/2. ;
157 }
158 
159 inline G4ThreeVector
160 G4TwistTrapParallelSide::
161 SurfacePoint( G4double phi, G4double u, G4bool isGlobal ) 
162 {
163   // function to calculate a point on the surface, given by parameters phi,u
164 
165   G4ThreeVector SurfPoint ( u*std::cos(phi) - Xcoef(phi)*std::sin(phi)
166                           + fdeltaX*phi/fPhiTwist,
167                             u*std::sin(phi) + Xcoef(phi)*std::cos(phi)
168                           + fdeltaY*phi/fPhiTwist,
169                             2*fDz*phi/fPhiTwist  );
170   if (isGlobal) { return (fRot * SurfPoint + fTrans); }
171   return SurfPoint;
172 }
173 
174 
175 inline
176 G4double G4TwistTrapParallelSide::GetBoundaryMin(G4double phi)
177 {
178   return  -(fPhiTwist*(fDx2 + fDx4 - fDy2plus1*fTAlph)
179           + 2*fDx4minus2*phi - 2*fDy2minus1*fTAlph*phi)/(2.*fPhiTwist) ;
180 }
181 
182 inline
183 G4double G4TwistTrapParallelSide::GetBoundaryMax(G4double phi)
184 {
185   return (fDx2 + fDx4 + fDy2plus1*fTAlph)/ 2.
186        + ((fDx4minus2 + fDy2minus1*fTAlph)*phi)/fPhiTwist ;
187 }
188 
189 inline
190 G4double  G4TwistTrapParallelSide::GetSurfaceArea()
191 {
192   return 2*fDx4plus2*fDz ;
193 }
194 
195 inline
196 G4ThreeVector G4TwistTrapParallelSide::NormAng( G4double phi, G4double u ) 
197 {
198   // function to calculate the norm at a given point on the surface
199   // replace a1-d1
200 
201   G4ThreeVector nvec(-2*fDz*std::sin(phi) ,
202                       2*fDz*std::cos(phi) ,
203                      -(fDy2minus1 + fPhiTwist*u + fdeltaY*std::cos(phi)
204                      -fdeltaX*std::sin(phi))) ;
205   return nvec.unit();
206 }
207 
208 #endif
209