Changeset 1347 for trunk/source/processes/hadronic/models/de_excitation/photon_evaporation/src/G4E1SingleProbability1.cc
- Timestamp:
- Dec 22, 2010, 3:52:27 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/de_excitation/photon_evaporation/src/G4E1SingleProbability1.cc
r819 r1347 24 24 // ******************************************************************** 25 25 // 26 // $Id: G4E1SingleProbability1.cc,v 1.5 2010/11/17 16:50:53 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-04-ref-00 $ 26 28 // 27 29 // Class G4E1SingleProbability1.cc … … 31 33 #include "G4ConstantLevelDensityParameter.hh" 32 34 #include "Randomize.hh" 35 #include "G4Pow.hh" 33 36 34 37 // Constructors and operators 35 38 // 36 39 37 G4E1SingleProbability1::G4E1SingleProbability1(const G4E1SingleProbability1& 38 ) : G4VEmissionProbability() 39 { 40 G4E1SingleProbability1::G4E1SingleProbability1() 41 {} 40 42 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 } 43 G4E1SingleProbability1::~G4E1SingleProbability1() 44 {} 68 45 69 46 // Calculate the emission probability … … 71 48 72 49 G4double G4E1SingleProbability1::EmissionProbDensity(const G4Fragment& frag, 73 constG4double exciteE)50 G4double exciteE) 74 51 { 75 52 … … 84 61 G4double theProb = 0.0; 85 62 86 const G4double Afrag = frag.GetA();87 const G4double Zfrag = frag.GetZ();88 constG4double Uexcite = frag.GetExcitationEnergy();63 G4int Afrag = frag.GetA_asInt(); 64 G4int Zfrag = frag.GetZ_asInt(); 65 G4double Uexcite = frag.GetExcitationEnergy(); 89 66 90 67 if( (Uexcite-exciteE) < 0.0 || exciteE < 0 || Uexcite <= 0) return theProb; … … 95 72 96 73 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); 100 75 101 76 G4double levelDensBef = std::exp(2.0*std::sqrt(aLevelDensityParam*Uexcite)); … … 109 84 G4double sigma0 = 2.5 * Afrag * millibarn; // millibarns 110 85 111 G4double Egdp = (40.3 / std::pow(Afrag,0.2) )*MeV;86 G4double Egdp = (40.3 / G4Pow::GetInstance()->powZ(Afrag,0.2) )*MeV; 112 87 G4double GammaR = 0.30 * Egdp; 113 88 114 G4double normC = 1.0 / ((pi * hbarc)*(pi * hbarc));89 const G4double normC = 1.0 / ((pi * hbarc)*(pi * hbarc)); 115 90 116 91 // CD … … 144 119 145 120 G4double G4E1SingleProbability1::EmissionProbability(const G4Fragment& frag, 146 constG4double exciteE)121 G4double exciteE) 147 122 { 148 123 … … 178 153 179 154 G4double G4E1SingleProbability1::EmissionIntegration(const G4Fragment& frag, 180 constG4double ,181 const G4double lowLim, constG4double upLim,182 constG4int numIters)155 G4double , 156 G4double lowLim, G4double upLim, 157 G4int numIters) 183 158 184 159 { … … 187 162 188 163 G4double x; 189 G4double root3 = 1.0/std::sqrt(3.0);164 const G4double root3 = 1.0/std::sqrt(3.0); 190 165 191 166 G4double Step = (upLim-lowLim)/(2.0*numIters); … … 212 187 } 213 188 214 G4E1SingleProbability1::~G4E1SingleProbability1() {}215 189 216 190
Note: See TracChangeset
for help on using the changeset viewer.