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 G4ProtonCoulombBarrier 00028 // Coded strictly according to Furihata's GEM paper 00029 // 00030 00031 #ifndef G4ProtonGEMCoulombBarrier_h 00032 #define G4ProtonGEMCoulombBarrier_h 1 00033 00034 #include "G4GEMCoulombBarrier.hh" 00035 #include "globals.hh" 00036 00037 class G4ProtonGEMCoulombBarrier : public G4GEMCoulombBarrier 00038 { 00039 public: 00040 G4ProtonGEMCoulombBarrier() : G4GEMCoulombBarrier(1,1) {} 00041 ~G4ProtonGEMCoulombBarrier() {} 00042 00043 private: 00044 G4ProtonGEMCoulombBarrier(const G4ProtonGEMCoulombBarrier & right); 00045 00046 const G4ProtonGEMCoulombBarrier & operator=(const G4ProtonGEMCoulombBarrier & right); 00047 G4bool operator==(const G4ProtonGEMCoulombBarrier & right) const; 00048 G4bool operator!=(const G4ProtonGEMCoulombBarrier & right) const; 00049 00050 private: 00051 G4double BarrierPenetrationFactor(const G4double aZ) const { 00052 // Data comes from 00053 // Dostrovsky, Fraenkel and Friedlander 00054 // Physical Review, vol 116, num. 3 1959 00055 // (JMQ 190709: according to notes added on proof) 00056 //dataK = {{20, 0.51}, {30, 0.60}, {40, 0.66}, {50, 0.68}}; 00057 // 00058 G4double K = 1.0; 00059 if (aZ >= 50){ 00060 K=0.68; 00061 } else if (aZ <= 20) { 00062 K=0.51; 00063 } 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; 00064 return K; 00065 } 00066 }; 00067 00068 #endif