G4INCL::IntersectionFactory::anonymous_namespace{G4INCLIntersection.hh} Namespace Reference

Functions
Intersection	getTrajectoryIntersection (const ThreeVector &x0, const ThreeVector &v, const G4double r, const G4bool earliest)

Function Documentation

getTrajectoryIntersection()

Intersection G4INCL::IntersectionFactory::anonymous_namespace{G4INCLIntersection.hh}::getTrajectoryIntersection	(	const ThreeVector &	x0,
		const ThreeVector &	v,
		const G4double	r,
		const G4bool	earliest
	)

Definition at line 104 of file G4INCLIntersection.hh.

                                                                                                                                   {
        const G4double scalarVelocity = v.mag();
        ThreeVector velocityUnitVector = v / scalarVelocity;
 
        ThreeVector positionTransverse = x0 - velocityUnitVector * x0.dot(velocityUnitVector);
        const G4double impactParameter = positionTransverse.mag();
 
        const G4double r2 = r*r;
        G4double distanceZ2 = r2 - impactParameter * impactParameter;
        if(distanceZ2 < 0.0)
          return Intersection(false, 0.0, ThreeVector());
 
        const G4double distanceZ = std::sqrt(distanceZ2);
        const ThreeVector position = positionTransverse + velocityUnitVector * (earliest ? -distanceZ : distanceZ);
        const G4double time = (position-x0).dot(velocityUnitVector)/scalarVelocity;
        return Intersection(true, time, position);
      }

References G4INCL::ThreeVector::dot(), and G4INCL::ThreeVector::mag().

Referenced by G4INCL::IntersectionFactory::getEarlierTrajectoryIntersection(), G4INCL::IntersectionFactory::getLaterTrajectoryIntersection(), and G4INCL::IntersectionFactory::getTrajectoryIntersections().