#include <G4GEMChannel.hh>

Inheritance diagram for G4GEMChannel:

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

G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)

virtual void	Dump () const

virtual G4Fragment *	EmittedFragment (G4Fragment *theNucleus)

	G4GEMChannel (G4int theA, G4int theZ, const G4String &aName, G4GEMProbability *aEmissionStrategy)

virtual G4double	GetEmissionProbability (G4Fragment *theNucleus)

virtual G4double	GetLifeTime (G4Fragment *theNucleus)

virtual void	Initialise ()

virtual void	RDMForced (G4bool)

virtual void	SetICM (G4bool)

void	SetLevelDensityParameter (G4VLevelDensityParameter *aLevelDensity)

void	SetOPTxs (G4int opt)

void	UseSICB (G4bool use)

virtual	~G4GEMChannel ()

Protected Attributes
G4int	OPTxs

G4bool	useSICB

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

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

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

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

G4double	SampleKineticEnergy (const G4Fragment &fragment)

Detailed Description

Definition at line 47 of file G4GEMChannel.hh.

Constructor & Destructor Documentation

◆ G4GEMChannel() [1/2]

G4GEMChannel::G4GEMChannel	(	G4int	theA,
		G4int	theZ,
		const G4String &	aName,
		G4GEMProbability *	aEmissionStrategy
	)

explicit

Definition at line 49 of file G4GEMChannel.cc.

                                                                 :
  G4VEvaporationChannel(aName),
  A(theA),
  Z(theZ),
  EmissionProbability(0.0),
  MaximalKineticEnergy(-CLHEP::GeV),
  theEvaporationProbabilityPtr(aEmissionStrategy),
  secID(-1)
{ 
  theCoulombBarrierPtr = new G4GEMCoulombBarrier(theA, theZ);
  theEvaporationProbabilityPtr->SetCoulomBarrier(theCoulombBarrierPtr);
  theLevelDensityPtr = new G4EvaporationLevelDensityParameter;
  MyOwnLevelDensity = true;
  EvaporatedMass = G4NucleiProperties::GetNuclearMass(A, Z);
  ResidualMass = CoulombBarrier = 0.0;
  fG4pow = G4Pow::GetInstance(); 
  ResidualZ = ResidualA = 0;
  fNucData = G4NuclearLevelData::GetInstance();
  secID = G4PhysicsModelCatalog::GetModelID("model_G4GEMChannel");
}

References A, CoulombBarrier, EvaporatedMass, fG4pow, fNucData, G4Pow::GetInstance(), G4NuclearLevelData::GetInstance(), G4PhysicsModelCatalog::GetModelID(), G4NucleiProperties::GetNuclearMass(), MyOwnLevelDensity, ResidualA, ResidualMass, ResidualZ, secID, G4GEMProbability::SetCoulomBarrier(), theCoulombBarrierPtr, theEvaporationProbabilityPtr, theLevelDensityPtr, and Z.

◆ ~G4GEMChannel()

G4GEMChannel::~G4GEMChannel ( )

virtual

Definition at line 71 of file G4GEMChannel.cc.

{
  if (MyOwnLevelDensity) { delete theLevelDensityPtr; }
  delete theCoulombBarrierPtr;
}

References MyOwnLevelDensity, theCoulombBarrierPtr, and theLevelDensityPtr.

◆ G4GEMChannel() [2/2]

G4GEMChannel::G4GEMChannel ( const G4GEMChannel & right )

privatedelete

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

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 254 of file G4GEMChannel.cc.

{
  theEvaporationProbabilityPtr->Dump();
}

References G4GEMProbability::Dump(), and theEvaporationProbabilityPtr.

◆ EmittedFragment()

G4Fragment * G4GEMChannel::EmittedFragment ( G4Fragment * theNucleus )

virtual

Reimplemented from G4VEvaporationChannel.

Definition at line 130 of file G4GEMChannel.cc.

{
  G4Fragment* evFragment = 0;
  G4double evEnergy = SampleKineticEnergy(*theNucleus) + EvaporatedMass;
 
  G4ThreeVector momentum = G4RandomDirection()*
    std::sqrt((evEnergy - EvaporatedMass)*(evEnergy + EvaporatedMass));
  
  G4LorentzVector EvaporatedMomentum(momentum, evEnergy);
  G4LorentzVector ResidualMomentum = theNucleus->GetMomentum();
  EvaporatedMomentum.boost(ResidualMomentum.boostVector());
  
  evFragment = new G4Fragment(A, Z, EvaporatedMomentum);
  if ( evFragment != nullptr ) { evFragment->SetCreatorModelID(secID); }
  ResidualMomentum -= EvaporatedMomentum;
  theNucleus->SetZandA_asInt(ResidualZ, ResidualA);
  theNucleus->SetMomentum(ResidualMomentum);
  theNucleus->SetCreatorModelID(secID);
  
  return evFragment; 
} 

References A, CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), EvaporatedMass, G4RandomDirection(), G4Fragment::GetMomentum(), ResidualA, ResidualZ, SampleKineticEnergy(), secID, G4Fragment::SetCreatorModelID(), G4Fragment::SetMomentum(), G4Fragment::SetZandA_asInt(), and Z.

◆ GetEmissionProbability()

G4double G4GEMChannel::GetEmissionProbability ( G4Fragment * theNucleus )

virtual

Implements G4VEvaporationChannel.

Definition at line 77 of file G4GEMChannel.cc.

{
  G4int anA = fragment->GetA_asInt();
  G4int aZ  = fragment->GetZ_asInt();
  ResidualA = anA - A;
  ResidualZ = aZ - Z;
  /*
  G4cout << "G4GEMChannel: Z= " << Z << "  A= " << A 
     << " FragmentZ= " << aZ << " FragmentA= " << anA
     << " Zres= " << ResidualZ << " Ares= " << ResidualA 
     << G4endl; 
  */
  // We only take into account channels which are physically allowed
  EmissionProbability = 0.0;
 
  // Only channels which are physically allowed are taken into account 
  if (ResidualA >= ResidualZ && ResidualZ >= 0 && ResidualA >= A) {
  
    //Effective excitation energy
    G4double ExEnergy = fragment->GetExcitationEnergy()
      - fNucData->GetPairingCorrection(aZ, anA);
    if(ExEnergy > 0.0) { 
      ResidualMass = G4NucleiProperties::GetNuclearMass(ResidualA, ResidualZ);
      G4double FragmentMass = fragment->GetGroundStateMass();
      G4double Etot = FragmentMass + ExEnergy;
      // Coulomb Barrier calculation
      CoulombBarrier = 
    theCoulombBarrierPtr->GetCoulombBarrier(ResidualA,ResidualZ,ExEnergy);
      /*  
      G4cout << "Eexc(MeV)= " << ExEnergy/MeV 
         << " CoulBarrier(MeV)= " << CoulombBarrier/MeV << G4endl;
      */
      if(Etot > ResidualMass + EvaporatedMass + CoulombBarrier) {
  
    // Maximal Kinetic Energy
    MaximalKineticEnergy = ((Etot-ResidualMass)*(Etot+ResidualMass) 
      + EvaporatedMass*EvaporatedMass)/(2.0*Etot) 
      - EvaporatedMass - CoulombBarrier;
 
    //G4cout << "CBarrier(MeV)= " << CoulombBarrier/MeV << G4endl;
 
    if (MaximalKineticEnergy > 0.0) { 
      // Total emission probability for this channel
      EmissionProbability = theEvaporationProbabilityPtr->
        EmissionProbability(*fragment, MaximalKineticEnergy);
    }
      }
    }
  }   
  //G4cout << "Prob= " << EmissionProbability << G4endl;
  return EmissionProbability;
}

References A, CoulombBarrier, EmissionProbability, EvaporatedMass, fNucData, G4Fragment::GetA_asInt(), G4VCoulombBarrier::GetCoulombBarrier(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetGroundStateMass(), G4NucleiProperties::GetNuclearMass(), G4NuclearLevelData::GetPairingCorrection(), G4Fragment::GetZ_asInt(), MaximalKineticEnergy, ResidualA, ResidualMass, ResidualZ, theCoulombBarrierPtr, theEvaporationProbabilityPtr, 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 G4GEMChannel::operator!= ( const G4GEMChannel & right ) const

privatedelete

◆ operator=()

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

privatedelete

◆ operator==()

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

privatedelete

◆ RDMForced()

void G4VEvaporationChannel::RDMForced ( G4bool )

virtualinherited

Reimplemented in G4PhotonEvaporation.

Definition at line 58 of file G4VEvaporationChannel.cc.

59{}

◆ SampleKineticEnergy()

G4double G4GEMChannel::SampleKineticEnergy ( const G4Fragment & fragment )

private

Definition at line 152 of file G4GEMChannel.cc.

{
  G4double U = fragment.GetExcitationEnergy();
  
  G4double Alpha = theEvaporationProbabilityPtr->CalcAlphaParam(fragment);
  G4double Beta = theEvaporationProbabilityPtr->CalcBetaParam(fragment);
 
  //                       ***RESIDUAL***
  //JMQ (September 2009) the following quantities  refer to the RESIDUAL:
  G4double delta0 = fNucData->GetPairingCorrection(ResidualZ,ResidualA);
  G4double Ux = (2.5 + 150.0/ResidualA)*MeV;
  G4double Ex = Ux + delta0;
  G4double InitialLevelDensity;
  //                    ***end RESIDUAL ***
  
  //                       ***PARENT***
  //JMQ (September 2009) the following quantities   refer to the PARENT:
  
  G4double deltaCN = fNucData->GetPairingCorrection(fragment.GetZ_asInt(),
                            fragment.GetA_asInt());
  G4double aCN = theLevelDensityPtr->LevelDensityParameter(fragment.GetA_asInt(),
                               fragment.GetZ_asInt(),
                               U-deltaCN);   
  G4double UxCN = (2.5 + 150.0/G4double(fragment.GetA_asInt()))*MeV;
  G4double ExCN = UxCN + deltaCN;
  G4double TCN = 1.0/(std::sqrt(aCN/UxCN) - 1.5/UxCN);
  //                       ***end PARENT***
  
  //JMQ quantities calculated for  CN in InitialLevelDensity
  if ( U < ExCN ) 
    {
      G4double E0CN = ExCN - TCN*(G4Log(TCN/MeV) - G4Log(aCN*MeV)/4.0 
                  - 1.25*G4Log(UxCN/MeV) + 2.0*std::sqrt(aCN*UxCN));
      InitialLevelDensity = (pi/12.0)*G4Exp((U-E0CN)/TCN)/TCN;
    } 
  else 
    {
      G4double x  = U-deltaCN;
      G4double x1 = std::sqrt(aCN*x);
      InitialLevelDensity = (pi/12.0)*G4Exp(2*x1)/(x*std::sqrt(x1));
    }
  
  G4double Spin = theEvaporationProbabilityPtr->GetSpin();
  //JMQ  BIG BUG fixed: hbarc instead of hbar_Planck !!!!
  //     it was also fixed in total probability 
  G4double gg = (2.0*Spin+1.0)*EvaporatedMass/(pi2* hbarc*hbarc);
  //
  //JMQ  fix on Rb and geometrical cross sections according to Furihata's paper 
  //                      (JAERI-Data/Code 2001-105, p6)
  G4double Rb = 0.0; 
  if (A > 4) 
    {
      G4double Ad = fG4pow->Z13(ResidualA);
      G4double Aj = fG4pow->Z13(A);
      Rb = (1.12*(Aj + Ad) - 0.86*((Aj+Ad)/(Aj*Ad))+2.85)*fermi;
    }
  else if (A>1)
    {
      G4double Ad = fG4pow->Z13(ResidualA);
      G4double Aj = fG4pow->Z13(A);
      Rb=1.5*(Aj+Ad)*fermi;
    }
  else 
    {
      G4double Ad = fG4pow->Z13(ResidualA);
      Rb = 1.5*Ad*fermi;
    }
  G4double GeometricalXS = pi*Rb*Rb; 
    
  G4double ConstantFactor = gg*GeometricalXS*Alpha*pi/(InitialLevelDensity*12);
  //JMQ : this is the assimptotic maximal kinetic energy of the ejectile 
  //      (obtained by energy-momentum conservation). 
  //      In general smaller than U-theSeparationEnergy 
  G4double theEnergy = MaximalKineticEnergy + CoulombBarrier;
  G4double KineticEnergy;
  G4double Probability;
 
  for(G4int i=0; i<100; ++i) {
    KineticEnergy =  CoulombBarrier + G4UniformRand()*(MaximalKineticEnergy);
    G4double edelta = theEnergy-KineticEnergy-delta0;
    Probability = ConstantFactor*(KineticEnergy + Beta);
    G4double a = 
      theLevelDensityPtr->LevelDensityParameter(ResidualA,ResidualZ,edelta);
    G4double T = 1.0/(std::sqrt(a/Ux) - 1.5/Ux);
    //JMQ fix in units
    
    if (theEnergy - KineticEnergy < Ex) {
      G4double E0 = Ex - T*(G4Log(T) - G4Log(a)*0.25
                - 1.25*G4Log(Ux) + 2.0*std::sqrt(a*Ux));
      Probability *= G4Exp((theEnergy-KineticEnergy-E0)/T)/T;
    } else {
      G4double e2 = edelta*edelta;
      Probability *= 
    G4Exp(2*std::sqrt(a*edelta) - 0.25*G4Log(a*edelta*e2*e2));
    }
    if(EmissionProbability*G4UniformRand() <= Probability) { break; }
  }
    
  return KineticEnergy;
} 

References A, Alpha, G4GEMProbability::CalcAlphaParam(), G4GEMProbability::CalcBetaParam(), CoulombBarrier, e2, EmissionProbability, EvaporatedMass, fermi, fG4pow, fNucData, G4Exp(), G4Log(), G4UniformRand, G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), G4NuclearLevelData::GetPairingCorrection(), G4GEMProbability::GetSpin(), G4Fragment::GetZ_asInt(), source.hepunit::hbarc, G4VLevelDensityParameter::LevelDensityParameter(), MaximalKineticEnergy, MeV, pi, pi2, ResidualA, ResidualZ, theEvaporationProbabilityPtr, theLevelDensityPtr, and G4Pow::Z13().

Referenced by EmittedFragment().

◆ SetICM()

void G4VEvaporationChannel::SetICM ( G4bool )

virtualinherited

Reimplemented in G4PhotonEvaporation.

Definition at line 55 of file G4VEvaporationChannel.cc.

56{}

Referenced by G4NeutronRadCapture::InitialiseModel().

◆ SetLevelDensityParameter()

void G4GEMChannel::SetLevelDensityParameter ( G4VLevelDensityParameter * aLevelDensity )

inline

Definition at line 63 of file G4GEMChannel.hh.

  {
    if (MyOwnLevelDensity) { delete theLevelDensityPtr; }
    theLevelDensityPtr = aLevelDensity;
    MyOwnLevelDensity = false;
  }

References MyOwnLevelDensity, and theLevelDensityPtr.

◆ 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 G4GEMChannel::A

private

Definition at line 85 of file G4GEMChannel.hh.

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

◆ CoulombBarrier

G4double G4GEMChannel::CoulombBarrier

private

Definition at line 98 of file G4GEMChannel.hh.

Referenced by G4GEMChannel(), GetEmissionProbability(), and SampleKineticEnergy().

◆ EmissionProbability

G4double G4GEMChannel::EmissionProbability

private

Definition at line 99 of file G4GEMChannel.hh.

Referenced by GetEmissionProbability(), and SampleKineticEnergy().

◆ EvaporatedMass

G4double G4GEMChannel::EvaporatedMass

private

Definition at line 96 of file G4GEMChannel.hh.

Referenced by EmittedFragment(), G4GEMChannel(), GetEmissionProbability(), and SampleKineticEnergy().

◆ fG4pow

G4Pow* G4GEMChannel::fG4pow

private

Definition at line 104 of file G4GEMChannel.hh.

Referenced by G4GEMChannel(), and SampleKineticEnergy().

◆ fNucData

G4NuclearLevelData* G4GEMChannel::fNucData

private

Definition at line 116 of file G4GEMChannel.hh.

Referenced by G4GEMChannel(), GetEmissionProbability(), and SampleKineticEnergy().

◆ MaximalKineticEnergy

G4double G4GEMChannel::MaximalKineticEnergy

private

Definition at line 102 of file G4GEMChannel.hh.

Referenced by GetEmissionProbability(), and SampleKineticEnergy().

◆ MyOwnLevelDensity

G4bool G4GEMChannel::MyOwnLevelDensity

private

Definition at line 110 of file G4GEMChannel.hh.

Referenced by G4GEMChannel(), SetLevelDensityParameter(), and ~G4GEMChannel().

◆ OPTxs

G4int G4VEvaporationChannel::OPTxs

protectedinherited

Definition at line 93 of file G4VEvaporationChannel.hh.

Referenced by G4EvaporationChannel::GetEmissionProbability().

◆ ResidualA

G4int G4GEMChannel::ResidualA

private

Definition at line 91 of file G4GEMChannel.hh.

Referenced by EmittedFragment(), G4GEMChannel(), GetEmissionProbability(), and SampleKineticEnergy().

◆ ResidualMass

G4double G4GEMChannel::ResidualMass

private

Definition at line 97 of file G4GEMChannel.hh.

Referenced by G4GEMChannel(), and GetEmissionProbability().

◆ ResidualZ

G4int G4GEMChannel::ResidualZ

private

Definition at line 94 of file G4GEMChannel.hh.

Referenced by EmittedFragment(), G4GEMChannel(), GetEmissionProbability(), and SampleKineticEnergy().

◆ secID

G4int G4GEMChannel::secID

private

Definition at line 119 of file G4GEMChannel.hh.

Referenced by EmittedFragment(), and G4GEMChannel().

◆ theCoulombBarrierPtr

G4VCoulombBarrier* G4GEMChannel::theCoulombBarrierPtr

private

Definition at line 114 of file G4GEMChannel.hh.

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

◆ theEvaporationProbabilityPtr

G4GEMProbability* G4GEMChannel::theEvaporationProbabilityPtr

private

Definition at line 107 of file G4GEMChannel.hh.

Referenced by Dump(), G4GEMChannel(), GetEmissionProbability(), and SampleKineticEnergy().

◆ theLevelDensityPtr

G4VLevelDensityParameter* G4GEMChannel::theLevelDensityPtr

private

Definition at line 111 of file G4GEMChannel.hh.

Referenced by G4GEMChannel(), SampleKineticEnergy(), SetLevelDensityParameter(), and ~G4GEMChannel().

◆ useSICB

G4bool G4VEvaporationChannel::useSICB

protectedinherited

Definition at line 94 of file G4VEvaporationChannel.hh.

◆ Z

G4int G4GEMChannel::Z

private

Definition at line 88 of file G4GEMChannel.hh.

Referenced by EmittedFragment(), G4GEMChannel(), 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/G4GEMChannel.hh
source/processes/hadronic/models/de_excitation/gem_evaporation/src/G4GEMChannel.cc

Private Attributes
G4int	A

G4double	CoulombBarrier

G4double	EmissionProbability

G4double	EvaporatedMass

G4Pow *	fG4pow

G4NuclearLevelData *	fNucData

G4double	MaximalKineticEnergy

G4bool	MyOwnLevelDensity

G4int	ResidualA

G4double	ResidualMass

G4int	ResidualZ

G4int	secID

G4VCoulombBarrier *	theCoulombBarrierPtr

G4GEMProbability *	theEvaporationProbabilityPtr

G4VLevelDensityParameter *	theLevelDensityPtr

G4int	Z

Public Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4GEMChannel() [1/2]

◆ ~G4GEMChannel()

◆ G4GEMChannel() [2/2]

Member Function Documentation

◆ BreakUpChain()

◆ BreakUpFragment()

◆ Dump()

◆ EmittedFragment()

◆ GetEmissionProbability()

◆ GetLifeTime()

◆ Initialise()

◆ operator!=()

◆ operator=()

◆ operator==()

◆ RDMForced()

◆ SampleKineticEnergy()

◆ SetICM()

◆ SetLevelDensityParameter()

◆ SetOPTxs()

◆ UseSICB()

Field Documentation

◆ A

◆ CoulombBarrier

◆ EmissionProbability

◆ EvaporatedMass

◆ fG4pow

◆ fNucData

◆ MaximalKineticEnergy

◆ MyOwnLevelDensity

◆ OPTxs

◆ ResidualA

◆ ResidualMass

◆ ResidualZ

◆ secID

◆ theCoulombBarrierPtr

◆ theEvaporationProbabilityPtr

◆ theLevelDensityPtr

◆ useSICB

◆ Z