#include <G4BEChargedChannel.hh>

Inheritance diagram for G4BEChargedChannel:

Public Member Functions
	G4BEChargedChannel ()

virtual	~G4BEChargedChannel ()

virtual void	calculateProbability ()

virtual G4DynamicParticle *	emit ()=0

virtual G4double	coulombFactor ()=0

G4double	coulombFactorForProton ()

G4double	qmFactorForProton ()

G4double	qmFactorForAlpha ()

G4double	sampleKineticEnergy ()

Public Member Functions inherited from G4BertiniEvaporationChannel
	G4BertiniEvaporationChannel ()

virtual	~G4BertiniEvaporationChannel ()

void	setVerboseLevel (G4int verbose)

void	setNucleusA (G4int inputA)

void	setNucleusZ (G4int inputZ)

G4int	getParticleA ()

G4int	getParticleZ ()

void	setExcitationEnergy (G4double inputE)

void	setQ (G4double inputQ)

void	setPairingCorrection (G4int isCorrection)

G4double	getLevelDensityParameter ()

G4String	getName ()

virtual G4double	getProbability ()

virtual void	setProbability (G4double newProb)

virtual void	calculateProbability ()=0

virtual G4double	qmFactor ()

virtual G4double	getQ ()

virtual G4double	getCoulomb ()

virtual G4double	getThresh ()

virtual G4DynamicParticle *	emit ()=0

Protected Attributes
G4double	A

G4double	spin

Protected Attributes inherited from G4BertiniEvaporationChannel
G4String	name

G4int	verboseLevel

G4int	nucleusA

G4int	nucleusZ

G4int	particleA

G4int	particleZ

G4double	exmass

G4double	emissionProbability

G4double	rho

G4double	correction

G4double	excitationEnergy

G4int	spin

Additional Inherited Members
Protected Member Functions inherited from G4BertiniEvaporationChannel
G4double	Q (G4double a, G4double z)

G4double	pairingEnergyProtons (G4int A)

G4double	pairingEnergyNeutrons (G4int N)

G4double	cameron (G4double a, G4double z)

G4double	cameronShellCorrectionP (G4int Z)

G4double	cameronShellCorrectionN (G4int N)

void	isotropicCosines (G4double &, G4double &, G4double &)

Detailed Description

Definition at line 38 of file G4BEChargedChannel.hh.

Constructor & Destructor Documentation

◆ G4BEChargedChannel()

G4BEChargedChannel::G4BEChargedChannel ( )

Definition at line 35 of file G4BEChargedChannel.cc.

{
    verboseLevel = 0;
}

◆ ~G4BEChargedChannel()

G4BEChargedChannel::~G4BEChargedChannel ( )

virtual

Definition at line 41 of file G4BEChargedChannel.cc.

42{

43}

Member Function Documentation

◆ calculateProbability()

void G4BEChargedChannel::calculateProbability ( )

virtual

Implements G4BertiniEvaporationChannel.

Definition at line 46 of file G4BEChargedChannel.cc.

{
  G4int residualZ = nucleusZ - particleZ;
  G4int residualA = nucleusA - particleA;
 
//    Check if nucleus is too small, if this evaporation channel
//    leads to an impossible residual nucleus or if there is no enough
//    energy.
  if ( nucleusA < 2.0 * particleA || 
       nucleusZ < 2.0 * particleZ ||
       residualA <= residualZ || 
       excitationEnergy - getThresh() - correction < 0 )
    {
      if ( verboseLevel >= 6 )
    G4cout << "G4BEChargedChannel : calculateProbability for " << getName() << " = 0 " << G4endl;
      emissionProbability = 0;
      return;
    }
 
  // In HETC88 s-s0 was used in std::exp( s ), in which s0 was either 50 or
  // max(s_i), where i goes over all channels.
 
  G4double levelParam = getLevelDensityParameter();
  G4double slevel = 2 * std::sqrt( levelParam  * ( excitationEnergy - getThresh() - correction ) );
  G4double constant = A / 2 * ( 2 * spin + 1 ) * ( 1 + coulombFactor() );
  G4double eye1     = ( slevel*slevel - 3 * slevel + 3 ) / ( 4 * levelParam*levelParam ) * std::exp( slevel );
 
  emissionProbability = constant * std::pow( G4double(residualA), 0.6666666 ) * eye1;
 
  if ( verboseLevel >= 6 )
    G4cout << "G4BEChargedChannel : calculateProbability for " << getName() << G4endl
       << "                    res A = " << residualA << G4endl 
       << "                    res Z = " << residualZ << G4endl 
       << "                 c*(c_i+1) = "<< constant << G4endl
       << "                  qmfactor = "<< qmFactor() << G4endl
       << "             coulombfactor = "<< coulombFactor() << G4endl
       << "                        E = " << excitationEnergy << G4endl
       << "               correction = " << correction << G4endl
       << "                     eye1 = " << eye1 << G4endl
       << "               levelParam = " << levelParam << G4endl
       << "                   thresh = " << getThresh() << G4endl
       << "                        s = " << s << G4endl
       << "              probability = " << emissionProbability << G4endl;
 
  return;
}

◆ coulombFactor()

virtual G4double G4BEChargedChannel::coulombFactor ( )

pure virtual

Implemented in G4BEDeuteronChannel, G4BEHe3Channel, G4BEHe4Channel, G4BEProtonChannel, and G4BETritonChannel.

Referenced by calculateProbability().

◆ coulombFactorForProton()

G4double G4BEChargedChannel::coulombFactorForProton ( )

Definition at line 128 of file G4BEChargedChannel.cc.

{
  //  Coefficient c_p:s for empirical cross section formula are
  //  defined with the proton constant.  See Dostrovsky, Phys. Rev.,
  //  vol. 116, 1959.
  G4double t[7] = { 0.08 , 0 , -0.06 , -0.1 , -0.1 , -0.1 , -0.1 };
  G4int Z = nucleusZ - particleZ;
 
  if ( Z >= 70.0 ) return t[6];
  if ( Z <= 10.0 ) return t[0];
  
  // Linear interpolation
  G4int   n = G4int( 0.1 * Z + 1.0 ); 
  G4float x = ( 10 * n - Z ) * 0.1; 
  G4double ret_val =  x * t[n - 2] + ( 1.0 - x ) * t[n-1];
 
  return ret_val;
}

Referenced by G4BEDeuteronChannel::coulombFactor(), G4BEProtonChannel::coulombFactor(), and G4BETritonChannel::coulombFactor().

◆ emit()

virtual G4DynamicParticle * G4BEChargedChannel::emit ( )

pure virtual

Implements G4BertiniEvaporationChannel.

Implemented in G4BEDeuteronChannel, G4BEHe3Channel, G4BEHe4Channel, G4BEProtonChannel, and G4BETritonChannel.

◆ qmFactorForAlpha()

G4double G4BEChargedChannel::qmFactorForAlpha ( )

Definition at line 166 of file G4BEChargedChannel.cc.

{
//  Coefficient k_alpha for empirical cross section formula presented
//  in Dostrovsky, Phys. Rev., vol. 116, 1959
 
  G4double t[7] = {  0.77, 0.81, 0.85, 0.89, 0.93, 0.97,  1.00 };
  G4int Z = nucleusZ - particleZ;
 
  if ( Z >= 70.0 ) return t[6];
  if ( Z <= 10.0 ) return t[0];
 
  // Linear interpolation
  G4int   n = G4int( 0.1 * Z + 1.0 ); 
  G4float x = ( 10 * n - Z ) * 0.1; 
  return x * t[n - 2] + ( 1.0 - x ) * t[n-1];
}

Referenced by G4BEHe3Channel::qmFactor(), and G4BEHe4Channel::qmFactor().

◆ qmFactorForProton()

G4double G4BEChargedChannel::qmFactorForProton ( )

Definition at line 148 of file G4BEChargedChannel.cc.

{
  //  Coefficient k_p for empirical cross section formula are defined
  //  with the proton constant.  See Dostrovsky, Phys. Rev., vol. 116,
  //  1959
  G4double t[7] = {  0.36, 0.51, 0.60, 0.66, 0.68, 0.69, 0.69 };
  G4int Z = nucleusZ - particleZ;
 
  if ( Z >= 70.0 ) return t[6];
  if ( Z <= 10.0 ) return t[0];
 
  // Linear interpolation
  G4int   n = G4int( 0.1 * Z + 1.0 ); 
  G4float x = ( 10 * n - Z ) * 0.1; 
  return x * t[n - 2] + ( 1.0 - x ) * t[n-1];
}

Referenced by G4BEDeuteronChannel::qmFactor(), G4BEProtonChannel::qmFactor(), and G4BETritonChannel::qmFactor().

◆ sampleKineticEnergy()

G4double G4BEChargedChannel::sampleKineticEnergy ( )

Definition at line 94 of file G4BEChargedChannel.cc.

{
  G4double levelParam;
  levelParam = getLevelDensityParameter();
  
  const G4double xMax  = excitationEnergy - getThresh() - correction; // maximum number
  const G4double xProb = ( - 1 + std::sqrt ( 1 + 4 * levelParam * xMax ) ) / ( 2 * levelParam ); // most probable value
  const G4double maxProb = xProb * std::exp ( 2 * std::sqrt ( levelParam * ( xMax - xProb ) ) ); // maximum value of P(x)
 
  // Sample x according to density function P(x) with rejection method
  G4double r1;
  G4double r2;
  G4int koe=0;
  do
    {
      r1 = G4UniformRand() * xMax;
      r2 = G4UniformRand() * maxProb;
      koe++;
    }
  while (  r1 * std::exp ( 2 * std::sqrt ( levelParam * ( xMax - r1 ) ) )  < r2 );
 
//  G4cout << "Q ch " << koe << G4endl;
  G4double kineticEnergy = r1 + getCoulomb(); // add coulomb potential;
 
  if ( verboseLevel >= 10 )
    G4cout << " G4BENeutronChannel : sampleKineticEnergy() " << G4endl
       << "       kinetic n e = " << kineticEnergy << G4endl
       << "        levelParam = " << levelParam << G4endl
       << "             thresh= " << getThresh() << G4endl;
 
  return kineticEnergy;
}

Referenced by G4BEDeuteronChannel::emit(), G4BEHe3Channel::emit(), G4BEHe4Channel::emit(), G4BEProtonChannel::emit(), and G4BETritonChannel::emit().

Member Data Documentation

◆ A

G4double G4BEChargedChannel::A

protected

Definition at line 53 of file G4BEChargedChannel.hh.

Referenced by calculateProbability(), G4BEDeuteronChannel::G4BEDeuteronChannel(), G4BEHe3Channel::G4BEHe3Channel(), G4BEHe4Channel::G4BEHe4Channel(), G4BEProtonChannel::G4BEProtonChannel(), and G4BETritonChannel::G4BETritonChannel().

◆ spin

G4double G4BEChargedChannel::spin