ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4InuclSpecialFunctions.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4InuclSpecialFunctions.hh
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 // 20100114 M. Kelsey -- Remove G4CascadeMomentum, use G4LorentzVector directly
28 // 20100319 M. Kelsey -- Add optional mass argument to generateWithFixedTheta;
29 // define new generateWithRandomAngles, encapsulating code; define
30 // cbrt() cube-root function (in math.h, but not in <math>!)
31 // 20100412 M. Kelsey -- Modify paraMaker[Truncated] to take buffer as argument
32 // 20100914 M. Kelsey -- Migrate to integer A and Z. Discard unused binding
33 // energy functions
34 // 20120608 M. Kelsey -- Fix variable-name "shadowing" compiler warnings.
35 // 20130308 M. Kelsey -- New function to encapsulate INUCL power expansion
36 // 20130808 M. Kelsey -- Convert paraMaker[Truncated] to class object, to
37 // allow thread isolation
38 // 20130829 M. Kelsey -- Address Coverity #52158, 52161 for copy actions.
39 // 20150619 M. Kelsey -- Overload G4cbrt for case of integer arguments
40 
41 #ifndef G4INUCL_SPECIAL_FUNC_HH
42 #define G4INUCL_SPECIAL_FUNC_HH
43 
44 #include "globals.hh"
45 #include "G4LorentzVector.hh"
46 #include <algorithm>
47 #include <vector>
48 
49 template <int N> class G4CascadeInterpolator;
50 
51 
52 namespace G4InuclSpecialFunctions {
54 
55  // NOTE: Used only by G4Fissioner
57 
58  G4double FermiEnergy(G4int A, G4int Z, G4int ntype);
59 
60  G4double getAL(G4int A);
61 
63 
65 
66  G4double G4cbrt(G4double x); // Can't use "cbrt" name, clashes with <math.h>
67  G4double G4cbrt(G4int n); // Use G4Pow::powN() here for speedup
68 
69  G4double inuclRndm(); // Wrapper for G4UniformRand()
70 
71  G4double randomInuclPowers(G4double ekin, // Power series in Ekin, S
72  const G4double (&coeff)[4][4]);
73 
74  G4double randomGauss(G4double sigma); // Gaussian distribution
75 
77 
78  std::pair<G4double, G4double> randomCOS_SIN();
79 
81 
82  // Optional mass argument will be used to fill G4LorentzVector correctly
84  G4double mass=0.);
85 
87 
88  // This is an object to be instantiated by client code
89  class paraMaker {
90  public:
91  paraMaker(G4int verbose=0);
92  ~paraMaker();
93 
94  // NOTE: Passing Z as double here, to be used as interpolation argument
95  void getParams(G4double Z, std::pair<std::vector<G4double>,
96  std::vector<G4double> >& parms);
97 
98  void getTruncated(G4double Z, std::pair<G4double, G4double>& parms);
99 
100  private:
103 
104  // No copy actions
105  paraMaker(const paraMaker& right);
107  };
108 }
109 
110 
111 #endif