F02ElectricFieldSetup Class Reference

#include <F02ElectricFieldSetup.hh>

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

Protected Member Functions
G4FieldManager *	GetGlobalFieldManager ()

Detailed Description

A class for control of the Electric Field of the detector.

The field for this case is uniform. It is simply a 'setup' class that creates the field and necessary other parts

Definition at line 57 of file F02ElectricFieldSetup.hh.

Constructor & Destructor Documentation

F02ElectricFieldSetup::F02ElectricFieldSetup

(

G4ThreeVector

fieldVector

)

Definition at line 93 of file F02ElectricFieldSetup.cc.

References GetGlobalFieldManager(), and UpdateField().

  : fFieldManager(0),
    fChordFinder(0),
    fEquation(0),
    fEMfield(0),
    fElFieldValue(),
    fStepper(0),
    fIntgrDriver(0),
    fStepperType(4),    // ClassicalRK4 -- the default stepper
    fMinStep(0.010*mm)  // minimal step of 10 microns
{
  fEMfield = new G4UniformElectricField(fieldVector);
  fEquation = new G4EqMagElectricField(fEMfield);

  fFieldManager = GetGlobalFieldManager();
  fFieldMessenger = new F02FieldMessenger(this);

  UpdateField();
}

F02ElectricFieldSetup::F02ElectricFieldSetup

(

)

Definition at line 71 of file F02ElectricFieldSetup.cc.

References python.hepunit::cm, GetGlobalFieldManager(), python.hepunit::kilovolt, and UpdateField().

 : fFieldManager(0),
   fChordFinder(0),
   fEquation(0),
   fEMfield(0),
   fElFieldValue(),
   fStepper(0),
   fIntgrDriver(0),
   fStepperType(4),    // ClassicalRK4 -- the default stepper
   fMinStep(0.010*mm)  // minimal step of 10 microns
{
  fEMfield = new G4UniformElectricField(
                   G4ThreeVector(0.0,100000.0*kilovolt/cm,0.0));
  fEquation = new G4EqMagElectricField(fEMfield);

  fFieldManager = GetGlobalFieldManager();
  fFieldMessenger = new F02FieldMessenger(this);
  UpdateField();
}

F02ElectricFieldSetup::~F02ElectricFieldSetup ( )

virtual

Definition at line 115 of file F02ElectricFieldSetup.cc.

{
  delete fChordFinder;
  delete fStepper;
  delete fEquation;
  delete fEMfield;
  delete fFieldMessenger;
}

Member Function Documentation

G4ThreeVector F02ElectricFieldSetup::GetConstantFieldValue

(

)

G4FieldManager * F02ElectricFieldSetup::GetGlobalFieldManager ( )

protected

Definition at line 243 of file F02ElectricFieldSetup.cc.

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

Referenced by F02ElectricFieldSetup(), and SetFieldValue().

{
//  Utility method

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

void F02ElectricFieldSetup::SetFieldValue

(

G4ThreeVector

fieldVector

)

Definition at line 218 of file F02ElectricFieldSetup.cc.

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

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

{
  if (fEMfield) delete fEMfield;

  // Set the value of the Global Field value to fieldVector

  // Find the Field Manager for the global field
  G4FieldManager* fieldMgr= GetGlobalFieldManager();

  if (fieldVector != G4ThreeVector(0.,0.,0.))
  {
    fEMfield = new G4UniformElectricField(fieldVector);
  }
  else
  {
    // If the new field's value is Zero, then it is best to
    //  insure that it is not used for propagation.
    fEMfield = 0;
  }
  fieldMgr->SetDetectorField(fEMfield);
  fEquation->SetFieldObj(fEMfield);  // must now point to the new field
}

void F02ElectricFieldSetup::SetFieldValue

(

G4double

fieldValue

)

Definition at line 207 of file F02ElectricFieldSetup.cc.

References SetFieldValue().

{
  // Set the value of the Global Field to fieldValue along Z

  G4ThreeVector fieldVector( 0.0, 0.0, fieldValue );

  SetFieldValue( fieldVector );
}

void F02ElectricFieldSetup::SetMinStep ( G4double  s )

inline

Definition at line 69 of file F02ElectricFieldSetup.hh.

Referenced by F02FieldMessenger::SetNewValue().

{ fMinStep = s ; }

void F02ElectricFieldSetup::SetStepper

(

)

Definition at line 150 of file F02ElectricFieldSetup.cc.

References G4cout, and G4endl.

Referenced by UpdateField().

{
// Set stepper according to the stepper type

  G4int nvar = 8;

  if (fStepper) delete fStepper;

  switch ( fStepperType )
  {
    case 0:
      fStepper = new G4ExplicitEuler( fEquation, nvar );
      G4cout<<"G4ExplicitEuler is calledS"<<G4endl;
      break;
    case 1:
      fStepper = new G4ImplicitEuler( fEquation, nvar );
      G4cout<<"G4ImplicitEuler is called"<<G4endl;
      break;
    case 2:
      fStepper = new G4SimpleRunge( fEquation, nvar );
      G4cout<<"G4SimpleRunge is called"<<G4endl;
      break;
    case 3:
      fStepper = new G4SimpleHeum( fEquation, nvar );
      G4cout<<"G4SimpleHeum is called"<<G4endl;
      break;
    case 4:
      fStepper = new G4ClassicalRK4( fEquation, nvar );
      G4cout<<"G4ClassicalRK4 (default) is called"<<G4endl;
      break;
    case 5:
      fStepper = new G4CashKarpRKF45( fEquation, nvar );
      G4cout<<"G4CashKarpRKF45 is called"<<G4endl;
      break;
    case 6:
      fStepper = 0; // new G4RKG3_Stepper( fEquation, nvar );
      G4cout<<"G4RKG3_Stepper is not currently working for Electric Field"
            <<G4endl;
      break;
    case 7:
      fStepper = 0; // new G4HelixExplicitEuler( fEquation );
      G4cout<<"G4HelixExplicitEuler is not valid for Electric Field"<<G4endl;
      break;
    case 8:
      fStepper = 0; // new G4HelixImplicitEuler( fEquation );
      G4cout<<"G4HelixImplicitEuler is not valid for Electric Field"<<G4endl;
      break;
    case 9:
      fStepper = 0; // new G4HelixSimpleRunge( fEquation );
      G4cout<<"G4HelixSimpleRunge is not valid for Electric Field"<<G4endl;
      break;
    default: fStepper = 0;
  }
}

void F02ElectricFieldSetup::SetStepperType ( G4int  i )

inline

Definition at line 65 of file F02ElectricFieldSetup.hh.

Referenced by F02FieldMessenger::SetNewValue().

{ fStepperType = i ; }

void F02ElectricFieldSetup::UpdateField

(

)

Definition at line 126 of file F02ElectricFieldSetup.cc.

References G4cout, G4endl, G4MagIntegratorStepper::GetNumberOfVariables(), python.hepunit::mm, G4FieldManager::SetChordFinder(), G4FieldManager::SetDetectorField(), and SetStepper().

Referenced by F02ElectricFieldSetup(), and F02FieldMessenger::SetNewValue().

{
// Register this field to 'global' Field Manager and
// Create Stepper and Chord Finder with predefined type, minstep (resp.)

  SetStepper();

  G4cout<<"The minimal step is equal to "<<fMinStep/mm<<" mm"<<G4endl;

  fFieldManager->SetDetectorField(fEMfield);

  if (fChordFinder) delete fChordFinder;

  fIntgrDriver = new G4MagInt_Driver(fMinStep,
                                     fStepper,
                                     fStepper->GetNumberOfVariables());

  fChordFinder = new G4ChordFinder(fIntgrDriver);

  fFieldManager->SetChordFinder(fChordFinder);
}

---

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

• F02ElectricFieldSetup.hh
• F02ElectricFieldSetup.cc