G4NuclearAbrasionGeometry Class Reference

#include <G4NuclearAbrasionGeometry.hh>

Public Member Functions
	G4NuclearAbrasionGeometry (G4double AP, G4double AT, G4double r)
	~G4NuclearAbrasionGeometry ()
void	SetPeripheralThreshold (G4double)
G4double	GetPeripheralThreshold ()
G4double	F ()
G4double	P ()
G4double	GetExcitationEnergyOfProjectile ()
G4double	GetExcitationEnergyOfTarget ()

Detailed Description

Definition at line 67 of file G4NuclearAbrasionGeometry.hh.

Constructor & Destructor Documentation

G4NuclearAbrasionGeometry::G4NuclearAbrasionGeometry	(	G4double	AP,
		G4double	AT,
		G4double	r
	)

Definition at line 73 of file G4NuclearAbrasionGeometry.cc.

References G4WilsonRadius::GetWilsonRadius(), and python.hepunit::MeV.

Referenced by GetExcitationEnergyOfTarget().

{
//
//
// Initialise variables for interaction geometry.
//
  G4WilsonRadius aR;
  AP = AP1;
  AT = AT1;
  rP = aR.GetWilsonRadius(AP);
  rT = aR.GetWilsonRadius(AT);
  r  = r1;
  n  = rP / (rP + rT);
  b  = r / (rP + rT);
  m  = rT / rP;
  Q  = (1.0 - b)/n;
  S  = Q * Q;
  T  = S * Q;
  R  = std::sqrt(m*n);
  U  = 1.0/m - 2.0;
//
//
// Initialise the threshold radius-ratio at which interactions are considered
// peripheral or central.
//  
  rth = 2.0/3.0;
  B   = 10.0 * MeV;
}

G4NuclearAbrasionGeometry::~G4NuclearAbrasionGeometry

(

)

Definition at line 104 of file G4NuclearAbrasionGeometry.cc.

{;}

Member Function Documentation

G4double G4NuclearAbrasionGeometry::F

(

)

Definition at line 147 of file G4NuclearAbrasionGeometry.cc.

Referenced by G4WilsonAbrasionModel::ApplyYourself(), GetExcitationEnergyOfProjectile(), and GetExcitationEnergyOfTarget().

{
//
//
// Initialise the value for F, then determine the actual value depending upon
// whether the projectile is larger or smaller than the target and these radii
// in relation to the impact parameter.
//
  G4double valueF = 0.0;

  if (rT > rP)
  {
    if (rT-rP<=r && r<=rT+rP) valueF = 0.75*R*S - 0.125*(3.0*R-1.0)*T;
    else                      valueF = 1.0;
  }
  else
  {
    if (rP-rT<=r && r<=rP+rT) valueF = 0.75*R*S - 0.125*(3.0*std::sqrt(n/m)-
      (1.0-std::pow(1.0-m*m,3.0/2.0))*std::sqrt(1.0-std::pow(1.0-m,2.0))/std::pow(m,3.0))*T;
    else                      valueF = (1.0-std::pow(1.0-m*m,3.0/2.0))*std::sqrt(1.0-b*b/n/n);
  }

  if (!(valueF <= 1.0 && valueF>= 0.0))
  {
    if (valueF > 1.0) valueF = 1.0;
    else         valueF = 0.0;
  }
  return valueF;
}

G4double G4NuclearAbrasionGeometry::GetExcitationEnergyOfProjectile

(

)

Definition at line 178 of file G4NuclearAbrasionGeometry.cc.

References F(), python.hepunit::fermi, python.hepunit::MeV, P(), and python.hepunit::pi.

Referenced by G4WilsonAbrasionModel::ApplyYourself().

{
  G4double F1 = F();
  G4double P1 = P();
  G4double Es = 0.0;

  Es = 0.95 * MeV * 4.0 * pi * rP*rP/fermi/fermi *
       (1.0+P1-std::pow(1.0-F1,2.0/3.0));
//  if (rT < rP && r < rP-rT)
  if ((r-rP)/rT < rth)
  {
    G4double omega = 0.0;
    if      (AP < 12.0)  omega = 1500.0;
    else if (AP <= 16.0) omega = 1500.0 - 320.0*(AP-12.0);
    Es *= 1.0 + F1*(5.0+omega*F1*F1);
  }
  
  if (Es < 0.0) 
    Es = 0.0;
  else if (Es > B * AP)
    Es = B * AP;
  return Es;
}

G4double G4NuclearAbrasionGeometry::GetExcitationEnergyOfTarget

(

)

Definition at line 203 of file G4NuclearAbrasionGeometry.cc.

References F(), python.hepunit::fermi, G4NuclearAbrasionGeometry(), python.hepunit::MeV, P(), and python.hepunit::pi.

Referenced by G4WilsonAbrasionModel::ApplyYourself().

{
  // This member function declares a new G4NuclearAbrasionGeometry object 
  // but with the projectile and target exchanged to determine the values
  // for F and P.  Determination of the excess surface area and excitation
  // energy is as above.

  G4NuclearAbrasionGeometry* revAbrasionGeometry =
    new G4NuclearAbrasionGeometry(AT, AP, r);
  G4double F1 = revAbrasionGeometry->F();
  G4double P1 = revAbrasionGeometry->P();
  G4double Es = 0.0;

  Es = 0.95 * MeV * 4.0 * pi * rT*rT/fermi/fermi *
       (1.0+P1-std::pow(1.0-F1,2.0/3.0));

//  if (rP < rT && r < rT-rP)
  if ((r-rT)/rP < rth) {
    G4double omega = 0.0;
    if      (AT < 12.0)  omega = 1500.0;
    else if (AT <= 16.0) omega = 1500.0 - 320.0*(AT-12.0);
    Es *= 1.0 + F1*(5.0+omega*F1*F1);
  }
  
  if (Es < 0.0)
    Es = 0.0;
  else if (Es > B * AT)
    Es = B * AT;

  delete revAbrasionGeometry;

  return Es;
}

G4double G4NuclearAbrasionGeometry::GetPeripheralThreshold

(

)

Definition at line 112 of file G4NuclearAbrasionGeometry.cc.

  {return rth;}

G4double G4NuclearAbrasionGeometry::P

(

)

Definition at line 116 of file G4NuclearAbrasionGeometry.cc.

Referenced by GetExcitationEnergyOfProjectile(), and GetExcitationEnergyOfTarget().

{
//
//
// Initialise the value for P, then determine the actual value depending upon
// whether the projectile is larger or smaller than the target and these radii
// in relation to the impact parameter.
//
  G4double valueP = 0.0;

  if (rT > rP)
  {
    if (rT-rP<=r && r<=rT+rP) valueP = 0.125*R*U*S - 0.125*(0.5*R*U+1.0)*T;
    else                      valueP = -1.0;
  }
  else
  {
    if (rP-rT<=r && r<=rP+rT) valueP = 0.125*R*U*S - 0.125*(0.5*std::sqrt(n/m)*U-
      (std::sqrt(1.0-m*m)/n - 1.0)*std::sqrt((2.0-m)/std::pow(m,5.0)))*T;
    else                      valueP = (std::sqrt(1.0-m*m)/n-1.0)*std::sqrt(1.0-b*b/n/n);
  }

  if (!(valueP <= 1.0 && valueP>= -1.0))
  {
    if (valueP > 1.0) valueP =  1.0;
    else         valueP = -1.0;
  }
  return valueP;
}

void G4NuclearAbrasionGeometry::SetPeripheralThreshold

(

G4double

rth1

)

Definition at line 108 of file G4NuclearAbrasionGeometry.cc.

  {if (rth1 > 0.0 && rth1 <= 1.0) rth = rth1;}

---

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

• G4NuclearAbrasionGeometry.hh
• G4NuclearAbrasionGeometry.cc