G4UnstableFragmentBreakUp Class Reference

#include <G4UnstableFragmentBreakUp.hh>

Inheritance diagram for G4UnstableFragmentBreakUp:

Collaboration diagram for G4UnstableFragmentBreakUp:

Public Member Functions
	G4UnstableFragmentBreakUp ()
virtual	~G4UnstableFragmentBreakUp ()
virtual G4bool	BreakUpChain (G4FragmentVector *, G4Fragment *) final
virtual G4double	GetEmissionProbability (G4Fragment *fragment) final
Public Member Functions inherited from G4VEvaporationChannel
	G4VEvaporationChannel (const G4String &aName="")
virtual	~G4VEvaporationChannel ()
virtual void	Initialise ()
virtual G4double	GetLifeTime (G4Fragment *theNucleus)
virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus)
G4FragmentVector *	BreakUpFragment (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)

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

explicit

Definition at line 58 of file G4UnstableFragmentBreakUp.cc.

{ 
  fLevelData = G4NuclearLevelData::GetInstance();
  for(G4int i=0; i<6; ++i) {
    masses[i] = G4NucleiProperties::GetNuclearMass(Afr[i], Zfr[i]);
  }
}

Here is the call graph for this function:

G4UnstableFragmentBreakUp::~G4UnstableFragmentBreakUp ( )

virtual

Definition at line 66 of file G4UnstableFragmentBreakUp.cc.

{}

Member Function Documentation

G4bool G4UnstableFragmentBreakUp::BreakUpChain ( G4FragmentVector *  results, G4Fragment *  nucleus  )

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 69 of file G4UnstableFragmentBreakUp.cc.

{
  //G4cout << "G4UnstableFragmentBreakUp::EmittedFragment" << G4endl;
  G4Fragment* frag = nullptr;
  
  G4int Z = nucleus->GetZ_asInt();
  G4int A = nucleus->GetA_asInt();

  // if the isotope is in the database it is not exotic
  // so, cannot be handled by this class
  if(fLevelData->GetLevelManager(Z, A)) { return false; }

  G4LorentzVector lv = nucleus->GetMomentum();
  G4double time = nucleus->GetCreationTime();

  G4double mass1(0.0), mass2(0.0);
 
  G4double mass = lv.mag();
  G4bool   done = false;

  G4int Amax = A;
  G4int Zres = Z;
  G4int Ares = A;

  for(G4int k=0; k<Amax; ++k) {
    // look for the decay channel with normal masses
    // without Coulomb barrier and paring corrections

    // 1 - recoil, 2 - emitted light ion 

    //G4cout << "Unstable decay #" << k << " Z= " << Z << " A= " << A << G4endl;

    G4double ekin = 0.0;
    G4int index = -1;
    for(G4int i=0; i<6; ++i) {
      Zres = Z - Zfr[i];
      Ares = A - Afr[i];
      for(G4int j=0; j<6; ++j) {
    if(Zres == Zfr[j] && Ares == Afr[j]) {
      /*
      G4cout << "i= " << i << " j= " << j << " Zres= " << Zres
         << " Ares= " << Ares << " dm= " << mass - masses[i] - masses[j]
         << G4endl;
      */ 
      if(mass >= masses[i] + masses[j]) {
        mass2 = masses[i];
        mass1 = masses[j];
        ekin  = 0.0;
        done = true;
        break;
      }
    }
      }
      if(done) {
    index = i;
        break;
      }
      if(Zres >= 0 && Ares >= Zres && Ares > 0) {
    G4double m = G4NucleiProperties::GetNuclearMass(Ares, Zres);
    G4double e = mass - m - masses[i];
    if(e > ekin) {
      mass2 = masses[i];
      mass1 = m;
      ekin  = e;
      index = i;
      if(fLevelData->GetLevelManager(Z, A)) {
        done = true;
            break;
      } 
    }
      }
    }

    // no decay channel - assume that primary mass is biased
    // only energy will be conserved
    if(index < 0) {
      for(G4int i=0; i<6; ++i) {
    Zres = Z - Zfr[i];
    Ares = A - Afr[i];
    if(Zres >= 0 && Ares >= Zres && Ares > 0) {
          
      G4double m = proton_mass_c2*Zres + neutron_mass_c2*(Ares - Zres);
      G4double e = mass - m - masses[i];
      if(e > ekin) {
        mass2 = masses[i];
        mass1 = m;
        ekin  = e;
        index = i;
      }
    }
      }
    }
    if(index < 0) {

      // further decay impossible
      // if no one decay sampled do not update primary 
      if(0 == k) { return false; }
      else       { break; }
    }

    // useful to left max excitation for the residual
    mass1 += ekin;

    // compute energy of light fragment
    G4double e2 = 0.5*((mass - mass1)*(mass + mass1) + mass2*mass2)/mass;
    e2 = std::max(e2, mass2);
    G4double mom = std::sqrt((e2 - mass2)*(e2 + mass2));

    // sample decay
    G4ThreeVector bst  = lv.boostVector();
    G4ThreeVector v = G4RandomDirection();
    G4LorentzVector mom2 = G4LorentzVector(v*mom, e2);
    mom2.boost(bst);  
    frag = new G4Fragment(Afr[index], Zfr[index], mom2);
    frag->SetCreationTime(time);
    results->push_back(frag);

    // residual
    lv -= mom2;
    Z  -= Zfr[index];
    A  -= Afr[index];
    
    mass = lv.mag();

    if(done) { break; }
  }

  // updated primary
  nucleus->SetZandA_asInt(Z, A);
  nucleus->SetMomentum(lv);
  return true;
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4UnstableFragmentBreakUp::GetEmissionProbability ( G4Fragment *  fragment )

finalvirtual

Implements G4VEvaporationChannel.

Definition at line 203 of file G4UnstableFragmentBreakUp.cc.

{
  return 0.0;
}

---

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

• source/geant4.10.03.p03/source/processes/hadronic/models/de_excitation/evaporation/include/G4UnstableFragmentBreakUp.hh
• source/geant4.10.03.p03/source/processes/hadronic/models/de_excitation/evaporation/src/G4UnstableFragmentBreakUp.cc