Add QA plots

Examples/JpsiGamKsKlKK/Jamfile

@@ -13,3 +13,4 @@ lib JpsiGamKlKsKK :
 exe MJpsiGamKlKsKKApp : MJpsiGamKlKsKKApp.cc JpsiGamKlKsKK : ;
 exe MJpsiGamKsKlKKParserTestApp : JpsiGamKsKlKKParserTestApp.cc JpsiGamKlKsKK : ;
 
+

Examples/JpsiGamKsKlKK/JpsiGamKsKlKKData.hh

@@ -11,6 +11,7 @@
 
 struct enumJpsiGamKsKlKKData{
   enum {V4_Psi=0, V4_KsKlKpKm_HeliPsi, V4_KsKl_HeliPsi,V4_KpKm_HeliPsi, V4_gamma_HeliPsi,
+	V4_Ks_HeliPsi, V4_Kl_HeliPsi, V4_Kp_HeliPsi, V4_Km_HeliPsi,
 	V4_KsKl_HeliKsKlKpKm,V4_KpKm_HeliKsKlKpKm ,V4_Ks_HeliKsKl,V4_Kp_HeliKpKm,
 	n4Vecs, Df_Psi=n4Vecs, Df_Spin0,Df_Spin2, Df_KsKl,Df_KpKm, nDfuncts};
   
@@ -18,6 +19,7 @@ struct enumJpsiGamKsKlKKData{
   {
     static std::string fitName[enumJpsiGamKsKlKKData::nDfuncts]
       ={"Psi", "KsKlKpKm_HeliPsi", "KsKl_HeliPsi","KpKm_HeliPsi","gamma_HeliPsi",
+	"Ks_HeliPsi", "Kl_HeliPsi", "Kp_HeliPsi", "Km_HeliPsi",
 	"KsKl_HeliKsKlKpKm","KpKm_HeliKsKlKpKm" ,"Ks_HeliKsKl","Kp_HeliKpKm",
 	"Df_Psi","Df_Spin0","Df_Spin2", "Df_KsKl","Df_KpKm"};
     if (t<0 || t>=enumJpsiGamKsKlKKData::nDfuncts) assert(0);

Examples/JpsiGamKsKlKK/JpsiGamKsKlKKEventList.cc

@@ -56,7 +56,7 @@ void JpsiGamKsKlKKEventList::read4Vecs(EventList& evtList, std::vector<EvtData*>
      Vector4<float>  V4_Ks_Lab( ks   );
      Vector4<float>  V4_KpKm_Lab( kp+km   );
      Vector4<float>  V4_Kp_Lab( kp   );
-     
+          
      Vector4<float>  V4_KsKlKpKm_HeliPsi( ks+kl+kp+km   );
      V4_KsKlKpKm_HeliPsi.Boost(V4_psi);
      
@@ -67,15 +67,18 @@ void JpsiGamKsKlKKEventList::read4Vecs(EventList& evtList, std::vector<EvtData*>
      Vector4<float>  V4_KsKl_HeliPsi( ks+kl  );
      V4_KsKl_HeliPsi.Boost(V4_psi);
      
-//      Vector4<float>  V4_KsKl_HeliKsKlKpKm( ks+kl   );
-//      V4_KsKl_HeliKsKlKpKm.Boost( ks+kl+kp+km  );
-
-     Vector4<float>  V4_KsKl_HeliKsKlKpKm=helicityVec(V4_all_Lab, V4_KsKlKpKm_Lab, V4_KsKl_Lab);
+     Vector4<float>  V4_Kl_HeliPsi( kl   );
+     V4_Kl_HeliPsi.Boost( V4_psi );
+     Vector4<float>  V4_Ks_HeliPsi( ks   );
+     V4_Ks_HeliPsi.Boost( V4_psi );
+     Vector4<float>  V4_Kp_HeliPsi( kp   );
+     V4_Kp_HeliPsi.Boost( V4_psi );
+     Vector4<float>  V4_Km_HeliPsi( km   );
+     V4_Km_HeliPsi.Boost( V4_psi );
      
-//      Vector4<float>  V4_KpKm_HeliKsKlKpKm( kp+km   );
-//      V4_KpKm_HeliKsKlKpKm.Boost( ks+kl+kp+km  );
-     Vector4<float>  V4_KpKm_HeliKsKlKpKm=helicityVec(V4_all_Lab, V4_KsKlKpKm_Lab, V4_KpKm_Lab);     
      
+     Vector4<float>  V4_KsKl_HeliKsKlKpKm=helicityVec(V4_all_Lab, V4_KsKlKpKm_Lab, V4_KsKl_Lab);
+     Vector4<float>  V4_KpKm_HeliKsKlKpKm=helicityVec(V4_all_Lab, V4_KsKlKpKm_Lab, V4_KpKm_Lab);     
      Vector4<float>  V4_Ks_HeliKsKl=helicityVec(V4_KsKlKpKm_Lab, V4_KsKl_Lab, V4_Ks_Lab);     
      Vector4<float>  V4_Kp_HeliKpKm=helicityVec(V4_KsKlKpKm_Lab, V4_KpKm_Lab, V4_Kp_Lab);
 
@@ -84,10 +87,13 @@ void JpsiGamKsKlKKEventList::read4Vecs(EventList& evtList, std::vector<EvtData*>
      evtData->FourVecs[enumJpsiGamKsKlKKData::V4_Psi] = V4_psi;
      evtData->FourVecs[enumJpsiGamKsKlKKData::V4_KsKlKpKm_HeliPsi] = V4_KsKlKpKm_HeliPsi;
      
-//      evtData->FourVecs[enumJpsiGamKsKlKKData::V4_KsKlKpKm_HeliPsi] = V4_KsKlKpKm_HeliPsi;
      evtData->FourVecs[enumJpsiGamKsKlKKData::V4_KsKl_HeliPsi] = V4_KsKl_HeliPsi;
      evtData->FourVecs[enumJpsiGamKsKlKKData::V4_KpKm_HeliPsi] = V4_KpKm_HeliPsi;
      evtData->FourVecs[enumJpsiGamKsKlKKData::V4_gamma_HeliPsi] = V4_gamma_HeliPsi;
+     evtData->FourVecs[enumJpsiGamKsKlKKData::V4_Ks_HeliPsi] = V4_Ks_HeliPsi;
+     evtData->FourVecs[enumJpsiGamKsKlKKData::V4_Kl_HeliPsi] = V4_Kl_HeliPsi;
+     evtData->FourVecs[enumJpsiGamKsKlKKData::V4_Kp_HeliPsi] = V4_Kp_HeliPsi;
+     evtData->FourVecs[enumJpsiGamKsKlKKData::V4_Km_HeliPsi] = V4_Km_HeliPsi;
      
      evtData->FourVecs[enumJpsiGamKsKlKKData::V4_KsKl_HeliKsKlKpKm] = V4_KsKl_HeliKsKlKpKm;
      evtData->FourVecs[enumJpsiGamKsKlKKData::V4_KpKm_HeliKsKlKpKm] = V4_KpKm_HeliKsKlKpKm;

PwaUtils/KinUtils.hh

@@ -42,4 +42,6 @@ T costDecHeli(const Vector4<T>& parent,const Vector4<T>& resonance,
 }
 
 
+
+
 #endif