BesMagneticField Class Reference

#include <BesMagneticField.hh>

List of all members.

Public Member Functions
	BesMagneticField ()
	~BesMagneticField ()
void	GetFieldValue (const double Point[3], double *Bfield) const
void	SetStepperType (G4int i)
void	SetStepper ()
void	SetMinStep (G4double s)
void	InitialiseAll ()
void	CreateStepperAndChordFinder ()
void	SetDeltaOneStep (double newvalue)
void	SetDeltaIntersection (double newvalue)
void	SetMinimumEpsilonStep (double newvalue)
void	SetMaximumEpsilonStep (double newvalue)
Protected Attributes
G4FieldManager *	fFieldManager
G4ChordFinder *	fChordFinder
G4Mag_UsualEqRhs *	fEquation
G4MagIntegratorStepper *	fStepper
G4int	fStepperType
G4double	fMinStep
BesMagneticFieldMessenger *	fFieldMessenger
IMagneticFieldSvc *	m_pIMF

---

Detailed Description

Definition at line 56 of file BesMagneticField.hh.

---

Constructor & Destructor Documentation

BesMagneticField::BesMagneticField

(

)

Definition at line 83 of file BesMagneticField.cc.

References InitialiseAll(), and m_pIMF.

  : fChordFinder(0), fStepper(0),m_pIMF(0)
{
  ISvcLocator* svcLocator = Gaudi::svcLocator();
  StatusCode sc = svcLocator->service("MagneticFieldSvc",m_pIMF);
  if(sc!=StatusCode::SUCCESS) {
    G4cout<< "Unable to open Magnetic field service"<<G4endl;
  }
  InitialiseAll();
}

BesMagneticField::~BesMagneticField

(

)

Definition at line 94 of file BesMagneticField.cc.

References fChordFinder, fEquation, fFieldMessenger, and fStepper.

{
  if(fFieldMessenger) delete fFieldMessenger;
  if(fChordFinder)   delete fChordFinder;
  if(fEquation)      delete fEquation;
  if(fStepper)       delete fStepper;
}

---

Member Function Documentation

void BesMagneticField::CreateStepperAndChordFinder

(

)

Definition at line 149 of file BesMagneticField.cc.

References fChordFinder, fFieldManager, fMinStep, fStepper, and SetStepper().

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

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

  fFieldManager->SetDetectorField(this );

  if(fChordFinder) delete fChordFinder;

   fChordFinder = new G4ChordFinder(this , fMinStep,fStepper);
 
   fChordFinder->SetDeltaChord(0.25*mm);
   fFieldManager->SetChordFinder( fChordFinder );
   fFieldManager->SetDeltaOneStep(1.0e-2*mm);
   fFieldManager->SetDeltaIntersection(1.0e-3*mm);
   fFieldManager->SetMinimumEpsilonStep(5.0e-5);
   fFieldManager->SetMaximumEpsilonStep(1.0e-3);

   G4PropagatorInField* fieldPropagator
                 = G4TransportationManager::GetTransportationManager()
                      ->GetPropagatorInField();
   G4cout<<"LargestAcceptableStep is "<<fieldPropagator->GetLargestAcceptableStep()/m<<G4endl; 
//read some values
    G4cout<<"field has created"<<G4endl;
    G4cout<<"fDelta_One_Step_Value is "<<fFieldManager->GetDeltaOneStep()<<G4endl;
    G4cout<<"fDelta_Intersection_Val is "<<fFieldManager->GetDeltaIntersection()<<G4endl;
    G4cout<<"fEpsilonMin is "<<fFieldManager->GetMinimumEpsilonStep()<<G4endl;
    G4cout<<"fEpsilonMax is "<< fFieldManager->GetMaximumEpsilonStep()<<G4endl;
    return;
}

void BesMagneticField::GetFieldValue	(	const double	Point[3],
		double *	Bfield
	)			const

Definition at line 103 of file BesMagneticField.cc.

References IMagneticFieldSvc::fieldVector(), ReadBoostRoot::GetField(), m_pIMF, IMagneticFieldSvc::uniFieldVector(), and x.

{
        double x=Point[0];
        double y=Point[1];
        double z=Point[2];

        HepPoint3D r(x,y,z);
        HepVector3D b;

        if(ReadBoostRoot::GetField()==2)
          m_pIMF->fieldVector(r,b);
        else
          m_pIMF->uniFieldVector(r,b);

        Bfield[0]=b.x();
        Bfield[1]=b.y();
        Bfield[2]=b.z();
        
//caogf debug   
//      ofstream haha("field_out.dat",ios_base::app);
//      haha<<x/mm<<" "<<y/mm<<" "<<z/mm<<" "<<Bfield[0]/tesla<<" "<<Bfield[1]/tesla<<" "<<Bfield[2]/tesla<<G4endl;
//      haha.close();
} 

void BesMagneticField::InitialiseAll

(

)

Definition at line 130 of file BesMagneticField.cc.

References CreateStepperAndChordFinder(), fEquation, fFieldManager, fFieldMessenger, fMinStep, and fStepperType.

Referenced by BesMagneticField().

{
    
  fFieldMessenger=new BesMagneticFieldMessenger(this); 
  fEquation = new G4Mag_UsualEqRhs(this); 
 
  fMinStep     = 0.01*mm ; // minimal step of 1 mm is default

  fStepperType =4;      // ClassicalRK4 is default stepper
  fFieldManager = G4TransportationManager::GetTransportationManager()
                                         ->GetFieldManager();
  G4cout<<"before CreateStepperAndChordFinder"<<G4endl;
  CreateStepperAndChordFinder();
}

void BesMagneticField::SetDeltaIntersection

(

double

newvalue

)

Definition at line 243 of file BesMagneticField.cc.

References fFieldManager.

Referenced by BesMagneticFieldMessenger::SetNewValue().

{
        fFieldManager = G4TransportationManager::GetTransportationManager()->GetFieldManager();
        fFieldManager->SetDeltaIntersection(newvalue);
}

void BesMagneticField::SetDeltaOneStep

(

double

newvalue

)

Definition at line 237 of file BesMagneticField.cc.

References fFieldManager.

Referenced by BesMagneticFieldMessenger::SetNewValue().

{
        fFieldManager = G4TransportationManager::GetTransportationManager()
                                                         ->GetFieldManager();
        fFieldManager->SetDeltaOneStep(newvalue);
}

void BesMagneticField::SetMaximumEpsilonStep

(

double

newvalue

)

Definition at line 253 of file BesMagneticField.cc.

References fFieldManager.

Referenced by BesMagneticFieldMessenger::SetNewValue().

{
        fFieldManager =G4TransportationManager::GetTransportationManager()->GetFieldManager();
        fFieldManager->SetMaximumEpsilonStep(newvalue);
}

void BesMagneticField::SetMinimumEpsilonStep

(

double

newvalue

)

Definition at line 248 of file BesMagneticField.cc.

References fFieldManager.

Referenced by BesMagneticFieldMessenger::SetNewValue().

{
        fFieldManager = G4TransportationManager::GetTransportationManager()->GetFieldManager();
        fFieldManager->SetMinimumEpsilonStep(newvalue);
}

void BesMagneticField::SetMinStep

(

G4double

s

)

[inline]

Definition at line 71 of file BesMagneticField.hh.

References fMinStep.

Referenced by BesMagneticFieldMessenger::SetNewValue().

{ fMinStep = s ; }

void BesMagneticField::SetStepper

(

)

Definition at line 185 of file BesMagneticField.cc.

References fEquation, fStepper, and fStepperType.

Referenced by CreateStepperAndChordFinder().

{
  if(fStepper) delete fStepper;

  switch ( fStepperType ) 
  {
    case 0:  
      fStepper = new G4ExplicitEuler( fEquation ); 
      G4cout<<"G4ExplicitEuler is called"<<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;
  }
  return; 
}

void BesMagneticField::SetStepperType

(

G4int

i

)

[inline]

Definition at line 67 of file BesMagneticField.hh.

References fStepperType.

Referenced by BesMagneticFieldMessenger::SetNewValue().

{ fStepperType = i ; }

---

Member Data Documentation

G4ChordFinder* BesMagneticField::fChordFinder [protected]

Definition at line 85 of file BesMagneticField.hh.

Referenced by CreateStepperAndChordFinder(), and ~BesMagneticField().

G4Mag_UsualEqRhs* BesMagneticField::fEquation [protected]

Definition at line 86 of file BesMagneticField.hh.

Referenced by InitialiseAll(), SetStepper(), and ~BesMagneticField().

G4FieldManager* BesMagneticField::fFieldManager [protected]

Definition at line 84 of file BesMagneticField.hh.

Referenced by CreateStepperAndChordFinder(), InitialiseAll(), SetDeltaIntersection(), SetDeltaOneStep(), SetMaximumEpsilonStep(), and SetMinimumEpsilonStep().

BesMagneticFieldMessenger* BesMagneticField::fFieldMessenger [protected]

Definition at line 94 of file BesMagneticField.hh.

Referenced by InitialiseAll(), and ~BesMagneticField().

G4double BesMagneticField::fMinStep [protected]

Definition at line 92 of file BesMagneticField.hh.

Referenced by CreateStepperAndChordFinder(), InitialiseAll(), and SetMinStep().

G4MagIntegratorStepper* BesMagneticField::fStepper [protected]

Definition at line 89 of file BesMagneticField.hh.

Referenced by CreateStepperAndChordFinder(), SetStepper(), and ~BesMagneticField().

G4int BesMagneticField::fStepperType [protected]

Definition at line 90 of file BesMagneticField.hh.

Referenced by InitialiseAll(), SetStepper(), and SetStepperType().

IMagneticFieldSvc* BesMagneticField::m_pIMF [protected]

Definition at line 95 of file BesMagneticField.hh.

Referenced by BesMagneticField(), and GetFieldValue().

---