G4ReactionProduct.cc File Reference

#include "G4ReactionProduct.hh"

Go to the source code of this file.

Functions
G4ReactionProduct	operator+ (const G4ReactionProduct &p1, const G4ReactionProduct &p2)
G4ReactionProduct	operator- (const G4ReactionProduct &p1, const G4ReactionProduct &p2)

Variables
G4ThreadLocal G4Allocator < G4ReactionProduct > *	aRPAllocator = 0

Function Documentation

G4ReactionProduct operator+	(	const G4ReactionProduct &	p1,
		const G4ReactionProduct &	p2
	)

Definition at line 235 of file G4ReactionProduct.cc.

References G4ReactionProduct::HasInitialStateParton(), G4ReactionProduct::SetFormationTime(), G4ReactionProduct::SetMass(), G4ReactionProduct::SetMomentum(), G4ReactionProduct::SetPositionInNucleus(), G4ReactionProduct::SetTotalEnergy(), test::x, CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), z, and CLHEP::Hep3Vector::z().

  {
    G4double totEnergy = p1.totalEnergy + p2.totalEnergy;
    G4double x = p1.momentum.x() + p2.momentum.x();
    G4double y = p1.momentum.y() + p2.momentum.y();
    G4double z = p1.momentum.z() + p2.momentum.z();
    G4double newMass = totEnergy*totEnergy - ( x*x + y*y + z*z );
    if( newMass < 0.0 )
      newMass = -1. * std::sqrt( -newMass );
    else
      newMass = std::sqrt( newMass );
    G4ReactionProduct result;
    result.SetMass( newMass );
    result.SetMomentum( x, y, z );
    result.SetTotalEnergy( totEnergy );
    result.SetPositionInNucleus( 0.0, 0.0, 0.0 );
    result.SetFormationTime(0.0);
    result.HasInitialStateParton(false);
    return result;
  }

G4ReactionProduct operator-	(	const G4ReactionProduct &	p1,
		const G4ReactionProduct &	p2
	)

Definition at line 257 of file G4ReactionProduct.cc.

References G4ReactionProduct::HasInitialStateParton(), G4ReactionProduct::SetFormationTime(), G4ReactionProduct::SetMass(), G4ReactionProduct::SetMomentum(), G4ReactionProduct::SetPositionInNucleus(), G4ReactionProduct::SetTotalEnergy(), test::x, CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), z, and CLHEP::Hep3Vector::z().

  {
    G4double totEnergy = p1.totalEnergy - p2.totalEnergy;
    G4double x = p1.momentum.x() - p2.momentum.x();
    G4double y = p1.momentum.y() - p2.momentum.y();
    G4double z = p1.momentum.z() - p2.momentum.z();
    G4double newMass = totEnergy*totEnergy - ( x*x + y*y + z*z );
    if( newMass < 0.0 )
      newMass = -1. * std::sqrt( -newMass );
    else
      newMass = std::sqrt( newMass );
    G4ReactionProduct result;
    result.SetMass( newMass );
    result.SetMomentum( x, y, z );
    result.SetTotalEnergy( totEnergy );
    result.SetPositionInNucleus( 0.0, 0.0, 0.0 );
    result.SetFormationTime(0.0);
    result.HasInitialStateParton(false);
    return result;
  }

Variable Documentation

G4ThreadLocal G4Allocator<G4ReactionProduct>* aRPAllocator = 0

Definition at line 33 of file G4ReactionProduct.cc.

Referenced by G4ReactionProduct::operator delete(), and G4ReactionProduct::operator new().