Changeset 1347 for trunk/source/processes/hadronic/models/de_excitation/photon_evaporation/src/G4ContinuumGammaDeexcitation.cc
- Timestamp:
- Dec 22, 2010, 3:52:27 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/de_excitation/photon_evaporation/src/G4ContinuumGammaDeexcitation.cc
r1340 r1347 24 24 // ******************************************************************** 25 25 // 26 // $Id: G4ContinuumGammaDeexcitation.cc,v 1. 7 2010/04/30 16:08:03vnivanch Exp $27 // GEANT4 tag $Name: geant4-09-0 3-ref-09$26 // $Id: G4ContinuumGammaDeexcitation.cc,v 1.8 2010/11/17 19:17:17 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-04-ref-00 $ 28 28 // 29 29 // ------------------------------------------------------------------- … … 67 67 G4ContinuumGammaDeexcitation::G4ContinuumGammaDeexcitation() 68 68 : _nucleusZ(0), _nucleusA(0), _levelManager(0) 69 { } 70 69 {} 71 70 72 71 G4ContinuumGammaDeexcitation::~G4ContinuumGammaDeexcitation() 73 { } 74 72 {} 75 73 76 74 G4VGammaTransition* G4ContinuumGammaDeexcitation::CreateTransition() 77 75 { 78 76 G4Fragment* nucleus = GetNucleus(); 79 G4int Z = static_cast<G4int>(nucleus->GetZ());80 G4int A = static_cast<G4int>(nucleus->GetA());77 G4int Z = nucleus->GetZ_asInt(); 78 G4int A = nucleus->GetA_asInt(); 81 79 G4double excitation = nucleus->GetExcitationEnergy(); 82 80 … … 91 89 G4cout << "G4ContinuumGammaDeexcitation::CreateTransition - Created" << G4endl; 92 90 } 93 G4VGammaTransition* gt = new G4ContinuumGammaTransition(_levelManager,Z,A,excitation,_verbose ); 91 G4VGammaTransition* gt = 92 new G4ContinuumGammaTransition(_levelManager,Z,A,excitation,_verbose ); 94 93 95 94 return gt; … … 99 98 G4bool G4ContinuumGammaDeexcitation::CanDoTransition() 100 99 { 101 //JMQ: far too small, creating sometimes continuum gammas instead of the right discrete ones 102 // (when excitation energy is slightly over maximum discrete energy): changed 100 //JMQ: far too small, creating sometimes continuum gammas instead 101 // of the right discrete ones (when excitation energy is slightly 102 // over maximum discrete energy): changed 103 103 // G4double tolerance = 10*eV; 104 104 const G4double tolerance = CLHEP::keV; … … 116 116 G4double excitation = nucleus->GetExcitationEnergy(); 117 117 118 // G4int A = (G4int)nucleus->GetA(); 119 // G4int Z = (G4int)nucleus->GetZ(); 120 if (_nucleusA<2 || _nucleusZ<3) 118 if (_nucleusZ < 2 || _nucleusA < 3) 121 119 { 122 120 if (_verbose > 1) {
Note: See TracChangeset
for help on using the changeset viewer.