G4INCLNuclearPotentialConstant.cc

Go to the documentation of this file.
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 // INCL++ intra-nuclear cascade model
00027 // Pekka Kaitaniemi, CEA and Helsinki Institute of Physics
00028 // Davide Mancusi, CEA
00029 // Alain Boudard, CEA
00030 // Sylvie Leray, CEA
00031 // Joseph Cugnon, University of Liege
00032 //
00033 #define INCLXX_IN_GEANT4_MODE 1
00034 
00035 #include "globals.hh"
00036 
00046 #include "G4INCLNuclearPotentialConstant.hh"
00047 #include "G4INCLParticleTable.hh"
00048 
00049 namespace G4INCL {
00050 
00051   namespace NuclearPotential {
00052 
00053     // Constructors
00054     NuclearPotentialConstant::NuclearPotentialConstant(const G4int A, const G4int Z, const G4bool aPionPotential)
00055       : INuclearPotential(A, Z, aPionPotential)
00056     {
00057       initialize();
00058     }
00059 
00060     // Destructor
00061     NuclearPotentialConstant::~NuclearPotentialConstant() {
00062     }
00063 
00064     void NuclearPotentialConstant::initialize() {
00065       const G4double mp = ParticleTable::getINCLMass(Proton);
00066       const G4double mn = ParticleTable::getINCLMass(Neutron);
00067 
00068       G4double theFermiMomentum;
00069       if(theA<ParticleTable::clusterTableASize && theZ<ParticleTable::clusterTableZSize)
00070         // Use momentum RMS from tables to define the Fermi momentum for light
00071         // nuclei
00072         theFermiMomentum = Math::sqrtFiveThirds * ParticleTable::getMomentumRMS(theA,theZ);
00073       else
00074         theFermiMomentum = PhysicalConstants::Pf;
00075 
00076       fermiMomentum[Proton] = theFermiMomentum;
00077       const G4double theProtonFermiEnergy = std::sqrt(theFermiMomentum*theFermiMomentum + mp*mp) - mp;
00078       fermiEnergy[Proton] = theProtonFermiEnergy;
00079 
00080       fermiMomentum[Neutron] = theFermiMomentum;
00081       const G4double theNeutronFermiEnergy = std::sqrt(theFermiMomentum*theFermiMomentum + mn*mn) - mn;
00082       fermiEnergy[Neutron] = theNeutronFermiEnergy;
00083 
00084       fermiEnergy[DeltaPlusPlus] = fermiEnergy.find(Proton)->second;
00085       fermiEnergy[DeltaPlus] = fermiEnergy.find(Proton)->second;
00086       fermiEnergy[DeltaZero] = fermiEnergy.find(Neutron)->second;
00087       fermiEnergy[DeltaMinus] = fermiEnergy.find(Neutron)->second;
00088 
00089       const G4double theAverageSeparationEnergy = 0.5*(ParticleTable::getSeparationEnergy(Proton,theA,theZ)+ParticleTable::getSeparationEnergy(Neutron,theA,theZ));
00090       separationEnergy[Proton] = theAverageSeparationEnergy;
00091       separationEnergy[Neutron] = theAverageSeparationEnergy;
00092 
00093       // Use separation energies from the ParticleTable
00094       vNucleon = 0.5*(theProtonFermiEnergy + theNeutronFermiEnergy) + theAverageSeparationEnergy;
00095       vDelta = vNucleon;
00096       separationEnergy[DeltaPlusPlus] = vDelta - fermiEnergy.find(DeltaPlusPlus)->second;
00097       separationEnergy[DeltaPlus] = vDelta - fermiEnergy.find(DeltaPlus)->second;
00098       separationEnergy[DeltaZero] = vDelta - fermiEnergy.find(DeltaZero)->second;
00099       separationEnergy[DeltaMinus] = vDelta - fermiEnergy.find(DeltaMinus)->second;
00100     }
00101 
00102     G4double NuclearPotentialConstant::computePotentialEnergy(const Particle *particle) const {
00103 
00104       switch( particle->getType() )
00105       {
00106         case Proton:
00107         case Neutron:
00108           return vNucleon;
00109           break;
00110 
00111         case PiPlus:
00112         case PiZero:
00113         case PiMinus:
00114           return computePionPotentialEnergy(particle);
00115           break;
00116 
00117         case DeltaPlusPlus:
00118         case DeltaPlus:
00119         case DeltaZero:
00120         case DeltaMinus:
00121           return vDelta;
00122           break;
00123         case UnknownParticle:
00124           ERROR("Trying to compute potential energy of an unknown particle.");
00125           return 0.0;
00126           break;
00127         default:
00128           ERROR("Trying to compute potential energy of a malformed particle.");
00129           return 0.0;
00130           break;
00131       }
00132     }
00133 
00134   }
00135 }
00136 

Generated on Mon May 27 17:48:36 2013 for Geant4 by  doxygen 1.4.7