113 G4double cosLambda = std::cos( sclambda );
114 G4double sinLambda = std::sin( sclambda );
115 G4double sinPhi = std::sin( scphi );
116 G4double cosPhi = std::cos( scphi );
122 G4ThreeVector vTN( cosLambda*cosPhi, cosLambda*sinPhi,sinLambda );
127 if(
iverbose >= 4)
G4cout <<
" SC2SD: vTN " << vTN <<
" vUN " << vUN <<
" vVN " << vVN <<
G4endl;
141 G4double T1R = 1. / ( vTN * vectorU );
171 G4cout <<
" UI " << UI <<
" VI " << VI <<
" vectorU " << vectorU <<
G4endl;
177 transfM[1][3] = -UI*( VK*cosz-UK*sinz)*T3R;
178 transfM[1][4] = -VI*( VK*cosz-UK*sinz)*T3R;
179 transfM[2][3] = UI*( VJ*cosz-UJ*sinz)*T3R;
180 transfM[2][4] = VI*( VJ*cosz-UJ*sinz)*T3R;
186 transfM[1][1] = -UK*T2R;
187 transfM[1][2] = VK*cosLambda*T2R;
188 transfM[2][1] = UJ*T2R;
189 transfM[2][2] = -VJ*cosLambda*T2R;
190 transfM[3][3] = VK*T1R;
191 transfM[3][4] = -UK*T1R;
192 transfM[4][3] = -VJ*T1R;
193 transfM[4][4] = UJ*T1R;
227 std::ios::fmtflags oldFlags = out.flags();
228 out.setf(std::ios::fixed,std::ios::floatfield);
CLHEP::Hep3Vector G4ThreeVector
static const G4double pos
std::ostream & operator<<(std::ostream &, const BasicVector3D< float > &)
BasicVector3D< T > cross(const BasicVector3D< T > &v) const
G4double GetLambda() const
G4ErrorSymMatrix G4ErrorTrajErr
G4ErrorMatrix BuildErrorMatrix(G4ErrorFreeTrajState &tpSC, const G4Vector3D &vecV, const G4Vector3D &vecW)
virtual void GetFieldValue(const double Point[4], double *fieldArr) const =0
G4Vector3D GetVectorW() const
G4ErrorSurfaceTrajParam fTrajParam
G4FieldManager * GetFieldManager() const
void DumpPosMomError(std::ostream &out=G4cout) const
G4ErrorFreeTrajParam GetParameters() const
const G4Field * GetDetectorField() const
G4ErrorSymMatrix similarity(const G4ErrorMatrix &m1) const
static G4TransportationManager * GetTransportationManager()
static constexpr double c_light
G4ErrorTrajErr GetError() const
G4Vector3D GetVectorV() const
static constexpr double cm
G4GLOB_DLL std::ostream G4cout
G4ErrorSurfaceTrajState(const G4String &partType, const G4Point3D &pos, const G4Vector3D &mom, const G4Plane3D &plane, const G4ErrorTrajErr &errmat=G4ErrorTrajErr(5, 0))
virtual void Dump(std::ostream &out=G4cout) const
static constexpr double GeV
static constexpr double km
static constexpr double pi
static constexpr double tesla
static G4ErrorPropagatorData * GetErrorPropagatorData()