G4DriverReporter Class Reference

#include <G4DriverReporter.hh>

Static Public Member Functions
static void	PrintStat_Aux (const G4FieldTrack &aFieldTrack, G4double requestStep, G4double actualStep, G4int subStepNo, G4double subStepSize, G4double dotVelocities)
static void	PrintStatus (const G4double *StartArr, G4double xstart, const G4double *CurrentArr, G4double xcurrent, G4double requestStep, unsigned int subStepNo, unsigned int noIntegrationVariables)
static void	PrintStatus (const G4FieldTrack &StartFT, const G4FieldTrack &CurrentFT, G4double requestStep, unsigned int subStepNo)

Detailed Description

Definition at line 41 of file G4DriverReporter.hh.

Member Function Documentation

PrintStat_Aux()

void G4DriverReporter::PrintStat_Aux ( const G4FieldTrack &  aFieldTrack, G4double  requestStep, G4double  actualStep, G4int  subStepNo, G4double  subStepSize, G4double  dotVelocities  )

static

Definition at line 136 of file G4DriverReporter.cc.

{
    const G4ThreeVector Position = aFieldTrack.GetPosition();
    const G4ThreeVector UnitVelocity = aFieldTrack.GetMomentumDir();
 
    G4int oldprec= G4cout.precision(noPrecision);
    
    if( subStepNo >= 0)
    {
       G4cout << std::setw( 5) << subStepNo << " ";
    }
    else
    {
       G4cout << std::setw( 5) << "Start" << " ";
    }
    G4double curveLen= aFieldTrack.GetCurveLength();
    G4cout << std::setw( 7) << curveLen;
    // G4cout.precision(noPrecision);
    G4cout << std::setw( prec9) << Position.x() << " "
           << std::setw( prec9) << Position.y() << " "
           << std::setw( prec9) << Position.z() << " "
           << std::setw( prec8) << UnitVelocity.x() << " "
           << std::setw( prec8) << UnitVelocity.y() << " "
           << std::setw( prec8) << UnitVelocity.z() << " ";
    G4cout.precision(3);
    G4double unitMagDif = UnitVelocity.mag2()-1.0;
    if( std::fabs( unitMagDif ) < 1.0e-15 ) { unitMagDif = 0.0; }        
    G4cout << std::setw( 8) << unitMagDif << " ";
    G4cout.precision(6);
    G4cout << std::setw(10) << dotVeloc_StartCurr << " ";
    G4cout.precision(oldprec);
    G4cout << std::setw( prec7) << aFieldTrack.GetKineticEnergy();
    G4cout << std::setw(12) << step_len << " ";
 
    static G4ThreadLocal G4double oldCurveLength = 0.0;
    static G4ThreadLocal G4double oldSubStepLength = 0.0;
    static G4ThreadLocal G4int oldSubStepNo = -1;
 
    G4double subStep_len = 0.0;
    if( curveLen > oldCurveLength )
    {
      subStep_len= curveLen - oldCurveLength;
    }
    else if (subStepNo == oldSubStepNo)
    {
      subStep_len= oldSubStepLength;
    }
    oldCurveLength= curveLen;
    oldSubStepLength= subStep_len;
 
    G4cout << std::setw(12) << subStep_len << " "; 
    G4cout << std::setw(12) << subStepSize << " "; 
    if( requestStep != -1.0 )
    {
      G4cout << std::setw( prec9) << requestStep << " ";
    }
    else
    {
       G4cout << std::setw( prec9) << " InitialStep " << " ";
    }
    G4cout << G4endl;
}

References G4cout, G4endl, G4ThreadLocal, G4FieldTrack::GetCurveLength(), G4FieldTrack::GetKineticEnergy(), G4FieldTrack::GetMomentumDir(), G4FieldTrack::GetPosition(), CLHEP::Hep3Vector::mag2(), noPrecision, prec7, prec8, prec9, CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

Referenced by PrintStatus().

PrintStatus() [1/2]

void G4DriverReporter::PrintStatus ( const G4double *  StartArr, G4double  xstart, const G4double *  CurrentArr, G4double  xcurrent, G4double  requestStep, unsigned int  subStepNo, unsigned int  noIntegrationVariables  )

static

Definition at line 38 of file G4DriverReporter.cc.

{
   G4FieldTrack  StartFT(G4ThreeVector(0,0,0),
                 G4ThreeVector(0,0,0), 0., 0., 0., 0. );
   G4FieldTrack  CurrentFT (StartFT);
 
   StartFT.LoadFromArray( StartArr, noIntegrationVariables); 
   StartFT.SetCurveLength( xstart);
   CurrentFT.LoadFromArray( CurrentArr, noIntegrationVariables); 
   CurrentFT.SetCurveLength( xcurrent );
 
   PrintStatus(StartFT, CurrentFT, requestStep, subStepNo ); 
}

References G4FieldTrack::LoadFromArray(), PrintStatus(), and G4FieldTrack::SetCurveLength().

Referenced by G4MagInt_Driver::AccurateAdvance(), G4OldMagIntDriver::AccurateAdvance(), and PrintStatus().

PrintStatus() [2/2]

void G4DriverReporter::PrintStatus ( const G4FieldTrack &  StartFT, const G4FieldTrack &  CurrentFT, G4double  requestStep, unsigned int  subStepNo  )

static

Definition at line 69 of file G4DriverReporter.cc.

{
    G4int verboseLevel= 2; // fVerboseLevel;
    G4int oldPrec= G4cout.precision(noPrecision);
    // G4cout.setf(ios_base::fixed,ios_base::floatfield);
        
    const G4ThreeVector StartPosition=       StartFT.GetPosition();
    const G4ThreeVector StartUnitVelocity=   StartFT.GetMomentumDir();
    const G4ThreeVector CurrentPosition=     CurrentFT.GetPosition();
    const G4ThreeVector CurrentUnitVelocity= CurrentFT.GetMomentumDir();
 
    G4double  DotStartCurrentVeloc= StartUnitVelocity.dot(CurrentUnitVelocity);
 
    G4double step_len= CurrentFT.GetCurveLength() - StartFT.GetCurveLength();
    G4double subStepSize = step_len;
     
    if( (subStepNo <= 1) || (verboseLevel > 3) )
    {
       subStepNo = - subStepNo;        // To allow printing banner
 
       G4cout << "------------------------------------------------------------------"
              << G4endl;
       G4cout << std::setw( 6)  << " " << std::setw( 25)
              << " G4DriverReporter: Current Position  and  Direction" << " "
              << G4endl; 
       G4cout << std::setw( 5) << "Step#" << " "
              << std::setw( prec7) << "s-curve" << " "
              << std::setw( prec9) << "X(mm)" << " "
              << std::setw( prec9) << "Y(mm)" << " "  
              << std::setw( prec9) << "Z(mm)" << " "
              << std::setw( prec8) << " N_x " << " "
              << std::setw( prec8) << " N_y " << " "
              << std::setw( prec8) << " N_z " << " "
              << std::setw( 6) << " N^2-1 " << " "
              << std::setw(10) << " N(0).N " << " "
              << std::setw( 7) << "KinEner " << " "
              << std::setw(12) << "Track-l" << " "   // Add the Sub-step ??
              << std::setw(12) << "Step-len" << " " 
              << std::setw(12) << "Step-len" << " " 
              << std::setw( 9) << "ReqStep" << " "  
              << G4endl;
    }
 
    G4cout.precision(noPrecision);
    
    if( (subStepNo <= 0) )
    {
      PrintStat_Aux( StartFT,  requestStep, 0., 
                       0,        0.0,         1.0);
    }
 
    // if( verboseLevel <= 3 )
    {
      G4cout.precision(noPrecision);
      PrintStat_Aux( CurrentFT, requestStep, step_len, 
                     subStepNo, subStepSize, DotStartCurrentVeloc );
    }
    G4cout << "------------------------------------------------------------------"
           << G4endl;
    G4cout.precision(oldPrec);
}

References CLHEP::Hep3Vector::dot(), G4cout, G4endl, G4FieldTrack::GetCurveLength(), G4FieldTrack::GetMomentumDir(), G4FieldTrack::GetPosition(), noPrecision, prec7, prec8, prec9, and PrintStat_Aux().

---

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

• source/geometry/magneticfield/include/G4DriverReporter.hh
• source/geometry/magneticfield/src/G4DriverReporter.cc