Geant4 9.6.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
Public Member Functions | Friends | List of all members
G4ErrorSurfaceTrajState Class Reference

#include <G4ErrorSurfaceTrajState.hh>

+ Inheritance diagram for G4ErrorSurfaceTrajState:

Public Member Functions

 G4ErrorSurfaceTrajState (const G4String &partType, const G4Point3D &pos, const G4Vector3D &mom, const G4Plane3D &plane, const G4ErrorTrajErr &errmat=G4ErrorTrajErr(5, 0))
 
 G4ErrorSurfaceTrajState (const G4String &partType, const G4Point3D &pos, const G4Vector3D &mom, const G4Vector3D &vecV, const G4Vector3D &vecW, const G4ErrorTrajErr &errmat=G4ErrorTrajErr(5, 0))
 
 G4ErrorSurfaceTrajState (G4ErrorFreeTrajState &tpSC, const G4Plane3D &plane)
 
 G4ErrorSurfaceTrajState (G4ErrorFreeTrajState &tpSC, const G4Vector3D &vecV, const G4Vector3D &vecW)
 
 ~G4ErrorSurfaceTrajState ()
 
void BuildErrorMatrix (G4ErrorFreeTrajState &tpSC, const G4Vector3D &vecV, const G4Vector3D &vecW)
 
virtual void Dump (std::ostream &out=G4cout) const
 
G4ErrorSurfaceTrajParam GetParameters () const
 
void SetParameters (const G4Point3D &pos, const G4Vector3D &mom, const G4Vector3D &vecV, const G4Vector3D &vecW)
 
void SetParameters (const G4Point3D &pos, const G4Vector3D &mom, const G4Plane3D &plane)
 
G4Vector3D GetVectorV () const
 
G4Vector3D GetVectorW () const
 
virtual void SetPosition (const G4Point3D pos)
 
virtual void SetMomentum (const G4Vector3D &mom)
 
- Public Member Functions inherited from G4ErrorTrajState
 G4ErrorTrajState ()
 
 G4ErrorTrajState (const G4String &partType, const G4Point3D &pos, const G4Vector3D &mom, const G4ErrorTrajErr &errmat=G4ErrorTrajErr(5, 0))
 
virtual ~G4ErrorTrajState ()
 
void SetData (const G4String &partType, const G4Point3D &pos, const G4Vector3D &mom)
 
void BuildCharge ()
 
virtual G4int PropagateError (const G4Track *)
 
virtual G4int Update (const G4Track *)
 
void UpdatePosMom (const G4Point3D &pos, const G4Vector3D &mom)
 
void DumpPosMomError (std::ostream &out=G4cout) const
 
virtual void Dump (std::ostream &out=G4cout) const =0
 
const G4String & GetParticleType () const
 
void SetParticleType (const G4String &partType)
 
G4Point3D GetPosition () const
 
virtual void SetPosition (const G4Point3D pos)
 
G4Vector3D GetMomentum () const
 
virtual void SetMomentum (const G4Vector3D &mom)
 
G4ErrorTrajErr GetError () const
 
virtual void SetError (G4ErrorTrajErr em)
 
G4Track * GetG4Track () const
 
void SetG4Track (G4Track *trk)
 
G4double GetCharge () const
 
void SetCharge (G4double ch)
 
virtual G4eTSType GetTSType () const
 

Friends

std::ostream & operator<< (std::ostream &, const G4ErrorSurfaceTrajState &ts)
 

Additional Inherited Members

- Protected Attributes inherited from G4ErrorTrajState
G4String fParticleType
 
G4Point3D fPosition
 
G4Vector3D fMomentum
 
G4double fCharge
 
G4ErrorTrajErr fError
 
G4eTSType theTSType
 
G4Track * theG4Track
 
G4int iverbose
 

Detailed Description

Definition at line 56 of file G4ErrorSurfaceTrajState.hh.

Constructor & Destructor Documentation

◆ G4ErrorSurfaceTrajState() [1/4]

G4ErrorSurfaceTrajState::G4ErrorSurfaceTrajState ( const G4String &  partType,
const G4Point3D &  pos,
const G4Vector3D &  mom,
const G4Plane3D &  plane,
const G4ErrorTrajErr &  errmat = G4ErrorTrajErr(5,0) 
)

Definition at line 60 of file G4ErrorSurfaceTrajState.cc.

64 : G4ErrorTrajState( partType, pos, mom, errmat )
65{
66 Init();
67 fTrajParam = G4ErrorSurfaceTrajParam( pos, mom, plane );
68
69}

◆ G4ErrorSurfaceTrajState() [2/4]

G4ErrorSurfaceTrajState::G4ErrorSurfaceTrajState ( const G4String &  partType,
const G4Point3D &  pos,
const G4Vector3D &  mom,
const G4Vector3D &  vecV,
const G4Vector3D &  vecW,
const G4ErrorTrajErr &  errmat = G4ErrorTrajErr(5,0) 
)

Definition at line 48 of file G4ErrorSurfaceTrajState.cc.

52 : G4ErrorTrajState( partType, pos, mom, errmat )
53{
54 Init();
55 fTrajParam = G4ErrorSurfaceTrajParam( pos, mom, vecU, vecV );
56}

◆ G4ErrorSurfaceTrajState() [3/4]

G4ErrorSurfaceTrajState::G4ErrorSurfaceTrajState ( G4ErrorFreeTrajState &  tpSC,
const G4Plane3D &  plane 
)

Definition at line 73 of file G4ErrorSurfaceTrajState.cc.

75 : G4ErrorTrajState( tpSC.GetParticleType(), tpSC.GetPosition(),
76 tpSC.GetMomentum() )
77{
78 // fParticleType = tpSC.GetParticleType();
79 // fPosition = tpSC.GetPosition();
80 // fMomentum = tpSC.GetMomentum();
81 fTrajParam = G4ErrorSurfaceTrajParam( fPosition, fMomentum, plane );
82 Init();
83
84 //----- Get the error matrix in SC coordinates
85 BuildErrorMatrix( tpSC, GetVectorV(), GetVectorW() );
86}
G4ErrorSurfaceTrajState::BuildErrorMatrix
void BuildErrorMatrix(G4ErrorFreeTrajState &tpSC, const G4Vector3D &vecV, const G4Vector3D &vecW)
Definition: G4ErrorSurfaceTrajState.cc:105
G4ErrorTrajState::GetMomentum
G4Vector3D GetMomentum() const
Definition: G4ErrorTrajState.hh:99
G4ErrorTrajState::GetParticleType
const G4String & GetParticleType() const
Definition: G4ErrorTrajState.hh:89
G4ErrorTrajState::GetPosition
G4Point3D GetPosition() const
Definition: G4ErrorTrajState.hh:94

◆ G4ErrorSurfaceTrajState() [4/4]

G4ErrorSurfaceTrajState::G4ErrorSurfaceTrajState ( G4ErrorFreeTrajState &  tpSC,
const G4Vector3D &  vecV,
const G4Vector3D &  vecW 
)

Definition at line 90 of file G4ErrorSurfaceTrajState.cc.

93 : G4ErrorTrajState( tpSC.GetParticleType(), tpSC.GetPosition(),
94 tpSC.GetMomentum() )
95{
96 fTrajParam = G4ErrorSurfaceTrajParam( fPosition, fMomentum, vecU, vecV );
97 theTSType = G4eTS_OS;
98 //----- Get the error matrix in SC coordinates
99 BuildErrorMatrix( tpSC, vecU, vecV );
100}
G4eTS_OS
@ G4eTS_OS
Definition: G4ErrorTrajState.hh:47

◆ ~G4ErrorSurfaceTrajState()

G4ErrorSurfaceTrajState::~G4ErrorSurfaceTrajState ( )
inline

Definition at line 81 of file G4ErrorSurfaceTrajState.hh.

81{}

Member Function Documentation

◆ BuildErrorMatrix()

void G4ErrorSurfaceTrajState::BuildErrorMatrix ( G4ErrorFreeTrajState &  tpSC,
const G4Vector3D &  vecV,
const G4Vector3D &  vecW 
)

Definition at line 104 of file G4ErrorSurfaceTrajState.cc.

107{
108 G4double sclambda = tpSC.GetParameters().GetLambda();
109 G4double scphi = tpSC.GetParameters().GetPhi();
110 if( G4ErrorPropagatorData::GetErrorPropagatorData()->GetMode() == G4ErrorMode_PropBackwards ){
111 sclambda *= -1;
112 scphi += CLHEP::pi;
113 }
114 G4double cosLambda = std::cos( sclambda );
115 G4double sinLambda = std::sin( sclambda );
116 G4double sinPhi = std::sin( scphi );
117 G4double cosPhi = std::cos( scphi );
118
119#ifdef G4EVERBOSE
120 if( iverbose >= 4) G4cout << " PM " << fMomentum.mag() << " pLambda " << sclambda << " pPhi " << scphi << G4endl;
121#endif
122
123 G4ThreeVector vTN( cosLambda*cosPhi, cosLambda*sinPhi,sinLambda );
124 G4ThreeVector vUN( -sinPhi, cosPhi, 0. );
125 G4ThreeVector vVN( -vTN.z()*vUN.y(),vTN.z()*vUN.x(), cosLambda );
126
127#ifdef G4EVERBOSE
128 if( iverbose >= 4) G4cout << " SC2SD: vTN " << vTN << " vUN " << vUN << " vVN " << vVN << G4endl;
129#endif
130 G4double UJ = vUN*GetVectorV();
131 G4double UK = vUN*GetVectorW();
132 G4double VJ = vVN*GetVectorV();
133 G4double VK = vVN*GetVectorW();
134
135
136 //--- Get transformation first
137 G4ErrorMatrix transfM(5, 5, 0 );
138 //--- Get magnetic field
139 const G4Field* field = G4TransportationManager::GetTransportationManager()->GetFieldManager()->GetDetectorField();
140
141 G4Vector3D vectorU = GetVectorV().cross( GetVectorW() );
142
143#ifdef G4EVERBOSE
144 if( iverbose >= 4) G4cout << "vectors " << vectorU << " " << GetVectorV() << " " << GetVectorW() << G4endl;
145#endif
146 G4double T1R = 1. / ( vTN * vectorU );
147
148 if( fCharge != 0 && field ) {
149 G4double pos[3]; pos[0] = fPosition.x()*cm; pos[1] = fPosition.y()*cm; pos[2] = fPosition.z()*cm;
150 G4double Hd[3];
151 field->GetFieldValue( pos, Hd );
152 G4ThreeVector H = G4ThreeVector( Hd[0], Hd[1], Hd[2] ) / tesla *10.; //in kilogauus
153 G4double magH = H.mag();
154 G4double invP = 1./(fMomentum.mag()/GeV);
155 G4double magHM = magH * invP;
156 if( magH != 0. ) {
157 G4double magHM2 = fCharge / magH;
158 G4double Q = -magHM * c_light/ (km/ns);
159#ifdef G4EVERBOSE
160 if( iverbose >= 4) G4cout << GeV << " Q " << Q << " magHM " << magHM << " c_light/(km/ns) " << c_light/(km/ns) << G4endl;
161#endif
162
163 G4double sinz = -H*vUN * magHM2;
164 G4double cosz = H*vVN * magHM2;
165 G4double T3R = Q * std::pow(T1R,3);
166 G4double UI = vUN * vectorU;
167 G4double VI = vVN * vectorU;
168#ifdef G4EVERBOSE
169 if( iverbose >= 4) {
170 G4cout << " T1R " << T1R << " T3R " << T3R << G4endl;
171 G4cout << " UI " << UI << " VI " << VI << " vectorU " << vectorU << G4endl;
172 G4cout << " UJ " << UJ << " VJ " << VJ << G4endl;
173 G4cout << " UK " << UK << " VK " << VK << G4endl;
174 }
175#endif
176
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;
181 }
182 }
183
184 G4double T2R = T1R * T1R;
185 transfM[0][0] = 1.;
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;
194
195#ifdef G4EVERBOSE
196 if( iverbose >= 4) G4cout << " SC2SD transf matrix A= " << transfM << G4endl;
197#endif
198 fError = G4ErrorTrajErr( tpSC.GetError().similarity( transfM ) );
199
200#ifdef G4EVERBOSE
201 if( iverbose >= 1) G4cout << "G4EP: error matrix SC2SD " << fError << G4endl;
202 if( iverbose >= 4) G4cout << "G4ErrorSurfaceTrajState from SC " << *this << G4endl;
203#endif
204
205}
G4ErrorMode_PropBackwards
@ G4ErrorMode_PropBackwards
Definition: G4ErrorPropagatorData.hh:48
G4ErrorTrajErr
G4ErrorSymMatrix G4ErrorTrajErr
Definition: G4ErrorTrajErr.hh:44
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:42
G4double
double G4double
Definition: G4Types.hh:64
G4endl
#define G4endl
Definition: G4ios.hh:52
G4cout
G4DLLIMPORT std::ostream G4cout
CLHEP::Hep3Vector::mag
double mag() const
G4ErrorFreeTrajParam::GetLambda
G4double GetLambda() const
Definition: G4ErrorFreeTrajParam.hh:73
G4ErrorFreeTrajState::GetParameters
G4ErrorFreeTrajParam GetParameters() const
Definition: G4ErrorFreeTrajState.hh:108
G4ErrorPropagatorData::GetErrorPropagatorData
static G4ErrorPropagatorData * GetErrorPropagatorData()
Definition: G4ErrorPropagatorData.cc:53
G4ErrorSymMatrix::similarity
G4ErrorSymMatrix similarity(const G4ErrorMatrix &m1) const
Definition: G4ErrorSymMatrix.cc:589
G4ErrorTrajState::fError
G4ErrorTrajErr fError
Definition: G4ErrorTrajState.hh:128
G4ErrorTrajState::GetError
G4ErrorTrajErr GetError() const
Definition: G4ErrorTrajState.hh:104
G4FieldManager::GetDetectorField
const G4Field * GetDetectorField() const
G4Field::GetFieldValue
virtual void GetFieldValue(const double Point[4], double *fieldArr) const =0
G4TransportationManager::GetTransportationManager
static G4TransportationManager * GetTransportationManager()
Definition: G4TransportationManager.cc:95
G4TransportationManager::GetFieldManager
G4FieldManager * GetFieldManager() const
HepGeom::BasicVector3D::cross
BasicVector3D< T > cross(const BasicVector3D< T > &v) const
Definition: BasicVector3D.h:275
ns
#define ns
Definition: xmlparse.cc:597

Referenced by G4ErrorSurfaceTrajState().

◆ Dump()

void G4ErrorSurfaceTrajState::Dump ( std::ostream &  out = G4cout) const
virtual

Implements G4ErrorTrajState.

Definition at line 217 of file G4ErrorSurfaceTrajState.cc.

218{
219 out << *this;
220}

◆ GetParameters()

G4ErrorSurfaceTrajParam G4ErrorSurfaceTrajState::GetParameters ( ) const
inline

Definition at line 94 of file G4ErrorSurfaceTrajState.hh.

95 { return fTrajParam; }

Referenced by G4ErrorFreeTrajState::G4ErrorFreeTrajState().

◆ GetVectorV()

G4Vector3D G4ErrorSurfaceTrajState::GetVectorV ( ) const
inline

Definition at line 112 of file G4ErrorSurfaceTrajState.hh.

113 { return fTrajParam.GetVectorV(); }

Referenced by BuildErrorMatrix(), G4ErrorFreeTrajState::G4ErrorFreeTrajState(), G4ErrorSurfaceTrajState(), SetMomentum(), and SetPosition().

◆ GetVectorW()

G4Vector3D G4ErrorSurfaceTrajState::GetVectorW ( ) const
inline

Definition at line 115 of file G4ErrorSurfaceTrajState.hh.

116 { return fTrajParam.GetVectorW(); }

Referenced by BuildErrorMatrix(), G4ErrorFreeTrajState::G4ErrorFreeTrajState(), G4ErrorSurfaceTrajState(), SetMomentum(), and SetPosition().

◆ SetMomentum()

virtual void G4ErrorSurfaceTrajState::SetMomentum ( const G4Vector3D &  mom)
inlinevirtual

Reimplemented from G4ErrorTrajState.

Definition at line 121 of file G4ErrorSurfaceTrajState.hh.

122 { SetParameters( fPosition, mom, GetVectorV(), GetVectorW() ); }
G4ErrorSurfaceTrajState::SetParameters
void SetParameters(const G4Point3D &pos, const G4Vector3D &mom, const G4Vector3D &vecV, const G4Vector3D &vecW)
Definition: G4ErrorSurfaceTrajState.hh:96

◆ SetParameters() [1/2]

void G4ErrorSurfaceTrajState::SetParameters ( const G4Point3D &  pos,
const G4Vector3D &  mom,
const G4Plane3D &  plane 
)
inline

Definition at line 104 of file G4ErrorSurfaceTrajState.hh.

106 {
107 fPosition = pos;
108 fMomentum = mom;
109 fTrajParam.SetParameters( pos, mom, plane );
110 }
G4ErrorSurfaceTrajParam::SetParameters
void SetParameters(const G4Point3D &pos, const G4Vector3D &mom, const G4Vector3D &vecV, const G4Vector3D &vecW)
Definition: G4ErrorSurfaceTrajParam.cc:83

◆ SetParameters() [2/2]

void G4ErrorSurfaceTrajState::SetParameters ( const G4Point3D &  pos,
const G4Vector3D &  mom,
const G4Vector3D &  vecV,
const G4Vector3D &  vecW 
)
inline

Definition at line 96 of file G4ErrorSurfaceTrajState.hh.

98 {
99 fPosition = pos;
100 fMomentum = mom;
101 fTrajParam.SetParameters( pos, mom, vecV, vecW );
102 }

Referenced by SetMomentum(), and SetPosition().

◆ SetPosition()

virtual void G4ErrorSurfaceTrajState::SetPosition ( const G4Point3D  pos)
inlinevirtual

Reimplemented from G4ErrorTrajState.

Definition at line 118 of file G4ErrorSurfaceTrajState.hh.

119 { SetParameters( pos, fMomentum, GetVectorV(), GetVectorW() ); }

Friends And Related Function Documentation

◆ operator<<

std::ostream & operator<< ( std::ostream &  out,
const G4ErrorSurfaceTrajState &  ts 
)
friend

Definition at line 224 of file G4ErrorSurfaceTrajState.cc.

225{
226 std::ios::fmtflags oldFlags = out.flags();
227 out.setf(std::ios::fixed,std::ios::floatfield);
228
229 ts.DumpPosMomError( out );
230
231 out << " G4ErrorSurfaceTrajState: Params: " << ts.fTrajParam << G4endl;
232 out.flags(oldFlags);
233 return out;
234}
G4ErrorTrajState::DumpPosMomError
void DumpPosMomError(std::ostream &out=G4cout) const
Definition: G4ErrorTrajState.cc:104
The documentation for this class was generated from the following files: