few modifications in IsobarTensorDecay.cc III

PwaUtils/IsobarTensorDecay.cc

@@ -139,8 +139,8 @@ void IsobarTensorDecay::fillWignerDs(std::map<std::string, Vector4<double> >& fs
   DebugMsg << name() << endmsg;  
   std::vector< boost::shared_ptr<const JPCLS> >::iterator itJPCLS;
   for(itJPCLS=theJPCLSAmps.begin(); itJPCLS!=theJPCLSAmps.end(); ++itJPCLS){
-    // (*itJPCLS)->print(std::cout);
-    // std::cout << std::endl;
+    //    (*itJPCLS)->print(std::cout);
+    //    std::cout << std::endl;
     Spin L=(*itJPCLS)->L;
     Spin S=(*itJPCLS)->S;
     int s1s2S=spinDaughter1+spinDaughter2+(*itJPCLS)->S;
@@ -216,7 +216,10 @@ void IsobarTensorDecay::fillWignerDs(std::map<std::string, Vector4<double> >& fs
             Tensor<complex<double> > result;
             if(add_lctForTensor){
               //              result = epsilonMotherProject | ( (_metricTensor | (_lctTensor | mother_4Vec )) | (orbTensor | chi12));
-              result = (epsilonMotherProject % mother_4Vec) | (orbTensor % chi12);
+	      result = epsilonMotherProject | ( ( (_lctTensor * mother_4Vec ) | orbTensor) | chi12);
+
+
+	      //              result = (epsilonMotherProject % mother_4Vec) | (orbTensor % chi12);
             }
             else{
               if( orbTensor.Rank()+chi12.Rank() == epsilonMotherProject.Rank()) result = epsilonMotherProject | (orbTensor % chi12);

pbarpUtils/pbarpTensorLh.cc

@@ -39,6 +39,7 @@
 #include "PwaUtils/FitParamsBase.hh"
 #include "PwaUtils/XdecAmpRegistry.hh"
 #include "Particle/Particle.hh"
+#include "Particle/ParticleTable.hh"
 #include "ErrLogger/ErrLogger.hh"
 
 #include <boost/bind.hpp>
@@ -46,12 +47,6 @@
 #include <boost/numeric/ublas/io.hpp>
 
 
-// pbarpTensorLh::pbarpTensorLh(boost::shared_ptr<const EvtDataBaseList> theEvtList) :
-//   pbarpBaseLh(theEvtList)
-// {
-//   initialize();
-// }
-
 pbarpTensorLh::pbarpTensorLh() :
   pbarpBaseLh()
 {
@@ -89,7 +84,27 @@ void  pbarpTensorLh::initialize(){
   fillMap(jpclsTripletp1States, _decAmps, _decAmpsTripletp1);
 
   std::vector< boost::shared_ptr<const JPCLS> > jpclsTripletm1States=_pbarpReactionPtr->jpclsTripletm1States();
-  fillMap(jpclsTripletm1States, _decAmps, _decAmpsTripletm1);  
+  fillMap(jpclsTripletm1States, _decAmps, _decAmpsTripletm1); 
+
+  double pbarMass = _absEnv->particleTable()->particle("antiproton")->mass();
+  double pMass = _absEnv->particleTable()->particle("proton")->mass();
+  double pbarMom = pbarpEnv::instance()->pbarMomentum();
+
+  Vector4<double> pbar4Vec(sqrt(pbarMass*pbarMass+pbarMom*pbarMom), 0.,0., pbarMom);
+  Vector4<double> p4Vec(pMass, 0.,0.,0.);
+  Vector4<double> allVec=pbar4Vec+p4Vec;
+
+  std::vector< boost::shared_ptr<const JPCLS> > jpclsStatesAll=_pbarpReactionPtr->jpclsStates();
+  std::vector< boost::shared_ptr<const JPCLS> >::const_iterator itJPCLS;
+  for(itJPCLS=jpclsStatesAll.begin(); itJPCLS!=jpclsStatesAll.end(); ++itJPCLS){
+    std::shared_ptr<OrbitalTensor> currentTensorPtr(new OrbitalTensor((*itJPCLS)->L));
+    currentTensorPtr->SetP4(pbar4Vec,p4Vec); 
+    _orbTensorMap[(*itJPCLS)->L]=currentTensorPtr;
+
+    std::shared_ptr<PolVector> currentPolVecPtr(new PolVector((*itJPCLS)->J));
+    currentPolVecPtr->SetP4(allVec, allVec.M());
+    _polVectorMap[(*itJPCLS)->J]=currentPolVecPtr; 
+    }
 }
 
 

pbarpUtils/pbarpTensorLh.hh

@@ -33,6 +33,7 @@
 #include <complex>
 #include <boost/shared_ptr.hpp>
 #include <boost/function.hpp>
+#include <memory>
 
 #include "qft++/topincludes/relativistic-quantum-mechanics.hh"
 
@@ -69,4 +70,7 @@ protected:
 private:
 
   void initialize();
+ 
+  std::map< Spin, std::shared_ptr<OrbitalTensor> > _orbTensorMap;
+  std::map< Spin, std::shared_ptr<PolVector> > _polVectorMap;
 };