G4ProjectileFragmentCrossSection Class Reference

#include <G4ProjectileFragmentCrossSection.hh>

Collaboration diagram for G4ProjectileFragmentCrossSection:

Public Member Functions
	G4ProjectileFragmentCrossSection ()
G4double	doit (G4double Ap, G4double Zp, G4double At, G4double Zt, G4double A, G4double Z)
void	testMe ()

Private Attributes
G4double	p_S [3]
G4double	p_P [3]
G4double	p_Delta [5]
G4double	p_R [3]
G4double	p_Un [2]
G4double	p_Up [4]
G4double	p_mn [3]
G4double	p_mp [3]
G4double	corr_d [3]
G4double	corr_r [3]
G4double	corr_y [3]

Detailed Description

Definition at line 39 of file G4ProjectileFragmentCrossSection.hh.

Constructor & Destructor Documentation

G4ProjectileFragmentCrossSection()

G4ProjectileFragmentCrossSection::G4ProjectileFragmentCrossSection ( )

inline

Definition at line 42 of file G4ProjectileFragmentCrossSection.hh.

  {
        p_S[1] = -2.38;                  // scale factor for xsect in barn    
        p_S[2] = 0.27;         

        p_P[1] = -2.5840E+00;            // slope of mass yield curve         
        p_P[2] = -7.5700E-03;

        p_Delta[1] = -1.0870E+00;        // centroid rel. to beta-stability
        p_Delta[2] = +3.0470E-02;
        p_Delta[3] = +2.1353E-04; 
        p_Delta[4] = +7.1350E+01; 

        p_R[1] = +0.885E+00;             // width parameter R
        p_R[2] = -9.8160E-03;

        p_Un[1] = 1.65;                  // slope par. n-rich ride of Z distr.

        p_Up[1] = 1.7880;                // slope par. p-rich ride of Z distr.  
        p_Up[2] = +4.7210E-03;   
        p_Up[3] = -1.3030E-05;   

        p_mn[1]  = 0.400;                // memory effect n-rich projectiles
        p_mn[2]  = 0.600;        

        p_mp[1] = -10.25;                // memory effect p-rich projectiles
        p_mp[2] = +10.1; 

        corr_d[1] = -25.0;       // correction close to proj.: centroid dzp
        corr_d[2] = 0.800;       
        corr_r[1] = +20.0;       // correction close to proj.: width R
        corr_r[2] = 0.820;
        corr_y[1] = 200.0;       // correction close to proj.: Yield_a
        corr_y[2] = 0.90;        
  }

Member Function Documentation

doit()

G4double G4ProjectileFragmentCrossSection::doit ( G4double  Ap, G4double  Zp, G4double  At, G4double  Zt, G4double  A, G4double  Z  )

inline

Definition at line 78 of file G4ProjectileFragmentCrossSection.hh.

  {
//  calculate mass yield
        G4double Ap13 = G4Pow::GetInstance()->powA(Ap, 1./3.);
        G4double At13 = G4Pow::GetInstance()->powA(At, 1./3.);
        G4double S = p_S[2] * (At13 + Ap13 + p_S[1]);
//  cout << "debug0 "<<S<<" "<<At13<<" "<<Ap13<<" "<<p_S[1]<<" "<<p_S[2]<<endl;
        G4double p    = G4Exp(p_P[2]*Ap + p_P[1]);
        G4double yield_a = p * S * G4Exp(-p * (Ap - A));
  cout << "debug1 "<<yield_a<<endl;
//   modification close to projectile
        G4double f_mod_y=1.0;
        if (A/Ap > corr_y[2])
    {
          f_mod_y=corr_y[1]*G4Pow::GetInstance()->powN(A/Ap-corr_y[2], 2) + 1.0;
        }
        yield_a= yield_a * f_mod_y;
  cout << "debug1 "<<yield_a<<endl;

//   calculate maximum of charge dispersion zprob
        G4double zbeta = A/(1.98+0.0155*G4Pow::GetInstance()->powA(A, (2./3.)));
        G4double zbeta_p = Ap/(1.98+0.0155*G4Pow::GetInstance()->powA(Ap, (2./3.)));
        G4double delta;
    if(A > p_Delta[4]) 
    {
          delta = p_Delta[1] + p_Delta[2]*A;
        }
    else
    {
          delta = p_Delta[3]*A*A;
        }

//   modification close to projectile
        G4double f_mod=1.0;
        if(A/Ap > corr_d[2]) 
    {
          f_mod = corr_d[1]*G4Pow::GetInstance()->powN(A/Ap-corr_d[2], 2) + 1.0;
        }
        delta = delta*f_mod;
        G4double zprob = zbeta+delta;

//   correction for proton- and neutron-rich projectiles
        G4double  dq;
    if((Zp-zbeta_p)>0) 
    {
          dq = G4Exp(p_mp[1] + G4double(A)/G4double(Ap)*p_mp[2]);
      cout << "dq "<<A<<" "<<Ap<<" "<<p_mp[1]
      <<" "<<p_mp[2]<<" "<<dq<<" "<<p_mp[1] + A/Ap*p_mp[2]<<endl;
    }
        else                       
        {
      dq = p_mn[1]*G4Pow::GetInstance()->powN(A/Ap, 2) + p_mn[2]*G4Pow::GetInstance()->powN(A/Ap, 4);
        }
        zprob = zprob + dq * (Zp-zbeta_p);

//   small corr. since Xe-129 and Pb-208 are not on Z_beta line
        zprob = zprob + 0.0020*A;
    cout <<"zprob "<<A<<" "<<dq<<" "<<Zp<<" "<<zbeta_p
         <<" "<<zbeta<<" "<<delta<<endl;

//  calculate width parameter R
        G4double r = G4Exp(p_R[1] + p_R[2]*A);

//  modification close to projectile
        f_mod=1.0;
        if (A/Ap > corr_r[2]) 
    {
          f_mod = corr_r[1]*Ap*G4Pow::GetInstance()->powN(A/Ap-corr_r[2], 4)+1.0;
        }
        r = r*f_mod;

//   change width according to dev. from beta-stability
        if ((Zp-zbeta_p) < 0.0) 
    { 
          r=r*(1.0-0.0833*std::abs(Zp-zbeta_p));
        }

//   calculate slope parameters u_n, u_p
        G4double u_n = p_Un[1];
        G4double u_p = p_Up[1] + p_Up[2]*A + p_Up[3]*A*A;

//   calculate charge dispersion
        G4double expo, fract;
    if((zprob-Z) > 0) 
    {
//     neutron-rich
          expo = -r*G4Pow::GetInstance()->powA(std::abs(zprob-Z), u_n);
          fract   =  G4Exp(expo)*std::sqrt(r/3.14159);
        }
    else
    {
//     proton-rich
          expo = -r*G4Pow::GetInstance()->powA(std::abs(zprob-Z), u_p);
          fract   =  G4Exp(expo)*std::sqrt(r/3.14159);
 cout << "1 "<<expo<<" "<<r<<" "<<zprob<<" "<<Z<<" "<<u_p<<endl;
//     go to exponential slope
          G4double dfdz = 1.2 + 0.647*G4Pow::GetInstance()->powA(A/2.,0.3);
          G4double z_exp = zprob + dfdz * G4Log(10.) / (2.*r);
          if( Z>z_exp ) 
      {
            expo = -r*G4Pow::GetInstance()->powA(std::abs(zprob-z_exp), u_p);
            fract   =  G4Exp(expo)*std::sqrt(r/3.14159)
                      / G4Pow::GetInstance()->powA(G4Pow::GetInstance()->powA(10, dfdz), Z-z_exp);
          }
        }

        cout << "debug "<<fract<<" "<<yield_a<<endl;
    G4double epaxv2=fract*yield_a;
    return epaxv2;
  }

Here is the call graph for this function:

Here is the caller graph for this function:

testMe()

void G4ProjectileFragmentCrossSection::testMe ( )

inline

Definition at line 189 of file G4ProjectileFragmentCrossSection.hh.

  {
    G4ProjectileFragmentCrossSection i;
    cout << i.doit(58, 28, 9, 4, 49, 28) << endl;
 // Sigma = 9.800163E-13 b
  }

Here is the call graph for this function:

Member Data Documentation

corr_d

G4double G4ProjectileFragmentCrossSection::corr_d[3]

private

Definition at line 204 of file G4ProjectileFragmentCrossSection.hh.

corr_r

G4double G4ProjectileFragmentCrossSection::corr_r[3]

private

Definition at line 205 of file G4ProjectileFragmentCrossSection.hh.

corr_y

G4double G4ProjectileFragmentCrossSection::corr_y[3]

private

Definition at line 206 of file G4ProjectileFragmentCrossSection.hh.

p_Delta

G4double G4ProjectileFragmentCrossSection::p_Delta[5]

private

Definition at line 198 of file G4ProjectileFragmentCrossSection.hh.

p_mn

G4double G4ProjectileFragmentCrossSection::p_mn[3]

private

Definition at line 202 of file G4ProjectileFragmentCrossSection.hh.

p_mp

G4double G4ProjectileFragmentCrossSection::p_mp[3]

private

Definition at line 203 of file G4ProjectileFragmentCrossSection.hh.

p_P

G4double G4ProjectileFragmentCrossSection::p_P[3]

private

Definition at line 197 of file G4ProjectileFragmentCrossSection.hh.

p_R

G4double G4ProjectileFragmentCrossSection::p_R[3]

private

Definition at line 199 of file G4ProjectileFragmentCrossSection.hh.

p_S

G4double G4ProjectileFragmentCrossSection::p_S[3]

private

Definition at line 196 of file G4ProjectileFragmentCrossSection.hh.

p_Un

G4double G4ProjectileFragmentCrossSection::p_Un[2]

private

Definition at line 200 of file G4ProjectileFragmentCrossSection.hh.

p_Up

G4double G4ProjectileFragmentCrossSection::p_Up[4]

private

Definition at line 201 of file G4ProjectileFragmentCrossSection.hh.

---

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

• G4ProjectileFragmentCrossSection.hh