G4KalbachCrossSection Class Reference

#include <G4KalbachCrossSection.hh>

Collaboration diagram for G4KalbachCrossSection:

Static Public Member Functions
static G4double	ComputePowerParameter (G4int resA, G4int idx)
static G4double	ComputeCrossSection (G4double K, G4double resA13, G4double amu1, G4int idx, G4int Z, G4int A, G4int resZ, G4int resA)

Detailed Description

Definition at line 77 of file G4KalbachCrossSection.hh.

Member Function Documentation

ComputeCrossSection()

static G4double G4KalbachCrossSection::ComputeCrossSection ( G4double  K, G4double  resA13, G4double  amu1, G4int  idx, G4int  Z, G4int  A, G4int  resZ, G4int  resA  )

inlinestatic

Definition at line 86 of file G4KalbachCrossSection.hh.

  {    
    G4double sig = 0.0;
    G4double signor = 1.0; 
    G4double lambda, mu, nu;
    G4double ec  = 0.5;
    if(0 < Z) {
      //JMQ 13.02.2009 tuning for improving cluster emission ddxs 
      //               (spallation benchmark) 
      G4double xx = 1.7;
      if(1 == A) { xx = 1.5; }
      ec = 1.44 * Z * resZ / (xx*resA13 + paramK[idx][10]);
    }
    G4double ecsq = ec*ec;
    G4double elab = K * (A + resA) / G4double(resA);    
    
    if(idx == 0) { // parameterization for neutron

      if(resA < 40)       { signor =0.7 + resA*0.0075; }
      else if(resA > 210) { signor = 1. + (resA-210)*0.004; }
      lambda = paramK[idx][3]/resA13 + paramK[idx][4];
      mu = paramK[idx][5]*resA13 + paramK[idx][6]*resA13*resA13;
      // JMQ 20.11.2008 very low energy behaviour corrected 
      //                (problem for A (apprx.)>60) fix for avoiding  
      //                neutron xs going to zero at very low energies
      nu = std::abs(paramK[idx][7]*resA13*resA + paramK[idx][8]*resA13*resA13 
            + paramK[idx][9]);

    } else { // parameterization for charged 
      // proton correction
      if(idx == 1) {
    if (resA <= 60)      { signor = 0.92; }
    else if (resA < 100) { signor = 0.8 + resA*0.002; }
      }
      lambda = paramK[idx][3]*resA + paramK[idx][4];
      mu = paramK[idx][5]*amu1;
      nu = amu1* (paramK[idx][7] + paramK[idx][8]*ec + paramK[idx][9]*ecsq);
    }

    // threashold cross section
    if(elab <= ec) {
      G4double p = paramK[idx][0];
      if(0 < Z) { p += paramK[idx][1]/ec + paramK[idx][2]/ecsq; }
      G4double a = -2*p*ec + lambda - nu/ecsq;
      G4double b = p*ecsq + mu + 2*nu/ec;
      G4double ecut;
      G4double det = a*a - 4*p*b;
      if (det > 0.0) { ecut = (std::sqrt(det) - a)/(p + p); }
      else           { ecut = a/(p + p); }

      // If ecut>0, sig=0 at elab=ecut
      if(elab > ecut) { 
    sig = (p*elab*elab + a*elab + b)*signor; 

    // extra proton correction
        if(1 == idx) {
      // c and w are for global correction factor for 
      // they are scaled down for light targets where ec is low.
      G4double cc = std::min(3.15, ec*0.5);
      G4double signor2 = (ec - elab - cc) *3.15/ (0.7*cc);
      sig /= (1. + G4Exp(signor2));
    }
      }

      // high energy cross section
    } else {
      // etest is the energy above which the rxn cross section is
      // compared with the geometrical limit and the max taken.

      // neutron parameters
      G4double etest  = 32.;
      G4double xnulam = 1.0;
      // parameters for charged
      if(0 < Z) {
        etest = 0.0;
    xnulam = nu / lambda;
    if(xnulam > spill)       { xnulam= 0.0; }
    else if (xnulam >= flow) { 
      if(1 == idx) { etest = std::sqrt(xnulam) + 7.; }
      else         { etest = 1.2 *std::sqrt(xnulam); }
    }
      }
      // ** For xnulam.gt.0, sig reaches a maximum at sqrt(xnulam).
      sig = (lambda*elab + mu + nu/elab)*signor;
      if (xnulam >= flow && elab >= etest) {
    G4double geom = std::sqrt(A*K);
    geom = 1.23*resA13 + paramK[idx][10] + 4.573/geom;
    geom = 31.416 * geom * geom;
    sig = std::max(sig, geom);
      }
    }
    sig = std::max(sig, 0.0);
    return sig;
  }

Here is the call graph for this function:

Here is the caller graph for this function:

ComputePowerParameter()

static G4double G4KalbachCrossSection::ComputePowerParameter ( G4int  resA, G4int  idx  )

inlinestatic

Definition at line 81 of file G4KalbachCrossSection.hh.

  {
    return G4Pow::GetInstance()->powZ(resA, paramK[idx][6]);
  }

Here is the call graph for this function:

Here is the caller graph for this function:

---

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

• G4KalbachCrossSection.hh