Ignore:
Timestamp:
Dec 22, 2010, 3:52:27 PM (13 years ago)
Author:
garnier
Message:

geant4 tag 9.4

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/source/processes/hadronic/models/de_excitation/photon_evaporation/src/G4E1SingleProbability1.cc

    r819 r1347  
    2424// ********************************************************************
    2525//
     26// $Id: G4E1SingleProbability1.cc,v 1.5 2010/11/17 16:50:53 vnivanch Exp $
     27// GEANT4 tag $Name: geant4-09-04-ref-00 $
    2628//
    2729//  Class G4E1SingleProbability1.cc
    … …  
    3133#include "G4ConstantLevelDensityParameter.hh"
    3234#include "Randomize.hh"
     35#include "G4Pow.hh"
    3336
    3437// Constructors and operators
    3538//
    3639
    37 G4E1SingleProbability1::G4E1SingleProbability1(const G4E1SingleProbability1&
    38                                                ) : G4VEmissionProbability()
    39 {
     40G4E1SingleProbability1::G4E1SingleProbability1()
     41{}
    4042
    41   throw G4HadronicException(__FILE__, __LINE__, "G4E1SingleProbability1::copy_constructor meant to not be accessible");
    42 
    43 }
    44 
    45 const G4E1SingleProbability1& G4E1SingleProbability1::
    46 operator=(const G4E1SingleProbability1& )
    47 {
    48 
    49   throw G4HadronicException(__FILE__, __LINE__, "G4E1SingleProbability1::operator= meant to not be accessible");
    50   return *this;
    51 }
    52 
    53 G4bool G4E1SingleProbability1::operator==(const G4E1SingleProbability1&
    54                                           ) const
    55 {
    56 
    57   return false;
    58 
    59 }
    60 
    61 G4bool G4E1SingleProbability1::operator!=(const G4E1SingleProbability1& )
    62 const
    63 {
    64 
    65   return true;
    66 
    67 }
     43G4E1SingleProbability1::~G4E1SingleProbability1()
     44{}
    6845
    6946// Calculate the emission probability
    … …  
    7148
    7249G4double G4E1SingleProbability1::EmissionProbDensity(const G4Fragment& frag,
    73                                                      const G4double exciteE)
     50                                                     G4double exciteE)
    7451{
    7552
    … …  
    8461  G4double theProb = 0.0;
    8562
    86   const G4double Afrag = frag.GetA();
    87   const G4double Zfrag = frag.GetZ();
    88   const G4double Uexcite = frag.GetExcitationEnergy();
     63  G4int Afrag = frag.GetA_asInt();
     64  G4int Zfrag = frag.GetZ_asInt();
     65  G4double Uexcite = frag.GetExcitationEnergy();
    8966
    9067  if( (Uexcite-exciteE) < 0.0 || exciteE < 0 || Uexcite <= 0) return theProb;
    … …  
    9572
    9673  G4ConstantLevelDensityParameter a;
    97   G4double aLevelDensityParam = a.LevelDensityParameter(static_cast<G4int>(Afrag),
    98                                                         static_cast<G4int>(Zfrag),
    99                                                         Uexcite);
     74  G4double aLevelDensityParam = a.LevelDensityParameter(Afrag,Zfrag,Uexcite);
    10075
    10176  G4double levelDensBef = std::exp(2.0*std::sqrt(aLevelDensityParam*Uexcite));
    … …  
    10984  G4double sigma0 = 2.5 * Afrag * millibarn;  // millibarns
    11085
    111   G4double Egdp = (40.3 / std::pow(Afrag,0.2) )*MeV;
     86  G4double Egdp = (40.3 / G4Pow::GetInstance()->powZ(Afrag,0.2) )*MeV;
    11287  G4double GammaR = 0.30 * Egdp;
    11388 
    114   G4double normC = 1.0 / ((pi * hbarc)*(pi * hbarc));
     89  const G4double normC = 1.0 / ((pi * hbarc)*(pi * hbarc));
    11590
    11691  // CD
    … …  
    144119
    145120G4double G4E1SingleProbability1::EmissionProbability(const G4Fragment& frag,
    146                                                      const G4double exciteE)
     121                                                     G4double exciteE)
    147122{
    148123
    … …  
    178153
    179154G4double G4E1SingleProbability1::EmissionIntegration(const G4Fragment& frag,
    180                              const G4double ,
    181                              const G4double lowLim, const G4double upLim,
    182                              const G4int numIters)
     155                                                    G4double ,
     156                                                     G4double lowLim, G4double upLim,
     157                                                    G4int numIters)
    183158
    184159{
    … …  
    187162
    188163  G4double x;
    189   G4double root3 = 1.0/std::sqrt(3.0);
     164  const G4double root3 = 1.0/std::sqrt(3.0);
    190165
    191166  G4double Step = (upLim-lowLim)/(2.0*numIters);
    … …  
    212187}
    213188
    214 G4E1SingleProbability1::~G4E1SingleProbability1() {}
    215189
    216190
Note: See TracChangeset for help on using the changeset viewer.