G4ChatterjeeCrossSection Class Reference

#include <G4ChatterjeeCrossSection.hh>

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

Detailed Description

Definition at line 40 of file G4ChatterjeeCrossSection.hh.

Member Function Documentation

ComputeCrossSection()

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

static

Definition at line 56 of file G4ChatterjeeCrossSection.cc.

{
  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);
  //G4cout << "Kc= " << Kc << " A= " << resA << " A13= " << resA13 << " cb= " << cb
  //     << " idx= " << idx << " Z= " << Z << " amu1= " << amu1 << " sig= " << sig <<G4endl;
  return sig;
}

References emax, G4INCL::Math::max(), G4INCL::Math::min(), paramC, and Z.

Referenced by G4PreCompoundFragment::CrossSection(), and G4EvaporationProbability::CrossSection().

---

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

• source/processes/hadronic/models/de_excitation/util/include/G4ChatterjeeCrossSection.hh
• source/processes/hadronic/models/de_excitation/util/src/G4ChatterjeeCrossSection.cc