#include <G4ParticleHPKallbachMannSyst.hh>

Public Member Functions
G4double	A (G4double anEnergy)

	G4ParticleHPKallbachMannSyst (G4double aCompoundFraction, G4double anIncidentEnergy, G4double anIncidentMass, G4double aProductEnergy, G4double aProductMass, G4double aResidualMass, G4int aResidualA, G4int aResidualZ, G4double aTargetMass, G4int aTargetA, G4int aTargetZ, G4int aProjectileA, G4int aProjectileZ, G4int aProductA, G4int aProductZ)

G4double	GetKallbachZero (G4double anEnergy)

G4double	Kallbach (G4double cosTh, G4double anEnergy)

G4double	Sample (G4double anEnergy)

G4double	SeparationEnergy (G4int Ac, G4int Nc, G4int AA, G4int ZA, G4int Abinding, G4int Zbinding)

	~G4ParticleHPKallbachMannSyst ()

Detailed Description

Definition at line 36 of file G4ParticleHPKallbachMannSyst.hh.

Constructor & Destructor Documentation

◆ G4ParticleHPKallbachMannSyst()

G4ParticleHPKallbachMannSyst::G4ParticleHPKallbachMannSyst	(	G4double	aCompoundFraction,
		G4double	anIncidentEnergy,
		G4double	anIncidentMass,
		G4double	aProductEnergy,
		G4double	aProductMass,
		G4double	aResidualMass,
		G4int	aResidualA,
		G4int	aResidualZ,
		G4double	aTargetMass,
		G4int	aTargetA,
		G4int	aTargetZ,
		G4int	aProjectileA,
		G4int	aProjectileZ,
		G4int	aProductA,
		G4int	aProductZ
	)

inline

Definition at line 40 of file G4ParticleHPKallbachMannSyst.hh.

  {
    theCompoundFraction = aCompoundFraction;
    theIncidentEnergy = anIncidentEnergy;
    theIncidentMass = anIncidentMass;
    theProductEnergy = aProductEnergy;
    theProductMass = aProductMass;
    theResidualMass = aResidualMass;
    theResidualA = aResidualA;
    theResidualZ = aResidualZ;
    theTargetMass = aTargetMass;
    theTargetA = aTargetA;
    theTargetZ = aTargetZ;
    theProjectileA=aProjectileA;
    theProjectileZ=aProjectileZ;
    theProductA=aProductA;
    theProductZ=aProductZ;
  }

References theCompoundFraction, theIncidentEnergy, theIncidentMass, theProductA, theProductEnergy, theProductMass, theProductZ, theProjectileA, theProjectileZ, theResidualA, theResidualMass, theResidualZ, theTargetA, theTargetMass, and theTargetZ.

◆ ~G4ParticleHPKallbachMannSyst()

G4ParticleHPKallbachMannSyst::~G4ParticleHPKallbachMannSyst ( )

inline

Definition at line 64 of file G4ParticleHPKallbachMannSyst.hh.

64{};

Member Function Documentation

◆ A()

G4double G4ParticleHPKallbachMannSyst::A ( G4double anEnergy )

Definition at line 100 of file G4ParticleHPKallbachMannSyst.cc.

{ 
  G4double result;
  G4double C1 = 0.04/MeV;
  G4double C2 = 1.8E-6/(MeV*MeV*MeV);
  G4double C3 = 6.7E-7/(MeV*MeV*MeV*MeV);
  
  G4double epsa = anEnergy*theTargetMass/(theTargetMass+theIncidentMass);
  G4int Ac = theTargetA+theProjectileA; 
  G4int Nc = Ac - theTargetZ-theProjectileZ;
  G4int AA = theTargetA;
  G4int ZA = theTargetZ;
  G4double ea = epsa+SeparationEnergy(Ac, Nc, AA, ZA,theProjectileA,theProjectileZ);
  G4double Et1 = 130*MeV;
  G4double R1 = std::min(ea, Et1);
  // theProductEnergy is still in CMS!!!
  G4double epsb = theProductEnergy*(theProductMass+theResidualMass)/theResidualMass;
  G4int AB = theResidualA;
  G4int ZB = theResidualZ;
  G4double eb = epsb+SeparationEnergy(Ac, Nc, AB, ZB,theProductA, theProductZ);
  G4double X1 = R1*eb/ea;
  G4double Et3 = 41*MeV;
  G4double R3 = std::min(ea, Et3);
  G4double X3 = R3*eb/ea;
 
  G4double Ma=1;
  G4double mb=1;
  if(theProjectileA==1 || (theProjectileZ==1 && theProjectileA==2)){Ma=1;}//neutron,proton,deuteron
  else if(theProjectileA==4 && theProjectileZ==2){Ma=0;}//alpha
  else if(theProjectileA==3 && (theProjectileZ==1 || theProjectileZ==2)){Ma=0.5;}//tritum,He3 : set intermediate value
  else
  {
    throw G4HadronicException(__FILE__, __LINE__, "Severe error in the sampling of Kallbach-Mann Systematics");
  }
  if(theProductA==1 && theProductZ==0){mb=1./2.;}//neutron
  else if(theProductA==4 && theProductZ==2){mb=2;}//alpha
  else{mb=1;}
 
  result = C1*X1 + C2*G4Pow::GetInstance()->powN(X1, 3) + C3*Ma*mb*G4Pow::GetInstance()->powN(X3, 4);
  return result;
 
}

References anonymous_namespace{paraMaker.cc}::AA, C1, C2, C3, G4Pow::GetInstance(), MeV, G4INCL::Math::min(), G4Pow::powN(), SeparationEnergy(), theIncidentMass, theProductA, theProductEnergy, theProductMass, theProductZ, theProjectileA, theProjectileZ, theResidualA, theResidualMass, theResidualZ, theTargetA, theTargetMass, and theTargetZ.

Referenced by GetKallbachZero(), Kallbach(), and mcscore.MCParticle::printout().

◆ GetKallbachZero()

G4double G4ParticleHPKallbachMannSyst::GetKallbachZero ( G4double anEnergy )

Definition at line 86 of file G4ParticleHPKallbachMannSyst.cc.

{
  G4double result;
  //delta 2.0e-16 in not good.
  //delta 4.0e-16 is OK 
  //safety factor of 2
  G4double delta = 8.0e-16;
  if ( std::abs (theCompoundFraction - 1 ) < delta ) { 
     theCompoundFraction = 1.0-delta;   
  } 
  result = 0.5 * (1./A(anEnergy)) * G4Log((1-theCompoundFraction)/(1+theCompoundFraction));
  return result;
}

References A(), G4Log(), and theCompoundFraction.

Referenced by Sample().

◆ Kallbach()

G4double G4ParticleHPKallbachMannSyst::Kallbach	(	G4double	cosTh,
		G4double	anEnergy
	)

Definition at line 76 of file G4ParticleHPKallbachMannSyst.cc.

{
  // Kallbach-Mann systematics without normalization.
  G4double result;
  G4double theX = A(anEnergy)*cosTh;
  result = 0.5*(G4Exp( theX)*(1+theCompoundFraction)
               +G4Exp(-theX)*(1-theCompoundFraction));
  return result;
}

References A(), G4Exp(), and theCompoundFraction.

Referenced by Sample().

◆ Sample()

G4double G4ParticleHPKallbachMannSyst::Sample ( G4double anEnergy )

Definition at line 44 of file G4ParticleHPKallbachMannSyst.cc.

{
  G4double result=0.;
  
  G4double zero = GetKallbachZero(anEnergy);
  if(zero>1) zero=1.;
  if(zero<-1)zero=-1.;
  G4double max = Kallbach(zero, anEnergy);
  G4double upper = Kallbach(1., anEnergy);
  G4double lower = Kallbach(-1., anEnergy);
  if(upper>max) max=upper;
  if(lower>max) max=lower;
  G4double value, random;
 
  G4int icounter=0;
  G4int icounter_max=1024;
  do
  {
    icounter++;
    if ( icounter > icounter_max ) {
       G4cout << "Loop-counter exceeded the threshold value at " << __LINE__ << "th line of " << __FILE__ << "." << G4endl;
       break;
    }
    result = 2.*G4UniformRand()-1;  
    value = Kallbach(result, anEnergy)/max;
    random = G4UniformRand();
  }
  while(random>value); // Loop checking, 11.05.2015, T. Koi
  
  return result;
}

References G4cout, G4endl, G4UniformRand, GetKallbachZero(), Kallbach(), G4INCL::Math::max(), and anonymous_namespace{G4CascadeDeexciteBase.cc}::zero.

Referenced by G4ParticleHPContAngularPar::Sample().

◆ SeparationEnergy()

G4double G4ParticleHPKallbachMannSyst::SeparationEnergy	(	G4int	Ac,
		G4int	Nc,
		G4int	AA,
		G4int	ZA,
		G4int	Abinding,
		G4int	Zbinding
	)

Definition at line 143 of file G4ParticleHPKallbachMannSyst.cc.

{
  G4double result;
  G4int NA = AA-ZA;
  G4int Zc = Ac-Nc;
  result = 15.68*(Ac-AA);
  result += -28.07*((Nc-Zc)*(Nc-Zc)/(G4double)Ac - (NA-ZA)*(NA-ZA)/(G4double)AA);
  result += -18.56*(G4Pow::GetInstance()->A23(G4double(Ac)) - G4Pow::GetInstance()->A23(G4double(AA)));
  result +=  33.22*((Nc-Zc)*(Nc-Zc)/G4Pow::GetInstance()->powA(G4double(Ac), 4./3.) - (NA-ZA)*(NA-ZA)/G4Pow::GetInstance()->powA(G4double(AA), 4./3.));
  result += -0.717*(Zc*Zc/G4Pow::GetInstance()->A13(G4double(Ac))-ZA*ZA/G4Pow::GetInstance()->A13(G4double(AA)));
  result +=  1.211*(Zc*Zc/(G4double)Ac-ZA*ZA/(G4double)AA);
  G4double totalBinding(0);
  if(Zbinding==0&&Abinding==1) totalBinding=0;
  if(Zbinding==1&&Abinding==1) totalBinding=0;
  if(Zbinding==1&&Abinding==2) totalBinding=2.224596;
  if(Zbinding==1&&Abinding==3) totalBinding=8.481798;
  if(Zbinding==2&&Abinding==3) totalBinding=7.718043;
  if(Zbinding==2&&Abinding==4) totalBinding=28.29566;
  result += -totalBinding;
  result *= MeV;
  return result;
}

References G4Pow::A13(), G4Pow::A23(), anonymous_namespace{paraMaker.cc}::AA, G4Pow::GetInstance(), MeV, and G4Pow::powA().

Referenced by A().

Field Documentation

◆ theCompoundFraction

G4double G4ParticleHPKallbachMannSyst::theCompoundFraction

private

Definition at line 78 of file G4ParticleHPKallbachMannSyst.hh.

Referenced by G4ParticleHPKallbachMannSyst(), GetKallbachZero(), and Kallbach().

◆ theIncidentEnergy

G4double G4ParticleHPKallbachMannSyst::theIncidentEnergy

private

Definition at line 79 of file G4ParticleHPKallbachMannSyst.hh.

Referenced by G4ParticleHPKallbachMannSyst().

◆ theIncidentMass

G4double G4ParticleHPKallbachMannSyst::theIncidentMass

private

Definition at line 80 of file G4ParticleHPKallbachMannSyst.hh.

Referenced by A(), and G4ParticleHPKallbachMannSyst().

◆ theProductA

G4int G4ParticleHPKallbachMannSyst::theProductA

private

Definition at line 91 of file G4ParticleHPKallbachMannSyst.hh.

Referenced by A(), and G4ParticleHPKallbachMannSyst().

◆ theProductEnergy

G4double G4ParticleHPKallbachMannSyst::theProductEnergy

private

Definition at line 81 of file G4ParticleHPKallbachMannSyst.hh.

Referenced by A(), and G4ParticleHPKallbachMannSyst().

◆ theProductMass

G4double G4ParticleHPKallbachMannSyst::theProductMass

private

Definition at line 82 of file G4ParticleHPKallbachMannSyst.hh.

Referenced by A(), and G4ParticleHPKallbachMannSyst().

◆ theProductZ

G4int G4ParticleHPKallbachMannSyst::theProductZ

private

Definition at line 92 of file G4ParticleHPKallbachMannSyst.hh.

Referenced by A(), and G4ParticleHPKallbachMannSyst().

◆ theProjectileA

G4int G4ParticleHPKallbachMannSyst::theProjectileA

private

Definition at line 89 of file G4ParticleHPKallbachMannSyst.hh.

Referenced by A(), and G4ParticleHPKallbachMannSyst().

◆ theProjectileZ

G4int G4ParticleHPKallbachMannSyst::theProjectileZ

private

Definition at line 90 of file G4ParticleHPKallbachMannSyst.hh.

Referenced by A(), and G4ParticleHPKallbachMannSyst().

◆ theResidualA

G4int G4ParticleHPKallbachMannSyst::theResidualA

private

Definition at line 85 of file G4ParticleHPKallbachMannSyst.hh.

Referenced by A(), and G4ParticleHPKallbachMannSyst().

◆ theResidualMass

G4double G4ParticleHPKallbachMannSyst::theResidualMass

private

Definition at line 83 of file G4ParticleHPKallbachMannSyst.hh.

Referenced by A(), and G4ParticleHPKallbachMannSyst().

◆ theResidualZ

G4int G4ParticleHPKallbachMannSyst::theResidualZ

private

Definition at line 86 of file G4ParticleHPKallbachMannSyst.hh.

Referenced by A(), and G4ParticleHPKallbachMannSyst().

◆ theTargetA

G4int G4ParticleHPKallbachMannSyst::theTargetA

private

Definition at line 87 of file G4ParticleHPKallbachMannSyst.hh.

Referenced by A(), and G4ParticleHPKallbachMannSyst().

◆ theTargetMass

G4double G4ParticleHPKallbachMannSyst::theTargetMass

private

Definition at line 84 of file G4ParticleHPKallbachMannSyst.hh.

Referenced by A(), and G4ParticleHPKallbachMannSyst().

◆ theTargetZ

G4int G4ParticleHPKallbachMannSyst::theTargetZ

private

Definition at line 88 of file G4ParticleHPKallbachMannSyst.hh.

Referenced by A(), and G4ParticleHPKallbachMannSyst().

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

source/processes/hadronic/models/particle_hp/include/G4ParticleHPKallbachMannSyst.hh
source/processes/hadronic/models/particle_hp/src/G4ParticleHPKallbachMannSyst.cc

Private Attributes
G4double	theCompoundFraction

G4double	theIncidentEnergy

G4double	theIncidentMass

G4int	theProductA

G4double	theProductEnergy

G4double	theProductMass

G4int	theProductZ

G4int	theProjectileA

G4int	theProjectileZ

G4int	theResidualA

G4double	theResidualMass

G4int	theResidualZ

G4int	theTargetA

G4double	theTargetMass

G4int	theTargetZ

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4ParticleHPKallbachMannSyst()

◆ ~G4ParticleHPKallbachMannSyst()

Member Function Documentation

◆ A()

◆ GetKallbachZero()

◆ Kallbach()

◆ Sample()

◆ SeparationEnergy()

Field Documentation

◆ theCompoundFraction

◆ theIncidentEnergy

◆ theIncidentMass

◆ theProductA

◆ theProductEnergy

◆ theProductMass

◆ theProductZ

◆ theProjectileA

◆ theProjectileZ

◆ theResidualA

◆ theResidualMass

◆ theResidualZ

◆ theTargetA

◆ theTargetMass

◆ theTargetZ