From 1512d0b58f0d9fb5bdc24cffd1438c4b6c830ba9 Mon Sep 17 00:00:00 2001
From: fysquare <freddy.flavigny@laposte.net>
Date: Thu, 14 May 2020 15:23:42 +0200
Subject: [PATCH] Fixing typo for QFS
---
NPLib/Physics/NPQFS.cxx | 6 +++---
NPLib/Physics/TReactionConditions.cxx | 2 +-
NPLib/Physics/TReactionConditions.h | 2 +-
3 files changed, 5 insertions(+), 5 deletions(-)
diff --git a/NPLib/Physics/NPQFS.cxx b/NPLib/Physics/NPQFS.cxx
index ed46fd18a..108d319c9 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 bd312dc2d..0a2fa85cc 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 dc4a40436..307af4ad2 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();}//!
--
GitLab