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

---

Detailed Description

Definition at line 58 of file G4HelixMixedStepper.hh.

---

Constructor & Destructor Documentation

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

Definition at line 57 of file G4HelixMixedStepper.cc.

References SetupStepper(), and SetVerbose().

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

G4HelixMixedStepper::~G4HelixMixedStepper

(

)

Definition at line 67 of file G4HelixMixedStepper.cc.

References PrintCalls().

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

---

Member Function Documentation

G4double G4HelixMixedStepper::DistChord

(

)

const [virtual]

Reimplemented from G4MagHelicalStepper.

Definition at line 145 of file G4HelixMixedStepper.cc.

References G4MagHelicalStepper::GetAngCurve(), G4MagHelicalStepper::GetRadHelix(), and G4INCL::Math::pi.

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

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

Implements G4MagHelicalStepper.

Definition at line 132 of file G4HelixMixedStepper.cc.

References G4MagHelicalStepper::AdvanceHelix().

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

    
               
}  

G4int G4HelixMixedStepper::IntegratorOrder

(

)

const [inline, virtual]

Implements G4MagIntegratorStepper.

Definition at line 93 of file G4HelixMixedStepper.hh.

{ return 4; }

void G4HelixMixedStepper::PrintCalls

(

)

Definition at line 172 of file G4HelixMixedStepper.cc.

References G4cout, and G4endl.

Referenced by ~G4HelixMixedStepper().

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

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

Definition at line 180 of file G4HelixMixedStepper.cc.

References G4cout, and G4endl.

Referenced by G4HelixMixedStepper().

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

void G4HelixMixedStepper::SetVerbose

(

G4int

newvalue

)

[inline]

Definition at line 86 of file G4HelixMixedStepper.hh.

Referenced by G4HelixMixedStepper().

{fVerbose=newvalue;}

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

Reimplemented from G4MagHelicalStepper.

---

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

• G4HelixMixedStepper.hh
• G4HelixMixedStepper.cc

---