G4ParticleHPKallbachMannSyst Class Reference

#include <G4ParticleHPKallbachMannSyst.hh>

Collaboration diagram for G4ParticleHPKallbachMannSyst:

Public Member Functions
	G4ParticleHPKallbachMannSyst (G4double aCompoundFraction, G4double anIncidentEnergy, G4double anIncidentMass, G4double aProductEnergy, G4double aProductMass, G4double aResidualMass, G4int aResidualA, G4int aResidualZ, G4double aTargetMass, G4int aTargetA, G4int aTargetZ)
	~G4ParticleHPKallbachMannSyst ()
G4double	Sample (G4double anEnergy)
G4double	Kallbach (G4double cosTh, G4double anEnergy)
G4double	GetKallbachZero (G4double anEnergy)
G4double	A (G4double anEnergy)
G4double	SeparationEnergy (G4int Ac, G4int Nc, G4int AA, G4int ZA)

Private Attributes
G4double	theCompoundFraction
G4double	theIncidentEnergy
G4double	theIncidentMass
G4double	theProductEnergy
G4double	theProductMass
G4double	theResidualMass
G4double	theTargetMass
G4int	theResidualA
G4int	theResidualZ
G4int	theTargetA
G4int	theTargetZ

Detailed Description

Definition at line 34 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  )

inline

Definition at line 38 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;
  }

~G4ParticleHPKallbachMannSyst()

G4ParticleHPKallbachMannSyst::~G4ParticleHPKallbachMannSyst ( )

inline

Definition at line 57 of file G4ParticleHPKallbachMannSyst.hh.

{};

Here is the call graph for this function:

Member Function Documentation

A()

G4double G4ParticleHPKallbachMannSyst::A

(

G4double

anEnergy

)

Definition at line 96 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+1; 
    G4int Nc = Ac - theTargetZ;
    G4int AA = theTargetA;
    G4int ZA = theTargetZ;
    G4double ea = epsa+SeparationEnergy(Ac, Nc, AA, ZA);
    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 );
  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(0);
  G4int productA = theTargetA+1-theResidualA;
  G4int productZ = theTargetZ-theResidualZ;
  if(productZ==0)
  {
    mb = 0.5;
  }
  else if(productZ==1)
  {
    mb = 1;
  }
  else if(productZ==2)
  {
    mb = 2;
    if(productA==3) mb=1;
  }
  else
  {
    throw G4HadronicException(__FILE__, __LINE__, "Severe error in the sampling of Kallbach-Mann Systematics");
  }
  
  //result = C1*X1 + C2*G4Pow::GetInstance()->powA(X1, 3.) + C3*Ma*mb*G4Pow::GetInstance()->powA(X3, 4.);
  result = C1*X1 + C2*G4Pow::GetInstance()->powN(X1, 3) + C3*Ma*mb*G4Pow::GetInstance()->powN(X3, 4);
  return result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetKallbachZero()

G4double G4ParticleHPKallbachMannSyst::GetKallbachZero

(

G4double

anEnergy

)

Definition at line 84 of file G4ParticleHPKallbachMannSyst.cc.

{
  G4double result;
  if ( theCompoundFraction == 1 ) 
  { 
     //G4cout << "080730b Adjust theCompoundFraction " << G4endl;
     theCompoundFraction *= (1-1.0e-15);   
  } 
  result = 0.5 * (1./A(anEnergy)) * G4Log((1-theCompoundFraction)/(1+theCompoundFraction));
  return result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Kallbach()

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

Definition at line 74 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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Sample()

G4double G4ParticleHPKallbachMannSyst::Sample

(

G4double

anEnergy

)

Definition at line 42 of file G4ParticleHPKallbachMannSyst.cc.

{
  G4double result;
  
  G4double zero = GetKallbachZero(anEnergy);
  if(zero>1) zero=1.;
  if(zero<-1)zero=-1.;
  G4double max = Kallbach(zero, anEnergy);
  double upper = Kallbach(1., anEnergy);
  double 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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

SeparationEnergy()

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

Definition at line 147 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)/Ac - (NA-ZA)*(NA-ZA)/AA);
  //result += -18.56*(std::pow(G4double(Ac), 2./3.) - std::pow(G4double(AA), 2./3.));
  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/std::pow(G4double(Ac),1./3.)-ZA*ZA/std::pow(G4double(AA),1./3.));
  result += -0.717*(Zc*Zc/G4Pow::GetInstance()->A13(G4double(Ac))-ZA*ZA/G4Pow::GetInstance()->A13(G4double(AA)));
  result +=  1.211*(Zc*Zc/Ac-ZA*ZA/AA);
  G4double totalBinding(0);
  G4int productA = theTargetA+1-theResidualA;
  G4int productZ = theTargetZ-theResidualZ;
  if(productZ==0&&productA==1) totalBinding=0;
  if(productZ==1&&productA==1) totalBinding=0;
  if(productZ==1&&productA==2) totalBinding=2.22;
  if(productZ==1&&productA==3) totalBinding=8.48;
  if(productZ==2&&productA==3) totalBinding=7.72;
  if(productZ==2&&productA==4) totalBinding=28.3;
  result += -totalBinding;
  result *= MeV;
  return result;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

theCompoundFraction

G4double G4ParticleHPKallbachMannSyst::theCompoundFraction

private

Definition at line 71 of file G4ParticleHPKallbachMannSyst.hh.

theIncidentEnergy

G4double G4ParticleHPKallbachMannSyst::theIncidentEnergy

private

Definition at line 72 of file G4ParticleHPKallbachMannSyst.hh.

theIncidentMass

G4double G4ParticleHPKallbachMannSyst::theIncidentMass

private

Definition at line 73 of file G4ParticleHPKallbachMannSyst.hh.

theProductEnergy

G4double G4ParticleHPKallbachMannSyst::theProductEnergy

private

Definition at line 74 of file G4ParticleHPKallbachMannSyst.hh.

theProductMass

G4double G4ParticleHPKallbachMannSyst::theProductMass

private

Definition at line 75 of file G4ParticleHPKallbachMannSyst.hh.

theResidualA

G4int G4ParticleHPKallbachMannSyst::theResidualA

private

Definition at line 78 of file G4ParticleHPKallbachMannSyst.hh.

theResidualMass

G4double G4ParticleHPKallbachMannSyst::theResidualMass

private

Definition at line 76 of file G4ParticleHPKallbachMannSyst.hh.

theResidualZ

G4int G4ParticleHPKallbachMannSyst::theResidualZ

private

Definition at line 79 of file G4ParticleHPKallbachMannSyst.hh.

theTargetA

G4int G4ParticleHPKallbachMannSyst::theTargetA

private

Definition at line 80 of file G4ParticleHPKallbachMannSyst.hh.

theTargetMass

G4double G4ParticleHPKallbachMannSyst::theTargetMass

private

Definition at line 77 of file G4ParticleHPKallbachMannSyst.hh.

theTargetZ

G4int G4ParticleHPKallbachMannSyst::theTargetZ

private

Definition at line 81 of file G4ParticleHPKallbachMannSyst.hh.

---

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

• G4ParticleHPKallbachMannSyst.hh
• G4ParticleHPKallbachMannSyst.cc