#include <G4GEMChannelVI.hh>

Inheritance diagram for G4GEMChannelVI:

Public Member Functions
virtual G4bool	BreakUpChain (G4FragmentVector theResult, G4Fragment theNucleus)

G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)

void	Dump () const final

G4Fragment *	EmittedFragment (G4Fragment *theNucleus) final

	G4GEMChannelVI (G4int theA, G4int theZ)

G4double	GetEmissionProbability (G4Fragment *theNucleus) final

virtual G4double	GetLifeTime (G4Fragment *theNucleus)

virtual void	Initialise ()

virtual void	RDMForced (G4bool)

virtual void	SetICM (G4bool)

void	SetOPTxs (G4int opt)

void	UseSICB (G4bool use)

	~G4GEMChannelVI () final

Protected Attributes
G4int	OPTxs

G4bool	useSICB

Private Member Functions
	G4GEMChannelVI (const G4GEMChannelVI &right)

G4bool	operator!= (const G4GEMChannelVI &right) const

const G4GEMChannelVI &	operator= (const G4GEMChannelVI &right)

G4bool	operator== (const G4GEMChannelVI &right) const

Private Attributes
G4int	A

const G4VCoulombBarrier *	cBarrier

G4double	evapMass

G4double	evapMass2

G4GEMProbabilityVI *	fProbability

G4int	fragA

G4int	fragZ

G4double	mass

const G4PairingCorrection *	pairingCorrection

G4int	resA

G4double	resMass

G4int	resZ

G4int	secID

G4int	Z

Detailed Description

Definition at line 40 of file G4GEMChannelVI.hh.

Constructor & Destructor Documentation

◆ G4GEMChannelVI() [1/2]

G4GEMChannelVI::G4GEMChannelVI	(	G4int	theA,
		G4int	theZ
	)

explicit

Definition at line 42 of file G4GEMChannelVI.cc.

  : A(theA), Z(theZ), secID(-1)
{ 
  G4NuclearLevelData* nData = G4NuclearLevelData::GetInstance();
  pairingCorrection = nData->GetPairingCorrection();
  const G4LevelManager* lManager = nullptr;
  if(A > 4) { lManager = nData->GetLevelManager(Z, A); }
  evapMass = G4NucleiProperties::GetNuclearMass(A, Z);
  evapMass2 = evapMass*evapMass;
 
  cBarrier = new G4CoulombBarrier(A, Z);
  fProbability = new G4GEMProbabilityVI(A, Z, lManager);
 
  resA = resZ = fragZ = fragA = 0;
  mass = resMass = 0.0;
  secID = G4PhysicsModelCatalog::GetModelID("model_G4GEMChannelVI");
}

References A, cBarrier, evapMass, evapMass2, fProbability, fragA, fragZ, G4NuclearLevelData::GetInstance(), G4NuclearLevelData::GetLevelManager(), G4PhysicsModelCatalog::GetModelID(), G4NucleiProperties::GetNuclearMass(), G4NuclearLevelData::GetPairingCorrection(), mass, pairingCorrection, resA, resMass, resZ, secID, and Z.

◆ ~G4GEMChannelVI()

G4GEMChannelVI::~G4GEMChannelVI ( )

final

Definition at line 60 of file G4GEMChannelVI.cc.

{
  delete cBarrier;
  delete fProbability;
}

References cBarrier, and fProbability.

◆ G4GEMChannelVI() [2/2]

G4GEMChannelVI::G4GEMChannelVI ( const G4GEMChannelVI & right )

private

Member Function Documentation

◆ BreakUpChain()

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

virtualinherited

Reimplemented in G4UnstableFragmentBreakUp, and G4PhotonEvaporation.

Definition at line 66 of file G4VEvaporationChannel.cc.

{
  return false;
}

Referenced by G4VEvaporationChannel::BreakUpFragment().

◆ 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 G4GEMChannelVI::Dump ( ) const

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 126 of file G4GEMChannelVI.cc.

127{}

◆ EmittedFragment()

G4Fragment * G4GEMChannelVI::EmittedFragment ( G4Fragment * theNucleus )

finalvirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 96 of file G4GEMChannelVI.cc.

{
  // assumed, that TotalProbability(...) was already called
  // if value iz zero no possiblity to sample final state
  G4Fragment* evFragment = nullptr;
  G4LorentzVector lv0 = theNucleus->GetMomentum();
  if(resA <= 4 || fProbability->GetProbability() == 0.0) {
    G4double ekin = 
      std::max(0.5*(mass*mass - resMass*resMass + evapMass2)/mass 
               - evapMass, 0.0);
    G4LorentzVector lv(std::sqrt(ekin*(ekin + 2.0*evapMass))
                       *G4RandomDirection(), ekin + evapMass);
    lv.boost(lv0.boostVector());
    evFragment = new G4Fragment(A, Z, lv);
    lv0 -= lv;
  } else {
    evFragment = fProbability->SampleEvaporationFragment();
    G4LorentzVector lv = evFragment->GetMomentum();
    lv.boost(lv0.boostVector());
    evFragment->SetMomentum(lv);
    lv0 -= lv; 
  }
  if(evFragment != nullptr) { evFragment->SetCreatorModelID(secID); }
  theNucleus->SetZandA_asInt(resZ, resA);
  theNucleus->SetMomentum(lv0);
  theNucleus->SetCreatorModelID(secID);
  
  return evFragment;  
} 

References A, CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), evapMass, evapMass2, fProbability, G4RandomDirection(), G4Fragment::GetMomentum(), mass, G4INCL::Math::max(), resA, resMass, resZ, G4GEMProbabilityVI::SampleEvaporationFragment(), secID, G4Fragment::SetCreatorModelID(), G4Fragment::SetMomentum(), G4Fragment::SetZandA_asInt(), and Z.

◆ GetEmissionProbability()

G4double G4GEMChannelVI::GetEmissionProbability ( G4Fragment * theNucleus )

finalvirtual

Implements G4VEvaporationChannel.

Definition at line 66 of file G4GEMChannelVI.cc.

{
  fProbability->ResetProbability();
  fragZ = fragment->GetZ_asInt();
  fragA = fragment->GetA_asInt();
  resZ = fragZ - Z;
  resA = fragA - A;
  if(resA < A || resA < resZ || resZ < 0 || (resA == A && resZ < Z)) { 
    return 0.0; 
  }
 
  const G4double exc = fragment->GetExcitationEnergy();
  const G4double delta0 = 
    std::max(pairingCorrection->GetPairingCorrection(fragA, fragZ),0.0);
  if(exc < delta0) { return 0.0; }
 
  resMass = G4NucleiProperties::GetNuclearMass(resA, resZ);
  const G4double fragM = fragment->GetGroundStateMass() + exc;
  const G4double CB = cBarrier->GetCoulombBarrier(resA, resZ, exc);
  
  const G4double delta1 = 
    std::max(0.0,pairingCorrection->GetPairingCorrection(resA,resZ));
  if(fragM <= resMass + CB + delta1) { return 0.0; }
 
  fProbability->SetDecayKinematics(resZ, resA, resMass, fragM);
  G4double prob = fProbability->ComputeTotalProbability(*fragment, CB);
  //G4cout<<"G4EvaporationChannel: probability= "<< prob <<G4endl;
  return prob;
}

References A, cBarrier, G4GEMProbabilityVI::ComputeTotalProbability(), fProbability, fragA, fragZ, G4Fragment::GetA_asInt(), G4VCoulombBarrier::GetCoulombBarrier(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetGroundStateMass(), G4NucleiProperties::GetNuclearMass(), G4PairingCorrection::GetPairingCorrection(), G4Fragment::GetZ_asInt(), G4INCL::Math::max(), pairingCorrection, resA, G4VEmissionProbability::ResetProbability(), resMass, resZ, G4VEmissionProbability::SetDecayKinematics(), and Z.

◆ 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 G4GEMChannelVI::operator!= ( const G4GEMChannelVI & right ) const

private

◆ operator=()

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

private

◆ operator==()

G4bool G4GEMChannelVI::operator== ( const G4GEMChannelVI & 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{}

◆ UseSICB()

void G4VEvaporationChannel::UseSICB ( G4bool use )

inlineinherited

Definition at line 116 of file G4VEvaporationChannel.hh.

117{}

Field Documentation

◆ A

G4int G4GEMChannelVI::A

private

Definition at line 65 of file G4GEMChannelVI.hh.

Referenced by EmittedFragment(), G4GEMChannelVI(), GetEmissionProbability(), and mcscore.MCParticle::printout().

◆ cBarrier

const G4VCoulombBarrier* G4GEMChannelVI::cBarrier

private

Definition at line 61 of file G4GEMChannelVI.hh.

Referenced by G4GEMChannelVI(), GetEmissionProbability(), and ~G4GEMChannelVI().

◆ evapMass

G4double G4GEMChannelVI::evapMass

private

Definition at line 75 of file G4GEMChannelVI.hh.

Referenced by EmittedFragment(), and G4GEMChannelVI().

◆ evapMass2

G4double G4GEMChannelVI::evapMass2

private

Definition at line 76 of file G4GEMChannelVI.hh.

Referenced by EmittedFragment(), and G4GEMChannelVI().

◆ fProbability

G4GEMProbabilityVI* G4GEMChannelVI::fProbability

private

Definition at line 63 of file G4GEMChannelVI.hh.

Referenced by EmittedFragment(), G4GEMChannelVI(), GetEmissionProbability(), and ~G4GEMChannelVI().

◆ fragA

G4int G4GEMChannelVI::fragA

private

Definition at line 69 of file G4GEMChannelVI.hh.

Referenced by G4GEMChannelVI(), and GetEmissionProbability().

◆ fragZ

G4int G4GEMChannelVI::fragZ

private

Definition at line 70 of file G4GEMChannelVI.hh.

Referenced by G4GEMChannelVI(), and GetEmissionProbability().

◆ mass

G4double G4GEMChannelVI::mass

private

Definition at line 73 of file G4GEMChannelVI.hh.

Referenced by EmittedFragment(), and G4GEMChannelVI().

◆ OPTxs

G4int G4VEvaporationChannel::OPTxs

protectedinherited

Definition at line 93 of file G4VEvaporationChannel.hh.

Referenced by G4EvaporationChannel::GetEmissionProbability().

◆ pairingCorrection

const G4PairingCorrection* G4GEMChannelVI::pairingCorrection

private

Definition at line 62 of file G4GEMChannelVI.hh.

Referenced by G4GEMChannelVI(), and GetEmissionProbability().

◆ resA

G4int G4GEMChannelVI::resA

private

Definition at line 67 of file G4GEMChannelVI.hh.

Referenced by EmittedFragment(), G4GEMChannelVI(), and GetEmissionProbability().

◆ resMass

G4double G4GEMChannelVI::resMass

private

Definition at line 74 of file G4GEMChannelVI.hh.

Referenced by EmittedFragment(), G4GEMChannelVI(), and GetEmissionProbability().

◆ resZ

G4int G4GEMChannelVI::resZ

private

Definition at line 68 of file G4GEMChannelVI.hh.

Referenced by EmittedFragment(), G4GEMChannelVI(), and GetEmissionProbability().

◆ secID

G4int G4GEMChannelVI::secID

private

Definition at line 71 of file G4GEMChannelVI.hh.

Referenced by EmittedFragment(), and G4GEMChannelVI().

◆ useSICB

G4bool G4VEvaporationChannel::useSICB

protectedinherited

Definition at line 94 of file G4VEvaporationChannel.hh.

◆ Z

G4int G4GEMChannelVI::Z

private

Definition at line 66 of file G4GEMChannelVI.hh.

Referenced by EmittedFragment(), G4GEMChannelVI(), GetEmissionProbability(), and mcscore.MCParticle::printout().

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

source/processes/hadronic/models/de_excitation/gem_evaporation/include/G4GEMChannelVI.hh
source/processes/hadronic/models/de_excitation/gem_evaporation/src/G4GEMChannelVI.cc

Public Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4GEMChannelVI() [1/2]

◆ ~G4GEMChannelVI()

◆ G4GEMChannelVI() [2/2]

Member Function Documentation

◆ BreakUpChain()

◆ BreakUpFragment()

◆ Dump()

◆ EmittedFragment()

◆ GetEmissionProbability()

◆ GetLifeTime()

◆ Initialise()

◆ operator!=()

◆ operator=()

◆ operator==()

◆ RDMForced()

◆ SetICM()

◆ SetOPTxs()

◆ UseSICB()

Field Documentation

◆ A

◆ cBarrier

◆ evapMass

◆ evapMass2

◆ fProbability

◆ fragA

◆ fragZ

◆ mass

◆ OPTxs

◆ pairingCorrection

◆ resA

◆ resMass

◆ resZ

◆ secID

◆ useSICB

◆ Z