G4HelixMixedStepper Class Reference

#include <G4HelixMixedStepper.hh>

Inheritance diagram for G4HelixMixedStepper:

Public Member Functions
	G4HelixMixedStepper (G4Mag_EqRhs *EqRhs, G4int fStepperNumber=0)
	~G4HelixMixedStepper ()
void	Stepper (const G4double y[], const G4double dydx[], G4double h, G4double yout[], G4double yerr[])
void	DumbStepper (const G4double y[], G4ThreeVector Bfld, G4double h, G4double yout[])
G4double	DistChord () const
void	SetVerbose (G4int newvalue)
void	PrintCalls ()
G4MagIntegratorStepper *	SetupStepper (G4Mag_EqRhs *EqRhs, G4int StepperName)
G4int	IntegratorOrder () const
Public Member Functions inherited from G4MagHelicalStepper
	G4MagHelicalStepper (G4Mag_EqRhs *EqRhs)
virtual	~G4MagHelicalStepper ()
virtual void	Stepper (const G4double y[], const G4double dydx[], G4double h, G4double yout[], G4double yerr[])
virtual void	DumbStepper (const G4double y[], G4ThreeVector Bfld, G4double h, G4double yout[])=0
G4double	DistChord () const
Public Member Functions inherited from G4MagIntegratorStepper
	G4MagIntegratorStepper (G4EquationOfMotion *Equation, G4int numIntegrationVariables, G4int numStateVariables=12)
virtual	~G4MagIntegratorStepper ()
virtual void	Stepper (const G4double y[], const G4double dydx[], G4double h, G4double yout[], G4double yerr[])=0
virtual G4double	DistChord () const =0
virtual void	ComputeRightHandSide (const G4double y[], G4double dydx[])
void	NormaliseTangentVector (G4double vec[6])
void	NormalisePolarizationVector (G4double vec[12])
void	RightHandSide (const double y[], double dydx[])
G4int	GetNumberOfVariables () const
G4int	GetNumberOfStateVariables () const
virtual G4int	IntegratorOrder () const =0
G4EquationOfMotion *	GetEquationOfMotion ()
void	SetEquationOfMotion (G4EquationOfMotion *newEquation)

Additional Inherited Members
Protected Member Functions inherited from G4MagHelicalStepper
void	LinearStep (const G4double yIn[], G4double h, G4double yHelix[]) const
void	AdvanceHelix (const G4double yIn[], G4ThreeVector Bfld, G4double h, G4double yHelix[], G4double yHelix2[]=0)
void	MagFieldEvaluate (const G4double y[], G4ThreeVector &Bfield)
G4double	GetInverseCurve (const G4double Momentum, const G4double Bmag)
void	SetAngCurve (const G4double Ang)
G4double	GetAngCurve () const
void	SetCurve (const G4double Curve)
G4double	GetCurve () const
void	SetRadHelix (const G4double Rad)
G4double	GetRadHelix () const

Detailed Description

Definition at line 58 of file G4HelixMixedStepper.hh.

Constructor & Destructor Documentation

G4HelixMixedStepper()

G4HelixMixedStepper::G4HelixMixedStepper	(	G4Mag_EqRhs *	EqRhs,
		G4int	fStepperNumber = 0
	)

Definition at line 57 of file G4HelixMixedStepper.cc.

  : G4MagHelicalStepper(EqRhs)
    
{
   SetVerbose(1); fNumCallsRK4=0; fNumCallsHelix=0;
   if(!fStepperNumber) fStepperNumber=4;
   fRK4Stepper =  SetupStepper(EqRhs, fStepperNumber);
}

~G4HelixMixedStepper()

G4HelixMixedStepper::~G4HelixMixedStepper

(

)

Definition at line 67 of file G4HelixMixedStepper.cc.

                                          {
     
     delete(fRK4Stepper);
     if (fVerbose>0){ PrintCalls();};
} 

Member Function Documentation

DistChord()

G4double G4HelixMixedStepper::DistChord ( ) const

virtual

Implements G4MagIntegratorStepper.

Definition at line 145 of file G4HelixMixedStepper.cc.

{
  // Implementation : must check whether h/R > 2 pi  !!
  //   If( h/R <  pi) use G4LineSection::DistLine
  //   Else           DistChord=R_helix
  //
  G4double distChord;
  G4double Ang_curve=GetAngCurve();
 
      
     if(Ang_curve<=pi){
       distChord=GetRadHelix()*(1-std::cos(0.5*Ang_curve));
     }
         else 
         if(Ang_curve<twopi){
           distChord=GetRadHelix()*(1+std::cos(0.5*(twopi-Ang_curve)));
         }
         else{
          distChord=2.*GetRadHelix();  
         }
 
   
 
  return distChord;
  
}

DumbStepper()

void G4HelixMixedStepper::DumbStepper ( const G4double  y[], G4ThreeVector  Bfld, G4double  h, G4double  yout[]  )

virtual

Implements G4MagHelicalStepper.

Definition at line 132 of file G4HelixMixedStepper.cc.

{
 
    
       AdvanceHelix(yIn, Bfld, h, yOut);
 
    
               
}  

IntegratorOrder()

G4int G4HelixMixedStepper::IntegratorOrder ( ) const

inlinevirtual

Implements G4MagIntegratorStepper.

Definition at line 93 of file G4HelixMixedStepper.hh.

{ return 4; }

PrintCalls()

void G4HelixMixedStepper::PrintCalls

(

)

Definition at line 172 of file G4HelixMixedStepper.cc.

{
  G4cout<<"In HelixMixedStepper::Number of calls to smallStepStepper = "<<fNumCallsRK4
        <<"  and Number of calls to Helix = "<<fNumCallsHelix<<G4endl;
}

Referenced by ~G4HelixMixedStepper().

SetupStepper()

G4MagIntegratorStepper * G4HelixMixedStepper::SetupStepper	(	G4Mag_EqRhs *	EqRhs,
		G4int	StepperName
	)

Definition at line 180 of file G4HelixMixedStepper.cc.

{
  G4MagIntegratorStepper* pStepper;
  if (fVerbose>0)G4cout<<"In G4HelixMixedStepper Stepper for small steps is "; 
  switch ( StepperNumber )
    {
     case 0: pStepper = new G4ExplicitEuler( pE ); if (fVerbose>0)G4cout<<"G4ExplicitEuler"<<G4endl; break;
     case 1: pStepper = new G4ImplicitEuler( pE ); if (fVerbose>0)G4cout<<"G4ImplicitEuler"<<G4endl; break;
     case 2: pStepper = new G4SimpleRunge( pE ); if (fVerbose>0)G4cout<<"G4SimpleRunge"<<G4endl; break;
     case 3: pStepper = new G4SimpleHeum( pE );  if (fVerbose>0)G4cout<<"G4SimpleHeum"<<G4endl;break;
     case 4: pStepper = new G4ClassicalRK4( pE ); if (fVerbose>0)G4cout<<"G4ClassicalRK4"<<G4endl; break;
     case 5: pStepper = new G4HelixExplicitEuler( pE ); if (fVerbose>0)G4cout<<"G4HelixExplicitEuler"<<G4endl; break;
     case 6: pStepper = new G4HelixImplicitEuler( pE ); if (fVerbose>0)G4cout<<"G4HelixImplicitEuler"<<G4endl; break;
     case 7: pStepper = new G4HelixSimpleRunge( pE ); if (fVerbose>0)G4cout<<"G4HelixSimpleRunge"<<G4endl; break;
     case 8: pStepper = new G4CashKarpRKF45( pE );    if (fVerbose>0)G4cout<<"G4CashKarpRKF45"<<G4endl; break;
     case 9: pStepper = new G4ExactHelixStepper( pE );  if (fVerbose>0)G4cout<<"G4ExactHelixStepper"<<G4endl;   break;
     case 10: pStepper = new G4RKG3_Stepper( pE );  if (fVerbose>0)G4cout<<"G4RKG3_Stepper"<<G4endl;   break;
      
      default: pStepper = new G4ClassicalRK4( pE );G4cout<<"Default G4ClassicalRK4"<<G4endl; break;
      
    }
  return pStepper;
}

Referenced by G4HelixMixedStepper().

SetVerbose()

void G4HelixMixedStepper::SetVerbose ( G4int  newvalue )

inline

Definition at line 86 of file G4HelixMixedStepper.hh.

{fVerbose=newvalue;}

Referenced by G4HelixMixedStepper().

Stepper()

void G4HelixMixedStepper::Stepper ( const G4double  y[], const G4double  dydx[], G4double  h, G4double  yout[], G4double  yerr[]  )

virtual

Reimplemented from G4MagHelicalStepper.

Definition at line 72 of file G4HelixMixedStepper.cc.

{
 
 //Estimation of the Stepping Angle
 
  G4ThreeVector Bfld;
  MagFieldEvaluate(yInput, Bfld); 
 
  G4double Bmag = Bfld.mag();
  const G4double *pIn = yInput+3;
  G4ThreeVector initVelocity= G4ThreeVector( pIn[0], pIn[1], pIn[2]);
  G4double      velocityVal = initVelocity.mag();
  G4double R_1;  
  G4double Ang_curve;
 
   R_1=std::abs(GetInverseCurve(velocityVal,Bmag));
   Ang_curve=R_1*Step;
   SetAngCurve(Ang_curve);
   SetCurve(std::abs(1/R_1));
   
 
   if(Ang_curve<0.33*pi){
     fNumCallsRK4++;   
     fRK4Stepper->Stepper(yInput,dydx,Step,yOut,yErr);
    
 
   }
    else{
      fNumCallsHelix++;
      const G4int nvar = 6 ;
      G4int i;
      G4double      yTemp[7], yIn[7] ;
      G4double yTemp2[7];
      G4ThreeVector  Bfld_midpoint;
    //  Saving yInput because yInput and yOut can be aliases for same array
        for(i=0;i<nvar;i++) yIn[i]=yInput[i];
 
      G4double h = Step * 0.5;
     // Do two half steps and full step
          AdvanceHelix(yIn,   Bfld,  h, yTemp,yTemp2);
          MagFieldEvaluate(yTemp, Bfld_midpoint) ;     
          AdvanceHelix(yTemp, Bfld_midpoint, h, yOut);
     // Error estimation
          for(i=0;i<nvar;i++) {
          yErr[i] = yOut[i] - yTemp2[i] ;
          
          }
    }
 
 
 
 
}

---

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

• G4HelixMixedStepper.hh
• G4HelixMixedStepper.cc