G4FieldTrack.cc File Reference

#include "G4FieldTrack.hh"

Go to the source code of this file.

Functions
std::ostream &	operator<< (std::ostream &os, const G4FieldTrack &SixVec)

Function Documentation

operator<<()

std::ostream & operator<<	(	std::ostream &	os,
		const G4FieldTrack &	SixVec
	)

Definition at line 33 of file G4FieldTrack.cc.

{
     const G4double* SixV = SixVec.SixVector;
     const G4int precPos= 9;   // For position
     const G4int precEp=  9;   // For Energy / momentum
     const G4int precLen= 12;  // For Length along track
     const G4int precSpin= 9;  // For polarisation
     const G4int precTime= 6;  // For time of flight
     const G4int oldpr= os.precision(precPos);
     os << " ( ";
     os << " X= " << SixV[0] << " " << SixV[1] << " "
                  << SixV[2] << " ";  // Position
     os.precision(precEp);     
     os << " P= " << SixV[3] << " " << SixV[4] << " "
                  << SixV[5] << " ";  // Momentum
     os << " Pmag= "
        << G4ThreeVector(SixV[3], SixV[4], SixV[5]).mag(); // mom magnitude
     os << " Ekin= " << SixVec.fKineticEnergy ;
     os.precision(precLen);
     os << " l= " << SixVec.GetCurveLength();
     os.precision(6);
     os << " m0= " <<   SixVec.fRestMass_c2;
     os << " (Pdir-1)= " <<  SixVec.fMomentumDir.mag()-1.0;
     if( SixVec.fLabTimeOfFlight > 0.0 )
     {
       os.precision(precTime);
     }
     else
     {
       os.precision(3);
     }
     os << " t_lab= "    << SixVec.fLabTimeOfFlight; 
     os << " t_proper= " << SixVec.fProperTimeOfFlight ;
     G4ThreeVector pol= SixVec.GetPolarization();
     if( pol.mag2() > 0.0 )
     {
        os.precision(precSpin);
        os << " PolV= " << pol; // SixVec.GetPolarization();
     }
     else
     {
        os << " PolV= (0,0,0) "; 
     }
     os << " ) ";
     os.precision(oldpr);
     return os;
}