#include <G4PreCompoundTriton.hh>

Inheritance diagram for G4PreCompoundTriton:

Public Member Functions
	G4PreCompoundTriton ()

virtual	~G4PreCompoundTriton ()

Public Member Functions inherited from G4PreCompoundIon
	G4PreCompoundIon (const G4ParticleDefinition , G4VCoulombBarrier aCoulombBarrier)

virtual	~G4PreCompoundIon ()

Public Member Functions inherited from G4PreCompoundFragment
	G4PreCompoundFragment (const G4ParticleDefinition , G4VCoulombBarrier aCoulombBarrier)

virtual	~G4PreCompoundFragment ()

G4double	CalcEmissionProbability (const G4Fragment &aFragment)

G4double	GetKineticEnergy (const G4Fragment &aFragment)

Public Member Functions inherited from G4VPreCompoundFragment
	G4VPreCompoundFragment (const G4ParticleDefinition , G4VCoulombBarrier aCoulombBarrier)

virtual	~G4VPreCompoundFragment ()

void	Initialize (const G4Fragment &aFragment)

G4ReactionProduct *	GetReactionProduct () const

G4int	GetA () const

G4int	GetZ () const

G4int	GetRestA () const

G4int	GetRestZ () const

G4double	ResidualA13 () const

G4double	GetCoulombBarrier () const

G4double	GetBindingEnergy () const

G4double	GetMaximalKineticEnergy () const

G4double	GetEnergyThreshold () const

G4double	GetEmissionProbability () const

G4double	GetNuclearMass () const

G4double	GetRestNuclearMass () const

G4double	GetReducedMass () const

const G4LorentzVector &	GetMomentum () const

void	SetMomentum (const G4LorentzVector &value)

const G4String	GetName () const

void	SetOPTxs (G4int)

void	UseSICB (G4bool)

Protected Member Functions
virtual G4double	GetRj (G4int NumberParticles, G4int NumberCharged)

virtual G4double	CrossSection (G4double ekin)

virtual G4double	FactorialFactor (G4int N, G4int P)

virtual G4double	CoalescenceFactor (G4int A)

virtual G4double	GetAlpha ()

G4double	GetOpt12 (G4double K)

G4double	GetOpt34 (G4double K)

Protected Member Functions inherited from G4PreCompoundIon
virtual G4double	ProbabilityDistributionFunction (G4double eKin, const G4Fragment &aFragment)

G4double	GetBeta ()

G4double	GetOpt0 (G4double ekin)

Protected Member Functions inherited from G4VPreCompoundFragment
G4bool	IsItPossible (const G4Fragment &aFragment) const

Additional Inherited Members
Protected Attributes inherited from G4VPreCompoundFragment
G4PreCompoundParameters *	theParameters

G4Pow *	g4pow

G4double	theEmissionProbability

G4double	theCoulombBarrier

G4int	OPTxs

G4bool	useSICB

Detailed Description

Definition at line 42 of file G4PreCompoundTriton.hh.

Constructor & Destructor Documentation

G4PreCompoundTriton::G4PreCompoundTriton ( )

Definition at line 47 of file G4PreCompoundTriton.cc.

References G4VPreCompoundFragment::GetA(), G4VPreCompoundFragment::GetRestA(), G4VPreCompoundFragment::GetRestZ(), G4VPreCompoundFragment::GetZ(), and G4VPreCompoundFragment::ResidualA13().

   : G4PreCompoundIon(G4Triton::Triton(), &theTritonCoulombBarrier)
 {
   ResidualA = GetRestA();
   ResidualZ = GetRestZ(); 
   theA = GetA();
   theZ = GetZ();
   ResidualAthrd = ResidualA13();
   FragmentAthrd = ResidualAthrd;
   FragmentA = theA + ResidualA;
 }

G4PreCompoundTriton::~G4PreCompoundTriton ( )

virtual

Definition at line 59 of file G4PreCompoundTriton.cc.

60 {}

Member Function Documentation

G4double G4PreCompoundTriton::CoalescenceFactor ( G4int A )

protectedvirtual

Implements G4PreCompoundIon.

Definition at line 67 of file G4PreCompoundTriton.cc.

 {
   return 243.0/G4double(A*A);
 }    

G4double G4PreCompoundTriton::CrossSection ( G4double ekin )

protectedvirtual

Implements G4PreCompoundIon.

Definition at line 90 of file G4PreCompoundTriton.cc.

References G4endl, G4VPreCompoundFragment::g4pow, G4VPreCompoundFragment::GetA(), G4PreCompoundIon::GetOpt0(), GetOpt12(), GetOpt34(), G4VPreCompoundFragment::GetRestA(), G4VPreCompoundFragment::GetRestZ(), G4VPreCompoundFragment::GetZ(), G4VPreCompoundFragment::OPTxs, G4VPreCompoundFragment::ResidualA13(), and G4Pow::Z13().

 {
   ResidualA = GetRestA();
   ResidualZ = GetRestZ(); 
   theA = GetA();
   theZ = GetZ();
   ResidualAthrd = ResidualA13();
   FragmentA = theA + ResidualA;
   FragmentAthrd = g4pow->Z13(FragmentA);
 
   if (OPTxs==0) { return GetOpt0( K); }
   else if( OPTxs==1 || OPTxs==2) { return GetOpt12( K); }
   else if (OPTxs==3 || OPTxs==4) { return GetOpt34( K); }
   else{
     std::ostringstream errOs;
     errOs << "BAD TRITON CROSS SECTION OPTION !!"  <<G4endl;
     throw G4HadronicException(__FILE__, __LINE__, errOs.str());
     return 0.;
   }
 }

G4double G4PreCompoundTriton::FactorialFactor	(	G4int	N,
		G4int	P
	)

protectedvirtual

Implements G4PreCompoundIon.

Definition at line 62 of file G4PreCompoundTriton.cc.

 {
   return G4double((N-3)*(P-2)*(N-2)*(P-1)*(N-1)*P)/6.0; 
 }

G4double G4PreCompoundTriton::GetAlpha ( )

protectedvirtual

Implements G4PreCompoundIon.

Definition at line 111 of file G4PreCompoundTriton.cc.

 {
   G4double C = 0.0;
   G4int aZ = theZ + ResidualZ;
   if (aZ >= 70) 
     {
       C = 0.10;
     } 
   else 
     {
       C = ((((0.15417e-06*aZ) - 0.29875e-04)*aZ + 0.21071e-02)*aZ - 0.66612e-01)*aZ + 0.98375; 
     }
  
   return 1.0 + C/3.0;
 }

G4double G4PreCompoundTriton::GetOpt12 ( G4double K )

protected

Definition at line 131 of file G4PreCompoundTriton.cc.

References G4VPreCompoundFragment::g4pow, G4INCL::Math::max(), python.hepunit::MeV, and G4Pow::powZ().

Referenced by CrossSection().

 {
   G4double Kc=K;
 
   // JMQ xsec is set constat above limit of validity
   if (K > 50*MeV) { Kc=50*MeV; }
 
   G4double landa ,mu ,nu ,p , Ec,q,r,ji,xs;
  
   G4double    p0 = -11.04;
   G4double    p1 = 619.1;
   G4double    p2 = -2147.;
   G4double    landa0 = -0.0426;
   G4double    landa1 = -10.33;
   G4double    mm0 = 601.9;
   G4double    mu1 = 0.37;
   G4double    nu0 = 583.0;
   G4double    nu1 = -546.2;
   G4double    nu2 = 1.718;  
   G4double    delta=1.2;            
 
   Ec = 1.44*theZ*ResidualZ/(1.5*ResidualAthrd+delta);
   p = p0 + p1/Ec + p2/(Ec*Ec);
   landa = landa0*ResidualA + landa1;
 
   G4double resmu1 = g4pow->powZ(ResidualA,mu1); 
   mu = mm0*resmu1;
   nu = resmu1*(nu0 + nu1*Ec + nu2*(Ec*Ec));
   q = landa - nu/(Ec*Ec) - 2*p*Ec;
   r = mu + 2*nu/Ec + p*(Ec*Ec);
   
   ji=std::max(Kc,Ec);
   if(Kc < Ec) { xs = p*Kc*Kc + q*Kc + r;}
   else {xs = p*(Kc - ji)*(Kc - ji) + landa*Kc + mu + nu*(2 - Kc/ji)/ji ;}
                  
   if (xs <0.0) {xs=0.0;}
               
   return xs;
 }

G4double G4PreCompoundTriton::GetOpt34 ( G4double K )

protected

Definition at line 172 of file G4PreCompoundTriton.cc.

References test::a, test::b, G4VPreCompoundFragment::g4pow, G4INCL::Math::max(), and G4Pow::powZ().

Referenced by CrossSection().

 {
   G4double landa, mu, nu, p , signor(1.),sig;
   G4double ec,ecsq,xnulam,etest(0.),a; 
   G4double b,ecut,cut,ecut2,geom,elab;
 
   G4double     flow = 1.e-18;
   G4double     spill= 1.e+18;
 
   G4double     p0 = -21.45;
   G4double     p1 = 484.7;
   G4double     p2 = -1608.;
   G4double     landa0 = 0.0186;
   G4double     landa1 = -8.90;
   G4double     mm0 = 686.3;
   G4double     mu1 = 0.325;
   G4double     nu0 = 368.9;
   G4double     nu1 = -522.2;
   G4double     nu2 = -4.998;  
   
   G4double      ra=0.80;
         
   //JMQ 13/02/09 increase of reduced radius to lower the barrier
   // ec = 1.44 * theZ * ResidualZ / (1.5*ResidualAthrd+ra);
   ec = 1.44 * theZ * ResidualZ / (1.7*ResidualAthrd+ra);
   ecsq = ec * ec;
   p = p0 + p1/ec + p2/ecsq;
   landa = landa0*ResidualA + landa1;
   a = g4pow->powZ(ResidualA,mu1);
   mu = mm0 * a;
   nu = a* (nu0+nu1*ec+nu2*ecsq);  
   xnulam = nu / landa;
   if (xnulam > spill) { xnulam=0.; }
   if (xnulam >= flow) { etest = 1.2 *std::sqrt(xnulam); }
  
   a = -2.*p*ec + landa - nu/ecsq;
   b = p*ecsq + mu + 2.*nu/ec;
   ecut = 0.;
   cut = a*a - 4.*p*b;
   if (cut > 0.) { ecut = std::sqrt(cut); }
   ecut = (ecut-a) / (p+p);
   ecut2 = ecut;
   //JMQ 290310 for avoiding unphysical increase below minimum (at ecut)
   // ecut<0 means that there is no cut with energy axis, i.e. xs is set 
   // to 0 bellow minimum
   //  if (cut < 0.) ecut2 = ecut - 2.;
   if (cut < 0.) { ecut2 = ecut; }
   elab = K * FragmentA / G4double(ResidualA);
   sig = 0.;
  
   if (elab <= ec) { //start for E<Ec
     if (elab > ecut2) { sig = (p*elab*elab+a*elab+b) * signor; }
   }           //end for E<Ec
   else {           //start for E>Ec
     sig = (landa*elab+mu+nu/elab) * signor;
     geom = 0.;
     if (xnulam < flow || elab < etest) { return sig; }
     geom = std::sqrt(theA*K);
     geom = 1.23*ResidualAthrd + ra + 4.573/geom;
     geom = 31.416 * geom * geom;
     sig = std::max(geom,sig);
   }           //end for E>Ec
   return sig;
 }

G4double G4PreCompoundTriton::GetRj	(	G4int	NumberParticles,
		G4int	NumberCharged
	)

protectedvirtual

Implements G4PreCompoundIon.

Definition at line 72 of file G4PreCompoundTriton.cc.

 {
   G4double rj = 0.0;
   if(nCharged >= 1 && (nParticles-nCharged) >= 2) {
     G4double denominator = 
       G4double(nParticles*(nParticles-1)*(nParticles-2));
     rj = G4double(3*nCharged*(nParticles-nCharged)*(nParticles-nCharged-1))
       /denominator; 
   }
   return rj;
 }

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

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation