#include <F01FieldSetup.hh>

Public Member Functions
	F01FieldSetup (G4ThreeVector)

	F01FieldSetup ()

virtual	~F01FieldSetup ()

void	SetStepperType (G4int i)

void	SetStepper ()

void	SetMinStep (G4double s)

void	InitialiseAll ()

void	CreateStepperAndChordFinder ()

void	SetFieldValue (G4ThreeVector fieldVector)

void	SetFieldValue (G4double fieldValue)

G4ThreeVector	GetConstantFieldValue ()

Protected Member Functions
G4FieldManager *	GetGlobalFieldManager ()

Protected Attributes
G4FieldManager *	fFieldManager

G4ChordFinder *	fChordFinder

G4Mag_UsualEqRhs *	fEquation

G4MagneticField *	fMagneticField

G4MagIntegratorStepper *	fStepper

G4int	fStepperType

G4double	fMinStep

F01FieldMessenger *	fFieldMessenger

Detailed Description

Definition at line 57 of file F01FieldSetup.hh.

Constructor & Destructor Documentation

F01FieldSetup::F01FieldSetup ( G4ThreeVector fieldVector )

Definition at line 67 of file F01FieldSetup.cc.

References G4cout, G4endl, and InitialiseAll().

  : fFieldManager(0),
    fChordFinder(0),
    fEquation(0),
    fMagneticField(new G4UniformMagField(fieldVector)),
    fStepper(0),
    fStepperType(0),
    fMinStep(0.),
    fFieldMessenger(0)
 {
   G4cout << " F01FieldSetup: magnetic field set to Uniform( "
          << fieldVector << " ) " << G4endl;
   InitialiseAll();
 }

F01FieldSetup::F01FieldSetup ( )

Definition at line 84 of file F01FieldSetup.cc.

References G4cout, G4endl, and InitialiseAll().

  : fFieldManager(0),
    fChordFinder(0),
    fEquation(0),
    fMagneticField(new G4UniformMagField(G4ThreeVector())),
    fStepper(0),
    fStepperType(0),
    fMinStep(0.),
    fFieldMessenger(0)
 {
   G4cout << " F01FieldSetup: magnetic field set to Uniform( 0.0, 0, 0 ) "
          << G4endl;
   InitialiseAll();
 }

F01FieldSetup::~F01FieldSetup ( )

virtual

Definition at line 118 of file F01FieldSetup.cc.

References fChordFinder, fFieldMessenger, fMagneticField, and fStepper.

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

Member Function Documentation

void F01FieldSetup::CreateStepperAndChordFinder ( )

Definition at line 128 of file F01FieldSetup.cc.

References fChordFinder, fFieldManager, fMagneticField, fMinStep, fStepper, G4cout, G4endl, python.hepunit::mm, G4FieldManager::SetChordFinder(), G4FieldManager::SetDetectorField(), and SetStepper().

Referenced by InitialiseAll(), and F01FieldMessenger::SetNewValue().

 {
   // Update field
 
   SetStepper();
   G4cout<<"The minimal step is equal to "<<fMinStep/mm<<" mm"<<G4endl;
 
   fFieldManager->SetDetectorField(fMagneticField );
 
   if (fChordFinder) delete fChordFinder;
 
   fChordFinder = new G4ChordFinder( fMagneticField, fMinStep,fStepper );
 
   fFieldManager->SetChordFinder( fChordFinder );
 }

G4ThreeVector F01FieldSetup::GetConstantFieldValue ( )

G4FieldManager * F01FieldSetup::GetGlobalFieldManager ( )

protected

Definition at line 262 of file F01FieldSetup.cc.

References G4TransportationManager::GetFieldManager(), and G4TransportationManager::GetTransportationManager().

Referenced by SetFieldValue().

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

void F01FieldSetup::InitialiseAll ( )

Definition at line 101 of file F01FieldSetup.cc.

References CreateStepperAndChordFinder(), fEquation, fFieldManager, fFieldMessenger, fMagneticField, fMinStep, fStepperType, G4TransportationManager::GetFieldManager(), G4TransportationManager::GetTransportationManager(), and python.hepunit::mm.

Referenced by F01FieldSetup().

 {
   fFieldMessenger = new F01FieldMessenger(this);
  
   fEquation = new G4Mag_UsualEqRhs(fMagneticField);
  
   fMinStep     = 1.0*mm; // minimal step of 1 mm is default
 
   fStepperType = 4;      // ClassicalRK4 is default stepper
 
   fFieldManager = G4TransportationManager::GetTransportationManager()
                     ->GetFieldManager();
   CreateStepperAndChordFinder();
 }

void F01FieldSetup::SetFieldValue ( G4ThreeVector fieldVector )

Definition at line 235 of file F01FieldSetup.cc.

References fEquation, fMagneticField, GetGlobalFieldManager(), G4FieldManager::SetDetectorField(), and G4EquationOfMotion::SetFieldObj().

Referenced by SetFieldValue(), and F01FieldMessenger::SetNewValue().

 {
   // Set the value of the Global Field
 
   if (fMagneticField) delete fMagneticField;
  
   if (fieldVector != G4ThreeVector(0.,0.,0.))
   {
     fMagneticField = new G4UniformMagField(fieldVector);
   }
   else
   {
     // If the new field's value is Zero, signal it as below
     // so that it is not used for propagation.
     fMagneticField = 0;
   }
 
   // Set this as the field of the global Field Manager
   GetGlobalFieldManager()->SetDetectorField(fMagneticField);
 
   // Now notify equation of new field
   fEquation->SetFieldObj( fMagneticField );
 
 }

void F01FieldSetup::SetFieldValue ( G4double fieldValue )

Definition at line 200 of file F01FieldSetup.cc.

References FatalException, fMagneticField, G4cout, G4endl, G4Exception(), python.hepunit::gauss, G4MagneticField::GetFieldValue(), and SetFieldValue().

 {
   // Set the value of the Global Field to fieldValue along Z
 
 #ifdef G4VERBOSE
   G4cout << "Setting Field strength to "
          << fieldStrength / gauss  << " Gauss."; // << G4endl;
 #endif
 
   G4ThreeVector fieldSetVec(0.0, 0.0, fieldStrength);
   this->SetFieldValue( fieldSetVec );
 
 #ifdef G4VERBOSE
   G4double fieldValue[6],  position[4];
   position[0] = position[1] = position[2] = position[3] = 0.0;
   if ( fieldStrength != 0.0 ) {
     fMagneticField->GetFieldValue( position, fieldValue);
     G4ThreeVector fieldVec(fieldValue[0], fieldValue[1], fieldValue[2]);
     // G4cout << " fMagneticField is now " << fMagneticField
     G4cout << " Magnetic field vector is "
            << fieldVec / gauss << " G " << G4endl;
   } else {
     if ( fMagneticField == 0 )
       G4cout << " Magnetic field pointer is null." << G4endl;
     else
       G4Exception("F01FieldSetup::SetFieldValue(double)",
                   "IncorrectForZeroField",
                   FatalException,
                   "fMagneticField ptr should be set to 0 for no field.");
   }
 #endif
 }

void F01FieldSetup::SetMinStep ( G4double s )

inline

Definition at line 69 of file F01FieldSetup.hh.

References fMinStep.

Referenced by F01FieldMessenger::SetNewValue().

69 { fMinStep = s; }

F01FieldSetup::fMinStep

G4double fMinStep

Definition: F01FieldSetup.hh:92

s

const XML_Char * s

Definition: include/expat.h:262

void F01FieldSetup::SetStepper ( )

Definition at line 146 of file F01FieldSetup.cc.

References fEquation, fStepper, fStepperType, G4cout, and G4endl.

Referenced by CreateStepperAndChordFinder().

 {
 // Set stepper according to the stepper type
 
   if (fStepper) delete fStepper;
 
   switch ( fStepperType )
   {
     case 0:
       fStepper = new G4ExplicitEuler( fEquation );
       G4cout<<"G4ExplicitEuler is calledS"<<G4endl;
       break;
     case 1:
       fStepper = new G4ImplicitEuler( fEquation );
       G4cout<<"G4ImplicitEuler is called"<<G4endl;
       break;
     case 2:
       fStepper = new G4SimpleRunge( fEquation );
       G4cout<<"G4SimpleRunge is called"<<G4endl;
       break;
     case 3:
       fStepper = new G4SimpleHeum( fEquation );
       G4cout<<"G4SimpleHeum is called"<<G4endl;
       break;
     case 4:
       fStepper = new G4ClassicalRK4( fEquation );
       G4cout<<"G4ClassicalRK4 (default) is called"<<G4endl;
       break;
     case 5:
       fStepper = new G4HelixExplicitEuler( fEquation );
       G4cout<<"G4HelixExplicitEuler is called"<<G4endl;
       break;
     case 6:
       fStepper = new G4HelixImplicitEuler( fEquation );
       G4cout<<"G4HelixImplicitEuler is called"<<G4endl;
       break;
     case 7:
       fStepper = new G4HelixSimpleRunge( fEquation );
       G4cout<<"G4HelixSimpleRunge is called"<<G4endl;
       break;
     case 8:
       fStepper = new G4CashKarpRKF45( fEquation );
       G4cout<<"G4CashKarpRKF45 is called"<<G4endl;
       break;
     case 9:
       fStepper = new G4RKG3_Stepper( fEquation );
       G4cout<<"G4RKG3_Stepper is called"<<G4endl;
       break;
     default: fStepper = 0;
   }
 }