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)
	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)
	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 InstallArea/include/BesSim/BesSim/BesMagneticField.hh.

Constructor & Destructor Documentation

BesMagneticField() [1/3]

BesMagneticField::BesMagneticField

(

)

Definition at line 83 of file bak_BesSim-00-04-14/src/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() [1/3]

BesMagneticField::~BesMagneticField

(

)

Definition at line 94 of file bak_BesSim-00-04-14/src/BesMagneticField.cc.

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

BesMagneticField() [2/3]

BesMagneticField::BesMagneticField

(

)

~BesMagneticField() [2/3]

BesMagneticField::~BesMagneticField

(

)

BesMagneticField() [3/3]

BesMagneticField::BesMagneticField

(

)

~BesMagneticField() [3/3]

BesMagneticField::~BesMagneticField

(

)

Member Function Documentation

CreateStepperAndChordFinder() [1/3]

void BesMagneticField::CreateStepperAndChordFinder

(

)

Definition at line 149 of file bak_BesSim-00-04-14/src/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().

CreateStepperAndChordFinder() [2/3]

void BesMagneticField::CreateStepperAndChordFinder

(

)

CreateStepperAndChordFinder() [3/3]

void BesMagneticField::CreateStepperAndChordFinder

(

)

GetFieldValue() [1/3]

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

Definition at line 103 of file bak_BesSim-00-04-14/src/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();
} 

GetFieldValue() [2/3]

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

GetFieldValue() [3/3]

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

InitialiseAll() [1/3]

void BesMagneticField::InitialiseAll

(

)

Definition at line 130 of file bak_BesSim-00-04-14/src/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().

InitialiseAll() [2/3]

void BesMagneticField::InitialiseAll

(

)

InitialiseAll() [3/3]

void BesMagneticField::InitialiseAll

(

)

SetDeltaIntersection() [1/3]

void BesMagneticField::SetDeltaIntersection

(

double

newvalue

)

Definition at line 243 of file bak_BesSim-00-04-14/src/BesMagneticField.cc.

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

Referenced by BesMagneticFieldMessenger::SetNewValue().

SetDeltaIntersection() [2/3]

void BesMagneticField::SetDeltaIntersection

(

double

newvalue

)

SetDeltaIntersection() [3/3]

void BesMagneticField::SetDeltaIntersection

(

double

newvalue

)

SetDeltaOneStep() [1/3]

void BesMagneticField::SetDeltaOneStep

(

double

newvalue

)

Definition at line 237 of file bak_BesSim-00-04-14/src/BesMagneticField.cc.

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

Referenced by BesMagneticFieldMessenger::SetNewValue().

SetDeltaOneStep() [2/3]

void BesMagneticField::SetDeltaOneStep

(

double

newvalue

)

SetDeltaOneStep() [3/3]

void BesMagneticField::SetDeltaOneStep

(

double

newvalue

)

SetMaximumEpsilonStep() [1/3]

void BesMagneticField::SetMaximumEpsilonStep

(

double

newvalue

)

Definition at line 253 of file bak_BesSim-00-04-14/src/BesMagneticField.cc.

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

Referenced by BesMagneticFieldMessenger::SetNewValue().

SetMaximumEpsilonStep() [2/3]

void BesMagneticField::SetMaximumEpsilonStep

(

double

newvalue

)

SetMaximumEpsilonStep() [3/3]

void BesMagneticField::SetMaximumEpsilonStep

(

double

newvalue

)

SetMinimumEpsilonStep() [1/3]

void BesMagneticField::SetMinimumEpsilonStep

(

double

newvalue

)

Definition at line 248 of file bak_BesSim-00-04-14/src/BesMagneticField.cc.

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

Referenced by BesMagneticFieldMessenger::SetNewValue().

SetMinimumEpsilonStep() [2/3]

void BesMagneticField::SetMinimumEpsilonStep

(

double

newvalue

)

SetMinimumEpsilonStep() [3/3]

void BesMagneticField::SetMinimumEpsilonStep

(

double

newvalue

)

SetMinStep() [1/3]

void BesMagneticField::SetMinStep ( G4double  s )

inline

Definition at line 71 of file InstallArea/include/BesSim/BesSim/BesMagneticField.hh.

{ fMinStep = s ; }

Referenced by BesMagneticFieldMessenger::SetNewValue().

SetMinStep() [2/3]

void BesMagneticField::SetMinStep ( G4double  s )

inline

Definition at line 71 of file Simulation/BOOST/BesSim/bak_BesSim-00-04-14/BesSim/BesMagneticField.hh.

{ fMinStep = s ; }

SetMinStep() [3/3]

void BesMagneticField::SetMinStep ( G4double  s )

inline

Definition at line 71 of file Simulation/BOOST/BesSim/BesSim-00-04-16/BesSim/BesMagneticField.hh.

{ fMinStep = s ; }

SetStepper() [1/3]

void BesMagneticField::SetStepper

(

)

Definition at line 185 of file bak_BesSim-00-04-14/src/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().

SetStepper() [2/3]

void BesMagneticField::SetStepper

(

)

SetStepper() [3/3]

void BesMagneticField::SetStepper

(

)

SetStepperType() [1/3]

void BesMagneticField::SetStepperType ( G4int  i )

inline

Definition at line 67 of file InstallArea/include/BesSim/BesSim/BesMagneticField.hh.

{ fStepperType = i ; }

Referenced by BesMagneticFieldMessenger::SetNewValue().

SetStepperType() [2/3]

void BesMagneticField::SetStepperType ( G4int  i )

inline

Definition at line 67 of file Simulation/BOOST/BesSim/bak_BesSim-00-04-14/BesSim/BesMagneticField.hh.

{ fStepperType = i ; }

SetStepperType() [3/3]

void BesMagneticField::SetStepperType ( G4int  i )

inline

Definition at line 67 of file Simulation/BOOST/BesSim/BesSim-00-04-16/BesSim/BesMagneticField.hh.

{ fStepperType = i ; }

Member Data Documentation

fChordFinder

G4ChordFinder * BesMagneticField::fChordFinder

protected

Definition at line 85 of file InstallArea/include/BesSim/BesSim/BesMagneticField.hh.

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

fEquation

G4Mag_UsualEqRhs * BesMagneticField::fEquation

protected

Definition at line 86 of file InstallArea/include/BesSim/BesSim/BesMagneticField.hh.

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

fFieldManager

G4FieldManager * BesMagneticField::fFieldManager

protected

Definition at line 84 of file InstallArea/include/BesSim/BesSim/BesMagneticField.hh.

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

fFieldMessenger

BesMagneticFieldMessenger * BesMagneticField::fFieldMessenger

protected

Definition at line 94 of file InstallArea/include/BesSim/BesSim/BesMagneticField.hh.

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

fMinStep

G4double BesMagneticField::fMinStep

protected

Definition at line 92 of file InstallArea/include/BesSim/BesSim/BesMagneticField.hh.

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

fStepper

G4MagIntegratorStepper * BesMagneticField::fStepper

protected

Definition at line 89 of file InstallArea/include/BesSim/BesSim/BesMagneticField.hh.

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

fStepperType

G4int BesMagneticField::fStepperType

protected

Definition at line 90 of file InstallArea/include/BesSim/BesSim/BesMagneticField.hh.

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

m_pIMF

IMagneticFieldSvc * BesMagneticField::m_pIMF

protected

Definition at line 95 of file InstallArea/include/BesSim/BesSim/BesMagneticField.hh.

Referenced by BesMagneticField(), and GetFieldValue().

---

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

• InstallArea/include/BesSim/BesSim/BesMagneticField.hh
• Simulation/BOOST/BesSim/bak_BesSim-00-04-14/BesSim/BesMagneticField.hh
• Simulation/BOOST/BesSim/BesSim-00-04-16/BesSim/BesMagneticField.hh
• bak_BesSim-00-04-14/src/BesMagneticField.cc
• BesSim-00-04-16/src/BesMagneticField.cc