diff --git a/NPLib/Physics/NPQFS.cxx b/NPLib/Physics/NPQFS.cxx
index ed46fd18a39860d83fb29933109ec9ccca4af2dd..108d319c9c4bb3bb41100e8a1baf71934662dd38 100644
--- a/NPLib/Physics/NPQFS.cxx
+++ b/NPLib/Physics/NPQFS.cxx
@@ -413,7 +413,7 @@ TGraph* QFS::GetTheta2VsTheta1(double AngleStep_CM){
   vector<double> vx;
   vector<double> vy;
   double theta1,phi1,E1,theta2,phi2,E2;
-  SetPhiCM(20.*TMath::Pi()/180.);
+  SetPhiCM(0.*TMath::Pi()/180.);
 
   for (double angle=0 ; angle <= 180 ; angle+=AngleStep_CM){
     SetThetaCM(angle*TMath::Pi()/180.);
@@ -433,7 +433,7 @@ TGraph* QFS::GetPhi2VsPhi1(double AngleStep_CM){
   vector<double> vx;
   vector<double> vy;
   double theta1,phi1,E1,theta2,phi2,E2;
-  SetThetaCM(20.*TMath::Pi()/180.);
+  SetThetaCM(0.*TMath::Pi()/180.);
 
   for (double angle=-180 ; angle <= 180 ; angle+=AngleStep_CM){
       SetPhiCM(angle*TMath::Pi()/180.);
@@ -536,7 +536,7 @@ void QFS::CalculateVariablesOld(){
     s = ma_off*ma_off + mT*mT + 2*mT*Ea_lab ; 
     fTotalEnergyImpulsionCM = TLorentzVector(0,0,0,sqrt(s));
     fEcm = sqrt(s) - m1 -m2;
-    if(fEcm<=0) { cout<<"ERROR Ecm negative =\t"<<fEcm<<endl;DumpFormatted(); return;}
+    if(fEcm<=0) { cout<<"ERROR Ecm negative =\t"<<fEcm<<endl;Dump(); return;}
 
     vector<double> theta1;
     vector<double> theta2;
diff --git a/NPLib/Physics/TReactionConditions.cxx b/NPLib/Physics/TReactionConditions.cxx
index bd312dc2dfd0d5be26b45b1a2b4e2023b7eddda7..0a2fa85cce65c6f9943fe709d5af60a67fa410b4 100644
--- a/NPLib/Physics/TReactionConditions.cxx
+++ b/NPLib/Physics/TReactionConditions.cxx
@@ -49,7 +49,7 @@ void TReactionConditions::Clear(){
     fRC_Vertex_Position_Z = -1;
     
     fRC_ThetaCM = -1;
-    fRC_Internal_Momentum = {-1, -1; -1};
+    fRC_Internal_Momentum = {-1, -1, -1};
     
     // emmitted particles
     fRC_Particle_Name.clear();
diff --git a/NPLib/Physics/TReactionConditions.h b/NPLib/Physics/TReactionConditions.h
index dc4a404361eea71237cced14655ec689568e63e0..307af4ad2143c808e047966cd0ec377c635601c8 100644
--- a/NPLib/Physics/TReactionConditions.h
+++ b/NPLib/Physics/TReactionConditions.h
@@ -120,7 +120,7 @@ public:
     double GetExcitation3         () const {return fRC_ExcitationEnergy3     ;}//!       
     double GetExcitation4         () const {return fRC_ExcitationEnergy4     ;}//!       
     double GetThetaCM             () const {return fRC_ThetaCM;}//!
-    TVector3 GetInternalMomentum  () const {return fRC_InternalMomentum;}//!
+    TVector3 GetInternalMomentum  () const {return fRC_Internal_Momentum;}//!
     
     // emmitted particles
     int GetParticleMultiplicity()                const {return fRC_Kinetic_Energy.size();}//!