Changeset 1315 for trunk/source/processes/hadronic/models/de_excitation/photon_evaporation/src/G4DiscreteGammaTransition.cc

Timestamp:

Jun 18, 2010, 11:42:07 AM (14 years ago)

Author:

garnier

Message:

update geant4-09-04-beta-cand-01 interfaces-V09-03-09 vis-V09-03-08

File:

• trunk/source/processes/hadronic/models/de_excitation/photon_evaporation/src/G4DiscreteGammaTransition.cc (modified) (9 diffs)

trunk/source/processes/hadronic/models/de_excitation/photon_evaporation/src/G4DiscreteGammaTransition.cc

@@ -37 +37 @@
 //
 //      Modifications:
+//        15 April 1999, Alessandro Brunengo (Alessandro.Brunengo@ge.infn.it)
+//              Added creation time evaluation for products of evaporation
+//
+//        21 Nov. 2001, Fan Lei (flei@space.qinetiq.com)
+//              i) added G4int _nucleusZ initialise it through the constructor
+//              ii) modified SelectGamma() to allow the generation of conversion electrons    
+//              iii) added #include G4AtomicShells.hh
+//      
 //        09 Sep. 2002, Fan Lei  (flei@space.qinetiq.com)
 //              Added renormalization to determine whether transition leads to
 //              electron or gamma in SelectGamma()
 //
-//        21 Nov. 2001, Fan Lei (flei@space.qinetiq.com)
-//              i) added G4int _nucleusZ initialise it through the constructor
-//              ii) modified SelectGamma() to allow the generation of conversion electrons      
-//              iii) added #include G4AtomicShells.hh
-//      
-//        15 April 1999, Alessandro Brunengo (Alessandro.Brunengo@ge.infn.it)
-//              Added creation time evaluation for products of evaporation
+//        19 April 2010, J. M. Quesada. 
+//              Corrections added for taking into account mismatch between tabulated 
+//              gamma energies and level energy differences (fake photons eliminated) 
+//
+//        9 May 2010, V.Ivanchenko
+//              Removed unphysical corretions of gamma energy; fixed default particle 
+//              as gamma; do not subtract bounding energy in case of electron emmision
 //      
 // -------------------------------------------------------------------
@@ -55 +63 @@
 #include "G4RandGeneralTmp.hh"
 #include "G4AtomicShells.hh"
+//JMQ: 
+#include "G4NuclearLevelStore.hh"
+#include "G4Pow.hh"
 
 G4DiscreteGammaTransition::G4DiscreteGammaTransition(const G4NuclearLevel& level): 
@@ -60 +71 @@
 { }
 
-G4DiscreteGammaTransition::G4DiscreteGammaTransition(const G4NuclearLevel& level, G4int Z): 
+//JMQ: now A is also needed in the constructor
+//G4DiscreteGammaTransition::G4DiscreteGammaTransition(const G4NuclearLevel& level, G4int Z): 
+G4DiscreteGammaTransition::G4DiscreteGammaTransition(const G4NuclearLevel& level, G4int Z, G4int A): 
   _nucleusZ(Z), _orbitE(-1), _bondE(0.), _aGamma(true), _icm(false), _gammaEnergy(0.), 
-  _level(level), _excitation(0.),  _gammaCreationTime(0.)
+  _level(level), _excitation(0.),  _gammaCreationTime(0.),_A(A),_Z(Z)
 {
   _verbose = 0;
+  //JMQ: added tolerence in the mismatch
+  _tolerance = CLHEP::keV;
 }
 
@@ -74 +89 @@
 void G4DiscreteGammaTransition::SelectGamma()
 {
-
+  // default gamma 
+  _aGamma = true;    
   _gammaEnergy = 0.;
-
+  
   G4int nGammas = _level.NumberOfGammas();
   if (nGammas > 0)
     {
       G4double random = G4UniformRand();
-
-      G4int iGamma = 0;
-      for(iGamma=0;iGamma < nGammas;iGamma++)
+      
+      G4int iGamma;
+      for(iGamma=0; iGamma<nGammas; ++iGamma)
         {
           if(random <= (_level.GammaCumulativeProbabilities())[iGamma])
-            break;
+            { break; }
         }
-
-
+      
       // Small correction due to the fact that there are mismatches between 
       // nominal level energies and emitted gamma energies
-
-      G4double eCorrection = _level.Energy() - _excitation;
-
+      
+      // 09.05.2010 VI : it is an error ?
+      G4double eCorrection = _level.Energy() - _excitation;      
       _gammaEnergy = (_level.GammaEnergies())[iGamma] - eCorrection;
-
+            
+      //JMQ: 
+      //1)If chosen gamma energy is close enough to excitation energy, the later
+      //  is used instead for gamma dacey to gs (it guarantees energy conservation)
+      //2)For energy conservation, level energy differences instead of  tabulated 
+      //  gamma energies must be used (origin of final fake photons)
+      
+      if(_excitation - _gammaEnergy < _tolerance)      
+        { 
+          _gammaEnergy =_excitation;
+        }
+      else
+        {                
+          _levelManager = G4NuclearLevelStore::GetInstance()->GetManager(_Z,_A);
+          _gammaEnergy = _excitation - 
+            _levelManager->NearestLevel(_excitation - _gammaEnergy)->Energy();
+        }
+      
       //  Warning: the following check is needed to avoid loops:
       //  Due essentially to missing nuclear levels in data files, it is
@@ -106 +138 @@
       //          but this change needs a more complex revision of actual design.
       //          I leave this for a later revision.
-
-      if (_gammaEnergy < _level.Energy()*10e-5) _gammaEnergy = _excitation;
+      
+      if (_gammaEnergy < _level.Energy()*10.e-5) _gammaEnergy = _excitation;
+
+      //G4cout << "G4DiscreteGammaTransition::SelectGamma: " << _gammaEnergy 
+      //             << " _icm: " << _icm << G4endl;
+
       // now decide whether Internal Coversion electron should be emitted instead
       if (_icm) {
@@ -144 +180 @@
               }
             }
-            if (_verbose > 0)
+            _bondE = G4AtomicShells::GetBindingEnergy(_nucleusZ, iShell);
+            if (_verbose > 0) {
               G4cout << "G4DiscreteGammaTransition: _nucleusZ = " <<_nucleusZ 
                      << " , iShell = " << iShell  
-                     << " , Shell binding energy = " << G4AtomicShells::GetBindingEnergy(_nucleusZ, iShell) / keV
+                     << " , Shell binding energy = " << _bondE/keV
                      << " keV " << G4endl;
-            _bondE = G4AtomicShells::GetBindingEnergy(_nucleusZ, iShell);
-            _gammaEnergy = _gammaEnergy - _bondE;
+            }
+
+            // 09.05.2010 VI : it is an error - cannot subtract bond energy from 
+            //                 transition energy here
+            //_gammaEnergy = _gammaEnergy - _bondE; 
+            //G4cout << "_gammaEnergy = " << _gammaEnergy << G4endl;
+
             _orbitE = iShell;     
             _aGamma = false ;   // emitted is not a gamma now 
           }
       }
-    
-      G4double tau = _level.HalfLife() / std::log(2.0);
-
-      G4double tMin = 0;
-      G4double tMax = 10.0 * tau;
+      
+      G4double tau = _level.HalfLife() / G4Pow::GetInstance()->logZ(2);
+
+      //09.05.2010 VI rewrite samling of decay time 
+      //              assuming ordinary exponential low
+      _gammaCreationTime = 0.;      
+      if(tau > 0.0) {  _gammaCreationTime = -tau*log(G4UniformRand()); }
+
+      //G4double tMin = 0;
+      //G4double tMax = 10.0 * tau;
       //  Original code, not very efficent
       //  G4int nBins = 200;
@@ -172 +219 @@
       //  G4RandGeneralTmp randGeneral(sampleArray, nBins);
       //G4double random = randGeneral.shoot();
-  
+      
       //_gammaCreationTime = tMin + (tMax - tMin) * random;
 
       // new code by Fan Lei
       //
+      /*
       if (tau != 0 ) 
       {
@@ -186 +234 @@
           //       << _gammaCreationTime/second << G4endl;
        } else { _gammaCreationTime=0.; }
+      */
     }
   return;
 }
 
-
-//G4bool G4DiscreteGammaTransition::IsAGamma()
-//{
-//  return _aGamma;
-//}
-
-
 G4double G4DiscreteGammaTransition::GetGammaEnergy()
 {
@@ -210 +252 @@
 {
   _excitation = energy;
-  return;
-}
-
-
-
-
-
-
+}
+
+
+
+
+
+