Changeset 962 for trunk/source/processes/hadronic/models/de_excitation/photon_evaporation/src/G4DiscreteGammaDeexcitation.cc
- Timestamp:
- Apr 6, 2009, 12:30:29 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/de_excitation/photon_evaporation/src/G4DiscreteGammaDeexcitation.cc
r819 r962 130 130 G4bool canDo = true; 131 131 132 if (_transition == 0) 133 { 134 canDo = false; 135 136 if (_verbose > 0) 137 G4cout 138 << "G4DiscreteGammaDeexcitation::CanDoTransition - Null transition " 139 << G4endl; 140 } 132 if (_transition == 0) { 133 canDo = false; 134 135 if (_verbose > 0) 136 G4cout 137 << "G4DiscreteGammaDeexcitation::CanDoTransition - Null transition " 138 << G4endl; 139 } 141 140 G4Fragment nucleus = GetNucleus(); 142 141 if (canDo) { … … 154 153 155 154 G4double excitation = nucleus.GetExcitationEnergy(); 155 156 156 if (canDo) { 157 if (excitation <= 0.) 158 { 159 canDo = false; 160 if (_verbose > 0) 161 G4cout 162 << "G4DiscreteGammaDeexcitation::CanDoTransition - Excitation <= 0" 163 << G4endl; 164 } 165 166 if (excitation > _levelManager->MaxLevelEnergy() + _tolerance) canDo = false; 167 if (excitation < _levelManager->MinLevelEnergy() - _tolerance) canDo = false; 168 // The following is a protection to avoid looping in case of elements with very low 169 // ensdf levels 170 if (excitation < _levelManager->MinLevelEnergy() * 0.9) canDo = false; 157 if (excitation <= 0.) { 158 canDo = false; 159 if (_verbose > 0) 160 G4cout 161 << "G4DiscreteGammaDeexcitation::CanDoTransition - Excitation <= 0" 162 << G4endl; 163 } else { 164 if (excitation > _levelManager->MaxLevelEnergy() + _tolerance) canDo = false; 165 if (excitation < _levelManager->MinLevelEnergy() - _tolerance) canDo = false; 166 // The following is a protection to avoid looping in case of elements with very low 167 // ensdf levels 168 if (excitation < _levelManager->MinLevelEnergy() * 0.9) canDo = false; 171 169 172 if (_verbose > 0) 173 { 170 if (_verbose > 0) { 174 171 G4cout << "G4DiscreteGammaDeexcitation::CanDoTransition - Excitation " 175 172 << excitation << ", Min-Max are " … … 177 174 << _levelManager->MaxLevelEnergy() << G4endl; 178 175 } 179 } 180 176 } 177 } 178 181 179 if (canDo) { 182 180 const G4NuclearLevel* level = _levelManager->NearestLevel(excitation); … … 186 184 canDo = false; 187 185 } 188 189 if (_verbose > 0) 190 { 191 G4cout << "G4DiscreteGammaDeexcitation::CanDoTransition - Halflife " 192 << level->HalfLife() << ", Calling from RDM " 193 << (_rdm ? " True " : " False ") << " Max-HL = " << _max_hl << G4endl; 194 } 195 } 196 if (_verbose > 0) 197 { 198 199 G4cout <<"G4DiscreteGammaDeexcitation::CanDoTransition - CanDo:" 200 << (canDo ? " True " : " False ") << G4endl; 201 } 202 203 186 if (_verbose > 0) { 187 G4cout << "G4DiscreteGammaDeexcitation::CanDoTransition - Halflife " 188 << level->HalfLife() << ", Calling from RDM " 189 << (_rdm ? " True " : " False ") << ", Max-HL = " << _max_hl << G4endl; 190 } 191 } 192 if (_verbose > 0) { 193 G4cout <<"G4DiscreteGammaDeexcitation::CanDoTransition - CanDo:" 194 << (canDo ? " True " : " False ") << G4endl; 195 } 196 204 197 return canDo; 205 198
Note: See TracChangeset
for help on using the changeset viewer.