G4FieldTrack Class Reference

#include <G4FieldTrack.hh>

Collaboration diagram for G4FieldTrack:

Public Types
enum	{ ncompSVEC = 12 }

Public Member Functions
	G4FieldTrack (const G4ThreeVector &pPosition, G4double LaboratoryTimeOfFlight, const G4ThreeVector &pMomentumDirection, G4double kineticEnergy, G4double restMass_c2, G4double charge, const G4ThreeVector &polarization, G4double magnetic_dipole_moment=0.0, G4double curve_length=0.0, G4double PDGspin=-1.0)
	G4FieldTrack (const G4FieldTrack &pFieldTrack)
	G4FieldTrack (char)
	~G4FieldTrack ()
void	UpdateState (const G4ThreeVector &pPosition, G4double LaboratoryTimeOfFlight, const G4ThreeVector &pMomentumDirection, G4double kineticEnergy)
void	UpdateFourMomentum (G4double kineticEnergy, const G4ThreeVector &momentumDirection)
void	SetChargeAndMoments (G4double charge, G4double magnetic_dipole_moment=DBL_MAX, G4double electric_dipole_moment=DBL_MAX, G4double magnetic_charge=DBL_MAX)
void	SetPDGSpin (G4double pdgSpin)
G4double	GetPDGSpin ()
	G4FieldTrack (const G4ThreeVector &pPosition, const G4ThreeVector &pMomentumDirection, G4double curve_length, G4double kineticEnergy, const G4double restMass_c2, G4double velocity, G4double LaboratoryTimeOfFlight=0.0, G4double ProperTimeOfFlight=0.0, const G4ThreeVector *pPolarization=0, G4double PDGspin=-1.0)
G4FieldTrack &	operator= (const G4FieldTrack &rStVec)
G4ThreeVector	GetMomentum () const
G4ThreeVector	GetPosition () const
const G4ThreeVector &	GetMomentumDir () const
G4ThreeVector	GetMomentumDirection () const
G4double	GetCurveLength () const
G4ThreeVector	GetPolarization () const
void	SetPolarization (const G4ThreeVector &vecPol)
G4double	GetLabTimeOfFlight () const
G4double	GetProperTimeOfFlight () const
G4double	GetKineticEnergy () const
G4double	GetCharge () const
G4double	GetRestMass () const
void	SetPosition (G4ThreeVector nPos)
void	SetMomentum (G4ThreeVector nMomDir)
void	SetMomentumDir (G4ThreeVector nMomDir)
void	SetRestMass (G4double Mass_c2)
void	SetCurveLength (G4double nCurve_s)
void	SetKineticEnergy (G4double nEnergy)
void	SetLabTimeOfFlight (G4double tofLab)
void	SetProperTimeOfFlight (G4double tofProper)
void	DumpToArray (G4double valArr[ncompSVEC]) const
void	LoadFromArray (const G4double valArr[ncompSVEC], G4int noVarsIntegrated)
void	InitialiseSpin (const G4ThreeVector &vecPolarization)
G4ThreeVector	GetSpin () const
void	SetSpin (G4ThreeVector vSpin)
const G4ChargeState *	GetChargeState () const

Private Member Functions
G4FieldTrack &	SetCurvePnt (const G4ThreeVector &pPosition, const G4ThreeVector &pMomentum, G4double s_curve)

Private Attributes
G4double	SixVector [6]
G4double	fDistanceAlongCurve
G4double	fKineticEnergy
G4double	fRestMass_c2
G4double	fLabTimeOfFlight
G4double	fProperTimeOfFlight
G4ThreeVector	fPolarization
G4ThreeVector	fMomentumDir
G4ChargeState	fChargeState

Friends
std::ostream &	operator<< (std::ostream &os, const G4FieldTrack &SixVec)

Detailed Description

Definition at line 54 of file G4FieldTrack.hh.

Member Enumeration Documentation

anonymous enum

anonymous enum

Enumerator
ncompSVEC

Definition at line 149 of file G4FieldTrack.hh.

{ ncompSVEC = 12 };

Constructor & Destructor Documentation

G4FieldTrack() [1/4]

G4FieldTrack::G4FieldTrack	(	const G4ThreeVector &	pPosition,
		G4double	LaboratoryTimeOfFlight,
		const G4ThreeVector &	pMomentumDirection,
		G4double	kineticEnergy,
		G4double	restMass_c2,
		G4double	charge,
		const G4ThreeVector &	polarization,
		G4double	magnetic_dipole_moment = 0.0,
		G4double	curve_length = 0.0,
		G4double	PDGspin = -1.0
	)

Definition at line 72 of file G4FieldTrack.cc.

:  fDistanceAlongCurve(curve_length),
   fKineticEnergy(kineticEnergy),
   fRestMass_c2(restMass_c2),
   fLabTimeOfFlight(LaboratoryTimeOfFlight), 
   fProperTimeOfFlight(0.),
   // fMomentumDir(pMomentumDirection),
   fChargeState(  charge, magnetic_dipole_moment, pdgSpin ) 
   // fChargeState(  charge, magnetic_dipole_moment ) , 
   // fPDGSpin( pdgSpin )
{
  UpdateFourMomentum( kineticEnergy, pMomentumDirection ); 
      // Sets momentum direction as well.

  SetPosition( pPosition ); 

  SetPolarization( vecPolarization ); 
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4FieldTrack() [2/4]

G4FieldTrack::G4FieldTrack

(

const G4FieldTrack &

pFieldTrack

)

G4FieldTrack() [3/4]

G4FieldTrack::G4FieldTrack

(

char

)

Definition at line 128 of file G4FieldTrack.cc.

  : fKineticEnergy(0.), fRestMass_c2(0.), fLabTimeOfFlight(0.),
    fProperTimeOfFlight(0.), fChargeState( DBL_MAX , DBL_MAX, -1 )
{
  G4ThreeVector Zero(0.0, 0.0, 0.0);
  SetCurvePnt( Zero, Zero, 0.0 );
  SetPolarization( Zero ); 
  // fInitialMomentumMag= 0.00; // Invalid
  // fLastMomentumMag= 0.0; 
}

Here is the call graph for this function:

~G4FieldTrack()

G4FieldTrack::~G4FieldTrack

(

)

G4FieldTrack() [4/4]

G4FieldTrack::G4FieldTrack	(	const G4ThreeVector &	pPosition,
		const G4ThreeVector &	pMomentumDirection,
		G4double	curve_length,
		G4double	kineticEnergy,
		const G4double	restMass_c2,
		G4double	velocity,
		G4double	LaboratoryTimeOfFlight = 0.0,
		G4double	ProperTimeOfFlight = 0.0,
		const G4ThreeVector *	pPolarization = 0,
		G4double	PDGspin = -1.0
	)

Definition at line 100 of file G4FieldTrack.cc.

 : fDistanceAlongCurve(curve_length),
   fKineticEnergy(kineticEnergy),
   fRestMass_c2(restMass_c2),
   fLabTimeOfFlight(pLaboratoryTimeOfFlight), 
   fProperTimeOfFlight(pProperTimeOfFlight),
   fChargeState( DBL_MAX, DBL_MAX, -1.0 ) //  charge not set 
{
  UpdateFourMomentum( kineticEnergy, pMomentumDirection ); 
      // Sets momentum direction as well.
    
  SetPosition( pPosition );    
  fChargeState.SetPDGSpin( pdgSpin );   

  G4ThreeVector PolarVec(0.0, 0.0, 0.0); 
  if( pPolarization )  { PolarVec= *pPolarization; }
  SetPolarization( PolarVec );
}

Here is the call graph for this function:

Member Function Documentation

DumpToArray()

void G4FieldTrack::DumpToArray ( G4double  valArr[ncompSVEC] ) const

inline

Here is the caller graph for this function:

GetCharge()

G4double G4FieldTrack::GetCharge ( ) const

inline

Here is the caller graph for this function:

GetChargeState()

const G4ChargeState* G4FieldTrack::GetChargeState ( ) const

inline

Definition at line 185 of file G4FieldTrack.hh.

{ return &fChargeState; } 

Here is the caller graph for this function:

GetCurveLength()

G4double G4FieldTrack::GetCurveLength ( ) const

inline

Here is the caller graph for this function:

GetKineticEnergy()

G4double G4FieldTrack::GetKineticEnergy ( ) const

inline

Here is the caller graph for this function:

GetLabTimeOfFlight()

G4double G4FieldTrack::GetLabTimeOfFlight ( ) const

inline

Here is the caller graph for this function:

GetMomentum()

G4ThreeVector G4FieldTrack::GetMomentum ( ) const

inline

Here is the caller graph for this function:

GetMomentumDir()

const G4ThreeVector& G4FieldTrack::GetMomentumDir ( ) const

inline

Here is the caller graph for this function:

GetMomentumDirection()

G4ThreeVector G4FieldTrack::GetMomentumDirection ( ) const

inline

Here is the caller graph for this function:

GetPDGSpin()

G4double G4FieldTrack::GetPDGSpin ( )

inline

Definition at line 94 of file G4FieldTrack.hh.

{ return fChargeState.GetPDGSpin(); } 

Here is the call graph for this function:

GetPolarization()

G4ThreeVector G4FieldTrack::GetPolarization ( ) const

inline

Here is the caller graph for this function:

GetPosition()

G4ThreeVector G4FieldTrack::GetPosition ( ) const

inline

Here is the caller graph for this function:

GetProperTimeOfFlight()

G4double G4FieldTrack::GetProperTimeOfFlight ( ) const

inline

Here is the caller graph for this function:

GetRestMass()

G4double G4FieldTrack::GetRestMass ( ) const

inline

Definition at line 126 of file G4FieldTrack.hh.

{ return fRestMass_c2; }

Here is the call graph for this function:

Here is the caller graph for this function:

GetSpin()

G4ThreeVector G4FieldTrack::GetSpin ( ) const

inline

Definition at line 161 of file G4FieldTrack.hh.

{ return GetPolarization(); } 

Here is the call graph for this function:

Here is the caller graph for this function:

InitialiseSpin()

void G4FieldTrack::InitialiseSpin ( const G4ThreeVector &  vecPolarization )

inline

Definition at line 159 of file G4FieldTrack.hh.

           { SetPolarization( vecPolarization );  } 

Here is the call graph for this function:

LoadFromArray()

void G4FieldTrack::LoadFromArray	(	const G4double	valArr[ncompSVEC],
		G4int	noVarsIntegrated
	)

Definition at line 166 of file G4FieldTrack.cc.

{
  G4int i;

  // Fill the variables not integrated with zero -- so it's clear !!
  G4double valArr[ncompSVEC];
  for( i=0; i<noVarsIntegrated; i++){
     valArr[i]= valArrIn[i];
  }
  for( i=noVarsIntegrated; i<ncompSVEC; i++) {
     valArr[i]= 0.0; 
  }

  SixVector[0]=valArr[0];
  SixVector[1]=valArr[1];
  SixVector[2]=valArr[2];
  SixVector[3]=valArr[3];
  SixVector[4]=valArr[4];
  SixVector[5]=valArr[5];

  G4ThreeVector Momentum(valArr[3],valArr[4],valArr[5]);

  G4double momentum_square= Momentum.mag2();
  fMomentumDir= Momentum.unit();

  fKineticEnergy = momentum_square / 
                   (std::sqrt(momentum_square+fRestMass_c2*fRestMass_c2)
                     + fRestMass_c2 ); 
  // The above equation is stable for small and large momenta

  // The following components may or may not be
  //    integrated over -- integration is optional
  // fKineticEnergy= valArr[6];

  fLabTimeOfFlight=valArr[7];
  fProperTimeOfFlight=valArr[8];
  G4ThreeVector  vecPolarization= G4ThreeVector(valArr[9],valArr[10],valArr[11]);
  SetPolarization( vecPolarization ); 

  // fMomentumDir=G4ThreeVector(valArr[13],valArr[14],valArr[15]);
  // fDistanceAlongCurve= valArr[]; 
}

Here is the call graph for this function:

Here is the caller graph for this function:

operator=()

G4FieldTrack& G4FieldTrack::operator= ( const G4FieldTrack &  rStVec )

inline

Here is the caller graph for this function:

SetChargeAndMoments()

void G4FieldTrack::SetChargeAndMoments	(	G4double	charge,
		G4double	magnetic_dipole_moment = DBL_MAX,
		G4double	electric_dipole_moment = DBL_MAX,
		G4double	magnetic_charge = DBL_MAX
	)

Definition at line 140 of file G4FieldTrack.cc.

{
  fChargeState.SetChargesAndMoments( charge,  
                                     magnetic_dipole_moment, 
                                     electric_dipole_moment,  
                                     magnetic_charge ); 

  // NOTE: Leaves Spin unchanged !
  // 
  // G4double pdgSpin= fChargeState.GetSpin();   // New Property of ChargeState (not well documented! )

  // IDEA: Improve the implementation using handles
  //   -- and handle to the old one (which can be shared by other copies) and
  //      must not be left to hang loose 
  // 
  // fpChargeState= new G4ChargeState(  charge, magnetic_dipole_moment, 
  //                 electric_dipole_moment, magnetic_charge  ); 
}

Here is the call graph for this function:

Here is the caller graph for this function:

SetCurveLength()

void G4FieldTrack::SetCurveLength ( G4double  nCurve_s )

inline

Here is the caller graph for this function:

SetCurvePnt()

G4FieldTrack& G4FieldTrack::SetCurvePnt ( const G4ThreeVector &  pPosition, const G4ThreeVector &  pMomentum, G4double  s_curve  )

inlineprivate

Here is the caller graph for this function:

SetKineticEnergy()

void G4FieldTrack::SetKineticEnergy ( G4double  nEnergy )

inline

Here is the caller graph for this function:

SetLabTimeOfFlight()

void G4FieldTrack::SetLabTimeOfFlight ( G4double  tofLab )

inline

Here is the caller graph for this function:

SetMomentum()

void G4FieldTrack::SetMomentum ( G4ThreeVector  nMomDir )

inline

Here is the caller graph for this function:

SetMomentumDir()

void G4FieldTrack::SetMomentumDir ( G4ThreeVector  nMomDir )

inline

Here is the caller graph for this function:

SetPDGSpin()

void G4FieldTrack::SetPDGSpin ( G4double  pdgSpin )

inline

Definition at line 93 of file G4FieldTrack.hh.

{ fChargeState.SetPDGSpin(pdgSpin); } 

Here is the call graph for this function:

SetPolarization()

void G4FieldTrack::SetPolarization ( const G4ThreeVector &  vecPol )

inline

Here is the caller graph for this function:

SetPosition()

void G4FieldTrack::SetPosition ( G4ThreeVector  nPos )

inline

Here is the caller graph for this function:

SetProperTimeOfFlight()

void G4FieldTrack::SetProperTimeOfFlight ( G4double  tofProper )

inline

Here is the caller graph for this function:

SetRestMass()

void G4FieldTrack::SetRestMass ( G4double  Mass_c2 )

inline

Definition at line 136 of file G4FieldTrack.hh.

{ fRestMass_c2= Mass_c2; }

Here is the call graph for this function:

SetSpin()

void G4FieldTrack::SetSpin ( G4ThreeVector  vSpin )

inline

Definition at line 162 of file G4FieldTrack.hh.

{ SetPolarization(vSpin); }

Here is the call graph for this function:

UpdateFourMomentum()

void G4FieldTrack::UpdateFourMomentum ( G4double  kineticEnergy, const G4ThreeVector &  momentumDirection  )

inline

Here is the caller graph for this function:

UpdateState()

void G4FieldTrack::UpdateState ( const G4ThreeVector &  pPosition, G4double  LaboratoryTimeOfFlight, const G4ThreeVector &  pMomentumDirection, G4double  kineticEnergy  )

inline

Friends And Related Function Documentation

operator<<

std::ostream& operator<< ( std::ostream &  os, const G4FieldTrack &  SixVec  )

friend

Definition at line 33 of file G4FieldTrack.cc.

{
     const G4double *SixV = SixVec.SixVector;
     const int precPos= 9;   // For position
     const int precEp=  9;   // For Energy / momentum
     const int precLen= 12;  // For Length along track
     const int precSpin= 9;  // For polarisation
     const int precTime= 6;  // For time of flight
     const int oldpr= os.precision(precPos);
     os << " ( ";
     os << " X= " << SixV[0] << " " << SixV[1] << " "
                  << SixV[2] << " ";  // Position
     os.precision(precEp);     
     os << " P= " << SixV[3] << " " << SixV[4] << " "
                  << SixV[5] << " ";  // Momentum
     os << " Pmag= "
        << G4ThreeVector(SixV[3], SixV[4], SixV[5]).mag(); // mom magnitude
     os << " Ekin= " << SixVec.fKineticEnergy ;
     os.precision(precLen);
     os << " l= " << SixVec.GetCurveLength();
     os.precision(6);
     os << " m0= " <<   SixVec.fRestMass_c2;
     os << " (Pdir-1)= " <<  SixVec.fMomentumDir.mag()-1.0;
     if( SixVec.fLabTimeOfFlight > 0.0 )   os.precision(precTime);
     else                                  os.precision(3); 
     os << " t_lab= "    << SixVec.fLabTimeOfFlight; 
     os << " t_proper= " << SixVec.fProperTimeOfFlight ;
     G4ThreeVector pol= SixVec.GetPolarization();
     if( pol.mag2() > 0.0 ){
        os.precision(precSpin);
        os << " PolV= " << pol; // SixVec.GetPolarization();
     }else{
        os << " PolV= (0,0,0) "; 
     }
     os << " ) ";
     os.precision(oldpr);
     return os;
}

Member Data Documentation

fChargeState

G4ChargeState G4FieldTrack::fChargeState

private

Definition at line 181 of file G4FieldTrack.hh.

fDistanceAlongCurve

G4double G4FieldTrack::fDistanceAlongCurve

private

Definition at line 171 of file G4FieldTrack.hh.

fKineticEnergy

G4double G4FieldTrack::fKineticEnergy

private

Definition at line 172 of file G4FieldTrack.hh.

fLabTimeOfFlight

G4double G4FieldTrack::fLabTimeOfFlight

private

Definition at line 174 of file G4FieldTrack.hh.

fMomentumDir

G4ThreeVector G4FieldTrack::fMomentumDir

private

Definition at line 177 of file G4FieldTrack.hh.

fPolarization

G4ThreeVector G4FieldTrack::fPolarization

private

Definition at line 176 of file G4FieldTrack.hh.

fProperTimeOfFlight

G4double G4FieldTrack::fProperTimeOfFlight

private

Definition at line 175 of file G4FieldTrack.hh.

fRestMass_c2

G4double G4FieldTrack::fRestMass_c2

private

Definition at line 173 of file G4FieldTrack.hh.

SixVector

G4double G4FieldTrack::SixVector[6]

private

Definition at line 170 of file G4FieldTrack.hh.

---

The documentation for this class was generated from the following files:

• G4FieldTrack.hh
• G4FieldTrack.cc