qft++/relativistic-quantum-mechanics/BlattWeisskopf.cc: fixed bug for L=1

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;
 

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