G4ChordFinderSaf Class Reference

#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	SetChargeMomentumMass (G4double pCharge, G4double pMomentum, G4double pMass)
void	SetIntegrationDriver (G4MagInt_Driver *IntegrationDriver)
G4MagInt_Driver *	GetIntegrationDriver ()
void	ResetStepEstimate ()
G4int	GetNoCalls ()
G4int	GetNoTrials ()
G4int	GetNoMaxTrials ()
virtual void	PrintStatistics ()
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)
virtual G4double	FindNextChord (const G4FieldTrack &yStart, G4double stepMax, G4FieldTrack &yEnd, G4double &dyErr, G4double epsStep, G4double *pNextStepForAccuracy, const G4ThreeVector latestSafetyOrigin, G4double latestSafetyRadius)
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() [1/2]

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

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::~G4ChordFinderSaf

(

)

Definition at line 57 of file G4ChordFinderSaf.cc.

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

Member Function Documentation

FindNextChord()

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.

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

PrintStatistics()

void G4ChordFinderSaf::PrintStatistics ( )

virtual

Reimplemented from G4ChordFinder.

Definition at line 64 of file G4ChordFinderSaf.cc.

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

Referenced by ~G4ChordFinderSaf().

---

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

• G4ChordFinderSaf.hh
• G4ChordFinderSaf.cc