#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 83 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 94 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 149 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 103 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 130 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 243 of file BesMagneticField.cc.

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

Referenced by BesMagneticFieldMessenger::SetNewValue().

◆ SetDeltaOneStep()

void BesMagneticField::SetDeltaOneStep ( double newvalue )

Definition at line 237 of file BesMagneticField.cc.

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

Referenced by BesMagneticFieldMessenger::SetNewValue().

◆ SetMaximumEpsilonStep()

void BesMagneticField::SetMaximumEpsilonStep ( double newvalue )

Definition at line 253 of file BesMagneticField.cc.

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

Referenced by BesMagneticFieldMessenger::SetNewValue().

◆ SetMinimumEpsilonStep()

void BesMagneticField::SetMinimumEpsilonStep ( double newvalue )

Definition at line 248 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.

71{ fMinStep = s ; }

s

XmlRpcServer s

Definition: HelloServer.cpp:11

Referenced by BesMagneticFieldMessenger::SetNewValue().

◆ SetStepper()

void BesMagneticField::SetStepper ( )

Definition at line 185 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.

67{ 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:

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ BesMagneticField()

◆ ~BesMagneticField()

Member Function Documentation

◆ CreateStepperAndChordFinder()

◆ GetFieldValue()

◆ InitialiseAll()

◆ SetDeltaIntersection()

◆ SetDeltaOneStep()

◆ SetMaximumEpsilonStep()

◆ SetMinimumEpsilonStep()

◆ SetMinStep()

◆ SetStepper()

◆ SetStepperType()

Member Data Documentation

◆ fChordFinder

◆ fEquation

◆ fFieldManager

◆ fFieldMessenger

◆ fMinStep

◆ fStepper

◆ fStepperType

◆ m_pIMF