SwitchingFieldManager Class Reference

G4FieldManager that sets tight stepping for muons; disables magnetic field for other particles. More...

#include <SwitchingFieldManager.h>

Inheritance diagram for SwitchingFieldManager:

Collaboration diagram for SwitchingFieldManager:

Public Member Functions
	SwitchingFieldManager (G4Field *field)
	~SwitchingFieldManager ()
void	ConfigureForTrack (const G4Track *) override
G4FieldManager *	Clone () const override

Private Attributes
G4Field *	m_savedField = nullptr
G4double	m_radiusXYmax = 1900
G4double	m_radiusXYmin = 1500
G4double	m_Zmax = 2000
G4double	m_offset = 0.0

Detailed Description

G4FieldManager that sets tight stepping for muons; disables magnetic field for other particles.

Authors

Evangelos Kourlitis, John Apostolakis

Date

2021-06-07

Definition at line 19 of file SwitchingFieldManager.h.

Constructor & Destructor Documentation

SwitchingFieldManager()

SwitchingFieldManager::SwitchingFieldManager

(

G4Field *

field

)

Definition at line 19 of file SwitchingFieldManager.cxx.

   :  G4FieldManager( detectorField ),
      m_savedField( detectorField )
{}

~SwitchingFieldManager()

SwitchingFieldManager::~SwitchingFieldManager

(

)

Definition at line 24 of file SwitchingFieldManager.cxx.

{}

Member Function Documentation

Clone()

G4FieldManager * SwitchingFieldManager::Clone ( ) const

override

Definition at line 55 of file SwitchingFieldManager.cxx.

{
   G4Field* cloneField = m_savedField->Clone();
   auto clone= new SwitchingFieldManager( cloneField );
 
   G4FieldManagerHelper::SetMinAndMaxEpsilonStep( clone, this->GetMinimumEpsilonStep(), this->GetMaximumEpsilonStep() );
   clone->SetDeltaOneStep(       this->GetDeltaOneStep() ); 
   clone->SetDeltaIntersection(  this->GetDeltaIntersection()  );
 
   clone->GetChordFinder()->SetDeltaChord( this->GetChordFinder()->GetDeltaChord() );
   // This lives in our chord finder ... 
   
   return clone;
}

ConfigureForTrack()

void SwitchingFieldManager::ConfigureForTrack ( const G4Track *  track )

override

Definition at line 29 of file SwitchingFieldManager.cxx.

{
   assert(  GetDetectorField() == m_savedField
         || GetDetectorField() == nullptr );
   // Must ensure that SetDetectorField() has NOT been called to establish
   //   a different field by an inside entity
   //   ... it would invalidate our assumptions !
 
   const bool isMuonTrack = (track->GetDefinition()==G4MuonPlus::Definition() ||
                             track->GetDefinition()==G4MuonMinus::Definition() ||
                             track->GetDefinition()==G4ChargedGeantino::ChargedGeantinoDefinition());
 
   const G4ThreeVector position = track->GetPosition();
   const G4double r2XY = position.x()*position.x() + position.y()*position.y();
   const G4double  Z = position.z();
   const G4double r2Min = (m_radiusXYmin + m_offset)*(m_radiusXYmin + m_offset);
   const G4double r2Max = (m_radiusXYmax - m_offset)*(m_radiusXYmax - m_offset);
 
   const bool inSide = (r2XY > r2Min) && (r2XY < r2Max) && (fabs(Z) < m_Zmax - m_offset);
   
   if ( inSide && DoesFieldExist() && !isMuonTrack )
      ChangeDetectorField( nullptr );
   else if ( !inSide && !DoesFieldExist() )
      ChangeDetectorField( m_savedField );
}

Member Data Documentation

m_offset

G4double SwitchingFieldManager::m_offset = 0.0

private

Definition at line 39 of file SwitchingFieldManager.h.

m_radiusXYmax

G4double SwitchingFieldManager::m_radiusXYmax = 1900

private

Definition at line 36 of file SwitchingFieldManager.h.

m_radiusXYmin

G4double SwitchingFieldManager::m_radiusXYmin = 1500

private

Definition at line 37 of file SwitchingFieldManager.h.

m_savedField

G4Field* SwitchingFieldManager::m_savedField = nullptr

private

Definition at line 34 of file SwitchingFieldManager.h.

m_Zmax

G4double SwitchingFieldManager::m_Zmax = 2000

private

Definition at line 38 of file SwitchingFieldManager.h.

---

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

• SwitchingFieldManager.h
• SwitchingFieldManager.cxx