G4PreCompoundIon.hh

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 // $Id$
00027 //
00028 // J. M. Quesada (August 2008).  
00029 // Based  on previous work by V. Lara
00030 //
00031 // Modified:
00032 // 20.08.2010 V.Ivanchenko added int Z and A and cleanup; added 
00033 //                        G4ParticleDefinition to constructor,
00034 //                        moved constructor and destructor to source,
00035 //                        added inline methods
00036 
00037 #ifndef G4PreCompoundIon_h
00038 #define G4PreCompoundIon_h 1
00039 
00040 #include "G4PreCompoundFragment.hh"
00041 
00042 class G4PreCompoundIon : public G4PreCompoundFragment
00043 {
00044 public:
00045 
00046   G4PreCompoundIon(const G4ParticleDefinition*,
00047                    G4VCoulombBarrier * aCoulombBarrier);
00048   
00049   virtual ~G4PreCompoundIon();
00050 
00051 protected:
00052           
00053   virtual G4double 
00054   ProbabilityDistributionFunction(G4double eKin, 
00055                                   const G4Fragment& aFragment);
00056 
00057   virtual G4double CrossSection(G4double ekin) = 0; 
00058 
00059   virtual G4double 
00060   GetRj(G4int NumberParticles, G4int NumberCharged) = 0; 
00061 
00062   virtual G4double FactorialFactor(G4int N, G4int P) = 0;
00063 
00064   virtual G4double CoalescenceFactor(G4int A) = 0; 
00065 
00066   virtual G4double GetAlpha() = 0;
00067 
00068   inline G4double GetBeta();
00069 
00070   inline G4double GetOpt0(G4double ekin);
00071 
00072 private:
00073 
00074   // default constructor
00075   G4PreCompoundIon();
00076   // operators
00077   G4PreCompoundIon(const G4PreCompoundIon &right);
00078   const G4PreCompoundIon& 
00079   operator= (const G4PreCompoundIon &right);
00080   G4int operator==(const G4PreCompoundIon &right) const;
00081   G4int operator!=(const G4PreCompoundIon &right) const;    
00082 
00083   G4double fact;
00084 };
00085 
00086 inline G4double G4PreCompoundIon::GetBeta()
00087 {
00088   return -GetCoulombBarrier();
00089 }
00090 
00091 // *********************** OPT=0 : Dostrovski's cross section  ***************
00092 inline G4double G4PreCompoundIon::GetOpt0(G4double K)
00093 {
00094   G4double r0 = theParameters->Getr0()*ResidualA13();
00095   // cross section is now given in mb (r0 is in mm) for the sake of consistency
00096   //with the rest of the options
00097   return 1.e+25*CLHEP::pi*r0*r0*ResidualA13()*GetAlpha()*(1.+GetBeta()/K);
00098 }
00099 
00100 #endif

Generated on Mon May 27 17:49:25 2013 for Geant4 by  doxygen 1.4.7