  gStyle->SetOptStat(110);
  gStyle->SetCanvasColor(0);
  gStyle->SetPalette(1);

  //   TCanvas* pad = new TCanvas();
  //   pad->Divide(2, 2);

  //   pad->cd(1);
  gKMinvHist = new TH1D("Statistics","Kaons",50,0.0,2.0);
  gKMinvHist->SetLineColor(2);
  gKMinvHist->GetXaxis()->SetTitle("Invariant Mass (GeV/c^{2})");
  gKMinvHist->GetYaxis()->SetTitle("Counts");
  gKMinvHist->SetFillColor(0);
  //   gKMinvHist->Draw();



    Int_t nTracks = fVSD->fTreeR->GetEntries();

    for(Int_t i = nCascades + nV0Tracks; i < nTracks; i++)
    {

       fVSD->fTreeR->GetEntry(i);
       TEveTrack* track1 = new TEveTrack(&fVSD->fR, trkProp);

       for(Int_t j = i; j < nTracks; j++)
       {

         fVSD->fTreeR->GetEntry(j);
         TEveTrack* track2 = new TEveTrack(&fVSD->fR, trkProp);

         Double_t px1 = track1->GetESDTrack->Px();
         Double_t py1 = track1->GetESDTrack->Py();
         Double_t pz1 = track1->GetESDTrack->Pz();
         Double_t mass1 = track1->GetESDTrack->M();
         Double_t E1 = TMath::Sqrt(px1*px1+py1*py1+pz1*pz1+mass1*mass1);

         Double_t px2 = track2->GetESDTrack->Px();
         Double_t py2 = track2->GetESDTrack->Py();
         Double_t pz2 = track2->GetESDTrack->Pz();
         Double_t mass2 = track2->GetESDTrack->M();
         Double_t E2 = TMath::Sqrt(px2*px2+py2*py2+pz2*pz2+mass2*mass2);

         Double_t minv = TMath::Sqrt((E1+E2)*(E1+E2)-(px1+px2)*(px1+px2)-(py1+py2)*(py1+py2)-(pz1+pz2)*(pz1+pz2));

       }

    }

    for(Int_t i = nCascades + nV0Tracks; i < nTracks; i++)
    {

       fVSD->fTreeR->GetEntry(i);
       TEveTrack* track1 = new TEveTrack(&fVSD->fR, trkProp);

       for(Int_t j = i; j < nTracks; j++)
       {

         fVSD->fTreeR->GetEntry(j);
         TEveTrack* track2 = new TEveTrack(&fVSD->fR, trkProp);

         if(track1->GetESDTrack->Sign() != track1->GetESDTrack->Sign())
         {

           Double_t px1 = track1->GetESDTrack->Px();
           Double_t py1 = track1->GetESDTrack->Py();
           Double_t pz1 = track1->GetESDTrack->Pz();
           Double_t mass1 = track1->GetESDTrack->M();
           Double_t E1 = TMath::Sqrt(px1*px1+py1*py1+pz1*pz1+mass1*mass1);

           Double_t px2 = track2->GetESDTrack->Px();
           Double_t py2 = track2->GetESDTrack->Py();
           Double_t pz2 = track2->GetESDTrack->Pz();
           Double_t mass2 = track2->GetESDTrack->M();
           Double_t E2 = TMath::Sqrt(px2*px2+py2*py2+pz2*pz2+mass2*mass2);

           Double_t minv = TMath::Sqrt((E1+E2)*(E1+E2)-(px1+px2)*(px1+px2)-(py1+py2)*(py1+py2)-(pz1+pz2)*(pz1+pz2));

         }

       }

    }
