Geant4 Cross Reference |
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 // 27 // 28 // 29 // Class description: 30 // 31 // Represents a trajectory state on a surface. 32 // It can be represented by the 5 variables: 33 // 1/p, v', w', v, w 34 // where v'=dv/du and w'=dw/du in an orthonormal coordinate system 35 // with axis u, v and w. 36 37 // History: 38 // 39 // - Created: P. Arce 40 // -------------------------------------------------------------------- 41 42 #ifndef G4ErrorSurfaceTrajState_hh 43 #define G4ErrorSurfaceTrajState_hh 44 45 #include "globals.hh" 46 47 #include "G4ErrorTrajState.hh" 48 #include "G4ErrorSurfaceTrajParam.hh" 49 #include "G4ErrorFreeTrajState.hh" 50 51 #include "G4Point3D.hh" 52 #include "G4Vector3D.hh" 53 #include "G4Plane3D.hh" 54 55 class G4ErrorSurfaceTrajState : public G4ErrorTrajState 56 { 57 public: // with description 58 G4ErrorSurfaceTrajState(const G4String& partType, const G4Point3D& pos, 59 const G4Vector3D& mom, const G4Plane3D& plane, 60 const G4ErrorTrajErr& errmat = G4ErrorTrajErr(5, 0)); 61 // Constructor by providing particle, position, momentum and 62 // G4Plane3D surface 63 64 G4ErrorSurfaceTrajState(const G4String& partType, const G4Point3D& pos, 65 const G4Vector3D& mom, const G4Vector3D& vecV, 66 const G4Vector3D& vecW, 67 const G4ErrorTrajErr& errmat = G4ErrorTrajErr(5, 0)); 68 // Constructor by providing particle, position, momentum and 69 // two vectors on surface 70 71 G4ErrorSurfaceTrajState(G4ErrorFreeTrajState& tpSC, const G4Plane3D& plane); 72 // Constructor by providing G4ErrorFreeTrajState and G4Plane3D surface 73 74 G4ErrorSurfaceTrajState(G4ErrorFreeTrajState& tpSC, const G4Vector3D& vecV, 75 const G4Vector3D& vecW, G4ErrorMatrix& transfM); 76 // Constructor by providing G4ErrorFreeTrajState and two vectors on surface 77 78 ~G4ErrorSurfaceTrajState() {} 79 G4ErrorMatrix BuildErrorMatrix(G4ErrorFreeTrajState& tpSC, 80 const G4Vector3D& vecV, 81 const G4Vector3D& vecW); 82 // Build the error matrix from a free state plus the vectors of the surface 83 84 virtual void Dump(std::ostream& out = G4cout) const; 85 86 friend std::ostream& operator<<(std::ostream&, 87 const G4ErrorSurfaceTrajState& ts); 88 89 // Set and Get methods 90 91 G4ErrorSurfaceTrajParam GetParameters() const { return fTrajParam; } 92 void SetParameters(const G4Point3D& pos, const G4Vector3D& mom, 93 const G4Vector3D& vecV, const G4Vector3D& vecW) 94 { 95 fPosition = pos; 96 fMomentum = mom; 97 fTrajParam.SetParameters(pos, mom, vecV, vecW); 98 } 99 100 void SetParameters(const G4Point3D& pos, const G4Vector3D& mom, 101 const G4Plane3D& plane) 102 { 103 fPosition = pos; 104 fMomentum = mom; 105 fTrajParam.SetParameters(pos, mom, plane); 106 } 107 108 G4Vector3D GetVectorV() const { return fTrajParam.GetVectorV(); } 109 110 G4Vector3D GetVectorW() const { return fTrajParam.GetVectorW(); } 111 112 virtual void SetPosition(const G4Point3D pos) 113 { 114 SetParameters(pos, fMomentum, GetVectorV(), GetVectorW()); 115 } 116 117 virtual void SetMomentum(const G4Vector3D& mom) 118 { 119 SetParameters(fPosition, mom, GetVectorV(), GetVectorW()); 120 } 121 122 private: 123 void Init(); 124 // Define Trajectory State type and build charge 125 126 private: 127 G4ErrorSurfaceTrajParam fTrajParam; 128 }; 129 130 #endif 131