Changeset 1347 for trunk/source/processes/hadronic/models/neutron_hp/src/G4NeutronHPDiscreteTwoBody.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/neutron_hp/src/G4NeutronHPDiscreteTwoBody.cc (modified) (7 diffs)

trunk/source/processes/hadronic/models/neutron_hp/src/G4NeutronHPDiscreteTwoBody.cc

@@ -30 +30 @@
 //080612 Bug fix contribution from Benoit Pirard and Laurent Desorgher (Univ. Bern) #2,3
 //080709 Bug fix Sampling Legendre expansion by T. Koi   
+//101110 Bug fix in MF=6, LAW=2 case; contribution from E. Mendoza, D. Cano-Ott (CIEMAT)
 //
 #include "G4NeutronHPDiscreteTwoBody.hh"
@@ -111 +112 @@
        for(i=0;i<theCoeff[it].GetNumberOfPoly(); i++)
        {
-         theStore.SetX(i, theCoeff[it].GetCoeff(i));
-         theStore.SetY(i, theCoeff[it].GetCoeff(i));
+         //101110
+         //theStore.SetX(i, theCoeff[it].GetCoeff(i));
+         //theStore.SetY(i, theCoeff[it].GetCoeff(i));
+         theStore.SetX(i/2, theCoeff[it].GetCoeff(i));
+         theStore.SetY(i/2, theCoeff[it].GetCoeff(i+1));
          i++;
        }
@@ -125 +129 @@
        for(i=0;i<theCoeff[it].GetNumberOfPoly(); i++)
        {
-         theStore.SetX(i, theCoeff[it].GetCoeff(i));
-         theStore.SetY(i, theCoeff[it].GetCoeff(i));
+         //101110
+         //theStore.SetX(i, theCoeff[it].GetCoeff(i));
+         //theStore.SetY(i, theCoeff[it].GetCoeff(i));
+         theStore.SetX(i/2, theCoeff[it].GetCoeff(i));
+         theStore.SetY(i/2, theCoeff[it].GetCoeff(i+1));
          i++;
        }
@@ -161 +168 @@
          for(i=0;i<theCoeff[it-1].GetNumberOfPoly(); i++)
          {
-           theBuff1.SetX(i, theCoeff[it-1].GetCoeff(i));
-           theBuff1.SetY(i, theCoeff[it-1].GetCoeff(i));
+           //101110
+           //theBuff1.SetX(i, theCoeff[it-1].GetCoeff(i));
+           //theBuff1.SetY(i, theCoeff[it-1].GetCoeff(i));
+           theBuff1.SetX(i/2, theCoeff[it-1].GetCoeff(i));
+           theBuff1.SetY(i/2, theCoeff[it-1].GetCoeff(i+1));
            i++;
          }
@@ -171 +181 @@
          for(i=0;i<theCoeff[it].GetNumberOfPoly(); i++)
          {
+           //theBuff2.SetX(i, theCoeff[it].GetCoeff(i));
+           //theBuff2.SetY(i, theCoeff[it].GetCoeff(i));
            theBuff2.SetX(i, theCoeff[it].GetCoeff(i));
-           theBuff2.SetY(i, theCoeff[it].GetCoeff(i));
+           theBuff2.SetY(i, theCoeff[it].GetCoeff(i+1));
            i++;
          }
@@ -215 +227 @@
          for(i=0;i<theCoeff[it-1].GetNumberOfPoly(); i++)
          {
-           theBuff1.SetX(i, theCoeff[it-1].GetCoeff(i));
-           theBuff1.SetY(i, theCoeff[it-1].GetCoeff(i));
+           //101110
+           //theBuff1.SetX(i, theCoeff[it-1].GetCoeff(i));
+           //theBuff1.SetY(i, theCoeff[it-1].GetCoeff(i));
+           theBuff1.SetX(i/2, theCoeff[it-1].GetCoeff(i));
+           theBuff1.SetY(i/2, theCoeff[it-1].GetCoeff(i+1));
            i++;
          }
@@ -226 +241 @@
          for(i=0;i<theCoeff[it].GetNumberOfPoly(); i++)
          {
-           theBuff2.SetX(i, theCoeff[it].GetCoeff(i));
-           theBuff2.SetY(i, theCoeff[it].GetCoeff(i));
+           //101110
+           //theBuff2.SetX(i, theCoeff[it].GetCoeff(i));
+           //theBuff2.SetY(i, theCoeff[it].GetCoeff(i));
+           theBuff2.SetX(i/2, theCoeff[it].GetCoeff(i));
+           theBuff2.SetY(i/2, theCoeff[it].GetCoeff(i+1));
            i++;
          }