G4FieldTrack Class Reference

#include <G4FieldTrack.hh>

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 (char)
	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=nullptr, G4double PDGspin=-1.0)
	~G4FieldTrack ()
	G4FieldTrack (const G4FieldTrack &pFieldTrack)
G4FieldTrack &	operator= (const G4FieldTrack &rStVec)
	G4FieldTrack (G4FieldTrack &&from) noexcept
G4FieldTrack &	operator= (G4FieldTrack &&from) noexcept
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 ()
G4ThreeVector	GetMomentum () const
G4ThreeVector	GetPosition () const
const G4ThreeVector &	GetMomentumDir () const
G4ThreeVector	GetMomentumDirection () const
G4double	GetCurveLength () const
G4ThreeVector	GetPolarization () const
void	SetPolarization (const G4ThreeVector &vecPol)
const G4ChargeState *	GetChargeState () const
G4double	GetLabTimeOfFlight () const
G4double	GetProperTimeOfFlight () const
G4double	GetKineticEnergy () const
G4double	GetCharge () const
G4double	GetRestMass () const
void	SetPosition (const G4ThreeVector &nPos)
void	SetMomentum (const G4ThreeVector &nMomDir)
void	SetMomentumDir (const 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 (const G4ThreeVector &vSpin)

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

Detailed Description

Definition at line 44 of file G4FieldTrack.hh.

Member Enumeration Documentation

anonymous enum

anonymous enum

Enumerator
ncompSVEC

Definition at line 144 of file G4FieldTrack.hh.

{ ncompSVEC = 12 };

Constructor & Destructor Documentation

G4FieldTrack() [1/5]

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 81 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 ); 
}

G4FieldTrack() [2/5]

G4FieldTrack::G4FieldTrack

(

char

)

Definition at line 136 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; 
}

G4FieldTrack() [3/5]

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 = nullptr,
		G4double	PDGspin = -1.0 )

Definition at line 108 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 != nullptr )  { PolarVec= *pPolarization; }
  SetPolarization( PolarVec );
}

~G4FieldTrack()

G4FieldTrack::~G4FieldTrack

(

)

G4FieldTrack() [4/5]

G4FieldTrack::G4FieldTrack ( const G4FieldTrack & pFieldTrack )

inline

G4FieldTrack() [5/5]

G4FieldTrack::G4FieldTrack ( G4FieldTrack && from )

inlinenoexcept

Member Function Documentation

DumpToArray()

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

inline

Referenced by G4BorisDriver::AccurateAdvance(), G4MagInt_Driver::AccurateAdvance(), G4OldMagIntDriver::AccurateAdvance(), G4BorisDriver::GetDerivatives(), G4MagInt_Driver::GetDerivatives(), G4OldMagIntDriver::GetDerivatives(), G4BorisDriver::QuickAdvance(), G4MagInt_Driver::QuickAdvance(), and G4OldMagIntDriver::QuickAdvance().

GetCharge()

G4double G4FieldTrack::GetCharge ( ) const

inline

Referenced by G4BorisDriver::AccurateAdvance(), G4PathFinder::ComputeStep(), G4BorisDriver::QuickAdvance(), and G4QSStepper< QSS >::reset().

GetChargeState()

const G4ChargeState * G4FieldTrack::GetChargeState ( ) const

inline

Referenced by G4PathFinder::DoNextCurvedStep().

GetCurveLength()

G4double G4FieldTrack::GetCurveLength ( ) const

inline

Referenced by G4BorisDriver::AccurateAdvance(), G4MagInt_Driver::AccurateAdvance(), G4OldMagIntDriver::AccurateAdvance(), G4ChordFinder::ApproxCurvePointS(), G4ChordFinder::ApproxCurvePointV(), G4VIntersectionLocator::CheckAndReEstimateEndpoint(), G4PropagatorInField::ComputeStep(), G4BrentLocator::EstimateIntersectionPoint(), G4MultiLevelLocator::EstimateIntersectionPoint(), G4SimpleLocator::EstimateIntersectionPoint(), G4LocatorChangeRecord::GetLength(), G4DriverReporter::PrintStat_Aux(), G4MagInt_Driver::PrintStat_Aux(), G4OldMagIntDriver::PrintStat_Aux(), G4DriverReporter::PrintStatus(), G4MagInt_Driver::PrintStatus(), G4OldMagIntDriver::PrintStatus(), G4PropagatorInField::printStatus(), G4VIntersectionLocator::printStatus(), G4BorisDriver::QuickAdvance(), G4MagInt_Driver::QuickAdvance(), G4OldMagIntDriver::QuickAdvance(), G4VIntersectionLocator::ReEstimateEndpoint(), and G4VIntersectionLocator::ReportReversedPoints().

GetKineticEnergy()

G4double G4FieldTrack::GetKineticEnergy ( ) const

inline

Referenced by G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), G4Transportation::AlongStepGetPhysicalInteractionLength(), G4DriverReporter::PrintStat_Aux(), G4MagInt_Driver::PrintStat_Aux(), and G4OldMagIntDriver::PrintStat_Aux().

GetLabTimeOfFlight()

G4double G4FieldTrack::GetLabTimeOfFlight ( ) const

inline

Referenced by G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), and G4Transportation::AlongStepGetPhysicalInteractionLength().

GetMomentum()

G4ThreeVector G4FieldTrack::GetMomentum ( ) const

inline

Referenced by G4PropagatorInField::ComputeStep(), G4PathFinder::DoNextCurvedStep(), G4PropagatorInField::printStatus(), G4VIntersectionLocator::printStatus(), G4QSStepper< QSS >::reset(), and G4FieldTrackUpdator::Update().

GetMomentumDir()

const G4ThreeVector & G4FieldTrack::GetMomentumDir ( ) const

inline

Referenced by G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), G4Transportation::AlongStepGetPhysicalInteractionLength(), G4PropagatorInField::ComputeStep(), G4BrentLocator::EstimateIntersectionPoint(), G4MultiLevelLocator::EstimateIntersectionPoint(), G4SimpleLocator::EstimateIntersectionPoint(), G4DriverReporter::PrintStat_Aux(), G4MagInt_Driver::PrintStat_Aux(), G4OldMagIntDriver::PrintStat_Aux(), G4DriverReporter::PrintStatus(), G4MagInt_Driver::PrintStatus(), G4OldMagIntDriver::PrintStatus(), G4PropagatorInField::printStatus(), and G4VIntersectionLocator::printStatus().

GetMomentumDirection()

G4ThreeVector G4FieldTrack::GetMomentumDirection ( ) const

inline

Referenced by G4ITPathFinder::ComputeStep(), G4PathFinder::ComputeStep(), G4ITPathFinder::DoNextLinearStep(), G4PathFinder::DoNextLinearStep(), G4BrentLocator::EstimateIntersectionPoint(), G4MultiLevelLocator::EstimateIntersectionPoint(), and G4SimpleLocator::EstimateIntersectionPoint().

GetPDGSpin()

G4double G4FieldTrack::GetPDGSpin ( )

inline

GetPolarization()

G4ThreeVector G4FieldTrack::GetPolarization ( ) const

inline

GetPosition()

G4ThreeVector G4FieldTrack::GetPosition ( ) const

inline

Referenced by G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), G4FastSimulationManagerProcess::AlongStepGetPhysicalInteractionLength(), G4ImportanceProcess::AlongStepGetPhysicalInteractionLength(), G4ParallelGeometriesLimiterProcess::AlongStepGetPhysicalInteractionLength(), G4ParallelWorldProcess::AlongStepGetPhysicalInteractionLength(), G4ParallelWorldScoringProcess::AlongStepGetPhysicalInteractionLength(), G4Transportation::AlongStepGetPhysicalInteractionLength(), G4WeightCutOffProcess::AlongStepGetPhysicalInteractionLength(), G4WeightWindowProcess::AlongStepGetPhysicalInteractionLength(), G4ChordFinder::ApproxCurvePointS(), G4ChordFinder::ApproxCurvePointV(), G4VIntersectionLocator::CheckAndReEstimateEndpoint(), G4ITPathFinder::ComputeStep(), G4PathFinder::ComputeStep(), G4PropagatorInField::ComputeStep(), G4ITPathFinder::DoNextCurvedStep(), G4PathFinder::DoNextCurvedStep(), G4ITPathFinder::DoNextLinearStep(), G4PathFinder::DoNextLinearStep(), G4BrentLocator::EstimateIntersectionPoint(), G4MultiLevelLocator::EstimateIntersectionPoint(), G4SimpleLocator::EstimateIntersectionPoint(), G4PathFinder::Locate(), G4DriverReporter::PrintStat_Aux(), G4MagInt_Driver::PrintStat_Aux(), G4OldMagIntDriver::PrintStat_Aux(), G4DriverReporter::PrintStatus(), G4MagInt_Driver::PrintStatus(), G4OldMagIntDriver::PrintStatus(), G4PropagatorInField::printStatus(), G4VIntersectionLocator::printStatus(), G4MagInt_Driver::QuickAdvance(), G4OldMagIntDriver::QuickAdvance(), G4VIntersectionLocator::ReEstimateEndpoint(), G4PathFinder::ReLocate(), and G4QSStepper< QSS >::reset().

GetProperTimeOfFlight()

G4double G4FieldTrack::GetProperTimeOfFlight ( ) const

inline

GetRestMass()

G4double G4FieldTrack::GetRestMass ( ) const

inline

Referenced by G4BorisDriver::AccurateAdvance(), G4PathFinder::DoNextCurvedStep(), G4MultiLevelLocator::EstimateIntersectionPoint(), G4BorisDriver::QuickAdvance(), and G4QSStepper< QSS >::reset().

GetSpin()

G4ThreeVector G4FieldTrack::GetSpin ( ) const

inline

Referenced by G4CoupledTransportation::AlongStepGetPhysicalInteractionLength(), and G4Transportation::AlongStepGetPhysicalInteractionLength().

InitialiseSpin()

void G4FieldTrack::InitialiseSpin ( const G4ThreeVector & vecPolarization )

inline

LoadFromArray()

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

Definition at line 175 of file G4FieldTrack.cc.

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

Referenced by G4BorisDriver::AccurateAdvance(), G4MagInt_Driver::AccurateAdvance(), G4OldMagIntDriver::AccurateAdvance(), G4DriverReporter::PrintStatus(), G4MagInt_Driver::PrintStatus(), G4OldMagIntDriver::PrintStatus(), G4BorisDriver::QuickAdvance(), G4MagInt_Driver::QuickAdvance(), and G4OldMagIntDriver::QuickAdvance().

operator=() [1/2]

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

inline

operator=() [2/2]

G4FieldTrack & G4FieldTrack::operator= ( G4FieldTrack && from )

inlinenoexcept

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 147 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  ); 
}

Referenced by G4FieldTrackUpdator::Update().

SetCurveLength()

void G4FieldTrack::SetCurveLength ( G4double nCurve_s )

inline

Referenced by G4BorisDriver::AccurateAdvance(), G4MagInt_Driver::AccurateAdvance(), G4OldMagIntDriver::AccurateAdvance(), G4DriverReporter::PrintStatus(), G4MagInt_Driver::PrintStatus(), G4OldMagIntDriver::PrintStatus(), G4BorisDriver::QuickAdvance(), G4MagInt_Driver::QuickAdvance(), and G4OldMagIntDriver::QuickAdvance().

SetKineticEnergy()

void G4FieldTrack::SetKineticEnergy ( G4double nEnergy )

inline

SetLabTimeOfFlight()

void G4FieldTrack::SetLabTimeOfFlight ( G4double tofLab )

inline

SetMomentum()

void G4FieldTrack::SetMomentum ( const G4ThreeVector & nMomDir )

inline

Referenced by G4MultiLevelLocator::EstimateIntersectionPoint().

SetMomentumDir()

void G4FieldTrack::SetMomentumDir ( const G4ThreeVector & nMomDir )

inline

SetPDGSpin()

void G4FieldTrack::SetPDGSpin ( G4double pdgSpin )

inline

Referenced by G4FieldTrackUpdator::Update().

SetPolarization()

void G4FieldTrack::SetPolarization ( const G4ThreeVector & vecPol )

inline

Referenced by G4FieldTrack(), G4FieldTrack(), G4FieldTrack(), and LoadFromArray().

SetPosition()

void G4FieldTrack::SetPosition ( const G4ThreeVector & nPos )

inline

Referenced by G4PathFinder::DoNextLinearStep(), G4BrentLocator::EstimateIntersectionPoint(), G4MultiLevelLocator::EstimateIntersectionPoint(), G4SimpleLocator::EstimateIntersectionPoint(), G4FieldTrack(), and G4FieldTrack().

SetProperTimeOfFlight()

void G4FieldTrack::SetProperTimeOfFlight ( G4double tofProper )

inline

Referenced by G4PathFinder::DoNextLinearStep(), and G4FieldTrackUpdator::Update().

SetRestMass()

void G4FieldTrack::SetRestMass ( G4double Mass_c2 )

inline

Referenced by G4FieldTrackUpdator::Update().

SetSpin()

void G4FieldTrack::SetSpin ( const G4ThreeVector & vSpin )

inline

Referenced by G4FieldTrackUpdator::Update().

UpdateFourMomentum()

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

inline

Referenced by G4FieldTrack(), and G4FieldTrack().

UpdateState()

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

inline

Referenced by G4FieldTrackUpdator::Update().

Friends And Related Symbol 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 G4int precPos= 9;   // For position
     const G4int precEp=  9;   // For Energy / momentum
     const G4int precLen= 12;  // For Length along track
     const G4int precSpin= 9;  // For polarisation
     const G4int precTime= 6;  // For time of flight
     const G4long 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;
}

---

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

• G4FieldTrack.hh
• G4FieldTrack.cc