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 42 of file G4GDecay3.hh.

Constructor & Destructor Documentation

G4GDecay3() [1/2]

G4GDecay3::G4GDecay3 ( )

inline

Definition at line 45 of file G4GDecay3.hh.

{;}

G4GDecay3() [2/2]

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

Definition at line 42 of file G4GDecay3.cc.

 : loopMax(100), mDaughter0(dMass0), mDaughter1(dMass1),
   mDaughter2(dMass2), pDaughter0(0.), pDaughter1(0.), pDaughter2(0.) 
{
  parentMass = std::max(pMass, mDaughter0 + mDaughter1 + mDaughter2 + CLHEP::keV);
}

~G4GDecay3()

G4GDecay3::~G4GDecay3 ( )

inline

Definition at line 48 of file G4GDecay3.hh.

{;}

Member Function Documentation

GetThreeBodyMomenta()

std::vector< G4ThreeVector > G4GDecay3::GetThreeBodyMomenta

(

)

Definition at line 96 of file G4GDecay3.cc.

{
 
  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 << "G4GDecay3::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