Ignore:
Timestamp:
Apr 6, 2009, 12:30:29 PM (15 years ago)
Author:
garnier
Message:

update processes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/source/processes/hadronic/models/de_excitation/photon_evaporation/src/G4DiscreteGammaDeexcitation.cc

    r819 r962  
    130130  G4bool canDo = true;
    131131   
    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  }
    141140  G4Fragment nucleus = GetNucleus();
    142141  if (canDo)  {
     
    154153
    155154  G4double excitation = nucleus.GetExcitationEnergy();
     155
    156156  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; 
    171169 
    172     if (_verbose > 0)
    173       {
     170      if (_verbose > 0) {
    174171        G4cout << "G4DiscreteGammaDeexcitation::CanDoTransition -  Excitation "
    175172               << excitation << ", Min-Max are "
     
    177174               << _levelManager->MaxLevelEnergy() << G4endl;
    178175      }
    179   }
    180 
     176    }
     177  }
     178 
    181179  if (canDo) {
    182180    const G4NuclearLevel* level = _levelManager->NearestLevel(excitation); 
     
    186184      canDo = false;
    187185    }
    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 
    204197  return canDo;
    205198     
Note: See TracChangeset for help on using the changeset viewer.