#include <G4UnstableFragmentBreakUp.hh>

Inheritance diagram for G4UnstableFragmentBreakUp:

Public Member Functions
	G4UnstableFragmentBreakUp ()

virtual	~G4UnstableFragmentBreakUp ()

virtual G4FragmentVector *	BreakUpFragment (G4Fragment *fragment)

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="Anonymous")

virtual	~G4VEvaporationChannel ()

virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus)

virtual G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)

virtual G4FragmentVector *	BreakUp (const G4Fragment &theNucleus)=0

virtual G4double	GetEmissionProbability (G4Fragment *theNucleus)=0

G4String	GetName () const

void	SetName (const G4String &aName)

void	SetOPTxs (G4int opt)

void	UseSICB (G4bool use)

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()

G4UnstableFragmentBreakUp::G4UnstableFragmentBreakUp ( )

Definition at line 57 of file G4UnstableFragmentBreakUp.cc.

  :verbose(0)
{ 
  fNistManager = G4NistManager::Instance();
  if(0 == Afr[0]) {
    G4int z[6] = {0, 1, 1, 1, 2, 2};
    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]);
    }
  }
}

◆ ~G4UnstableFragmentBreakUp()

G4UnstableFragmentBreakUp::~G4UnstableFragmentBreakUp ( )

virtual

Definition at line 72 of file G4UnstableFragmentBreakUp.cc.

73{}

Member Function Documentation

◆ BreakUp()

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

virtual

Implements G4VEvaporationChannel.

Definition at line 154 of file G4UnstableFragmentBreakUp.cc.

{
  return 0;
}

◆ BreakUpFragment()

G4FragmentVector * G4UnstableFragmentBreakUp::BreakUpFragment ( G4Fragment * fragment )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 80 of file G4UnstableFragmentBreakUp.cc.

{
  //G4cout << "G4UnstableFragmentBreakUp::BreakUpFragment" << G4endl;
  G4FragmentVector * theResult = new G4FragmentVector();
 
  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];
  }
 
  // updated fragment
  if( theResult->size() > 0) {
    nucleus->SetZandA_asInt(Z, A);
    nucleus->SetMomentum(lv);
  }
 
  return theResult;
}