#include <G4DecayWithSpin.hh>

Inheritance diagram for G4DecayWithSpin:

Collaboration diagram for G4DecayWithSpin:

Public Member Functions
	G4DecayWithSpin (const G4String &processName="DecayWithSpin")

virtual	~G4DecayWithSpin ()

Public Member Functions inherited from G4Decay
	G4Decay (const G4String &processName="Decay")

virtual	~G4Decay ()

virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

virtual G4bool	IsApplicable (const G4ParticleDefinition &)

virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &track, G4ForceCondition *condition)

virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)

virtual void	StartTracking (G4Track *)

virtual void	EndTracking ()

void	SetExtDecayer (G4VExtDecayer *)

const G4VExtDecayer *	GetExtDecayer () const

G4double	GetRemainderLifeTime () const

void	SetVerboseLevel (G4int value)

G4int	GetVerboseLevel () const

Public Member Functions inherited from G4VRestDiscreteProcess
	G4VRestDiscreteProcess (const G4String &, G4ProcessType aType=fNotDefined)

	G4VRestDiscreteProcess (G4VRestDiscreteProcess &)

virtual	~G4VRestDiscreteProcess ()

virtual G4double	AlongStepGetPhysicalInteractionLength (const G4Track &, G4double, G4double, G4double &, G4GPILSelection *)

virtual G4VParticleChange *	AlongStepDoIt (const G4Track &, const G4Step &)

Public Member Functions inherited from G4VProcess
	G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)

	G4VProcess (const G4VProcess &right)

virtual	~G4VProcess ()

G4int	operator== (const G4VProcess &right) const

G4int	operator!= (const G4VProcess &right) const

G4double	GetCurrentInteractionLength () const

void	SetPILfactor (G4double value)

G4double	GetPILfactor () const

G4double	AlongStepGPIL (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)

G4double	AtRestGPIL (const G4Track &track, G4ForceCondition *condition)

G4double	PostStepGPIL (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)

virtual void	PreparePhysicsTable (const G4ParticleDefinition &)

virtual G4bool	StorePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)

virtual G4bool	RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)

const G4String &	GetPhysicsTableFileName (const G4ParticleDefinition *, const G4String &directory, const G4String &tableName, G4bool ascii=false)

const G4String &	GetProcessName () const

G4ProcessType	GetProcessType () const

void	SetProcessType (G4ProcessType)

G4int	GetProcessSubType () const

void	SetProcessSubType (G4int)

virtual void	SetProcessManager (const G4ProcessManager *)

virtual const G4ProcessManager *	GetProcessManager ()

virtual void	ResetNumberOfInteractionLengthLeft ()

G4double	GetNumberOfInteractionLengthLeft () const

G4double	GetTotalNumberOfInteractionLengthTraversed () const

G4bool	isAtRestDoItIsEnabled () const

G4bool	isAlongStepDoItIsEnabled () const

G4bool	isPostStepDoItIsEnabled () const

virtual void	DumpInfo () const

void	SetVerboseLevel (G4int value)

G4int	GetVerboseLevel () const

virtual void	SetMasterProcess (G4VProcess *masterP)

const G4VProcess *	GetMasterProcess () const

virtual void	BuildWorkerPhysicsTable (const G4ParticleDefinition &part)

virtual void	PrepareWorkerPhysicsTable (const G4ParticleDefinition &)

Protected Member Functions
virtual G4VParticleChange *	PostStepDoIt (const G4Track &aTrack, const G4Step &aStep)

virtual G4VParticleChange *	AtRestDoIt (const G4Track &aTrack, const G4Step &aStep)

Protected Member Functions inherited from G4Decay
virtual G4VParticleChange *	DecayIt (const G4Track &aTrack, const G4Step &aStep)

virtual void	DaughterPolarization (const G4Track &aTrack, G4DecayProducts *products)

virtual G4double	GetMeanFreePath (const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition)

virtual G4double	GetMeanLifeTime (const G4Track &aTrack, G4ForceCondition *condition)

Protected Member Functions inherited from G4VProcess
void	SubtractNumberOfInteractionLengthLeft (G4double previousStepSize)

void	ClearNumberOfInteractionLengthLeft ()

Private Member Functions
G4ThreeVector	Spin_Precession (const G4Step &aStep, G4ThreeVector B, G4double deltatime)

Additional Inherited Members
Static Public Member Functions inherited from G4VProcess
static const G4String &	GetProcessTypeName (G4ProcessType)

Protected Attributes inherited from G4Decay
G4int	verboseLevel

const G4double	HighestValue

G4double	fRemainderLifeTime

G4ParticleChangeForDecay	fParticleChangeForDecay

G4VExtDecayer *	pExtDecayer

Protected Attributes inherited from G4VProcess
const G4ProcessManager *	aProcessManager

G4VParticleChange *	pParticleChange

G4ParticleChange	aParticleChange

G4double	theNumberOfInteractionLengthLeft

G4double	currentInteractionLength

G4double	theInitialNumberOfInteractionLength

G4String	theProcessName

G4String	thePhysicsTableFileName

G4ProcessType	theProcessType

G4int	theProcessSubType

G4double	thePILfactor

G4bool	enableAtRestDoIt

G4bool	enableAlongStepDoIt

G4bool	enablePostStepDoIt

G4int	verboseLevel

Detailed Description

Definition at line 43 of file G4DecayWithSpin.hh.

Constructor & Destructor Documentation

◆ G4DecayWithSpin()

G4DecayWithSpin::G4DecayWithSpin ( const G4String & processName = "DecayWithSpin" )

Definition at line 52 of file G4DecayWithSpin.cc.

                                                            :G4Decay(processName)
 {
   // set Process Sub Type   
   SetProcessSubType(static_cast<int>(DECAY_WithSpin));
 
 }

Here is the call graph for this function:

◆ ~G4DecayWithSpin()

G4DecayWithSpin::~G4DecayWithSpin ( )

virtual

Definition at line 59 of file G4DecayWithSpin.cc.

59 {}

Member Function Documentation

◆ AtRestDoIt()

G4VParticleChange * G4DecayWithSpin::AtRestDoIt	(	const G4Track &	aTrack,
		const G4Step &	aStep
	)

protectedvirtual

Reimplemented from G4Decay.

Definition at line 117 of file G4DecayWithSpin.cc.

 {
 
 // get particle
   const G4DynamicParticle* aParticle = aTrack.GetDynamicParticle();
   const G4ParticleDefinition* aParticleDef = aParticle->GetDefinition();
 
 // get parent_polarization
   G4ThreeVector parent_polarization = aParticle->GetPolarization();
 
   if(parent_polarization == G4ThreeVector(0,0,0))
   {
     // Generate random polarization direction
 
     G4double cost = 1. - 2.*G4UniformRand();
     G4double sint = std::sqrt((1.-cost)*(1.+cost));
 
     G4double phi = twopi*G4UniformRand();
     G4double sinp = std::sin(phi);
     G4double cosp = std::cos(phi);
 
     G4double px = sint*cosp;
     G4double py = sint*sinp;
     G4double pz = cost;
 
     parent_polarization.setX(px);
     parent_polarization.setY(py);
     parent_polarization.setZ(pz);
 
   }else{
 
     G4FieldManager* fieldMgr = aStep.GetTrack()->GetVolume()->
                                      GetLogicalVolume()->GetFieldManager();
 
     if (!fieldMgr) {
        G4TransportationManager *transportMgr =
                          G4TransportationManager::GetTransportationManager();
        G4PropagatorInField* fFieldPropagator = 
                                         transportMgr->GetPropagatorInField();
        if (fFieldPropagator) fieldMgr = 
                                   fFieldPropagator->GetCurrentFieldManager();
     }
 
     const G4Field* field = NULL;
     if(fieldMgr)field = fieldMgr->GetDetectorField();
 
     if (field) {
 
        G4double point[4];
        point[0] = (aStep.GetPostStepPoint()->GetPosition())[0];
        point[1] = (aStep.GetPostStepPoint()->GetPosition())[1];
        point[2] = (aStep.GetPostStepPoint()->GetPosition())[2];
        point[3] = aTrack.GetGlobalTime();
 
        G4double fieldValue[6];
        field -> GetFieldValue(point,fieldValue);
 
        G4ThreeVector B(fieldValue[0],fieldValue[1],fieldValue[2]);
 
        // Call the spin precession only for non-zero mag. field
        if (B.mag2() > 0.) parent_polarization = 
                                  Spin_Precession(aStep,B,fRemainderLifeTime);
 
     }
   }
 
 // decay table
   G4DecayTable *decaytable = aParticleDef->GetDecayTable();
   if (decaytable) {
     for (G4int ip=0; ip<decaytable->entries(); ip++){
       decaytable->GetDecayChannel(ip)->SetPolarization(parent_polarization);
     }
   } 
 
   G4ParticleChangeForDecay* pParticleChangeForDecay;
   pParticleChangeForDecay = (G4ParticleChangeForDecay*)G4Decay::DecayIt(aTrack,aStep);
 
   pParticleChangeForDecay->ProposePolarization(parent_polarization);
 
   return pParticleChangeForDecay;
 
 }

Here is the call graph for this function:

◆ PostStepDoIt()

G4VParticleChange * G4DecayWithSpin::PostStepDoIt	(	const G4Track &	aTrack,
		const G4Step &	aStep
	)

protectedvirtual

Reimplemented from G4Decay.

Definition at line 61 of file G4DecayWithSpin.cc.

 {
   if ( (aTrack.GetTrackStatus() == fStopButAlive ) ||
        (aTrack.GetTrackStatus() == fStopAndKill )   ){
     fParticleChangeForDecay.Initialize(aTrack);
     return &fParticleChangeForDecay;
   }
 
 // get particle
   const G4DynamicParticle* aParticle = aTrack.GetDynamicParticle();
   const G4ParticleDefinition* aParticleDef = aParticle->GetDefinition();
 
 // get parent_polarization
   G4ThreeVector parent_polarization = aParticle->GetPolarization();
 
   if(parent_polarization == G4ThreeVector(0,0,0))
   {
     // Generate random polarization direction
 
     G4double cost = 1. - 2.*G4UniformRand();
     G4double sint = std::sqrt((1.-cost)*(1.+cost));
 
     G4double phi = twopi*G4UniformRand();
     G4double sinp = std::sin(phi);
     G4double cosp = std::cos(phi);
 
     G4double px = sint*cosp;
     G4double py = sint*sinp;
     G4double pz = cost;
 
     parent_polarization.setX(px);
     parent_polarization.setY(py);
     parent_polarization.setZ(pz);
 
   }
 
 // decay table
   G4DecayTable *decaytable = aParticleDef->GetDecayTable();
   if (decaytable) {
     for (G4int ip=0; ip<decaytable->entries(); ip++){
       decaytable->GetDecayChannel(ip)->SetPolarization(parent_polarization);
     }
   } 
 
   G4ParticleChangeForDecay* pParticleChangeForDecay;
   pParticleChangeForDecay = (G4ParticleChangeForDecay*)G4Decay::DecayIt(aTrack,aStep);
 
   pParticleChangeForDecay->ProposePolarization(parent_polarization);
 
   //G4cout << parent_polarization.x() << ", "
   //       << parent_polarization.y() << ", "
   //       << parent_polarization.z() << G4endl;
 
   return pParticleChangeForDecay;
 }

Here is the call graph for this function:

◆ Spin_Precession()

G4ThreeVector G4DecayWithSpin::Spin_Precession	(	const G4Step &	aStep,
		G4ThreeVector	B,
		G4double	deltatime
	)

private

Definition at line 200 of file G4DecayWithSpin.cc.

 {
   G4double Bnorm = std::sqrt(sqr(B[0])  + sqr(B[1]) +sqr(B[2]) );
 
   G4double q = aStep.GetTrack()->GetDefinition()->GetPDGCharge();
   G4double a = 1.165922e-3;
   G4double s_omega = 8.5062e+7*rad/(s*kilogauss);
 
   G4double omega = -(q*s_omega)*(1.+a) * Bnorm;
 
   G4double rotationangle = deltatime * omega;
 
   G4Transform3D SpinRotation = G4Rotate3D(rotationangle,B.unit());
 
   G4Vector3D Spin = aStep.GetTrack() -> GetPolarization();
 
   G4Vector3D newSpin = SpinRotation * Spin;
 
 #ifdef G4VERBOSE
   if (GetVerboseLevel()>2) {
     G4double normspin = std::sqrt(Spin*Spin);
     G4double normnewspin = std::sqrt(newSpin*newSpin);
     //G4double cosalpha = Spin*newSpin/normspin/normnewspin;
     //G4double alpha = std::acos(cosalpha);
 
     G4cout << "AT REST::: PARAMETERS " << G4endl;
     G4cout << "Initial spin  : " << Spin  << G4endl;
     G4cout << "Delta time    : " << deltatime  << G4endl;
     G4cout << "Rotation angle: " << rotationangle/rad  << G4endl;
     G4cout << "New spin      : " << newSpin  << G4endl;
     G4cout << "Checked norms : " << normspin <<" " << normnewspin << G4endl;
   }
 #endif
 
   return newSpin;
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

Public Member Functions

Protected Member Functions

Private Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4DecayWithSpin()

◆ ~G4DecayWithSpin()

Member Function Documentation

◆ AtRestDoIt()

◆ PostStepDoIt()

◆ Spin_Precession()