00001 // 00002 // ******************************************************************** 00003 // * License and Disclaimer * 00004 // * * 00005 // * The Geant4 software is copyright of the Copyright Holders of * 00006 // * the Geant4 Collaboration. It is provided under the terms and * 00007 // * conditions of the Geant4 Software License, included in the file * 00008 // * LICENSE and available at http://cern.ch/geant4/license . These * 00009 // * include a list of copyright holders. * 00010 // * * 00011 // * Neither the authors of this software system, nor their employing * 00012 // * institutes,nor the agencies providing financial support for this * 00013 // * work make any representation or warranty, express or implied, * 00014 // * regarding this software system or assume any liability for its * 00015 // * use. Please see the license in the file LICENSE and URL above * 00016 // * for the full disclaimer and the limitation of liability. * 00017 // * * 00018 // * This code implementation is the result of the scientific and * 00019 // * technical work of the GEANT4 collaboration. * 00020 // * By using, copying, modifying or distributing the software (or * 00021 // * any work based on the software) you agree to acknowledge its * 00022 // * use in resulting scientific publications, and indicate your * 00023 // * acceptance of all terms of the Geant4 Software license. * 00024 // ******************************************************************** 00025 // 00026 // 00027 // J. M. Quesada (July 2009) based on G4TritonCoulombBarrier 00028 // Coded strictly according to Furihata's GEM paper 00029 // 00030 00031 #ifndef G4TritonGEMCoulombBarrier_h 00032 #define G4TritonGEMCoulombBarrier_h 1 00033 00034 #include "G4GEMCoulombBarrier.hh" 00035 #include "globals.hh" 00036 00037 class G4TritonGEMCoulombBarrier : public G4GEMCoulombBarrier 00038 { 00039 public: 00040 00041 G4TritonGEMCoulombBarrier() : G4GEMCoulombBarrier(3,1) {} 00042 ~G4TritonGEMCoulombBarrier() {} 00043 00044 private: 00045 G4TritonGEMCoulombBarrier(const G4TritonGEMCoulombBarrier & right); 00046 00047 const G4TritonGEMCoulombBarrier & operator=(const G4TritonGEMCoulombBarrier & right); 00048 G4bool operator==(const G4TritonGEMCoulombBarrier & right) const; 00049 G4bool operator!=(const G4TritonGEMCoulombBarrier & right) const; 00050 00051 private: 00052 00053 G4double BarrierPenetrationFactor(G4double aZ) const 00054 { 00055 // Data comes from 00056 // Dostrovsky, Fraenkel and Friedlander 00057 // Physical Review, vol 116, num. 3 1959 00058 // (JMQ 190709: according to notes added on proof) 00059 //dataK = {{20, 0.51}, {30, 0.60}, {40, 0.66}, {50, 0.68}}; 00060 // 00061 G4double K = 1.0; 00062 if (aZ >= 50){ 00063 K=0.68; 00064 } else if (aZ <= 20) { 00065 K=0.51; 00066 } else K=0.28445+0.0115956*aZ+0.000026329*aZ*aZ-2.18583*1e-6*aZ*aZ*aZ+3.7083*1e-9*aZ*aZ*aZ*aZ; 00067 return K+0.12; 00068 } 00069 }; 00070 00071 #endif