JpsiGamEtaPiPi: added 3 particle decay

Examples/JpsiGamEtaPiPi/JpsiGamEtaPiPiData.hh

@@ -15,7 +15,7 @@ struct enumJpsiGamEtaPiPiData{
 	  V4_Eta_HeliEtaPipPim, V4_Pip_HeliEtaPipPim, V4_Pim_HeliEtaPipPim, V4_EtaPip_HeliEtaPipPim, V4_EtaPim_HeliEtaPipPim, V4_PipPim_HeliEtaPipPim,
 	  V4_Pip_HeliPipPim, V4_Pim_HeliPipPim, V4_Pip_HeliEtaPip, V4_Pim_HeliEtaPim, V4_Eta_HeliEtaPip, V4_Eta_HeliEtaPim,
 	  V4_normPipPimDecHeliEtaPipPim, V4_normEtaPipDecHeliEtaPipPim,
-	  n4Vecs, Df_Psi=n4Vecs, Df_Spin0, Df_Spin2, Df_KsKl, Df_KpKm, nDfuncts
+	  n4Vecs, Df_Psi=n4Vecs, Df_etapipidec, nDfuncts
   };
   
   static const std::string& name(unsigned int t)
@@ -25,7 +25,7 @@ struct enumJpsiGamEtaPiPiData{
     	"gamma_HeliPsi", "Eta_HeliPsi", "Pip_HeliPsi", "Pim_HeliPsi",
     	"Eta_HeliEtaPipPim", "Pip_HeliEtaPipPim", "Pim_HeliEtaPipPim", "EtaPip_HeliEtaPipPim", "EtaPim_HeliEtaPipPim", "PipPim_HeliEtaPipPim",
     	"Pip_HeliPipPim", "Pim_HeliPipPim", "Pip_HeliEtaPip", "Pim_HeliEtaPim", "Eta_HeliEtaPip", "Eta_HeliEtaPim",
-    	"Df_Psi","Df_Spin0","Df_Spin2", "Df_KsKl","Df_KpKm"
+    	"Df_Psi","Df_etapipidec"
     };
     if (t<0 || t>=enumJpsiGamEtaPiPiData::nDfuncts) assert(0);
     return fitName[t];

Examples/JpsiGamEtaPiPi/JpsiGamEtaPiPiEventList.cc

@@ -147,6 +147,14 @@ void JpsiGamEtaPiPiEventList::read4Vecs(EventList& evtList, std::vector<EvtData*
 			}
 		}
 
+             for (Spin J_etapipi=0; J_etapipi<=2; J_etapipi++){
+               for (Spin lam_etapipi=-J_etapipi; lam_etapipi<=J_etapipi; lam_etapipi++){
+		evtData->WignerDs[enumJpsiGamEtaPiPiData::Df_etapipidec][J_etapipi][lam_etapipi][0]
+		=Wigner_D(V4_normPipPimDecHeliEtaPipPim.Phi(),V4_normPipPimDecHeliEtaPipPim.Theta(),0,J_etapipi,lam_etapipi,0);
+               }
+             }
+
+
 		/*		//WignerD function for 2+ -> phi phi
 		Spin jTensor =2;
 		for(Spin M=-jTensor; M<=jTensor; M++){

Examples/JpsiGamEtaPiPi/JpsiGamEtaPiPiProdLh.cc

@@ -58,9 +58,17 @@ double JpsiGamEtaPiPiProdLh::calcEvtIntensity(EvtData* theData, fitParams& thePa
     
   }
 
+//   std::map<int, map<Spin,map<Spin,map<Spin,complex<double> > > > > f1Amps;
   if(_f1Hyp){
     std::map< boost::shared_ptr<const JPCLS>, double, pawian::Collection::SharedPtrLess > PsiTof1GamMag=theParamVal.Mags[paramEnumJpsiGamEtaPiPi::PsiToF1Gamma];
     std::map< boost::shared_ptr<const JPCLS>, double, pawian::Collection::SharedPtrLess > PsiTof1GamPhi=theParamVal.Phis[paramEnumJpsiGamEtaPiPi::PsiToF1Gamma];
+
+	for (Spin helf1=-1; helf1<2; helf1++){
+	  JmpGmp+= psiToXGammaAmp(1, 1, helf1, 1, theData, PsiTof1GamMag, PsiTof1GamPhi); 
+          JmpGmm+= psiToXGammaAmp(1, 1, helf1, -1, theData, PsiTof1GamMag, PsiTof1GamPhi); 
+          JmmGmp+= psiToXGammaAmp(-1, 1, helf1, 1, theData, PsiTof1GamMag, PsiTof1GamPhi);
+          JmmGmm+= psiToXGammaAmp(-1, 1, helf1, -1, theData, PsiTof1GamMag, PsiTof1GamPhi);
+	}
 //     JmpGmp+=etaGammaAmp(1, 0, 1, theData, PsiToEtaGamMag, PsiToEtaGamPhi );
 //     JmpGmm+=etaGammaAmp(1, 0, -1, theData,  PsiToEtaGamMag, PsiToEtaGamPhi  );
 //     JmmGmp+=etaGammaAmp(-1, 0, 1, theData,  PsiToEtaGamMag, PsiToEtaGamPhi );
@@ -103,7 +111,7 @@ complex<double> JpsiGamEtaPiPiProdLh::psiToXGammaAmp(Spin Minit, Spin jX, Spin l
        *Clebsch(jX, lamX, 1, -lamGamma, PsiState->S, lambda  )
        *conj( theData->WignerDs[enumJpsiGamEtaPiPiData::Df_Psi][PsiState->J][Minit][lambda]  );
      
-     result+= amp;
+     result+= amp*conj(theData->WignerDs[enumJpsiGamEtaPiPiData::Df_etapipidec][jX][lamX][0]);
    }
 
    return result;