F03FieldSetup Class Reference

#include <F03FieldSetup.hh>

Collaboration diagram for F03FieldSetup:

Public Member Functions
	F03FieldSetup ()
virtual	~F03FieldSetup ()
void	SetStepperType (G4int i)
void	SetStepper ()
void	SetMinStep (G4double s)
void	SetFieldValue (G4ThreeVector fieldVector)
void	SetFieldValue (G4double fieldValue)
G4ThreeVector	GetConstantFieldValue ()
void	UpdateField ()
G4FieldManager *	GetLocalFieldManager ()

Protected Member Functions
G4FieldManager *	GetGlobalFieldManager ()

Protected Attributes
G4FieldManager *	fFieldManager
G4FieldManager *	fLocalFieldManager
G4ChordFinder *	fChordFinder
G4ChordFinder *	fLocalChordFinder
G4Mag_UsualEqRhs *	fEquation
G4Mag_UsualEqRhs *	fLocalEquation
G4MagneticField *	fMagneticField
G4MagneticField *	fLocalMagneticField
G4MagIntegratorStepper *	fStepper
G4MagIntegratorStepper *	fLocalStepper
G4int	fStepperType
G4double	fMinStep
F03FieldMessenger *	fFieldMessenger

Detailed Description

A class for setting up the Magnetic Field

It also creates the necessary classes to control accuracy of propagation. In this example

• There is a global field for most of the setup;
• A local field overides it for some volume(s) and it assumed to be uniform.

Definition at line 56 of file F03FieldSetup.hh.

Constructor & Destructor Documentation

F03FieldSetup::F03FieldSetup

(

)

Definition at line 66 of file F03FieldSetup.cc.

 : fFieldManager(0),
   fLocalFieldManager(0),
   fChordFinder(0),
   fLocalChordFinder(0),
   fEquation(0),
   fLocalEquation(0),
   fMagneticField(0),
   fLocalMagneticField(0),
   fStepper(0),
   fLocalStepper(0),
   fFieldMessenger(0)
{
  fMagneticField = new G4UniformMagField(G4ThreeVector(3.3*tesla,
                                                       0.0, // 0.5*tesla,
                                                       0.0));
  fLocalMagneticField = new G4UniformMagField(G4ThreeVector(3.3*tesla,
                                                            0.0, // 0.5*tesla,
                                                            0.0));

  fFieldMessenger = new F03FieldMessenger(this);
 
  fEquation = new G4Mag_UsualEqRhs(fMagneticField);
  fLocalEquation = new G4Mag_UsualEqRhs(fLocalMagneticField);
 
  fMinStep     = 0.25*mm ; // minimal step of 1 mm is default
  fStepperType = 4 ;       // ClassicalRK4 is default stepper

  fFieldManager = GetGlobalFieldManager();
  fLocalFieldManager = new G4FieldManager();

  UpdateField();
}

Here is the call graph for this function:

F03FieldSetup::~F03FieldSetup ( )

virtual

Definition at line 102 of file F03FieldSetup.cc.

{
  delete fMagneticField;
  delete fChordFinder;
  delete fStepper;
  delete fFieldMessenger;
}

Member Function Documentation

G4ThreeVector F03FieldSetup::GetConstantFieldValue

(

)

Definition at line 239 of file F03FieldSetup.cc.

{
  static G4double fieldValue[6],  position[4];
  position[0] = position[1] = position[2] = position[3] = 0.0;

  fMagneticField->GetFieldValue( position, fieldValue);
  G4ThreeVector fieldVec(fieldValue[0], fieldValue[1], fieldValue[2]);

  return fieldVec;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4FieldManager * F03FieldSetup::GetGlobalFieldManager ( )

protected

Definition at line 231 of file F03FieldSetup.cc.

{
  return G4TransportationManager::GetTransportationManager()
                                  ->GetFieldManager();
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4FieldManager* F03FieldSetup::GetLocalFieldManager ( )

inline

Definition at line 74 of file F03FieldSetup.hh.

{ return fLocalFieldManager;}

Here is the caller graph for this function:

void F03FieldSetup::SetFieldValue

(

G4ThreeVector

fieldVector

)

Definition at line 203 of file F03FieldSetup.cc.

{
  if(fMagneticField) delete fMagneticField;

  if(fieldVector != G4ThreeVector(0.,0.,0.))
  {
    fMagneticField = new  G4UniformMagField(fieldVector);
  }
  else
  {
    // If the new field's value is Zero, then
    // setting the pointer to zero ensures
    // that it is not used for propagation.
    fMagneticField = 0;
  }

  // Either
  //   - UpdateField() to reset all (ChordFinder, Equation);
  // UpdateField();
  //     or simply update the field manager & equation of motion
  //     with pointer to new field
  GetGlobalFieldManager()->SetDetectorField(fMagneticField);
  fEquation->SetFieldObj( fMagneticField );

}

Here is the call graph for this function:

Here is the caller graph for this function:

void F03FieldSetup::SetFieldValue

(

G4double

fieldValue

)

Definition at line 195 of file F03FieldSetup.cc.

{
  G4ThreeVector fieldSetVec(0.0, 0.0, fieldStrength);
  SetFieldValue( fieldSetVec );
}

Here is the call graph for this function:

void F03FieldSetup::SetMinStep ( G4double  s )

inline

Definition at line 66 of file F03FieldSetup.hh.

{ fMinStep = s; }

Here is the caller graph for this function:

void F03FieldSetup::SetStepper

(

)

Definition at line 133 of file F03FieldSetup.cc.

{
  if (fStepper) delete fStepper;

  switch ( fStepperType )
  {
    case 0:
      fStepper = new G4ExplicitEuler( fEquation );
      fLocalStepper = new G4ExplicitEuler( fLocalEquation );
      G4cout<<"G4ExplicitEuler is called"<<G4endl;
      break;
    case 1:
      fStepper = new G4ImplicitEuler( fEquation );
      fLocalStepper = new G4ImplicitEuler( fLocalEquation );
      G4cout<<"G4ImplicitEuler is called"<<G4endl;
      break;
    case 2:
      fStepper = new G4SimpleRunge( fEquation );
      fLocalStepper = new G4SimpleRunge( fLocalEquation );
      G4cout<<"G4SimpleRunge is called"<<G4endl;
      break;
    case 3:
      fStepper = new G4SimpleHeum( fEquation );
      fLocalStepper = new G4SimpleHeum( fLocalEquation );
      G4cout<<"G4SimpleHeum is called"<<G4endl;
      break;
    case 4:
      fStepper = new G4ClassicalRK4( fEquation );
      fLocalStepper = new G4ClassicalRK4( fLocalEquation );
      G4cout<<"G4ClassicalRK4 (default) is called"<<G4endl;
      break;
    case 5:
      fStepper = new G4HelixExplicitEuler( fEquation );
      fLocalStepper = new G4HelixExplicitEuler( fLocalEquation );
      G4cout<<"G4HelixExplicitEuler is called"<<G4endl;
      break;
    case 6:
      fStepper = new G4HelixImplicitEuler( fEquation );
      fLocalStepper = new G4HelixImplicitEuler( fLocalEquation );
      G4cout<<"G4HelixImplicitEuler is called"<<G4endl;
      break;
    case 7:
      fStepper = new G4HelixSimpleRunge( fEquation );
      fLocalStepper = new G4HelixSimpleRunge( fLocalEquation );
      G4cout<<"G4HelixSimpleRunge is called"<<G4endl;
      break;
    case 8:
      fStepper = new G4CashKarpRKF45( fEquation );
      fLocalStepper = new G4CashKarpRKF45( fLocalEquation );
      G4cout<<"G4CashKarpRKF45 is called"<<G4endl;
      break;
    case 9:
      fStepper = new G4RKG3_Stepper( fEquation );
      fLocalStepper = new G4RKG3_Stepper( fLocalEquation );
      G4cout<<"G4RKG3_Stepper is called"<<G4endl;
      break;
    default: fStepper = 0;
  }
}

Here is the caller graph for this function:

void F03FieldSetup::SetStepperType ( G4int  i )

inline

Definition at line 62 of file F03FieldSetup.hh.

{ fStepperType = i; }

Here is the caller graph for this function:

void F03FieldSetup::UpdateField

(

)

Definition at line 112 of file F03FieldSetup.cc.

{
  SetStepper();
  G4cout<<"The minimal step is equal to "<<fMinStep/mm<<" mm"<<G4endl;

  fFieldManager->SetDetectorField(fMagneticField);
  fLocalFieldManager->SetDetectorField(fLocalMagneticField);

  if (fChordFinder) delete fChordFinder;
  if (fLocalChordFinder) delete fLocalChordFinder;

  fChordFinder = new G4ChordFinder(fMagneticField, fMinStep, fStepper);
  fLocalChordFinder = new G4ChordFinder(fLocalMagneticField,
                                        fMinStep,fLocalStepper);

  fFieldManager->SetChordFinder(fChordFinder);
  fLocalFieldManager->SetChordFinder(fLocalChordFinder);
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

G4ChordFinder* F03FieldSetup::fChordFinder

protected

Definition at line 84 of file F03FieldSetup.hh.

G4Mag_UsualEqRhs* F03FieldSetup::fEquation

protected

Definition at line 86 of file F03FieldSetup.hh.

G4FieldManager* F03FieldSetup::fFieldManager

protected

Definition at line 82 of file F03FieldSetup.hh.

F03FieldMessenger* F03FieldSetup::fFieldMessenger

protected

Definition at line 97 of file F03FieldSetup.hh.

G4ChordFinder* F03FieldSetup::fLocalChordFinder

protected

Definition at line 85 of file F03FieldSetup.hh.

G4Mag_UsualEqRhs* F03FieldSetup::fLocalEquation

protected

Definition at line 87 of file F03FieldSetup.hh.

G4FieldManager* F03FieldSetup::fLocalFieldManager

protected

Definition at line 83 of file F03FieldSetup.hh.

G4MagneticField* F03FieldSetup::fLocalMagneticField

protected

Definition at line 89 of file F03FieldSetup.hh.

G4MagIntegratorStepper* F03FieldSetup::fLocalStepper

protected

Definition at line 92 of file F03FieldSetup.hh.

G4MagneticField* F03FieldSetup::fMagneticField

protected

Definition at line 88 of file F03FieldSetup.hh.

G4double F03FieldSetup::fMinStep

protected

Definition at line 95 of file F03FieldSetup.hh.

G4MagIntegratorStepper* F03FieldSetup::fStepper

protected

Definition at line 91 of file F03FieldSetup.hh.

G4int F03FieldSetup::fStepperType

protected

Definition at line 93 of file F03FieldSetup.hh.

---

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

• source/geant4.10.03.p03/examples/extended/field/field03/include/F03FieldSetup.hh
• source/geant4.10.03.p03/examples/extended/field/field03/src/F03FieldSetup.cc