G4ChatterjeeCrossSection Class Reference

#include <G4ChatterjeeCrossSection.hh>

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

Detailed Description

Definition at line 59 of file G4ChatterjeeCrossSection.hh.

Member Function Documentation

static G4double G4ChatterjeeCrossSection::ComputeCrossSection ( G4double  K, G4double  cb, G4double  resA13, G4double  amu1, G4int  idx, G4int  Z, G4int  resA  )

inlinestatic

Definition at line 68 of file G4ChatterjeeCrossSection.hh.

  {
    G4double sig;
    G4double Kc  = std::min(K, emax);

    // parameterisation for neutron
    if(0 == Z) {
      G4double landa = paramC[idx][3]/resA13 + paramC[idx][4];
      G4double mu = (paramC[idx][5] + paramC[idx][6]*resA13)*resA13;
      G4double nu = std::abs((paramC[idx][7]*resA + paramC[idx][8]*resA13)*resA13 
                 + paramC[idx][9]);
      sig = landa*Kc + mu + nu/Kc;

      // parameterisation for charged 
    } else {
      //JMQ 20.04.2015 1.5 F
      G4double ec = cb;
      //G4double ec = 1.44 * Z * resZ / (1.5*resA13 + paramC[idx][10]);
      G4double ecsq = ec*ec;
      G4double p = paramC[idx][0] + paramC[idx][1]/ec + paramC[idx][2]/(ecsq); 
      G4double landa = paramC[idx][3]*resA + paramC[idx][4];
      G4double mu = paramC[idx][5]*amu1;
      G4double nu = amu1* (paramC[idx][7] + paramC[idx][8]*ec + paramC[idx][9]*ecsq);
      G4double q = landa - nu/ecsq - 2*p*ec;
      G4double r = mu + 2*nu/ec + p*ecsq;
      G4double ji= std::max(Kc, ec);
      if(Kc < ec) { sig = p*Kc*Kc + q*Kc + r; }
      else { sig = p*(Kc - ji)*(Kc - ji) + landa*Kc + mu + nu*(2 - Kc/ji)/ji; }
    }
    sig = std::max(sig, 0.0);
    return sig;
  }

Here is the call graph for this function:

Here is the caller graph for this function:

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

inlinestatic

Definition at line 63 of file G4ChatterjeeCrossSection.hh.

  {
    return G4Pow::GetInstance()->powZ(resA, paramC[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:

• geant4.10.03.p01/source/processes/hadronic/models/de_excitation/util/include/G4ChatterjeeCrossSection.hh