Changeset 1347 for trunk/source/processes/hadronic/models/de_excitation/util/src/G4PairingCorrection.cc

Timestamp:

Dec 22, 2010, 3:52:27 PM (13 years ago)

Author:

garnier

Message:

geant4 tag 9.4

File:

• trunk/source/processes/hadronic/models/de_excitation/util/src/G4PairingCorrection.cc (modified) (3 diffs)

trunk/source/processes/hadronic/models/de_excitation/util/src/G4PairingCorrection.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-//
-// $Id: G4PairingCorrection.cc,v 1.6 2009/11/30 10:33:33 gcosmo Exp $
-// GEANT4 tag $Name: geant4-09-03-ref-09 $
+// $Id: G4PairingCorrection.cc,v 1.8 2010/11/15 12:41:58 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-04-ref-00 $
 //
 // Hadronic Process: Nuclear De-excitations
@@ -33 +32 @@
 
 #include "G4PairingCorrection.hh"
-
 
 G4PairingCorrection* G4PairingCorrection::theInstance = 0;
@@ -49 +47 @@
 G4PairingCorrection* G4PairingCorrection::GetInstance()
 {
-  static G4PairingCorrection theCorrections;
-  if (!theInstance)  { theInstance = &theCorrections; }
+  if (!theInstance)  { 
+    static G4PairingCorrection theCorrections;
+    theInstance = &theCorrections; 
+  }
   return theInstance;
 }   
+
+G4double G4PairingCorrection::GetPairingCorrection(G4int A, G4int Z) const
+{
+  G4double PCorrection = 0.0;
+  G4int N = A - Z;
+  if (theCookPairingCorrections->IsInTableThisN(N) &&
+      theCookPairingCorrections->IsInTableThisZ(Z)) 
+    PCorrection = theCookPairingCorrections->GetParingCorrection(A,Z);
+  else if (theCameronGilbertPairingCorrections->IsInTableThisN(N) &&
+           theCameronGilbertPairingCorrections->IsInTableThisZ(Z))
+    PCorrection = theCameronGilbertPairingCorrections->GetPairingCorrection(A,Z);
+  else {
+    const G4double PairingConstant = 12.0*MeV;
+    G4double Pair = (1 - Z + 2*(Z/2)) + (1 - N + 2*(N/2));
+    PCorrection = Pair*PairingConstant/std::sqrt(static_cast<G4double>(A));
+  }
+  return std::max(PCorrection,0.0);
+}
+
+
+G4double G4PairingCorrection::GetFissionPairingCorrection(G4int A, G4int Z) const 
+{
+  const G4double PairingConstant = 14.0*MeV;
+  G4int N = A - Z;
+  G4double Pair = (1 - Z + 2*(Z/2)) + (1 - N + 2*(N/2));
+  G4double PCorrection = Pair*PairingConstant/std::sqrt(static_cast<G4double>(A));
+  return PCorrection;
+}