Defines the physics process that avoids hit invocation in sensitive volumes. More...

#include <ExtMagFieldLimitProcess.h>

Inheritance diagram for ExtMagFieldLimitProcess:

Public Member Functions
	ExtMagFieldLimitProcess (const G4String &processName="ExtMagFieldLimitProcess")
	constructor

	~ExtMagFieldLimitProcess ()
	destructor

G4VParticleChange *	PostStepDoIt (const G4Track &, const G4Step &)
	Do nothing special after the particle has stepped.

G4double	PostStepGetPhysicalInteractionLength (const G4Track &, G4double, G4ForceCondition *)
	Returns the step length after each step.

G4double	GetMeanFreePath (const G4Track &, G4double, G4ForceCondition *)
	Returns the mean free path (always infinity!) after each step.

G4double	GetStepLimit () const
	Returns the step limit.

void	SetStepLimit (G4double val)
	Sets the step limit.

Protected Attributes
const G4Field *	m_field
	Stores the pointer to the magnetic field class.

G4double	m_stepLimit
	Stores the step limit.

Detailed Description

Defines the physics process that avoids hit invocation in sensitive volumes.

Definition at line 27 of file ExtMagFieldLimitProcess.h.

Constructor & Destructor Documentation

◆ ExtMagFieldLimitProcess()

ExtMagFieldLimitProcess ( const G4String & processName = "ExtMagFieldLimitProcess" )

explicit

constructor

Definition at line 22 of file ExtMagFieldLimitProcess.cc.

                                                                            :
  G4VDiscreteProcess(processName),
  m_stepLimit(kInfinity)   // user may change this with a geant4 UI command
{
  m_field = G4TransportationManager::GetTransportationManager()->GetFieldManager()->GetDetectorField();
  if (false) {
    G4Track aTrack;
    G4Step aStep;
    G4ForceCondition* condition = nullptr;
    GetMeanFreePath(aTrack, 0.0, condition);
    PostStepDoIt(aTrack, aStep);
    PostStepGetPhysicalInteractionLength(aTrack, 0.0, condition);
  }
 
}

◆ ~ExtMagFieldLimitProcess()

~ExtMagFieldLimitProcess ( )

destructor

Definition at line 38 of file ExtMagFieldLimitProcess.cc.

39{

40}

Member Function Documentation

◆ GetMeanFreePath()

G4double GetMeanFreePath	(	const G4Track &	,
		G4double	,
		G4ForceCondition *
	)

Returns the mean free path (always infinity!) after each step.

Definition at line 42 of file ExtMagFieldLimitProcess.cc.

{
  return kInfinity;
}

◆ GetStepLimit()

G4double GetStepLimit ( ) const

inline

Returns the step limit.

Definition at line 47 of file ExtMagFieldLimitProcess.h.

47{ return m_stepLimit; }

◆ PostStepDoIt()

G4VParticleChange * PostStepDoIt	(	const G4Track &	track,
		const G4Step &
	)

Do nothing special after the particle has stepped.

Definition at line 69 of file ExtMagFieldLimitProcess.cc.

{
  aParticleChange.Initialize(track);
  return &aParticleChange;
}

◆ PostStepGetPhysicalInteractionLength()

G4double PostStepGetPhysicalInteractionLength	(	const G4Track &	aTrack,
		G4double	,
		G4ForceCondition *	condition
	)

Returns the step length after each step.

Definition at line 48 of file ExtMagFieldLimitProcess.cc.

{
  *condition = NotForced;
  G4double stepLength = kInfinity;
  if (m_field != 0) {
    G4ThreeVector trkPosi = aTrack.GetPosition();
    G4double pos1[3] = { trkPosi.x(), trkPosi.y(), trkPosi.z() };
    G4double h1[3] = { 0.0, 0.0, 0.0 };
    m_field->GetFieldValue(pos1, h1);
    G4ThreeVector BVec(h1[0], h1[1], h1[2]);
    G4double pmag = aTrack.GetMomentum().mag();
    G4double BPerpMom = BVec.cross(aTrack.GetMomentum()).mag() / pmag;   // LEP
    if (BPerpMom != 0.0) {
      stepLength = m_stepLimit * pmag / BPerpMom;
    }
    B2DEBUG(300, "ExtMagFieldLimitProcess::PostStepGetPhysicalInteractionLength() stepLength "
            << stepLength << " B " << BPerpMom << " BVec " << BVec << " pmag " << pmag);
  }
  return stepLength;
}

◆ SetStepLimit()

void SetStepLimit ( G4double val )

inline

Sets the step limit.

Definition at line 50 of file ExtMagFieldLimitProcess.h.

50{ m_stepLimit = val; }

Member Data Documentation

◆ m_field

const G4Field* m_field

protected

Stores the pointer to the magnetic field class.

Definition at line 55 of file ExtMagFieldLimitProcess.h.

◆ m_stepLimit

G4double m_stepLimit