Geant4_10
G4AlphaGEMCoulombBarrier.hh
Go to the documentation of this file.
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 // J. M. Quesada (July 2009) based on G4AlphaCoulombBarrier
28 // Coded strictly according to Furihata's GEM paper
29 //
30 
31 
32 #ifndef G4AlphaGEMCoulombBarrier_h
33 #define G4AlphaGEMCoulombBarrier_h 1
34 
35 #include "G4GEMCoulombBarrier.hh"
36 #include "globals.hh"
37 
39 {
40 public:
43 
44 private:
46 
47  const G4AlphaGEMCoulombBarrier & operator=(const G4AlphaGEMCoulombBarrier & right);
48  G4bool operator==(const G4AlphaGEMCoulombBarrier & right) const;
49  G4bool operator!=(const G4AlphaGEMCoulombBarrier & right) const;
50 
51 public:
53  {
54  // Data comes from
55  // Dostrovsky, Fraenkel and Friedlander
56  // Physical Review, vol 116, num. 3 1959
57  // (JMQ 190709: according to notes added on proof)
58  // dataKa = {{20, 0.81}, {30, 0.85}, {40, 0.89}, {50, 0.93}};
59  //
60  G4double K = 1.0;
61  if (aZ >= 50){
62  K=0.93;
63  } else if (aZ <= 20) {
64  K=0.81;
65  } else K=0.729802+0.00402544*aZ-1.17276*1e-6*aZ*aZ+2.31248*1e-8*aZ*aZ*aZ-1.65177*1e-10*aZ*aZ*aZ*aZ;
66  return K;
67  }
68 };
69 
70 #endif
G4double BarrierPenetrationFactor(G4double aZ) const
bool G4bool
Definition: G4Types.hh:79
double G4double
Definition: G4Types.hh:76