#include <G4ChordFinderSaf.hh>

Inheritance diagram for G4ChordFinderSaf:

Public Member Functions
	G4ChordFinderSaf (G4MagInt_Driver *pIntegrationDriver)

	G4ChordFinderSaf (G4MagneticField theMagField, G4double stepMinimum, G4MagIntegratorStepper pItsStepper)

	~G4ChordFinderSaf ()

G4double	FindNextChord (const G4FieldTrack &yStart, G4double stepMax, G4FieldTrack &yEnd, G4double &dyErrPos, G4double epsStep, G4double *pStepForAccuracy, const G4ThreeVector latestSafetyOrigin, G4double latestSafetyRadius)

void	PrintStatistics ()

Public Member Functions inherited from G4ChordFinder
	G4ChordFinder (G4MagInt_Driver *pIntegrationDriver)

	G4ChordFinder (G4MagneticField itsMagField, G4double stepMinimum=1.0e-2, G4MagIntegratorStepper pItsStepper=0)

virtual	~G4ChordFinder ()

G4double	AdvanceChordLimited (G4FieldTrack &yCurrent, G4double stepInitial, G4double epsStep_Relative, const G4ThreeVector latestSafetyOrigin, G4double lasestSafetyRadius)

G4FieldTrack	ApproxCurvePointS (const G4FieldTrack &curveAPointVelocity, const G4FieldTrack &curveBPointVelocity, const G4FieldTrack &ApproxCurveV, const G4ThreeVector &currentEPoint, const G4ThreeVector &currentFPoint, const G4ThreeVector &PointG, G4bool first, G4double epsStep)

G4FieldTrack	ApproxCurvePointV (const G4FieldTrack &curveAPointVelocity, const G4FieldTrack &curveBPointVelocity, const G4ThreeVector &currentEPoint, G4double epsStep)

G4double	InvParabolic (const G4double xa, const G4double ya, const G4double xb, const G4double yb, const G4double xc, const G4double yc)

G4double	GetDeltaChord () const

void	SetDeltaChord (G4double newval)

void	SetIntegrationDriver (G4MagInt_Driver *IntegrationDriver)

G4MagInt_Driver *	GetIntegrationDriver ()

void	ResetStepEstimate ()

G4int	GetNoCalls ()

G4int	GetNoTrials ()

G4int	GetNoMaxTrials ()

G4int	SetVerbose (G4int newvalue=1)

void	SetFractions_Last_Next (G4double fractLast=0.90, G4double fractNext=0.95)

void	SetFirstFraction (G4double fractFirst)

void	TestChordPrint (G4int noTrials, G4int lastStepTrial, G4double dChordStep, G4double nextStepTrial)

G4double	GetFirstFraction ()

G4double	GetFractionLast ()

G4double	GetFractionNextEstimate ()

G4double	GetMultipleRadius ()

Additional Inherited Members
Protected Member Functions inherited from G4ChordFinder
void	AccumulateStatistics (G4int noTrials)

G4bool	AcceptableMissDist (G4double dChordStep) const

G4double	NewStep (G4double stepTrialOld, G4double dChordStep, G4double &stepEstimate_Unconstrained)

void	PrintDchordTrial (G4int noTrials, G4double stepTrial, G4double oldStepTrial, G4double dChordStep)

G4double	GetLastStepEstimateUnc ()

void	SetLastStepEstimateUnc (G4double stepEst)

Detailed Description

Definition at line 42 of file G4ChordFinderSaf.hh.

Constructor & Destructor Documentation

G4ChordFinderSaf::G4ChordFinderSaf ( G4MagInt_Driver * pIntegrationDriver )

Definition at line 33 of file G4ChordFinderSaf.cc.

   : G4ChordFinder(pIntegrationDriver)
 {
     // check the values and set the other parameters
     // fNoInitialRadBig=0;    fNoInitialRadSmall=0;  
     // fNoTrialsRadBig=0;     fNoTrialsRadSmall=0; 
    
 }

G4ChordFinderSaf::G4ChordFinderSaf	(	G4MagneticField *	theMagField,
		G4double	stepMinimum,
		G4MagIntegratorStepper *	pItsStepper
	)

Definition at line 44 of file G4ChordFinderSaf.cc.

   : G4ChordFinder( theMagField, stepMinimum, pItsStepper )
 {
   //  Let  G4ChordFinder create the  Driver, the Stepper and EqRhs ...
   // ...
 }

G4ChordFinderSaf::~G4ChordFinderSaf ( )

Definition at line 57 of file G4ChordFinderSaf.cc.

References PrintStatistics(), and G4ChordFinder::SetVerbose().

 {
   if( SetVerbose(0) ) { PrintStatistics();  } 
    // Set verbosity 0, so that will be called in base class again
 }

Member Function Documentation

G4double G4ChordFinderSaf::FindNextChord	(	const G4FieldTrack &	yStart,
		G4double	stepMax,
		G4FieldTrack &	yEnd,
		G4double &	dyErrPos,
		G4double	epsStep,
		G4double *	pStepForAccuracy,
		const G4ThreeVector	latestSafetyOrigin,
		G4double	latestSafetyRadius
	)

virtual

Reimplemented from G4ChordFinder.

Definition at line 118 of file G4ChordFinderSaf.cc.

References G4ChordFinder::AcceptableMissDist(), G4ChordFinder::AccumulateStatistics(), CalculatePointInside(), CalculatePointSafety(), G4MagInt_Driver::ComputeNewStepSize(), G4cout, G4endl, G4ChordFinder::GetFirstFraction(), G4ChordFinder::GetFractionLast(), G4ChordFinder::GetIntegrationDriver(), G4ChordFinder::GetLastStepEstimateUnc(), G4FieldTrack::GetPosition(), G4INCL::Math::max(), G4INCL::Math::min(), G4FieldTrack::ncompSVEC, G4ChordFinder::NewStep(), G4MagInt_Driver::QuickAdvance(), and G4ChordFinder::SetLastStepEstimateUnc().

 {
   // G4int       stepRKnumber=0;
   G4FieldTrack yCurrent=  yStart;  
   G4double    stepTrial, stepForAccuracy;
   G4double    dydx[G4FieldTrack::ncompSVEC]; 
 
   //  1.)  Try to "leap" to end of interval
   //  2.)  Evaluate if resulting chord gives d_chord that is good enough.
   //     2a.)  If d_chord is not good enough, find one that is.
   
   G4bool    validEndPoint= false;
   G4double  dChordStep, lastStepLength; //  stepOfLastGoodChord;
 
   GetIntegrationDriver()-> GetDerivatives( yCurrent, dydx )  ;
 
   G4int     noTrials=0;
   const G4double safetyFactor= GetFirstFraction(); // was 0.999
 
   // Figure out the starting safety
   G4double      startSafety=  
     CalculatePointSafety( latestSafetyOrigin, 
               latestSafetyRadius, 
               yCurrent.GetPosition() );
 
   G4double
     likelyGood = std::max( startSafety , 
                         safetyFactor *  GetLastStepEstimateUnc() );
 
   stepTrial  = std::min( stepMax,  likelyGood ); 
 
   G4MagInt_Driver *pIntgrDriver= G4ChordFinder::GetIntegrationDriver(); 
   G4double newStepEst_Uncons= 0.0; 
   do
   { 
      G4double stepForChord;  
      yCurrent = yStart;    // Always start from initial point
 
    //            ************
      pIntgrDriver->QuickAdvance( yCurrent, dydx, stepTrial, 
                                  dChordStep, dyErrPos);
      //            ************
 
      G4ThreeVector EndPointCand= yCurrent.GetPosition(); 
      G4bool  endPointInSafetySphere= 
        CalculatePointInside(latestSafetyOrigin, latestSafetyRadius, EndPointCand);
 
      // We check whether the criterion is met here.
      validEndPoint = AcceptableMissDist(dChordStep)
                        || endPointInSafetySphere;
                       //  && (dyErrPos < eps) ;
 
      lastStepLength = stepTrial; 
 
      // This method estimates to step size for a good chord.
      stepForChord = NewStep(stepTrial, dChordStep, newStepEst_Uncons );
 
      if( ! validEndPoint ) {
         if( stepTrial<=0.0 )
       stepTrial = stepForChord; 
         else if (stepForChord <= stepTrial) 
           // Reduce by a fraction, possibly up to 20% 
           stepTrial = std::min( stepForChord, 
                                 GetFractionLast() * stepTrial); 
         else
           stepTrial *= 0.1;
 
         // if(dbg) G4cerr<<"Dchord too big. Try new hstep="<<stepTrial<<G4endl;
      }
      // #ifdef  TEST_CHORD_PRINT
      // TestChordPrint( noTrials, lastStepLength, dChordStep, stepTrial );
      // #endif
 
      noTrials++; 
   }
   while( ! validEndPoint );   // End of do-while  RKD 
 
   AccumulateStatistics( noTrials );
 
   // Should we update newStepEst_Uncons for a 'long step' via safety ??
   if( newStepEst_Uncons > 0.0  ){ 
     SetLastStepEstimateUnc( newStepEst_Uncons );
   }
 
   // stepOfLastGoodChord = stepTrial;
 
   if( pStepForAccuracy ){ 
      // Calculate the step size required for accuracy, if it is needed
      G4double dyErr_relative = dyErrPos/(epsStep*lastStepLength);
      if( dyErr_relative > 1.0 ) {
         stepForAccuracy =
            pIntgrDriver->ComputeNewStepSize( dyErr_relative,
                                              lastStepLength );
      }else{
         stepForAccuracy = 0.0;   // Convention to show step was ok 
      }
      *pStepForAccuracy = stepForAccuracy;
   }
 
 #ifdef  TEST_CHORD_PRINT
   static int dbg=0;
   if( dbg ) 
     G4cout << "ChordF/FindNextChord:  NoTrials= " << noTrials 
            << " StepForGoodChord=" << std::setw(10) << stepTrial << G4endl;
 #endif
 
   yEnd=  yCurrent;  
   return stepTrial; 
 }

void G4ChordFinderSaf::PrintStatistics ( )

virtual

Reimplemented from G4ChordFinder.

Definition at line 64 of file G4ChordFinderSaf.cc.

References G4cout, G4endl, and G4ChordFinder::PrintStatistics().

Referenced by ~G4ChordFinderSaf().

 {
   // Print Statistics
   G4cout << "G4ChordFinderSaf statistics report: " << G4endl;
   G4ChordFinder::PrintStatistics();
 
 /*******************
   G4cout 
     << "  No radbig calls " << std::setw(10) << fNoInitialRadBig 
     << " trials " << std::setw(10) << fNoTrialsRadBig
     << "  - over " << std::setw(10) << fNoTrialsRadBig - fNoInitialRadBig
     << G4endl
     << "  No radsm  calls " << std::setw(10) << fNoInitialRadSmall 
     << " trials "  << std::setw(10) << fNoTrialsRadSmall
     << "  - over " << std::setw(10) << fNoTrialsRadSmall - fNoInitialRadSmall
     << G4endl;
   G4cout
     << "  *** Limiting stepTrial via if Delta_chord < R_curvature " 
     << "   for large to angle from Delta_chord / R_curv "  
     << "   and for small with multiple " << GetMultipleRadius()
     << G4endl; 
 ********************/
 }

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

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation