G4INCLNuclearPotentialIsospin.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 "G4INCLNuclearPotentialIsospin.hh"
00047 #include "G4INCLNuclearPotentialConstant.hh"
00048 #include "G4INCLParticleTable.hh"
00049 #include "G4INCLGlobals.hh"
00050 
00051 namespace G4INCL {
00052 
00053   namespace NuclearPotential {
00054 
00055     // Constructors
00056     NuclearPotentialIsospin::NuclearPotentialIsospin(const G4int A, const G4int Z, const G4bool aPionPotential)
00057       : INuclearPotential(A, Z, aPionPotential)
00058     {
00059       initialize();
00060     }
00061 
00062     // Destructor
00063     NuclearPotentialIsospin::~NuclearPotentialIsospin() {}
00064 
00065     void NuclearPotentialIsospin::initialize() {
00066       const G4double ZOverA = ((G4double) theZ) / ((G4double) theA);
00067 
00068       const G4double mp = ParticleTable::getINCLMass(Proton);
00069       const G4double mn = ParticleTable::getINCLMass(Neutron);
00070 
00071       G4double theFermiMomentum;
00072       if(theA<ParticleTable::clusterTableASize && theZ<ParticleTable::clusterTableZSize)
00073         // Use momentum RMS from tables to define the Fermi momentum for light
00074         // nuclei
00075         theFermiMomentum = Math::sqrtFiveThirds * ParticleTable::getMomentumRMS(theA,theZ);
00076       else
00077         theFermiMomentum = PhysicalConstants::Pf;
00078 
00079       fermiMomentum[Proton] = theFermiMomentum * Math::pow13(2.*ZOverA);
00080       const G4double theProtonFermiEnergy = std::sqrt(fermiMomentum[Proton]*fermiMomentum[Proton] + mp*mp) - mp;
00081       fermiEnergy[Proton] = theProtonFermiEnergy;
00082       // Use separation energies from the ParticleTable
00083       const G4double theProtonSeparationEnergy = ParticleTable::getSeparationEnergy(Proton,theA,theZ);
00084       separationEnergy[Proton] = theProtonSeparationEnergy;
00085       vProton = theProtonFermiEnergy + theProtonSeparationEnergy;
00086 
00087       fermiMomentum[Neutron] = theFermiMomentum * Math::pow13(2.*(1.-ZOverA));
00088       const G4double theNeutronFermiEnergy = std::sqrt(fermiMomentum[Neutron]*fermiMomentum[Neutron] + mn*mn) - mn;
00089       fermiEnergy[Neutron] = theNeutronFermiEnergy;
00090       // Use separation energies from the ParticleTable
00091       const G4double theNeutronSeparationEnergy = ParticleTable::getSeparationEnergy(Neutron,theA,theZ);
00092       separationEnergy[Neutron] = theNeutronSeparationEnergy;
00093       vNeutron = theNeutronFermiEnergy + theNeutronSeparationEnergy;
00094 
00095       vDeltaPlus = vProton;
00096       vDeltaZero = vNeutron;
00097       vDeltaPlusPlus = 2.*vDeltaPlus - vDeltaZero;
00098       vDeltaMinus = 2.*vDeltaZero - vDeltaPlus;
00099 
00100       separationEnergy[DeltaPlusPlus] = 2.*theProtonSeparationEnergy - theNeutronSeparationEnergy;
00101       separationEnergy[DeltaPlus] = theProtonSeparationEnergy;
00102       separationEnergy[DeltaZero] = theNeutronSeparationEnergy;
00103       separationEnergy[DeltaMinus] = 2.*theNeutronSeparationEnergy - theProtonSeparationEnergy;
00104 
00105       fermiEnergy[DeltaPlusPlus] = vDeltaPlusPlus - separationEnergy[DeltaPlusPlus];
00106       fermiEnergy[DeltaPlus] = vDeltaPlus - separationEnergy[DeltaPlus];
00107       fermiEnergy[DeltaZero] = vDeltaZero - separationEnergy[DeltaZero];
00108       fermiEnergy[DeltaMinus] = vDeltaMinus - separationEnergy[DeltaMinus];
00109     }
00110 
00111     G4double NuclearPotentialIsospin::computePotentialEnergy(const Particle *particle) const {
00112 
00113       switch( particle->getType() )
00114       {
00115         case Proton:
00116           return vProton;
00117           break;
00118         case Neutron:
00119           return vNeutron;
00120           break;
00121 
00122         case PiPlus:
00123         case PiZero:
00124         case PiMinus:
00125           return computePionPotentialEnergy(particle);
00126           break;
00127 
00128         case DeltaPlusPlus:
00129           return vDeltaPlusPlus;
00130           break;
00131         case DeltaPlus:
00132           return vDeltaPlus;
00133           break;
00134         case DeltaZero:
00135           return vDeltaZero;
00136           break;
00137         case DeltaMinus:
00138           return vDeltaMinus;
00139           break;
00140       case Composite:
00141         ERROR("No potential computed for particle of type Cluster.");
00142         return 0.0;
00143         break;
00144       case UnknownParticle:
00145         ERROR("Trying to compute potential energy for an unknown particle.");
00146         return 0.0;
00147         break;
00148       }
00149 
00150       ERROR("There is no potential for this type of particle.");
00151       return 0.0;
00152     }
00153 
00154   }
00155 }
00156 

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