diff --git a/Examples/Psi2STo2K2PiGam/viewHistograms.C b/Examples/Psi2STo2K2PiGam/viewHistograms.C
index 2a2b8756a1ccf6a7a8695c260def4b2fc0f541de..f4a76489cabfe2168e9ce9b3c03df2b75d18ae2a 100644
--- a/Examples/Psi2STo2K2PiGam/viewHistograms.C
+++ b/Examples/Psi2STo2K2PiGam/viewHistograms.C
@@ -19,11 +19,16 @@ float datainvmassK2pi2;
 float datainvmassK1pi1pi2;
 float datainvmassK2pi1pi2;
 float datainvmasspi1pi2;
+float datainvmassK1K2pi1;
+float datainvmassK1K2pi2;
+float datainvmassK1K2;
 
 float datacosthetapi1;
 float datacosthetapi2;
 float datacosthetapipiViaK1pi1pi2;
 float datacosthetapipiViaK2pi1pi2;
+float datacosthetaK1K2ViaK1K2pi1;
+float datacosthetaK1K2ViaK1K2pi2;
 
 std::vector<float> vecdatainvmassK1pi1;
 std::vector<float> vecdatainvmassK1pi2;
@@ -38,6 +43,12 @@ std::vector<float> vecdatainvmasspi1pi2;
 std::vector<float> vecdatacosthetapipiViaK1pi1pi2;
 std::vector<float> vecdatacosthetapipiViaK2pi1pi2;
 
+std::vector<float> vecdatainvmassK1K2pi1;
+std::vector<float> vecdatainvmassK1K2pi2;
+std::vector<float> vecdatainvmassK1K2;
+std::vector<float> vecdatacosthetaK1K2ViaK1K2pi1;
+std::vector<float> vecdatacosthetaK1K2ViaK1K2pi2;
+
 float datacosThetapi1;
 float datacosThetapi2;
 
diff --git a/qft++/relativistic-quantum-mechanics/BlattWeisskopf.cc b/qft++/relativistic-quantum-mechanics/BlattWeisskopf.cc
index 0d9082eae76307fc3ad23ec1002b6e22b58d4ecb..7d24affc175ea87a17fa00f6f527e7bd55268349 100644
--- a/qft++/relativistic-quantum-mechanics/BlattWeisskopf.cc
+++ b/qft++/relativistic-quantum-mechanics/BlattWeisskopf.cc
@@ -78,7 +78,7 @@ double BlattWeisskopf::compute(double p) const
 
     if(0 == _LL) result=1.;
     else
-      if(1 == _LL) result=sqrt(x/(1.0+x));
+      if(1 == _LL) result=sqrt(2.*x/(1.0+x));
       else
 	if(2 == _LL) result=sqrt((13.*x*x)/((x-3.)*(x-3.)+9.*x));
 	else