#include <G4UnstableFermiFragment.hh>

Inheritance diagram for G4UnstableFermiFragment:

Public Member Functions
	G4UnstableFermiFragment (G4int anA, G4int aZ, G4int Pol, G4double ExE)

virtual	~G4UnstableFermiFragment ()

virtual G4FragmentVector *	GetFragment (const G4LorentzVector &) const

Public Member Functions inherited from G4VFermiFragment
	G4VFermiFragment (G4int anA, G4int aZ, G4int Pol, G4double ExE)

virtual	~G4VFermiFragment ()

G4int	GetA (void) const

G4int	GetZ (void) const

G4int	GetPolarization (void) const

G4double	GetExcitationEnergy (void) const

G4double	GetFragmentMass (void) const

G4double	GetTotalEnergy (void) const

G4bool	IsStable () const

Protected Attributes
std::vector< G4double >	Masses

std::vector< G4int >	Charges

std::vector< G4int >	AtomNum

Protected Attributes inherited from G4VFermiFragment
G4bool	isStable

G4int	A

G4int	Z

G4int	Polarization

G4double	ExcitEnergy

G4double	fragmentMass

Detailed Description

Definition at line 40 of file G4UnstableFermiFragment.hh.

Constructor & Destructor Documentation

G4UnstableFermiFragment::G4UnstableFermiFragment	(	G4int	anA,
		G4int	aZ,
		G4int	Pol,
		G4double	ExE
	)

Definition at line 36 of file G4UnstableFermiFragment.cc.

References G4VFermiFragment::isStable.

   : G4VFermiFragment(anA,aZ,Pol,ExE)
 {
   isStable = false;
 } 

G4UnstableFermiFragment::~G4UnstableFermiFragment ( )

virtual

Definition at line 43 of file G4UnstableFermiFragment.cc.

44 {}

Member Function Documentation

G4FragmentVector * G4UnstableFermiFragment::GetFragment ( const G4LorentzVector & aMomentum ) const

virtual

Implements G4VFermiFragment.

Definition at line 47 of file G4UnstableFermiFragment.cc.

References AtomNum, CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), Charges, G4FermiPhaseSpaceDecay::Decay(), CLHEP::HepLorentzVector::m(), Masses, N, and test::v.

 {
   std::vector<G4LorentzVector*> * P = thePhaseSpace.Decay(aMomentum.m(), Masses);
 
   G4ThreeVector Beta = aMomentum.boostVector();
 
   G4FragmentVector * theResult = new G4FragmentVector;
   size_t N = P->size();
 
   for (size_t i=0; i<N; ++i)
     {
       G4LorentzVector* v = (*P)[i];
       v->boost(Beta);
       theResult->push_back(new G4Fragment(AtomNum[i],Charges[i],*v));
       
       delete v;
     }
   delete P;
   
   return theResult;
 }