| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 //
27 //
28 // V.Ivanchenko 26.02.2018
29 //
30
31 #include "G4ChatterjeeCrossSection.hh"
32 #include <CLHEP/Units/SystemOfUnits.h>
33 #include "G4Pow.hh"
34
35 static const G4double emax = 50*CLHEP::MeV;
36
37 // index: 0-neutron, 1-proton, 2-deuteron, 3-t
38 // parameters: p0, p1, p2, landa0, landa1, mu0
39
40 static const G4double paramC[6][11] = {
41 // n from Becchetti and Greenlees
42 { 0., 0., 0., 18.57, -22.93, 381.7,
43 // p from Bechetti and Greenlees
44 {15.72, 9.65, -449., 0.00437,-16.58, 244.7,
45 // d from Lohr and Haeberli
46 {-38.21,922.6,-2804.,-0.0323, -5.48, 336.1,
47 // t from Becchetti and Greenlees
48 {-11.04,619.1,-2147., 0.0426, -10.33, 601.9,
49 // 3he from Becchetti and Greenlees
50 {-3.06, 278.5,-1389.,-0.00535,-11.16, 555.5,
51 // alpha from huizenga and igo
52 {10.95, -85.2, 1146., 0.0643, -13.96, 781.2,
53 };
54
55 G4double
56 G4ChatterjeeCrossSection::ComputeCrossSection(
57
58 G4int idx, G4int Z, G4in
59 {
60 G4double sig;
61 G4double Kc = std::min(K, emax);
62
63 // parameterisation for neutron
64 if(0 == Z) {
65 G4double landa = paramC[idx][3]/resA13 + p
66 G4double mu = (paramC[idx][5] + paramC[idx
67 G4double nu = std::abs((paramC[idx][7]*res
68 + paramC[idx][9]);
69 sig = landa*Kc + mu + nu/Kc;
70
71 // parameterisation for charged
72 } else {
73 //JMQ 20.04.2015 1.5 F
74 G4double ec = cb;
75 //G4double ec = 1.44 * Z * resZ / (1.5*res
76 G4double ecsq = ec*ec;
77 G4double p = paramC[idx][0] + paramC[idx][
78 G4double landa = paramC[idx][3]*resA + par
79 G4double mu = paramC[idx][5]*amu1;
80 G4double nu = amu1* (paramC[idx][7] + para
81 G4double q = landa - nu/ecsq - 2*p*ec;
82 G4double r = mu + 2*nu/ec + p*ecsq;
83 G4double ji= std::max(Kc, ec);
84 if(Kc < ec) { sig = p*Kc*Kc + q*Kc + r; }
85 else { sig = p*(Kc - ji)*(Kc - ji) + landa
86 }
87 sig = std::max(sig, 0.0);
88 //G4cout << "Kc= " << Kc << " A= " << resA <
89 // << " idx= " << idx << " Z= " << Z << "
90 return sig;
91 }
92