#include <G4UnstableFragmentBreakUp.hh>

Inheritance diagram for G4UnstableFragmentBreakUp:

Public Member Functions
virtual G4bool	BreakUpChain (G4FragmentVector , G4Fragment ) final

G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)

virtual void	Dump () const

virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus)

	G4UnstableFragmentBreakUp ()

virtual G4double	GetEmissionProbability (G4Fragment *fragment) final

virtual G4double	GetLifeTime (G4Fragment *theNucleus)

virtual void	Initialise ()

virtual void	RDMForced (G4bool)

virtual void	SetICM (G4bool)

void	SetOPTxs (G4int opt)

void	SetVerbose (G4int val)

void	UseSICB (G4bool use)

virtual	~G4UnstableFragmentBreakUp ()

Protected Attributes
G4int	OPTxs

G4bool	useSICB

Private Member Functions
	G4UnstableFragmentBreakUp (const G4UnstableFragmentBreakUp &right)

G4bool	operator!= (const G4UnstableFragmentBreakUp &right) const

const G4UnstableFragmentBreakUp &	operator= (const G4UnstableFragmentBreakUp &right)

G4bool	operator== (const G4UnstableFragmentBreakUp &right) const

Private Attributes
G4NuclearLevelData *	fLevelData

G4int	fSecID

G4int	fVerbose

G4double	masses [6]

Static Private Attributes
static const G4int	Afr [6] = {1, 1, 2, 3, 3, 4}

static const G4int	Zfr [6] = {0, 1, 1, 1, 2, 2}

Detailed Description

Definition at line 56 of file G4UnstableFragmentBreakUp.hh.

Constructor & Destructor Documentation

◆ G4UnstableFragmentBreakUp() [1/2]

G4UnstableFragmentBreakUp::G4UnstableFragmentBreakUp ( )

explicit

Definition at line 59 of file G4UnstableFragmentBreakUp.cc.

                                                     : fVerbose(1), fSecID(-1)
{ 
  fLevelData = G4NuclearLevelData::GetInstance();
  for(G4int i=0; i<6; ++i) {
    masses[i] = G4NucleiProperties::GetNuclearMass(Afr[i], Zfr[i]);
  }
  fSecID = G4PhysicsModelCatalog::GetModelID("model_G4UnstableFragmentBreakUp");
}

References Afr, fLevelData, fSecID, G4NuclearLevelData::GetInstance(), G4PhysicsModelCatalog::GetModelID(), G4NucleiProperties::GetNuclearMass(), masses, and Zfr.

◆ ~G4UnstableFragmentBreakUp()

G4UnstableFragmentBreakUp::~G4UnstableFragmentBreakUp ( )

virtual

Definition at line 68 of file G4UnstableFragmentBreakUp.cc.

69{}

◆ G4UnstableFragmentBreakUp() [2/2]

G4UnstableFragmentBreakUp::G4UnstableFragmentBreakUp ( const G4UnstableFragmentBreakUp & right )

private

Member Function Documentation

◆ BreakUpChain()

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

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 71 of file G4UnstableFragmentBreakUp.cc.

{
  //G4cout << "G4UnstableFragmentBreakUp::EmittedFragment" << G4endl;
  G4Fragment* frag = nullptr;
  
  G4int Z = nucleus->GetZ_asInt();
  G4int A = nucleus->GetA_asInt();
 
  G4LorentzVector lv = nucleus->GetMomentum();
  G4double time = nucleus->GetCreationTime();
 
  G4double mass1(0.0), mass2(0.0);
 
  // look for the decay channel with normal masses
  // without Coulomb barrier and paring corrections
  // 1 - recoil, 2 - emitted light ion 
  if(fVerbose > 1) {
    G4cout << "#Unstable decay " << " Z= " << Z << " A= " << A 
       << " Eex(MeV)= " << nucleus->GetExcitationEnergy() << G4endl;
  }
  const G4double tolerance = 10*CLHEP::eV;
  const G4double dmlimit   = 0.2*CLHEP::MeV;
  G4double mass = lv.mag();
  G4double exca = -1000.0;
  G4bool isChannel = false;
  G4int idx = -1;
  for(G4int i=0; i<6; ++i) {
    G4int Zres = Z - Zfr[i];
    G4int Ares = A - Afr[i];
    if(Zres >= 0 && Ares >= Zres && Ares >= Afr[i]) {
      if(Ares <= 4) {
    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;
        */ 
        G4double delm = mass - masses[i] - masses[j];
        if(delm > exca) {
          mass2 = masses[i]; // emitted
          mass1 = masses[j]; // recoil
              exca  = delm; 
          idx = i;
              if(delm > 0.0) { isChannel = true; }
          break;
        }
      }
    }
      }
      if(isChannel) { break; }
      // no simple channel
      G4double mres = G4NucleiProperties::GetNuclearMass(Ares, Zres);
      G4double e = mass - mres - masses[i];
      // select excited state
      const G4LevelManager* lman = fLevelData->GetLevelManager(Zres, Ares);
      if(lman && e >= 0.0) {
    mass2 = masses[i];
    mass1 = mres + e*G4UniformRand();
    idx = i;
    isChannel = true;
    break;
      } 
      // if physical channel is not identified
      // check excitation energy
      if(e > exca) {
    mass2 = masses[i];
    mass1 = mres;
        if(e > 0.0) { mass1 += e; }
    exca  = e;
    idx   = i;
      }
    }
  }
  G4double massmin = mass1 + mass2;
  if(mass < massmin) {
    if(mass + dmlimit < massmin) { return false; }
    if(fVerbose > 1) {
      G4cout << "#Unstable decay correction: Z= " << Z << " A= " << A 
         << " idx= " << idx
         << " deltaM(MeV)= " << mass - massmin
         << G4endl;
    }      
    mass = massmin;
    G4double e = std::max(lv.e(), mass + tolerance);
    G4double mom = std::sqrt((e - mass)*(e + mass));
    G4ThreeVector dir = lv.vect().unit();
    lv.set(dir*mom, e);
  }
 
  // 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[idx], Zfr[idx], mom2);
  frag->SetCreationTime(time);
  frag->SetCreatorModelID(fSecID);
  results->push_back(frag);
 
  // residual
  lv -= mom2;
  Z  -= Zfr[idx];
  A  -= Afr[idx];
    
  nucleus->SetZandA_asInt(Z, A);
  nucleus->SetMomentum(lv);
  nucleus->SetCreatorModelID(fSecID);
  return true;
}

References A, Afr, CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), CLHEP::HepLorentzVector::e(), e2, CLHEP::eV, fLevelData, fSecID, fVerbose, G4cout, G4endl, G4RandomDirection(), G4UniformRand, G4Fragment::GetA_asInt(), G4Fragment::GetCreationTime(), G4Fragment::GetExcitationEnergy(), G4NuclearLevelData::GetLevelManager(), G4Fragment::GetMomentum(), G4NucleiProperties::GetNuclearMass(), G4Fragment::GetZ_asInt(), CLHEP::HepLorentzVector::mag(), masses, G4INCL::Math::max(), CLHEP::MeV, CLHEP::HepLorentzVector::set(), G4Fragment::SetCreationTime(), G4Fragment::SetCreatorModelID(), G4Fragment::SetMomentum(), G4Fragment::SetZandA_asInt(), CLHEP::Hep3Vector::unit(), CLHEP::HepLorentzVector::vect(), Z, and Zfr.

Referenced by G4Evaporation::BreakFragment().

◆ BreakUpFragment()

G4FragmentVector * G4VEvaporationChannel::BreakUpFragment ( G4Fragment * theNucleus )

inlineinherited

Definition at line 106 of file G4VEvaporationChannel.hh.

{
  G4FragmentVector* results = new G4FragmentVector();
  BreakUpChain(results, theNucleus);
  return results;
}

References G4VEvaporationChannel::BreakUpChain().

Referenced by G4NeutronRadCapture::ApplyYourself().

◆ Dump()

void G4VEvaporationChannel::Dump ( ) const

virtualinherited

Reimplemented in G4GEMChannel, and G4GEMChannelVI.

Definition at line 71 of file G4VEvaporationChannel.cc.

72{}

◆ EmittedFragment()

G4Fragment * G4VEvaporationChannel::EmittedFragment ( G4Fragment * theNucleus )

virtualinherited

Reimplemented in G4GEMChannel, G4GEMChannelVI, G4PhotonEvaporation, G4EvaporationChannel, and G4CompetitiveFission.

Definition at line 61 of file G4VEvaporationChannel.cc.

{
  return nullptr;
}

◆ GetEmissionProbability()

G4double G4UnstableFragmentBreakUp::GetEmissionProbability ( G4Fragment * fragment )

finalvirtual

Implements G4VEvaporationChannel.

Definition at line 188 of file G4UnstableFragmentBreakUp.cc.

{
  return 0.0;
}

◆ GetLifeTime()

G4double G4VEvaporationChannel::GetLifeTime ( G4Fragment * theNucleus )

virtualinherited

Definition at line 50 of file G4VEvaporationChannel.cc.

{
  return 0.0;
}

◆ Initialise()

void G4VEvaporationChannel::Initialise ( )

virtualinherited

Reimplemented in G4PhotonEvaporation, and G4EvaporationChannel.

Definition at line 47 of file G4VEvaporationChannel.cc.

48{}

Referenced by G4EvaporationChannel::Initialise(), and G4NeutronRadCapture::InitialiseModel().

◆ 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

◆ RDMForced()

void G4VEvaporationChannel::RDMForced ( G4bool )

virtualinherited

Reimplemented in G4PhotonEvaporation.

Definition at line 58 of file G4VEvaporationChannel.cc.

59{}

◆ SetICM()

void G4VEvaporationChannel::SetICM ( G4bool )

virtualinherited

Reimplemented in G4PhotonEvaporation.

Definition at line 55 of file G4VEvaporationChannel.cc.

56{}

Referenced by G4NeutronRadCapture::InitialiseModel().

◆ SetOPTxs()

void G4VEvaporationChannel::SetOPTxs ( G4int opt )

inlineinherited

Definition at line 113 of file G4VEvaporationChannel.hh.

114{}

◆ SetVerbose()

void G4UnstableFragmentBreakUp::SetVerbose ( G4int val )

inline

Definition at line 69 of file G4UnstableFragmentBreakUp.hh.

69{ fVerbose = val; }

References fVerbose.

Referenced by G4Evaporation::InitialiseChannels().

◆ UseSICB()

void G4VEvaporationChannel::UseSICB ( G4bool use )

inlineinherited

Definition at line 116 of file G4VEvaporationChannel.hh.

117{}

Field Documentation

◆ Afr

const G4int G4UnstableFragmentBreakUp::Afr = {1, 1, 2, 3, 3, 4}

staticprivate

Definition at line 80 of file G4UnstableFragmentBreakUp.hh.

Referenced by BreakUpChain(), and G4UnstableFragmentBreakUp().

◆ fLevelData

G4NuclearLevelData* G4UnstableFragmentBreakUp::fLevelData

private

Definition at line 83 of file G4UnstableFragmentBreakUp.hh.

Referenced by BreakUpChain(), and G4UnstableFragmentBreakUp().

◆ fSecID

G4int G4UnstableFragmentBreakUp::fSecID

private

Definition at line 86 of file G4UnstableFragmentBreakUp.hh.

Referenced by BreakUpChain(), and G4UnstableFragmentBreakUp().

◆ fVerbose

G4int G4UnstableFragmentBreakUp::fVerbose

private

Definition at line 84 of file G4UnstableFragmentBreakUp.hh.

Referenced by BreakUpChain(), and SetVerbose().

◆ masses

G4double G4UnstableFragmentBreakUp::masses[6]

private

Definition at line 81 of file G4UnstableFragmentBreakUp.hh.

Referenced by BreakUpChain(), and G4UnstableFragmentBreakUp().

◆ OPTxs

G4int G4VEvaporationChannel::OPTxs

protectedinherited

Definition at line 93 of file G4VEvaporationChannel.hh.

Referenced by G4EvaporationChannel::GetEmissionProbability().

◆ useSICB

G4bool G4VEvaporationChannel::useSICB

protectedinherited

Definition at line 94 of file G4VEvaporationChannel.hh.

◆ Zfr

const G4int G4UnstableFragmentBreakUp::Zfr = {0, 1, 1, 1, 2, 2}

staticprivate

Definition at line 79 of file G4UnstableFragmentBreakUp.hh.

Referenced by BreakUpChain(), and G4UnstableFragmentBreakUp().

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

source/processes/hadronic/models/de_excitation/evaporation/include/G4UnstableFragmentBreakUp.hh
source/processes/hadronic/models/de_excitation/evaporation/src/G4UnstableFragmentBreakUp.cc

Public Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4UnstableFragmentBreakUp() [1/2]

◆ ~G4UnstableFragmentBreakUp()

◆ G4UnstableFragmentBreakUp() [2/2]

Member Function Documentation

◆ BreakUpChain()

◆ BreakUpFragment()

◆ Dump()

◆ EmittedFragment()

◆ GetEmissionProbability()

◆ GetLifeTime()

◆ Initialise()

◆ operator!=()

◆ operator=()

◆ operator==()

◆ RDMForced()

◆ SetICM()

◆ SetOPTxs()

◆ SetVerbose()

◆ UseSICB()

Field Documentation

◆ Afr

◆ fLevelData

◆ fSecID

◆ fVerbose

◆ masses

◆ OPTxs

◆ useSICB

◆ Zfr