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

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

@@ -28 +28 @@
 // A prototype of the low energy neutron transport model.
 //
+// 100413 Fix bug in incidence energy by T. Koi  
+//
 #include "G4NeutronHPEnAngCorrelation.hh"
 #include "G4LorentzRotation.hh"
 #include "G4LorentzVector.hh"
+#include "G4RotationMatrix.hh"
 
 G4ReactionProduct * G4NeutronHPEnAngCorrelation::SampleOne(G4double anEnergy)
@@ -67 +70 @@
   {
     // simplify and double check @
-    G4ThreeVector the3Neutron = theNeutron.GetMomentum();
+    G4ThreeVector the3Neutron = theNeutron.GetMomentum(); //theNeutron has value in LAB
     G4double nEnergy = theNeutron.GetTotalEnergy();
-    G4ThreeVector the3Target = theTarget.GetMomentum();
+    G4ThreeVector the3Target = theTarget.GetMomentum();  //theTarget has value in LAB
     G4double tEnergy = theTarget.GetTotalEnergy();
     G4double totE = nEnergy+tEnergy;
@@ -80 +83 @@
     G4ReactionProduct aNeutron;
     aNeutron.Lorentz(theNeutron, theCMS);
-    anEnergy = aNeutron.GetKineticEnergy();
+    //TKDB 100413 
+    //ENDF-6 Formats Manual ENDF-102
+    //CHAPTER 6. FILE 6: PRODUCT ENERGY-ANGLE DISTRIBUTIONS
+    //LCT Reference system for secondary energy and angle (incident energy is always given in the LAB system)
+    //anEnergy = aNeutron.GetKineticEnergy();
+    anEnergy = theNeutron.GetKineticEnergy(); //should be same argumment of "anEnergy"
 
     G4LorentzVector Ptmp (aNeutron.GetMomentum(), aNeutron.GetTotalEnergy());
+
     toZ.rotateZ(-1*Ptmp.phi());
     toZ.rotateY(-1*Ptmp.theta());
-
   }
   theTotalMeanEnergy=0;
-  G4LorentzRotation toLab(toZ.inverse());
+  G4LorentzRotation toLab(toZ.inverse()); //toLab only change axis NOT to LAB system
   for(i=0; i<nProducts; i++)
   {
@@ -110 +118 @@
         it->operator[](ii)->SetMomentum(pTmp1.vect());
         it->operator[](ii)->SetTotalEnergy(pTmp1.e());
-        if(frameFlag==1) // target rest
+        if(frameFlag==1) // target rest //TK 100413 should be LAB?
         {
-          it->operator[](ii)->Lorentz(*(it->operator[](ii)), -1.*theTarget);
+          it->operator[](ii)->Lorentz(*(it->operator[](ii)), -1.*theTarget); //TK 100413 Is this really need?
         }
         else if(frameFlag==2) // CMS