Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer << 3 // * DISCLAIMER * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers * 6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which * 7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: * 8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license * 9 // * include a list of copyright holders. << 10 // * 9 // * * 11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file << 14 // * use. * 16 // * for the full disclaimer and the limitatio << 17 // * 15 // * * 18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the * 19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. * 20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work * 21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this * 22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. * 23 // * acceptance of all terms of the Geant4 Sof << 24 // ******************************************* 21 // ******************************************************************** 25 // 22 // 26 // 23 // >> 24 // $Id: G4FissionProbability.cc,v 1.4 2003/11/14 16:54:18 hpw Exp $ >> 25 // GEANT4 tag $Name: geant4-06-00 $ 27 // 26 // 28 // Hadronic Process: Nuclear De-excitations 27 // Hadronic Process: Nuclear De-excitations 29 // by V. Lara (Oct 1998) 28 // by V. Lara (Oct 1998) 30 // 29 // 31 // << 32 // J.M.Quesada (14 february 2009) bug fixed in << 33 // parenthesis << 34 30 35 31 36 #include "G4FissionProbability.hh" 32 #include "G4FissionProbability.hh" 37 #include "G4PhysicalConstants.hh" << 38 #include "G4NuclearLevelData.hh" << 39 #include "G4PairingCorrection.hh" 33 #include "G4PairingCorrection.hh" 40 #include "G4EvaporationLevelDensityParameter.h << 34 41 #include "G4FissionLevelDensityParameter.hh" << 35 42 #include "G4Exp.hh" << 36 43 << 37 G4FissionProbability::G4FissionProbability(const G4FissionProbability &) : G4VEmissionProbability() 44 G4FissionProbability::G4FissionProbability() : << 45 G4VEmissionProbability(0, 0), << 46 theEvapLDP(new G4EvaporationLevelDensityPara << 47 theFissLDP(new G4FissionLevelDensityParamete << 48 ownEvapLDP(true), << 49 ownFissLDP(true) << 50 { 38 { 51 fPairCorr = G4NuclearLevelData::GetInstance( << 39 throw G4HadronicException(__FILE__, __LINE__, "G4FissionProbability::copy_constructor meant to not be accessable"); 52 } 40 } 53 41 54 G4FissionProbability::~G4FissionProbability() << 42 >> 43 >> 44 >> 45 const G4FissionProbability & G4FissionProbability::operator=(const G4FissionProbability &) 55 { 46 { 56 if (ownEvapLDP) delete theEvapLDP; << 47 throw G4HadronicException(__FILE__, __LINE__, "G4FissionProbability::operator= meant to not be accessable"); 57 if (ownFissLDP) delete theFissLDP; << 48 return *this; 58 } 49 } 59 50 60 G4double << 51 61 G4FissionProbability::EmissionProbability(cons << 52 G4bool G4FissionProbability::operator==(const G4FissionProbability &) const 62 G4double MaximalKineticEnergy) << 63 // Compute integrated probability of fission << 64 { 53 { 65 if (MaximalKineticEnergy <= 0.0) { return 0. << 54 return false; 66 G4int A = fragment.GetA_asInt(); << 55 } 67 G4int Z = fragment.GetZ_asInt(); << 56 68 G4double U = fragment.GetExcitationEnergy(); << 57 G4bool G4FissionProbability::operator!=(const G4FissionProbability &) const >> 58 { >> 59 return true; >> 60 } >> 61 >> 62 >> 63 G4double G4FissionProbability::EmissionProbability(const G4Fragment & fragment, const G4double MaximalKineticEnergy) >> 64 // Compute integrated probability of fission channel >> 65 { >> 66 if (MaximalKineticEnergy <= 0.0) return 0.0; >> 67 G4double A = fragment.GetA(); >> 68 G4double Z = fragment.GetZ(); >> 69 G4double U = fragment.GetExcitationEnergy(); 69 70 70 G4double Ucompound = U - fPairCorr->GetPairi << 71 G4double Ucompound = U - G4PairingCorrection::GetInstance()->GetPairingCorrection(static_cast<G4int>(A), 71 G4double Ufission = U - fPairCorr->GetFissio << 72 static_cast<G4int>(Z)); 72 if(Ucompound < 0.0 || Ufission < 0.0) { retu << 73 G4double Ufission = U - G4PairingCorrection::GetInstance()->GetFissionPairingCorrection(static_cast<G4int>(A), >> 74 static_cast<G4int>(Z)); 73 75 74 G4double SystemEntropy = << 76 G4double SystemEntropy = 2.0*sqrt(theEvapLDP.LevelDensityParameter(static_cast<G4int>(A), 75 2.0*std::sqrt(theEvapLDP->LevelDensityPara << 77 static_cast<G4int>(Z), >> 78 Ucompound)*Ucompound); 76 79 77 G4double afission = theFissLDP->LevelDensity << 80 G4double afission = theFissLDP.LevelDensityParameter(static_cast<G4int>(A), >> 81 static_cast<G4int>(Z), >> 82 Ufission); 78 83 79 G4double Cf = 2.0*std::sqrt(afission*Maximal << 84 G4double Cf = 2.0*sqrt(afission*MaximalKineticEnergy); 80 G4double Exp1 = (SystemEntropy <= 160.0) ? G << 81 G4double Exp2 = (SystemEntropy-Cf <= 160.0) << 82 85 83 // JMQ 14/02/09 BUG fixed in fission probabi << 86 // G4double Q1 = 1.0 + (Cf - 1.0)*exp(Cf); 84 // at denominator) << 87 // G4double Q2 = 4.0*pi*afission*exp(SystemEntropy); 85 G4double probability = (Exp1 + (Cf-1.0)*Exp2 << 88 86 << 89 // G4double probability = Q1/Q2; 87 return probability; << 90 >> 91 >> 92 G4double Exp1 = 0.0; >> 93 if (SystemEntropy <= 160.0) Exp1 = exp(-SystemEntropy); >> 94 // @@@@@@@@@@@@@@@@@ hpw changed max to min - cannot notify vicente now since cern mail gave up on me... >> 95 G4double Exp2 = exp( std::min(700.0,Cf-SystemEntropy) ); >> 96 >> 97 G4double probability = (Exp1 + (1.0-Cf)*Exp2) / 4.0*pi*afission; >> 98 >> 99 return probability; 88 } 100 } 89 101 90 102