BesMagneticField Class Reference

#include <BesMagneticField.hh>

Inheritance diagram for BesMagneticField:

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::BesMagneticField

(

)

Definition at line 86 of file BesMagneticField.cc.

  : 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::~BesMagneticField

(

)

Definition at line 97 of file BesMagneticField.cc.

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

Member Function Documentation

CreateStepperAndChordFinder()

void BesMagneticField::CreateStepperAndChordFinder

(

)

Definition at line 152 of file BesMagneticField.cc.

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

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

GetFieldValue()

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

Definition at line 106 of file BesMagneticField.cc.

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

InitialiseAll()

void BesMagneticField::InitialiseAll

(

)

Definition at line 133 of file BesMagneticField.cc.

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

Referenced by BesMagneticField().

SetDeltaIntersection()

void BesMagneticField::SetDeltaIntersection

(

double

newvalue

)

Definition at line 246 of file BesMagneticField.cc.

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

Referenced by BesMagneticFieldMessenger::SetNewValue().

SetDeltaOneStep()

void BesMagneticField::SetDeltaOneStep

(

double

newvalue

)

Definition at line 240 of file BesMagneticField.cc.

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

Referenced by BesMagneticFieldMessenger::SetNewValue().

SetMaximumEpsilonStep()

void BesMagneticField::SetMaximumEpsilonStep

(

double

newvalue

)

Definition at line 256 of file BesMagneticField.cc.

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

Referenced by BesMagneticFieldMessenger::SetNewValue().

SetMinimumEpsilonStep()

void BesMagneticField::SetMinimumEpsilonStep

(

double

newvalue

)

Definition at line 251 of file BesMagneticField.cc.

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

Referenced by BesMagneticFieldMessenger::SetNewValue().

SetMinStep()

void BesMagneticField::SetMinStep ( G4double s )

inline

Definition at line 71 of file BesMagneticField.hh.

{ fMinStep = s ; }

Referenced by BesMagneticFieldMessenger::SetNewValue().

SetStepper()

void BesMagneticField::SetStepper

(

)

Definition at line 188 of file BesMagneticField.cc.

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

Referenced by CreateStepperAndChordFinder().

SetStepperType()

void BesMagneticField::SetStepperType ( G4int i )

inline

Definition at line 67 of file BesMagneticField.hh.

{ fStepperType = i ; }

Referenced by BesMagneticFieldMessenger::SetNewValue().

Member Data Documentation

fChordFinder

G4ChordFinder* BesMagneticField::fChordFinder

protected

Definition at line 85 of file BesMagneticField.hh.

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

fEquation

G4Mag_UsualEqRhs* BesMagneticField::fEquation

protected

Definition at line 86 of file BesMagneticField.hh.

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

fFieldManager

G4FieldManager* BesMagneticField::fFieldManager

protected

Definition at line 84 of file BesMagneticField.hh.

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

fFieldMessenger

BesMagneticFieldMessenger* BesMagneticField::fFieldMessenger

protected

Definition at line 94 of file BesMagneticField.hh.

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

fMinStep

G4double BesMagneticField::fMinStep

protected

Definition at line 92 of file BesMagneticField.hh.

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

fStepper

G4MagIntegratorStepper* BesMagneticField::fStepper

protected

Definition at line 89 of file BesMagneticField.hh.

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

fStepperType

G4int BesMagneticField::fStepperType

protected

Definition at line 90 of file BesMagneticField.hh.

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

m_pIMF

IMagneticFieldSvc* BesMagneticField::m_pIMF

protected

Definition at line 95 of file BesMagneticField.hh.

Referenced by BesMagneticField(), and GetFieldValue().

---

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

• BesMagneticField.hh
• BesMagneticField.cc