#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)

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/2]

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() [1/2]

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

◆ BesMagneticField() [2/2]

BesMagneticField::BesMagneticField ( )

◆ ~BesMagneticField() [2/2]

BesMagneticField::~BesMagneticField ( )

Member Function Documentation

◆ CreateStepperAndChordFinder() [1/2]

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().

◆ CreateStepperAndChordFinder() [2/2]

void BesMagneticField::CreateStepperAndChordFinder ( )

◆ GetFieldValue() [1/2]

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

◆ GetFieldValue() [2/2]

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

◆ InitialiseAll() [1/2]

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().

◆ InitialiseAll() [2/2]

void BesMagneticField::InitialiseAll ( )

◆ SetDeltaIntersection() [1/2]

void BesMagneticField::SetDeltaIntersection ( double newvalue )

Definition at line 243 of file BesMagneticField.cc.

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

Referenced by BesMagneticFieldMessenger::SetNewValue().

◆ SetDeltaIntersection() [2/2]

void BesMagneticField::SetDeltaIntersection ( double newvalue )

◆ SetDeltaOneStep() [1/2]

void BesMagneticField::SetDeltaOneStep ( double newvalue )

Definition at line 237 of file BesMagneticField.cc.

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

Referenced by BesMagneticFieldMessenger::SetNewValue().

◆ SetDeltaOneStep() [2/2]

void BesMagneticField::SetDeltaOneStep ( double newvalue )

◆ SetMaximumEpsilonStep() [1/2]

void BesMagneticField::SetMaximumEpsilonStep ( double newvalue )

Definition at line 253 of file BesMagneticField.cc.

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

Referenced by BesMagneticFieldMessenger::SetNewValue().

◆ SetMaximumEpsilonStep() [2/2]

void BesMagneticField::SetMaximumEpsilonStep ( double newvalue )

◆ SetMinimumEpsilonStep() [1/2]

void BesMagneticField::SetMinimumEpsilonStep ( double newvalue )

Definition at line 248 of file BesMagneticField.cc.

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

Referenced by BesMagneticFieldMessenger::SetNewValue().

◆ SetMinimumEpsilonStep() [2/2]

void BesMagneticField::SetMinimumEpsilonStep ( double newvalue )

◆ SetMinStep() [1/2]

void BesMagneticField::SetMinStep ( G4double s )

inline

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

71{ fMinStep = s ; }

s

XmlRpcServer s

Definition: HelloServer.cpp:11

Referenced by BesMagneticFieldMessenger::SetNewValue().

◆ SetMinStep() [2/2]

void BesMagneticField::SetMinStep ( G4double s )

inline

Definition at line 71 of file Simulation/BOOST/BesSim/BesSim-00-01-24/BesSim/BesMagneticField.hh.

71{ fMinStep = s ; }

◆ SetStepper() [1/2]

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