EDepSim::EMFieldSetup Class Reference

#include <EDepSimEMFieldSetup.hh>

Public Member Functions
	EMFieldSetup (G4ElectroMagneticField *field, G4FieldManager *m=0)
	EMFieldSetup (G4FieldManager *m=0)
	Create a zero field for the field manager. More...
virtual	~EMFieldSetup ()
void	SetStepperType (G4int i)
void	SetStepper ()
void	SetMinStep (G4double s)

Protected Member Functions
G4FieldManager *	GetGlobalFieldManager ()
void	UpdateField ()

Private Attributes
G4FieldManager *	fFieldManager
G4ChordFinder *	fChordFinder
G4EqMagElectricField *	fEquation
G4ElectroMagneticField *	fEMfield
G4MagIntegratorStepper *	fStepper
G4MagInt_Driver *	fIntgrDriver
G4int	fStepperType
G4double	fMinStep

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 58 of file EDepSimEMFieldSetup.hh.

Constructor & Destructor Documentation

EDepSim::EMFieldSetup::EMFieldSetup	(	G4ElectroMagneticField *	field,
		G4FieldManager *	m = 0
	)

Create a new field. The first argument is a general EM field (for example, G4UniformMagneticField, or G4UniformElectricField. The second argument is the field manager to setup.

Definition at line 90 of file EDepSimEMFieldSetup.cc.

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

    UpdateField();
}

EDepSim::EMFieldSetup::EMFieldSetup

(

G4FieldManager *

m = 0

)

Create a zero field for the field manager.

Definition at line 72 of file EDepSimEMFieldSetup.cc.

    : fFieldManager(fieldManager),
      fChordFinder(0),
      fEquation(0),
      fEMfield(0),
      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,500.0*volt/cm,0.0));
    fEquation = new G4EqMagElectricField(fEMfield);

    UpdateField();
}

EDepSim::EMFieldSetup::~EMFieldSetup ( )

virtual

Definition at line 108 of file EDepSimEMFieldSetup.cc.

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

Member Function Documentation

G4FieldManager * EDepSim::EMFieldSetup::GetGlobalFieldManager ( )

protected

Definition at line 201 of file EDepSimEMFieldSetup.cc.

{
//  Utility method

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

void EDepSim::EMFieldSetup::SetMinStep ( G4double  s )

inline

Definition at line 75 of file EDepSimEMFieldSetup.hh.

{ fMinStep = s ; }

void EDepSim::EMFieldSetup::SetStepper

(

)

Definition at line 144 of file EDepSimEMFieldSetup.cc.

{
// 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 EDepSim::EMFieldSetup::SetStepperType ( G4int  i )

inline

Definition at line 71 of file EDepSimEMFieldSetup.hh.

{ fStepperType = i ; }

void EDepSim::EMFieldSetup::UpdateField ( )

protected

Definition at line 118 of file EDepSimEMFieldSetup.cc.

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

    if (!fFieldManager) fFieldManager = GetGlobalFieldManager();

    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);
}

Member Data Documentation

G4ChordFinder* EDepSim::EMFieldSetup::fChordFinder

private

Definition at line 89 of file EDepSimEMFieldSetup.hh.

G4ElectroMagneticField* EDepSim::EMFieldSetup::fEMfield

private

Definition at line 93 of file EDepSimEMFieldSetup.hh.

G4EqMagElectricField* EDepSim::EMFieldSetup::fEquation

private

Definition at line 91 of file EDepSimEMFieldSetup.hh.

G4FieldManager* EDepSim::EMFieldSetup::fFieldManager

private

Definition at line 87 of file EDepSimEMFieldSetup.hh.

G4MagInt_Driver* EDepSim::EMFieldSetup::fIntgrDriver

private

Definition at line 96 of file EDepSimEMFieldSetup.hh.

G4double EDepSim::EMFieldSetup::fMinStep

private

Definition at line 100 of file EDepSimEMFieldSetup.hh.

G4MagIntegratorStepper* EDepSim::EMFieldSetup::fStepper

private

Definition at line 95 of file EDepSimEMFieldSetup.hh.

G4int EDepSim::EMFieldSetup::fStepperType

private

Definition at line 98 of file EDepSimEMFieldSetup.hh.

---

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

• edep-sim/src/EDepSimEMFieldSetup.hh
• edep-sim/src/EDepSimEMFieldSetup.cc