G4UnstableFragmentBreakUp Class Reference

#include <G4UnstableFragmentBreakUp.hh>

Inheritance diagram for G4UnstableFragmentBreakUp:

Collaboration diagram for G4UnstableFragmentBreakUp:

Public Member Functions
	G4UnstableFragmentBreakUp ()
virtual	~G4UnstableFragmentBreakUp ()
virtual G4FragmentVector *	BreakUpFragment (G4Fragment *fragment)
virtual G4bool	BreakUpChain (G4FragmentVector *theResult, G4Fragment *theNucleus)
virtual G4Fragment *	EmittedFragment (G4Fragment *fragment)
virtual G4FragmentVector *	BreakUp (const G4Fragment &fragment)
virtual G4double	GetEmissionProbability (G4Fragment *fragment)
void	SetVerboseLevel (G4int val)
Public Member Functions inherited from G4VEvaporationChannel
	G4VEvaporationChannel (const G4String &aName="")
virtual	~G4VEvaporationChannel ()
virtual void	Initialise ()
virtual G4double	GetLifeTime (G4Fragment *theNucleus)
virtual void	Dump () const
virtual void	SetICM (G4bool)
virtual void	RDMForced (G4bool)
virtual G4double	GetFinalLevelEnergy (G4int Z, G4int A, G4double energy)
virtual G4double	GetUpperLevelEnergy (G4int Z, G4int A)
G4double	GetMaxLevelEnergy (G4int Z, G4int A)
G4double	GetNearestLevelEnergy (G4int Z, G4int A, G4double energy)
void	SetPhotonEvaporation (G4VEvaporationChannel *p)
void	SetOPTxs (G4int opt)
void	UseSICB (G4bool use)

Private Member Functions
	G4UnstableFragmentBreakUp (const G4UnstableFragmentBreakUp &right)
const G4UnstableFragmentBreakUp &	operator= (const G4UnstableFragmentBreakUp &right)
G4bool	operator== (const G4UnstableFragmentBreakUp &right) const
G4bool	operator!= (const G4UnstableFragmentBreakUp &right) const

Private Attributes
G4int	verbose
G4int	Zfr [6]
G4int	Afr [6]
G4double	masses [6]
G4NistManager *	fNistManager

Additional Inherited Members
Protected Attributes inherited from G4VEvaporationChannel
G4int	OPTxs
G4bool	useSICB

Detailed Description

Definition at line 57 of file G4UnstableFragmentBreakUp.hh.

Constructor & Destructor Documentation

G4UnstableFragmentBreakUp() [1/2]

G4UnstableFragmentBreakUp::G4UnstableFragmentBreakUp

(

)

Definition at line 53 of file G4UnstableFragmentBreakUp.cc.

  :verbose(0)
{ 
  fNistManager = G4NistManager::Instance();
  const G4int z[6] = {0, 1, 1, 1, 2, 2};
  const G4int a[6] = {1, 1, 2, 3, 3, 4};
  for(G4int i=0; i<6; ++i) {
    Zfr[i] = z[i];
    Afr[i] = a[i];
    masses[i] = G4NucleiProperties::GetNuclearMass(a[i], z[i]);
  }
}

Here is the call graph for this function:

~G4UnstableFragmentBreakUp()

G4UnstableFragmentBreakUp::~G4UnstableFragmentBreakUp ( )

virtual

Definition at line 66 of file G4UnstableFragmentBreakUp.cc.

{}

G4UnstableFragmentBreakUp() [2/2]

G4UnstableFragmentBreakUp::G4UnstableFragmentBreakUp ( const G4UnstableFragmentBreakUp &  right )

private

Member Function Documentation

BreakUp()

G4FragmentVector * G4UnstableFragmentBreakUp::BreakUp ( const G4Fragment &  fragment )

virtual

Implements G4VEvaporationChannel.

Definition at line 153 of file G4UnstableFragmentBreakUp.cc.

{
  return 0;
}

BreakUpChain()

G4bool G4UnstableFragmentBreakUp::BreakUpChain ( G4FragmentVector *  theResult, G4Fragment *  theNucleus  )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 82 of file G4UnstableFragmentBreakUp.cc.

{
  //G4cout << "G4UnstableFragmentBreakUp::BreakUpChain" << G4endl;

  G4int Z = nucleus->GetZ_asInt();
  G4int A = nucleus->GetA_asInt();
  G4int Amax = A;
  G4LorentzVector lv = nucleus->GetMomentum();
  G4double time = nucleus->GetCreationTime();

  G4double deltaE, mass, mass1(0.0), mass2(0.0);
  G4int i, index;

  // Starts loop over evaporated particles, loop is limited by number
  // of nucleons
  for(G4int ia=0; ia<Amax; ++ia) {
 
    mass = lv.mag();
    deltaE = 0.0;
    index  = -1;
    for(i=0; i<6; ++i) {
      G4int Zres = Z - Zfr[i];
      G4int Ares = A - Afr[i];
      if(Zres >= 0 && Ares >= Zres && Ares > 0) {
    G4double m1 = G4NucleiProperties::GetNuclearMass(Ares, Zres);
    G4double de = mass - m1 - masses[i];
        if(de > deltaE) {
      mass1 = m1;
          mass2 = masses[i];
          deltaE= de;
          index = i;
    }
      }
    }

    // no decay channels
    if(index < 0) { break; }

    // compute energy of light fragment
    G4double e2 = 0.5*((mass - mass1)*(mass + mass1) + mass2*mass2)/mass;
    if(e2 < mass2) { break; }

    // sample decay
    G4ThreeVector bst  = lv.boostVector();

    G4double cosTheta = 1. - 2. * G4UniformRand(); 
    G4double sinTheta = std::sqrt(1. - cosTheta * cosTheta);
    G4double phi = twopi * G4UniformRand();
    G4double mom = std::sqrt((e2 - mass2)*(e2 + mass2));
    G4LorentzVector mom2(mom * sinTheta * std::cos(phi),
             mom * sinTheta * std::sin(phi),
             mom * cosTheta,
             e2);
    mom2.boost(bst);  
    G4Fragment* fr = new G4Fragment(Afr[index], Zfr[index], mom2);
    fr->SetCreationTime(time);
    theResult->push_back(fr);

    // residual
    lv -= mom2;
    Z  -= Zfr[index];
    A  -= Afr[index];
  }

  nucleus->SetZandA_asInt(Z, A);
  nucleus->SetMomentum(lv);
  theResult->push_back(nucleus);
  return false;
}

Here is the call graph for this function:

Here is the caller graph for this function:

BreakUpFragment()

G4FragmentVector * G4UnstableFragmentBreakUp::BreakUpFragment ( G4Fragment *  fragment )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 74 of file G4UnstableFragmentBreakUp.cc.

{
  //G4cout << "G4UnstableFragmentBreakUp::BreakUpFragment" << G4endl;
  G4FragmentVector * theResult = new G4FragmentVector();
  BreakUpChain(theResult, nucleus);
  return theResult;
}

Here is the call graph for this function:

EmittedFragment()

G4Fragment * G4UnstableFragmentBreakUp::EmittedFragment ( G4Fragment *  fragment )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 69 of file G4UnstableFragmentBreakUp.cc.

{
  return 0;
}

GetEmissionProbability()

G4double G4UnstableFragmentBreakUp::GetEmissionProbability ( G4Fragment *  fragment )

virtual

Implements G4VEvaporationChannel.

Definition at line 158 of file G4UnstableFragmentBreakUp.cc.

{
  return 0.0;
}

operator!=()

G4bool G4UnstableFragmentBreakUp::operator!= ( const G4UnstableFragmentBreakUp &  right ) const

private

operator=()

const G4UnstableFragmentBreakUp& G4UnstableFragmentBreakUp::operator= ( const G4UnstableFragmentBreakUp &  right )

private

operator==()

G4bool G4UnstableFragmentBreakUp::operator== ( const G4UnstableFragmentBreakUp &  right ) const

private

SetVerboseLevel()

void G4UnstableFragmentBreakUp::SetVerboseLevel ( G4int  val )

inline

Definition at line 100 of file G4UnstableFragmentBreakUp.hh.

{
  verbose = val;
}

Member Data Documentation

Afr

G4int G4UnstableFragmentBreakUp::Afr[6]

private

Definition at line 94 of file G4UnstableFragmentBreakUp.hh.

fNistManager

G4NistManager* G4UnstableFragmentBreakUp::fNistManager

private

Definition at line 97 of file G4UnstableFragmentBreakUp.hh.

masses

G4double G4UnstableFragmentBreakUp::masses[6]

private

Definition at line 95 of file G4UnstableFragmentBreakUp.hh.

verbose

G4int G4UnstableFragmentBreakUp::verbose

private

Definition at line 91 of file G4UnstableFragmentBreakUp.hh.

Zfr

G4int G4UnstableFragmentBreakUp::Zfr[6]

private

Definition at line 93 of file G4UnstableFragmentBreakUp.hh.

---

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

• G4UnstableFragmentBreakUp.hh
• G4UnstableFragmentBreakUp.cc