| Geant4 Cross Reference |
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 //
28 // Hadronic Process: Nuclear De-excitations
29 // by V. Lara
30 //
31 // Modified:
32 // 21.03.2013 V.Ivanchenko redesigned and cleaned up
33
34 #include "G4CookPairingCorrections.hh"
35 #include <CLHEP/Units/SystemOfUnits.h>
36
37 // Data comes from:
38 // J.L. Cook, H. Ferguson and A.R.de L. Musgrove, Aust. J. Phys., 20, 477(1967)
39
40
41 // P(Z) 68 values from Z = 28 to Z = 95
42 G4double G4CookPairingCorrections::PairingZTable[] = {
43 1.28, 0.26, 0.88, 0.19, 1.35, -0.05, 1.52, -0.09, 1.17, 0.04,
44 1.24, 0.29, 1.09, 0.26, 1.17, 0.23, 1.15, -0.08, 1.35, 0.34,
45 1.05, 0.28, 1.27, 0.00, 1.05, 0.00, 1.00, 0.09, 1.20, 0.20,
46 1.40, 0.93, 1.00, -0.20, 1.19, 0.09, 0.97, 0.00, 0.92, 0.11,
47 0.68, 0.05, 0.68, -0.22, 0.79, 0.09, 0.69, 0.01, 0.72, 0.00,
48 0.40, 0.16, 0.73, 0.00, 0.46, 0.17, 0.89, 0.00, 0.79, 0.00,
49 0.89, 0.00, 0.81, -0.06, 0.69, -0.20, 0.71, -0.12
50 };
51
52
53 // P(N) 118 values from N = 33 to N = 150
54 G4double G4CookPairingCorrections::PairingNTable[] = {
55 0.08, 1.41, -0.08, 1.50, -0.05, 2.24, -0.47, 1.43, -0.15, 1.44,
56 0.06, 1.56, 0.25, 1.57, -0.16, 1.46, 0.00, 0.93, 0.01, 0.62,
57 -0.50, 1.42, 0.13, 1.52, -0.65, 0.80, -0.08, 1.29, -0.47, 1.25,
58 -0.44, 0.97, 0.08, 1.65, -0.11, 1.26, -0.46, 1.06, 0.22, 1.55,
59 -0.07, 1.37, 0.10, 1.20, -0.27, 0.92, -0.35, 1.19, 0.00, 1.05,
60 -0.25, 1.61, -0.21, 0.90, -0.21, 0.74, -0.38, 0.72, -0.34, 0.92,
61 -0.26, 0.94, 0.01, 0.65, -0.36, 0.83, 0.11, 0.67, 0.05, 1.00,
62 0.51, 1.04, 0.33, 0.68, -0.27, 0.81, 0.09, 0.75, 0.17, 0.86,
63 0.14, 1.10, -0.22, 0.84, -0.47, 0.48, 0.02, 0.88, 0.24, 0.52,
64 0.27, 0.41, -0.05, 0.38, 0.15, 0.67, 0.00, 0.61, 0.00, 0.78,
65 0.00, 0.67, 0.00, 0.67, 0.00, 0.79, 0.00, 0.60, 0.04, 0.64,
66 -0.06, 0.45, 0.05, 0.26, -0.22, 0.39, 0.00, 0.39
67 };
68
69
70 G4CookPairingCorrections::G4CookPairingCorrections()
71 {
72 for(size_t i=0; i<ZTableSize; ++i) { PairingZTable[i] *= CLHEP::MeV; }
73 for(size_t i=0; i<NTableSize; ++i) { PairingNTable[i] *= CLHEP::MeV; }
74 }
75
76
77