diff --git a/Examples/Psi2STo2K2PiGam/Jamfile b/Examples/Psi2STo2K2PiGam/Jamfile
index 2c1c7c221e16866a0f6c63af2412fedfd20c6485..e6c2eb1ce6f9a041de6ff2ae65c35c471acc8d49 100644
--- a/Examples/Psi2STo2K2PiGam/Jamfile
+++ b/Examples/Psi2STo2K2PiGam/Jamfile
@@ -15,4 +15,4 @@ exe ParserTestApp : ParserTestApp.cc Psi2STo2K2PiGam : ;
 exe Mpsi2STo2K2PiGamTestApp : Mpsi2STo2K2PiGamTestApp.cc Psi2STo2K2PiGam : ;
 exe qaMixedSampleCacheApp : qaMixedSampleCacheApp.cc Psi2STo2K2PiGam : ;
 exe qaMixedSampleFastApp : qaMixedSampleFastApp.cc Psi2STo2K2PiGam : ;
-exe qaMixedSampleTestApp : qaMixedSampleTestApp.cc Psi2STo2K2PiGam : ;
+
diff --git a/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamEvtList.cc b/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamEvtList.cc
index 621734fa315d117d057d3b578feb21cd995f41ec..17bc022e478429e77fb9d57c1e1f02fcfced6269 100644
--- a/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamEvtList.cc
+++ b/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamEvtList.cc
@@ -309,6 +309,12 @@ Psi2STo2K2PiGamEvtData* Psi2STo2K2PiGamEvtList::fillEvtData(Event* anEvent, int
 	 thePsi2STo2K2PiGamEvtData->DfK1mTof0Km[1][lam1][lam2]=Wigner_D(PiPi_HeliKmPi0Pi0_4V.Phi(), PiPi_HeliKmPi0Pi0_4V.Theta(),0, 1, lam1,lam2);
        }
      }
+
+    for (Spin Kj=0; Kj<=3; Kj++){
+ 	 thePsi2STo2K2PiGamEvtData->DfKjpTof0Kp[Kj][0][0]=Wigner_D(PiPi_HeliKpPi0Pi0_4V.Phi(), PiPi_HeliKpPi0Pi0_4V.Theta(),0, Kj, 0,0);
+	 thePsi2STo2K2PiGamEvtData->DfKjmTof0Km[Kj][0][0]=Wigner_D(PiPi_HeliKmPi0Pi0_4V.Phi(), PiPi_HeliKmPi0Pi0_4V.Theta(),0, Kj, 0,0);
+     }
+
   return thePsi2STo2K2PiGamEvtData;
 }
 
diff --git a/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamHist.cc b/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamHist.cc
index 34f414ed9661ab3985d680205758916a6d6b5d32..9bb2bf0c85370eab67f514b861058e39a1313695 100644
--- a/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamHist.cc
+++ b/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamHist.cc
@@ -192,11 +192,11 @@ void Psi2STo2K2PiGamHist::initRootStuff(const std::string& fileName)
 
   initHistMap();
 
-  _dataTuple = new TNtuple("dataTuple","dataTuple","phiKKpipi:costhetaKKpipi:phiK1pi1pi2:costhetaK1pi1pi2:phiK2pi1pi2:costhetaK2pi1pi2:phiK1pi1:costhetaK1pi1:phiK1pi2:costhetaK1pi2:phiK2pi1:costhetaK2pi1:phiK2pi2:costhetaK2pi2:phipi1:costhetapi1:phipi2:costhetapi2:mk1pi1pi2:mk2pi1pi2:mk1pi1:mk1pi2:mk2pi1:mk2pi2:mpipi:costhetapipiViaK1pipi:costhetapipiViaK2pipi:phipipiViaK1pipi:costhetapiViapipi:phipiViapipi:weight:datatype");
-
-  _mcTuple = new TNtuple("mcTuple","mcTuple","phiKKpipi:costhetaKKpipi:phiK1pi1pi2:costhetaK1pi1pi2:phiK2pi1pi2:costhetaK2pi1pi2:phiK1pi1:costhetaK1pi1:phiK1pi2:costhetaK1pi2:phiK2pi1:costhetaK2pi1:phiK2pi2:costhetaK2pi2:phipi1:costhetapi1:phipi2:costhetapi2:mk1pi1pi2:mk2pi1pi2:mk1pi1:mk1pi2:mk2pi1:mk2pi2:mpipi:costhetapipiViaK1pipi:costhetapipiViaK2pipi:phipipiViaK1pipi:costhetapiViapipi:phipiViapipi:weight:datatype");
+  _dataTuple = new TNtuple("dataTuple","dataTuple","phiKKpipi:costhetaKKpipi:phiK1pi1pi2:costhetaK1pi1pi2:phiK2pi1pi2:costhetaK2pi1pi2:phiK1pi1:costhetaK1pi1:phiK1pi2:costhetaK1pi2:phiK2pi1:costhetaK2pi1:phiK2pi2:costhetaK2pi2:phipi1:costhetapi1:phipi2:costhetapi2:mk1k2pi1:mk1k2pi2:mk1k2:mk1pi1pi2:mk2pi1pi2:mk1pi1:mk1pi2:mk2pi1:mk2pi2:mpipi:costhetaK1ViaK1K2:costhetaKKViaK1K2pi1:costhetaKKViaK1K2pi2:costhetapipiViaK1pipi:costhetapipiViaK2pipi:phipipiViaK1pipi:costhetapiViapipi:phipiViapipi:weight:datatype");
 
+  //  _mcTuple = new TNtuple("mcTuple","mcTuple","phiKKpipi:costhetaKKpipi:phiK1pi1pi2:costhetaK1pi1pi2:phiK2pi1pi2:costhetaK2pi1pi2:phiK1pi1:costhetaK1pi1:phiK1pi2:costhetaK1pi2:phiK2pi1:costhetaK2pi1:phiK2pi2:costhetaK2pi2:phipi1:costhetapi1:phipi2:costhetapi2:mK1K2pi1:mK1K2pi2:mK1K2:mk1pi1pi2:mk2pi1pi2:mk1pi1:mk1pi2:mk2pi1:mk2pi2:mpipi:costhetapipiViaK1pipi:costhetapipiViaK2pipi:phipipiViaK1pipi:costhetapiViapipi:phipiViapipi:weight:datatype");
 
+  _mcTuple = new TNtuple("mcTuple","mcTuple","phiKKpipi:costhetaKKpipi:phiK1pi1pi2:costhetaK1pi1pi2:phiK2pi1pi2:costhetaK2pi1pi2:phiK1pi1:costhetaK1pi1:phiK1pi2:costhetaK1pi2:phiK2pi1:costhetaK2pi1:phiK2pi2:costhetaK2pi2:phipi1:costhetapi1:phipi2:costhetapi2:mk1k2pi1:mk1k2pi2:mk1k2:mk1pi1pi2:mk2pi1pi2:mk1pi1:mk1pi2:mk2pi1:mk2pi2:mpipi:costhetaK1ViaK1K2:costhetaKKFromK1K2pi1:costhetaKKFromK1K2pi2:costhetapipiViaK1pipi:costhetapipiViaK2pipi:phipipiViaK1pipi:costhetapiViapipi:phipiViapipi:weight:datatype");
 }
 
 void Psi2STo2K2PiGamHist::initHistMap(){
@@ -319,6 +319,9 @@ void Psi2STo2K2PiGamHist::writeNTuple(TNtuple* theTuple, const Psi2STo2K2PiGamEv
   float CosThetapi1             = -theData->Kp_HeliKpPi0_4V.CosTheta();
   float Phipi2                  = -theData->Kp_HeliKpPi1_4V.Phi();
   float CosThetapi2             = -theData->Kp_HeliKpPi1_4V.CosTheta();
+  float mK1K2pi1                = theData->KKPi0_HeliChic0_4V.M(); 
+  float mK1K2pi2                = theData->KKPi1_HeliChic0_4V.M();  
+  float mK1K2                   = theData->KpKm_HeliChic0_4V.M(); 
   float mK1pi1pi2               = theData->KpPiPi_HeliChic0_4V.M();
   float mK2pi1pi2               = theData->KmPiPi_HeliChic0_4V.M();
   float mK1pi1                  = theData->KpPi0_HeliChic0_4V.M();
@@ -327,6 +330,9 @@ void Psi2STo2K2PiGamHist::writeNTuple(TNtuple* theTuple, const Psi2STo2K2PiGamEv
   float mK2pi2                  = theData->KmPi1_HeliChic0_4V.M();
 
   float mpipi                   = theData->PiPi_HeliChic0_4V.M();
+  float CosThetaK1FromK1K2      = theData->Km_HeliKmKp_4V.CosTheta();
+  float CosThetaKKFromK1K2Pi1   = theData->KK_HeliKKPi0_4V.CosTheta(); 
+  float CosThetaKKFromK1K2Pi2   = theData->KK_HeliKKPi1_4V.CosTheta(); 
   float CosThetaPiPiFromK1PiPi  = theData->PiPi_HeliKpPi0Pi0_4V.CosTheta();
   float CosThetaPiPiFromK2PiPi  = theData->PiPi_HeliKmPi0Pi0_4V.CosTheta();
   float PhiPiPiFromK1PiPi  = theData->PiPi_HeliKpPi0Pi0_4V.Phi();
@@ -356,22 +362,27 @@ void Psi2STo2K2PiGamHist::writeNTuple(TNtuple* theTuple, const Psi2STo2K2PiGamEv
   fourVectors[15] = CosThetapi1;
   fourVectors[16] = Phipi2;
   fourVectors[17] = CosThetapi2;
-  fourVectors[18] = mK1pi1pi2;
-  fourVectors[19] = mK2pi1pi2;
-  fourVectors[20] = mK1pi1;
-  fourVectors[21] = mK1pi2;
-  fourVectors[22] = mK2pi1;
-  fourVectors[23] = mK2pi2;
-  
-  fourVectors[24] = mpipi;
-  fourVectors[25] = CosThetaPiPiFromK1PiPi;
-  fourVectors[26] = CosThetaPiPiFromK2PiPi;
-  fourVectors[27] = PhiPiPiFromK1PiPi;
-  fourVectors[28] = CosThetaPiFromPiPi;
-  fourVectors[29] = PhiPiFromPiPi;
-
-  fourVectors[30] = evtweight;
-  fourVectors[31] = datatype;
+  fourVectors[18] = mK1K2pi1;
+  fourVectors[19] = mK1K2pi2;
+  fourVectors[20] = mK1K2;
+  fourVectors[21] = mK1pi1pi2;
+  fourVectors[22] = mK2pi1pi2;
+  fourVectors[23] = mK1pi1;
+  fourVectors[24] = mK1pi2;
+  fourVectors[25] = mK2pi1;
+  fourVectors[26] = mK2pi2;
+  fourVectors[27] = mpipi;
+  fourVectors[28] = CosThetaK1FromK1K2;
+  fourVectors[29] = CosThetaKKFromK1K2Pi1;
+  fourVectors[30] = CosThetaKKFromK1K2Pi2;
+  fourVectors[31] = CosThetaPiPiFromK1PiPi;
+  fourVectors[32] = CosThetaPiPiFromK2PiPi;
+  fourVectors[33] = PhiPiPiFromK1PiPi;
+  fourVectors[34] = CosThetaPiFromPiPi;
+  fourVectors[35] = PhiPiFromPiPi;
+
+  fourVectors[36] = evtweight;
+  fourVectors[37] = datatype;
 
   //  cout << evtweight << endl;
 
diff --git a/Examples/Psi2STo2K2PiGam/viewHistograms.C b/Examples/Psi2STo2K2PiGam/viewHistograms.C
index f4a76489cabfe2168e9ce9b3c03df2b75d18ae2a..a9ea41cca697bc58a4f7a41e0a039ab5b30eab93 100644
--- a/Examples/Psi2STo2K2PiGam/viewHistograms.C
+++ b/Examples/Psi2STo2K2PiGam/viewHistograms.C
@@ -1,13 +1,13 @@
-#include <vector>
-#include <map>
+#include <vector.h>
+#include <map.h>
 #include <string.h>
 #include <math.h>
 #include <stdlib.h>
-#include <TFile>
+#include <TFile.h>
 #include <TH1F.h>
 #include <TH2F.h>
-#include <TNtuple>
-#include <TCanvas>
+#include <TNtuple.h>
+#include <TCanvas.h>
 #include <TF1.h>
 
 TNtuple* ntdata;
@@ -29,6 +29,14 @@ float datacosthetapipiViaK1pi1pi2;
 float datacosthetapipiViaK2pi1pi2;
 float datacosthetaK1K2ViaK1K2pi1;
 float datacosthetaK1K2ViaK1K2pi2;
+      
+float datacosthetaK1pi1ViaK1pi1pi2;
+float datacosthetaK1pi2ViaK1pi1pi2;
+float datacosthetaK2pi1ViaK2pi1pi2;
+float datacosthetaK2pi2ViaK2pi1pi2;
+
+float datacosthetaK1ViaK1K2;
+float datacosthetaK2ViaK1K2;
 
 std::vector<float> vecdatainvmassK1pi1;
 std::vector<float> vecdatainvmassK1pi2;
@@ -49,16 +57,30 @@ std::vector<float> vecdatainvmassK1K2;
 std::vector<float> vecdatacosthetaK1K2ViaK1K2pi1;
 std::vector<float> vecdatacosthetaK1K2ViaK1K2pi2;
 
+std::vector<float> vecdatacosthetaK1pi1ViaK1pi1pi2;
+std::vector<float> vecdatacosthetaK1pi2ViaK1pi1pi2;
+std::vector<float> vecdatacosthetaK2pi1ViaK2pi1pi2;
+std::vector<float> vecdatacosthetaK2pi2ViaK2pi1pi2;
+                 
+std::vector<float> vecdatacosthetaK1ViaK1K2;
+std::vector<float> vecdatacosthetaK2ViaK1K2;
+
 float datacosThetapi1;
 float datacosThetapi2;
 
 Int_t numOfEntriesData;
 
 TH2F* datainvMassKpivsKpi;
-TH1F* datainvMassKpiangular;
+TH1F* datainvMassKpiAngular;
 
 TH2F* datainvMassKpipivspipi;
-TH1F* datainvMassKpipiangular;
+TH1F* datainvMassKpipiAngular;
+
+TH2F* datainvMassKKpivsKK;
+TH1F* datainvMassKKpiAngular;
+
+TH2F* datainvMassKpivsKpipi;
+TH1F* datainvMassKpivsKpipiAngular;
 
 bool printToPDF = false;
 
@@ -89,6 +111,19 @@ void viewHistograms(TString fname="bin/gcc-4.1.2/debug/link-static/Psi2STo2K2PiG
   ntdata->SetBranchAddress("costhetapipiViaK1pipi", &datacosthetapipiViaK1pi1pi2);
   ntdata->SetBranchAddress("costhetapipiViaK2pipi", &datacosthetapipiViaK2pi1pi2);
 
+  ntdata->SetBranchAddress("mk1k2pi1", &datainvmassK1K2pi1);
+  ntdata->SetBranchAddress("mk1k2pi2", &datainvmassK1K2pi2);
+  ntdata->SetBranchAddress("mk1k2", &datainvmassK1K2);
+  ntdata->SetBranchAddress("costhetaKKViaK1K2pi1", &datacosthetaK1K2ViaK1K2pi1);
+  ntdata->SetBranchAddress("costhetaKKViaK1K2pi2", &datacosthetaK1K2ViaK1K2pi2);
+                                             
+  ntdata->SetBranchAddress("costhetaK1pi1", &datacosthetaK1pi1ViaK1pi1pi2);
+  ntdata->SetBranchAddress("costhetaK1pi2", &datacosthetaK1pi2ViaK1pi1pi2);
+  ntdata->SetBranchAddress("costhetaK2pi1", &datacosthetaK2pi1ViaK2pi1pi2);
+  ntdata->SetBranchAddress("costhetaK2pi2", &datacosthetaK2pi2ViaK2pi1pi2);
+
+  ntdata->SetBranchAddress("costhetaK1ViaK1K2", &datacosthetaK1ViaK1K2);
+
   for(int i = 0; i<numOfEntriesData; i++) {
     ntdata->GetEntry(i);
     vecdatainvmassK1pi1.push_back(datainvmassK1pi1);
@@ -103,6 +138,19 @@ void viewHistograms(TString fname="bin/gcc-4.1.2/debug/link-static/Psi2STo2K2PiG
     vecdatainvmasspi1pi2.push_back(datainvmasspi1pi2);
     vecdatacosthetapipiViaK1pi1pi2.push_back(datacosthetapipiViaK1pi1pi2);
     vecdatacosthetapipiViaK2pi1pi2.push_back(datacosthetapipiViaK2pi1pi2);
+    
+    vecdatainvmassK1K2pi1.push_back(datainvmassK1K2pi1);
+    vecdatainvmassK1K2pi2.push_back(datainvmassK1K2pi2);
+    vecdatainvmassK1K2.push_back(datainvmassK1K2);
+    vecdatacosthetaK1K2ViaK1K2pi1.push_back(datacosthetaK1K2ViaK1K2pi1);
+    vecdatacosthetaK1K2ViaK1K2pi2.push_back(datacosthetaK1K2ViaK1K2pi2);  
+
+    vecdatacosthetaK1pi1ViaK1pi1pi2.push_back(datacosthetaK1pi1ViaK1pi1pi2);
+    vecdatacosthetaK1pi2ViaK1pi1pi2.push_back(datacosthetaK1pi2ViaK1pi1pi2);
+    vecdatacosthetaK2pi1ViaK2pi1pi2.push_back(datacosthetaK2pi1ViaK2pi1pi2);
+    vecdatacosthetaK2pi2ViaK2pi1pi2.push_back(datacosthetaK2pi2ViaK2pi1pi2);
+    
+    vecdatacosthetaK1ViaK1K2.push_back(datacosthetaK1ViaK1K2);
   }
 
   std::vector<TH1F*> histVectData;
@@ -131,7 +179,8 @@ void viewHistograms(TString fname="bin/gcc-4.1.2/debug/link-static/Psi2STo2K2PiG
   histVectMc.push_back(cosK892FittedHist);
   histVectMc.push_back(cosK1430FittedHist);
   histVectMc.push_back(cosK1430ViaK892FittedHist);
-  
+
+  /*  
   TCanvas* cmain = new TCanvas("cmain","cmain",1400,600);
   cmain->Divide(4,2);
   for(int i=0; i<histVectData.size(); i++) {
@@ -198,12 +247,27 @@ void viewHistograms(TString fname="bin/gcc-4.1.2/debug/link-static/Psi2STo2K2PiG
     }
   }
 
+  */
+
   datainvMassKpivsKpi = new TH2F("datainvMassKpivsKpi", "datainvMassKpivsKpi", 48, 0.5, 2.9, 48, 0.5, 2.9);
   datainvMassKpiAngular = new TH1F("datainvMassKpiAngular", "datainvMassKpiAngular", 20, -1., 1.);
 
+  datainvMassKpivsKpi_2 = new TH2F("datainvMassKpivsKpi_2", "datainvMassKpivsKpi_2", 48, 0.5, 2.9, 48, 0.5, 2.9);
+  datainvMassKpiAngular_2 = new TH1F("datainvMassKpiAngular_2", "datainvMassKpiAngular_2", 20, -1., 1.);
+
   datainvMassKpipivspipi = new TH2F("datainvMassKpipivspipi", "datainvMassKpipivspipi", 46, 0.2, 2.5, 48, 0.8, 3.2);
   datainvMassKpipiAngular = new TH1F("datainvMassKpipiAngular", "datainvMassKpipiAngular", 20, -1., 1.);
 
+  datainvMassKKpivsKK = new TH2F("datainvMassKKpivsKK", "datainvMassKKpivsKK", 48, 0.9, 3.3, 48, 1.1, 3.5);
+  datainvMassKKpiAngular = new TH1F("datainvMassKKpivsKKAngular", "datainvMassKKpivsKKAngular", 20, -1., 1.);
+
+  datainvMassKpivsKpipi = new TH2F("datainvMassKpivsKpipi", "datainvMassKpivsKpipi",  48, 0.5, 2.9, 48, 0.5, 2.9);
+  datainvMassKpivsKpipiAngular = new TH1F("datainvMassKpivsKpipiAngular", "datainvMassKpivsKpipiAngular", 20, -1., 1.);
+
+  datainvMasspipivsKK = new TH2F("datainvMasspipivsKK", "datainvMasspipivsKK", 48, 0.9, 3.3, 46, 0.2, 2.5);
+  datainvMasspipivsKKAngular = new TH1F("datainvMasspipivsKKAngular", "datainvMasspipivsKKAngular", 20, -1., 1.);
+
+  
   for(i=0; i<numOfEntriesData; i++) {
     if( (vecdatainvmassK1pi1[i]>0.870 && vecdatainvmassK1pi1[i]<0.930) && (vecdatainvmassK2pi2[i]>1.8 && vecdatainvmassK2pi2[i]<2.3) )
       {
@@ -244,8 +308,8 @@ void viewHistograms(TString fname="bin/gcc-4.1.2/debug/link-static/Psi2STo2K2PiG
 //   fit1->SetLineWidth(3);
 //   datainvMassKpiAngular->Fit("fit1","LMRS","",-10.,10.);
 
-  float lower_limit = 2.3; // 2.3;
-  float upper_limit = 2.6; // 2.6;
+  float lower_limit = 2.50; // 2.3;
+  float upper_limit = 2.60; // 2.6;
     
   for(i=0; i<numOfEntriesData; i++) {
     if( (vecdatainvmasspi1pi2[i] > 0.94 && vecdatainvmasspi1pi2[i] < 1.01 ) && (vecdatainvmassK1pi1pi2[i] > lower_limit && vecdatainvmassK1pi1pi2[i] < upper_limit) )
@@ -273,6 +337,121 @@ void viewHistograms(TString fname="bin/gcc-4.1.2/debug/link-static/Psi2STo2K2PiG
   fit2->SetLineColor(4);
   fit2->SetLineWidth(3);
   datainvMassKpipiAngular->Fit("fit2","LMRS","",-1.,1.);
+  cangular2->Print("KpipivsKpi_2500-2600.pdf");
+
+
+
+  for(i=0; i<numOfEntriesData; i++) {
+    if( (vecdatainvmassK1K2[i] > 1.7 && vecdatainvmassK1K2[i] < 1.77 ) && (vecdatainvmassK1K2pi1[i] > 2.25 && vecdatainvmassK1K2pi1[i] < 2.42) )
+      {
+	datainvMassKKpivsKK->Fill(vecdatainvmassK1K2[i], vecdatainvmassK1K2pi1[i]);
+	datainvMassKKpiAngular->Fill(vecdatacosthetaK1K2ViaK1K2pi1[i]);
+      }
+    if( (vecdatainvmassK1K2[i] > 1.7 && vecdatainvmassK1K2[i] < 1.77 ) && (vecdatainvmassK1K2pi2[i] > 2.25 && vecdatainvmassK1K2pi2[i] < 2.42) )
+      {
+	datainvMassKKpivsKK->Fill(vecdatainvmassK1K2[i], vecdatainvmassK1K2pi2[i]);
+	datainvMassKKpiAngular->Fill(vecdatacosthetaK1K2ViaK1K2pi2[i]);
+      }
+  } 
+
+  TCanvas* cangular3 = new TCanvas("cangular3","cangular3",1200,600);
+  cangular3->Divide(2,1);
+  cangular3->cd(1);
+  datainvMassKKpivsKK->Draw("colz");
+  cangular3->cd(2);
+  datainvMassKKpiAngular->Draw();
+
+
+  
+  for(int i=0; i<numOfEntriesData; i++) {
+    if( (vecdatainvmassK1pi1[i] > 0.80 && vecdatainvmassK1pi1[i] < 0.95 ) && (vecdatainvmassK1pi1pi2[i] > 1.35 && vecdatainvmassK1pi1pi2[i] < 1.5) )
+      {
+	datainvMassKpivsKpipi->Fill(vecdatainvmassK1pi1[i], vecdatainvmassK1pi1pi2[i]);
+	datainvMassKpivsKpipiAngular->Fill(vecdatacosthetaK1pi1ViaK1pi1pi2[i]);
+      }
+    if( (vecdatainvmassK1pi2[i] > 0.80 && vecdatainvmassK1pi2[i] < 0.95 ) && (vecdatainvmassK1pi1pi2[i] > 1.35 && vecdatainvmassK1pi1pi2[i] < 1.5) )
+      {
+	datainvMassKpivsKpipi->Fill(vecdatainvmassK1pi2[i], vecdatainvmassK1pi1pi2[i]);
+	datainvMassKpivsKpipiAngular->Fill(vecdatacosthetaK1pi2ViaK1pi1pi2[i]);
+      }
+    if( (vecdatainvmassK2pi1[i] > 0.80 && vecdatainvmassK2pi1[i] < 0.95 ) && (vecdatainvmassK2pi1pi2[i] > 1.35 && vecdatainvmassK2pi1pi2[i] < 1.5) )
+      {
+	datainvMassKpivsKpipi->Fill(vecdatainvmassK2pi1[i], vecdatainvmassK2pi1pi2[i]);
+	datainvMassKpivsKpipiAngular->Fill(vecdatacosthetaK2pi1ViaK2pi1pi2[i]);
+      }
+    if( (vecdatainvmassK2pi2[i] > 0.80 && vecdatainvmassK2pi2[i] < 0.95 ) && (vecdatainvmassK2pi1pi2[i] > 1.35 && vecdatainvmassK2pi1pi2[i] < 1.5) )
+      {
+	datainvMassKpivsKpipi->Fill(vecdatainvmassK2pi2[i], vecdatainvmassK2pi1pi2[i]);
+	datainvMassKpivsKpipiAngular->Fill(vecdatacosthetaK2pi2ViaK2pi1pi2[i]);
+      }
+
+  } 
+  
+  
+  TCanvas* cangular4 = new TCanvas("cangular4","cangular4",1200,600);
+  cangular4->Divide(2,1);
+  cangular4->cd(1);
+  datainvMassKpivsKpipi->Draw("colz");
+  cangular4->cd(2);
+  datainvMassKpivsKpipiAngular->Draw();
+  
+
+
+  
+  for(i=0; i<numOfEntriesData; i++) {
+    if( !(vecdatainvmasspi1pi2[i]>0.9 && vecdatainvmasspi1pi2[i]<1.01)) {
+      if( (vecdatainvmassK1pi1[i]>1.35 && vecdatainvmassK1pi1[i]<1.5) && (vecdatainvmassK2pi2[i]>1.7 && vecdatainvmassK2pi2[i]<1.9) )
+	{
+	  datainvMassKpivsKpi_2->Fill(vecdatainvmassK1pi1[i],vecdatainvmassK2pi2[i]);
+	  datainvMassKpivsKpi_2->Fill(vecdatainvmassK2pi2[i],vecdatainvmassK1pi1[i]);
+	  datainvMassKpiAngular_2->Fill(vecdatacosthetapi2[i]);
+	}
+      if( (vecdatainvmassK1pi2[i]>1.35 && vecdatainvmassK1pi2[i]<1.5) && (vecdatainvmassK2pi1[i]>1.7 && vecdatainvmassK2pi1[i]<1.9) )
+	{
+	  datainvMassKpivsKpi_2->Fill(vecdatainvmassK1pi2[i],vecdatainvmassK2pi1[i]);
+	  datainvMassKpivsKpi_2->Fill(vecdatainvmassK2pi1[i],vecdatainvmassK1pi2[i]);
+	  datainvMassKpiAngular_2->Fill(vecdatacosthetapi1[i]);
+	}
+      if( (vecdatainvmassK2pi1[i]>1.35 && vecdatainvmassK2pi1[i]<1.5) && (vecdatainvmassK1pi2[i]>1.7 && vecdatainvmassK1pi2[i]<1.9) )    
+	{
+	  datainvMassKpivsKpi_2->Fill(vecdatainvmassK2pi1[i],vecdatainvmassK1pi2[i]);
+	  datainvMassKpivsKpi_2->Fill(vecdatainvmassK1pi2[i],vecdatainvmassK2pi1[i]);
+	  datainvMassKpiAngular_2->Fill(vecdatacosthetapi2[i]);
+	}
+      if( (vecdatainvmassK2pi2[i]>1.35 && vecdatainvmassK2pi2[i]<1.5) && (vecdatainvmassK1pi1[i]>1.7 && vecdatainvmassK1pi1[i]<1.9) )
+	{
+	  datainvMassKpivsKpi_2->Fill(vecdatainvmassK2pi2[i],vecdatainvmassK1pi1[i]);
+	  datainvMassKpivsKpi_2->Fill(vecdatainvmassK1pi1[i],vecdatainvmassK2pi2[i]);
+	  datainvMassKpiAngular_2->Fill(vecdatacosthetapi1[i]);
+	}
+    }
+  }       
+
+  TCanvas* cangular5 = new TCanvas("cangular5","cangular5",1200,600);
+  cangular5->Divide(2,1);
+  cangular5->cd(1);
+  datainvMassKpivsKpi_2->Draw("colz");
+  cangular5->cd(2);
+  datainvMassKpiAngular_2->Draw();
+  
+
+
+  
+  for(i=0; i<numOfEntriesData; i++) {
+    if( (vecdatainvmasspi1pi2[i]>0.94 && vecdatainvmasspi1pi2[i]<1.01) && (vecdatainvmassK1K2[i]>2.33 && vecdatainvmassK1K2[i]<2.42) )
+      {
+	datainvMasspipivsKK->Fill(vecdatainvmassK1K2[i],vecdatainvmasspi1pi2[i]);
+	datainvMasspipivsKKAngular->Fill(vecdatacosthetaK1ViaK1K2[i]);
+      }
+  }       
+
+  TCanvas* cangular6 = new TCanvas("cangular6","cangular6",1200,600);
+  cangular6->Divide(2,1);
+  cangular6->cd(1);
+  datainvMasspipivsKK->Draw("colz");
+  cangular6->cd(2);
+  datainvMasspipivsKKAngular->Draw();
+  
 
 }