G4GDecay3 Class Reference

#include <G4GDecay3.hh>

Public Member Functions
	G4GDecay3 ()
	G4GDecay3 (const G4double &pMass, const G4double &dMass0, const G4double &dMass1, const G4double &dMass2)
	~G4GDecay3 ()
std::vector< G4ThreeVector >	GetThreeBodyMomenta ()

Detailed Description

Definition at line 43 of file G4GDecay3.hh.

Constructor & Destructor Documentation

G4GDecay3::G4GDecay3 ( )

inline

Definition at line 46 of file G4GDecay3.hh.

{;}

G4GDecay3::G4GDecay3	(	const G4double &	pMass,
		const G4double &	dMass0,
		const G4double &	dMass1,
		const G4double &	dMass2
	)

Definition at line 41 of file G4GDecay3.cc.

 : loopMax(100), parentMass(pMass), mDaughter0(dMass0), mDaughter1(dMass1),
   mDaughter2(dMass2), pDaughter0(0.), pDaughter1(0.), pDaughter2(0.) {;}

G4GDecay3::~G4GDecay3 ( )

inline

Definition at line 49 of file G4GDecay3.hh.

{;}

Member Function Documentation

std::vector< G4ThreeVector > G4GDecay3::GetThreeBodyMomenta

(

)

Definition at line 98 of file G4GDecay3.cc.

References G4cerr, G4endl, G4UniformRand, CLHEP::Hep3Vector::setX(), CLHEP::Hep3Vector::setY(), CLHEP::Hep3Vector::setZ(), and python.hepunit::twopi.

{

  std::vector<G4ThreeVector> pVect;

  if (CalculateMomentumMagnitudes() ) {
    
    // Calculate directions
    G4double costheta = 2.*G4UniformRand()-1.;
    G4double sintheta = std::sqrt((1.0-costheta)*(1.0+costheta));
    G4double phi = twopi*G4UniformRand();
    G4double sinphi = std::sin(phi);
    G4double cosphi = std::cos(phi);
    G4ThreeVector direction0(sintheta*cosphi, sintheta*sinphi, costheta);

    G4double costhetan = (pDaughter1*pDaughter1 - pDaughter2*pDaughter2
                          - pDaughter0*pDaughter0)/(2.0*pDaughter2*pDaughter0);
    G4double sinthetan = std::sqrt((1.0-costhetan)*(1.0+costhetan));
    G4double phin = twopi*G4UniformRand();
    G4double sinphin = std::sin(phin);
    G4double cosphin = std::cos(phin);
    G4ThreeVector direction2;
    direction2.setX(sinthetan*cosphin*costheta*cosphi -
                    sinthetan*sinphin*sinphi + costhetan*sintheta*cosphi);
    direction2.setY(sinthetan*cosphin*costheta*sinphi +
                    sinthetan*sinphin*cosphi + costhetan*sintheta*sinphi);
    direction2.setZ(-sinthetan*cosphin*sintheta + costhetan*costheta);

    // Return momentum vectors
    pVect.push_back(pDaughter0*direction0);
    pVect.push_back(-direction0*pDaughter0 - direction2*pDaughter2);
    pVect.push_back(pDaughter2*direction2);

  } else {
    G4cerr << "G4GDecay::GetThreeBodyMomenta: " << loopMax
           << " or more loops in momentum magnitude calculation " << G4endl;
  }

  return pVect;
}

---

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

• G4GDecay3.hh
• G4GDecay3.cc