#include <G4FermiBreakUpVI.hh>

Inheritance diagram for G4FermiBreakUpVI:

Public Member Functions
void	BreakFragment (G4FragmentVector , G4Fragment theNucleus) final

	G4FermiBreakUpVI ()

void	Initialise () final

G4bool	IsApplicable (G4int ZZ, G4int AA, G4double etot) const final

void	SetVerbose (G4int val)

	~G4FermiBreakUpVI () final

Protected Attributes
G4int	verbose

Private Member Functions
	G4FermiBreakUpVI (const G4FermiBreakUpVI &right)=delete

void	InitialisePool ()

G4bool	operator!= (const G4FermiBreakUpVI &right) const =delete

const G4FermiBreakUpVI &	operator= (const G4FermiBreakUpVI &right)=delete

G4bool	operator== (const G4FermiBreakUpVI &right) const =delete

G4bool	SampleDecay ()

Private Attributes
G4int	A

G4ThreeVector	boostVector

G4double	elim

G4double	excitation

std::vector< const G4FermiFragment * >	frag

const G4FermiFragment *	frag1

const G4FermiFragment *	frag2

G4LorentzVector	lv0

std::vector< G4LorentzVector >	lvect

G4double	mass

G4int	maxA

G4int	maxZ

std::vector< G4double >	prob

CLHEP::HepRandomEngine *	rndmEngine

G4int	secID

G4int	spin

const G4FermiDecayProbability *	theDecay

G4double	tolerance

G4int	Z

Static Private Attributes
static G4FermiFragmentsPoolVI *	thePool = nullptr

Detailed Description

Definition at line 46 of file G4FermiBreakUpVI.hh.

Constructor & Destructor Documentation

◆ G4FermiBreakUpVI() [1/2]

G4FermiBreakUpVI::G4FermiBreakUpVI ( )

explicit

Definition at line 47 of file G4FermiBreakUpVI.cc.

  : theDecay(nullptr), rndmEngine(nullptr), maxZ(9), maxA(17), secID(-1)
{
  frag.reserve(10);
  lvect.reserve(10);
  Z = A = spin = 0;
  secID = G4PhysicsModelCatalog::GetModelID("model_G4FermiBreakUpVI");
  mass = elim = excitation = 0.0;
  tolerance = CLHEP::MeV;  
  frag1 = frag2 = nullptr;
  prob.resize(12,0.0);
  Initialise();
}

References A, elim, excitation, frag, frag1, frag2, G4PhysicsModelCatalog::GetModelID(), Initialise(), lvect, mass, CLHEP::MeV, prob, secID, spin, tolerance, and Z.

◆ ~G4FermiBreakUpVI()

G4FermiBreakUpVI::~G4FermiBreakUpVI ( )

final

Definition at line 61 of file G4FermiBreakUpVI.cc.

{
  if(G4Threading::IsMasterThread()) { 
    delete thePool;
    thePool = nullptr;
  }
}

References G4Threading::IsMasterThread(), and thePool.

◆ G4FermiBreakUpVI() [2/2]

G4FermiBreakUpVI::G4FermiBreakUpVI ( const G4FermiBreakUpVI & right )

privatedelete

Member Function Documentation

◆ BreakFragment()

void G4FermiBreakUpVI::BreakFragment	(	G4FragmentVector *	theResult,
		G4Fragment *	theNucleus
	)

finalvirtual

Implements G4VFermiBreakUp.

Definition at line 98 of file G4FermiBreakUpVI.cc.

{
  if(verbose > 1) {
    G4cout << "### G4FermiBreakUpVI::BreakFragment start new fragment " 
           << G4endl;
    G4cout << *theNucleus << G4endl;
  }
 
  // initial fragment
  Z = theNucleus->GetZ_asInt();
  A = theNucleus->GetA_asInt();
  excitation = theNucleus->GetExcitationEnergy();
  mass = theNucleus->GetGroundStateMass() + excitation;
  spin = -1;
 
  lv0 = theNucleus->GetMomentum();
  rndmEngine = G4Random::getTheEngine();
 
  // sample first decay of an initial state
  // if not possible to decay - exit
  if(!SampleDecay()) {
    return; 
  }
 
  G4double time = theNucleus->GetCreationTime();
  delete theNucleus;
 
  static const G4int imax = 100; 
 
  // loop over vector of Fermi fragments
  // vector may grow at each iteraction
  for(size_t i=0; i<frag.size(); ++i) {
    Z = frag[i]->GetZ();
    A = frag[i]->GetA();
    spin = frag[i]->GetSpin();
    mass = frag[i]->GetTotalEnergy();
    lv0 = lvect[i];
    if(verbose > 1) {
      G4cout << "# FermiFrag " << i << ".  Z= " << Z << " A= " << A 
         << " mass= " << mass << " exc= " 
         << frag[i]->GetExcitationEnergy() << G4endl;
    }
    // stable fragment
    if(!SampleDecay()) {
      if(verbose > 1) { G4cout << "   New G4Fragment" << G4endl; }
      G4Fragment* f = new G4Fragment(A, Z, lv0);
      f->SetSpin(0.5*spin);
      f->SetCreationTime(time);
      f->SetCreatorModelID(secID);
      theResult->push_back(f);
    }
    // limit the loop
    if(i == imax) {
      break;
    }
  }
  frag.clear();
  lvect.clear();
}

References A, excitation, frag, G4cout, G4endl, G4Fragment::GetA_asInt(), G4Fragment::GetCreationTime(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetGroundStateMass(), G4Fragment::GetMomentum(), G4Fragment::GetZ_asInt(), imax, lv0, lvect, mass, rndmEngine, SampleDecay(), secID, G4Fragment::SetCreationTime(), G4Fragment::SetCreatorModelID(), G4Fragment::SetSpin(), spin, G4VFermiBreakUp::verbose, and Z.

◆ Initialise()

void G4FermiBreakUpVI::Initialise ( )

finalvirtual

Implements G4VFermiBreakUp.

Definition at line 69 of file G4FermiBreakUpVI.cc.

{
  if(verbose > 1) {
    G4cout << "### G4FermiBreakUpVI::Initialise(): " << thePool << G4endl;
  }
  if(thePool == nullptr) { InitialisePool(); }
  theDecay = thePool->FermiDecayProbability();
  elim = thePool->GetEnergyLimit();
}

References elim, G4FermiFragmentsPoolVI::FermiDecayProbability(), G4cout, G4endl, G4FermiFragmentsPoolVI::GetEnergyLimit(), InitialisePool(), theDecay, thePool, and G4VFermiBreakUp::verbose.

Referenced by G4FermiBreakUpVI().

◆ InitialisePool()

void G4FermiBreakUpVI::InitialisePool ( )

private

Definition at line 79 of file G4FermiBreakUpVI.cc.

{
#ifdef G4MULTITHREADED
  G4MUTEXLOCK(&G4FermiBreakUpVI::FermiBreakUpVIMutex);
#endif
  if(thePool == nullptr) {
    thePool = new G4FermiFragmentsPoolVI();
  }
#ifdef G4MULTITHREADED
  G4MUTEXUNLOCK(&G4FermiBreakUpVI::FermiBreakUpVIMutex);
#endif
}

References G4MUTEXLOCK, G4MUTEXUNLOCK, and thePool.

Referenced by Initialise().

◆ IsApplicable()

G4bool G4FermiBreakUpVI::IsApplicable	(	G4int	ZZ,
		G4int	AA,
		G4double	etot
	)		const

finalvirtual

Implements G4VFermiBreakUp.

Definition at line 92 of file G4FermiBreakUpVI.cc.

{
  return (ZZ < maxZ && AA < maxA && AA > 0 && eexc <= elim 
      && thePool->HasChannels(ZZ, AA, eexc));
}

References anonymous_namespace{paraMaker.cc}::AA, elim, G4FermiFragmentsPoolVI::HasChannels(), maxZ, and thePool.

◆ operator!=()

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

privatedelete

◆ operator=()

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

privatedelete

◆ operator==()

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

privatedelete

◆ SampleDecay()

G4bool G4FermiBreakUpVI::SampleDecay ( )

private

Definition at line 159 of file G4FermiBreakUpVI.cc.

{
  const G4FermiChannels* chan = thePool->ClosestChannels(Z, A, mass);
  if(!chan) { return false; }
  size_t nn = chan->GetNumberOfChannels();
  if(verbose > 1) {
    G4cout << "== SampleDecay " << nn << " channels Eex= " 
       << chan->GetExcitation() << G4endl;
  }
  if(0 == nn) { return false; }
 
  const G4FermiPair* fpair = nullptr;
 
  // one unstable fragment
  if(1 == nn) {
    fpair = chan->GetPair(0);
 
    // more pairs
  } else {
    
    // in static probabilities may be used
    if(std::abs(excitation - chan->GetExcitation()) < tolerance) {
      fpair = chan->SamplePair(rndmEngine->flat());
      
    } else {
 
      // recompute probabilities
      const std::vector<const G4FermiPair*>& pvect = chan->GetChannels();
      if(nn > 12) { prob.resize(nn, 0.0); }
      G4double ptot = 0.0;
      if(verbose > 2) { 
    G4cout << "Start recompute probabilities" << G4endl;
      }
      for(size_t i=0; i<nn; ++i) {
    ptot += theDecay->ComputeProbability(Z, A, -1, mass,
                         pvect[i]->GetFragment1(),
                         pvect[i]->GetFragment2());
    prob[i] = ptot;
    if(verbose > 2) { 
      G4cout << i << ". " << prob[i]
         << " Z1= " << pvect[i]->GetFragment1()->GetZ()
         << " A1= " << pvect[i]->GetFragment1()->GetA()
         << " Z2= " << pvect[i]->GetFragment2()->GetZ()
         << " A2= " << pvect[i]->GetFragment2()->GetA()
         << G4endl;
    }
      }
      ptot *= rndmEngine->flat();
      for(size_t i=0; i<nn; ++i) {
        if(ptot <= prob[i] || i+1 == nn) {
          fpair = pvect[i];
      break;
    }
      }
    }
  }
  if(!fpair) { return false; }
 
  frag1 = fpair->GetFragment1();
  frag2 = fpair->GetFragment2();
  
  G4double mass1 = frag1->GetTotalEnergy();
  G4double mass2 = frag2->GetTotalEnergy();
  if(verbose > 2) {
    G4cout << " M= " << mass << " M1= " << mass1 << "  M2= " << mass2 
       << " Exc1= " << frag1->GetExcitationEnergy() 
       << " Exc2= " << frag2->GetExcitationEnergy() << G4endl;
  }
  // sample fragment1
  G4double e1 = 0.5*(mass*mass - mass2*mass2 + mass1*mass1)/mass;
  //G4cout << "ekin1= " << e1 - mass1 << G4endl;
  G4double p1(0.0);
  if(e1 > mass1) {
    p1 = std::sqrt((e1 - mass1)*(e1 + mass1));
  } else {
    e1 = mass1;
  }
  G4ThreeVector v = G4RandomDirection();
  G4LorentzVector lv1 = G4LorentzVector(v*p1, e1);
 
  // compute kinematics
  boostVector = lv0.boostVector();  
  lv1.boost(boostVector);
 
  lv0 -= lv1;
 
  G4double e2 = lv0.e();
  if(e2 < mass2) {
    lv0.set(0.,0.,0.,mass2);
  }
 
  frag.push_back(frag1);
  frag.push_back(frag2);
  lvect.push_back(lv1);
  lvect.push_back(lv0);
 
  return true;
}

References A, CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), boostVector, G4FermiFragmentsPoolVI::ClosestChannels(), G4FermiDecayProbability::ComputeProbability(), CLHEP::HepLorentzVector::e(), e1, e2, excitation, CLHEP::HepRandomEngine::flat(), frag, frag1, frag2, G4cout, G4endl, G4RandomDirection(), G4FermiChannels::GetChannels(), G4FermiChannels::GetExcitation(), G4FermiFragment::GetExcitationEnergy(), G4FermiPair::GetFragment1(), G4FermiPair::GetFragment2(), G4FermiChannels::GetNumberOfChannels(), G4FermiChannels::GetPair(), G4FermiFragment::GetTotalEnergy(), lv0, lvect, mass, G4InuclParticleNames::nn, prob, rndmEngine, G4FermiChannels::SamplePair(), CLHEP::HepLorentzVector::set(), theDecay, thePool, tolerance, G4VFermiBreakUp::verbose, and Z.

Referenced by BreakFragment().

◆ SetVerbose()

void G4VFermiBreakUp::SetVerbose ( G4int val )

inlineinherited

Definition at line 59 of file G4VFermiBreakUp.hh.

59{ verbose = val; }

References G4VFermiBreakUp::verbose.

Referenced by G4ExcitationHandler::SetParameters().

Field Documentation

◆ A

G4int G4FermiBreakUpVI::A

private

Definition at line 84 of file G4FermiBreakUpVI.hh.

Referenced by BreakFragment(), G4FermiBreakUpVI(), mcscore.MCParticle::printout(), and SampleDecay().

◆ boostVector

G4ThreeVector G4FermiBreakUpVI::boostVector

private

Definition at line 98 of file G4FermiBreakUpVI.hh.

Referenced by SampleDecay().

◆ elim

G4double G4FermiBreakUpVI::elim

private

Definition at line 92 of file G4FermiBreakUpVI.hh.

Referenced by G4FermiBreakUpVI(), Initialise(), and IsApplicable().

◆ excitation

G4double G4FermiBreakUpVI::excitation

private

Definition at line 90 of file G4FermiBreakUpVI.hh.

Referenced by BreakFragment(), G4FermiBreakUpVI(), and SampleDecay().

◆ frag

std::vector<const G4FermiFragment*> G4FermiBreakUpVI::frag

private

Definition at line 101 of file G4FermiBreakUpVI.hh.

Referenced by BreakFragment(), G4FermiBreakUpVI(), and SampleDecay().

◆ frag1

const G4FermiFragment* G4FermiBreakUpVI::frag1

private

Definition at line 94 of file G4FermiBreakUpVI.hh.

Referenced by G4FermiBreakUpVI(), and SampleDecay().

◆ frag2

const G4FermiFragment* G4FermiBreakUpVI::frag2

private

Definition at line 95 of file G4FermiBreakUpVI.hh.

Referenced by G4FermiBreakUpVI(), and SampleDecay().

◆ lv0

G4LorentzVector G4FermiBreakUpVI::lv0

private

Definition at line 97 of file G4FermiBreakUpVI.hh.

Referenced by BreakFragment(), and SampleDecay().

◆ lvect

std::vector<G4LorentzVector> G4FermiBreakUpVI::lvect

private

Definition at line 102 of file G4FermiBreakUpVI.hh.

Referenced by BreakFragment(), G4FermiBreakUpVI(), and SampleDecay().

◆ mass

G4double G4FermiBreakUpVI::mass

private

Definition at line 89 of file G4FermiBreakUpVI.hh.

Referenced by BreakFragment(), G4FermiBreakUpVI(), and SampleDecay().

◆ maxA

G4int G4FermiBreakUpVI::maxA

private

Definition at line 81 of file G4FermiBreakUpVI.hh.

◆ maxZ

G4int G4FermiBreakUpVI::maxZ

private

Definition at line 80 of file G4FermiBreakUpVI.hh.

Referenced by IsApplicable().

◆ prob

std::vector<G4double> G4FermiBreakUpVI::prob

private

Definition at line 100 of file G4FermiBreakUpVI.hh.

Referenced by G4FermiBreakUpVI(), and SampleDecay().

◆ rndmEngine

CLHEP::HepRandomEngine* G4FermiBreakUpVI::rndmEngine

private

Definition at line 78 of file G4FermiBreakUpVI.hh.

Referenced by BreakFragment(), and SampleDecay().

◆ secID

G4int G4FermiBreakUpVI::secID

private

Definition at line 87 of file G4FermiBreakUpVI.hh.

Referenced by BreakFragment(), and G4FermiBreakUpVI().

◆ spin

G4int G4FermiBreakUpVI::spin

private

Definition at line 85 of file G4FermiBreakUpVI.hh.

Referenced by BreakFragment(), and G4FermiBreakUpVI().

◆ theDecay

const G4FermiDecayProbability* G4FermiBreakUpVI::theDecay

private

Definition at line 76 of file G4FermiBreakUpVI.hh.

Referenced by Initialise(), and SampleDecay().

◆ thePool

G4FermiFragmentsPoolVI * G4FermiBreakUpVI::thePool = nullptr

staticprivate

Definition at line 75 of file G4FermiBreakUpVI.hh.

Referenced by Initialise(), InitialisePool(), IsApplicable(), SampleDecay(), and ~G4FermiBreakUpVI().

◆ tolerance

G4double G4FermiBreakUpVI::tolerance

private

Definition at line 91 of file G4FermiBreakUpVI.hh.

Referenced by G4FermiBreakUpVI(), and SampleDecay().

◆ verbose

G4int G4VFermiBreakUp::verbose

protectedinherited

Definition at line 63 of file G4VFermiBreakUp.hh.

Referenced by BreakFragment(), Initialise(), SampleDecay(), and G4VFermiBreakUp::SetVerbose().

◆ Z

G4int G4FermiBreakUpVI::Z

private

Definition at line 83 of file G4FermiBreakUpVI.hh.

Referenced by BreakFragment(), G4FermiBreakUpVI(), mcscore.MCParticle::printout(), and SampleDecay().

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

source/processes/hadronic/models/de_excitation/fermi_breakup/include/G4FermiBreakUpVI.hh
source/processes/hadronic/models/de_excitation/fermi_breakup/src/G4FermiBreakUpVI.cc

Public Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4FermiBreakUpVI() [1/2]

◆ ~G4FermiBreakUpVI()

◆ G4FermiBreakUpVI() [2/2]

Member Function Documentation

◆ BreakFragment()

◆ Initialise()

◆ InitialisePool()

◆ IsApplicable()

◆ operator!=()

◆ operator=()

◆ operator==()

◆ SampleDecay()

◆ SetVerbose()

Field Documentation

◆ A

◆ boostVector

◆ elim

◆ excitation

◆ frag

◆ frag1

◆ frag2

◆ lv0

◆ lvect

◆ mass

◆ maxA

◆ maxZ

◆ prob

◆ rndmEngine

◆ secID

◆ spin

◆ theDecay

◆ thePool

◆ tolerance

◆ verbose

◆ Z