void Test()
{

  float *vz  = new float[10];
  float *svz = new float[10];
  int   *id  = new int[10];
  
  cout<<"Before Sorting : " <<endl;
  for(int i=0;i<10;i++)
    {
      vz[i] = gRandom->Gaus(0,5);
      cout<<"vz : " <<vz[i]<<endl;
    }

  cout<<"\t"<<endl;
  TMath::Sort(10,vz,id,false);

  cout<<"After Sorting : " <<endl;
  for(int i=0;i<10;i++)
    {
      svz[i] = vz[id[i]];
      cout<<"Vz : " <<svz[i]<<endl;
    }
  cout<<"\t"<<endl;

  float CurrVz = -7.15398;
  Long64_t kId = TMath::BinarySearch(10,svz,CurrVz);
  cout<<"CurrVz : " <<CurrVz<<"\t"<<"closest vertex : "<<vz[id[kId]]<<"\t"<<kId<<endl;
  cout<<"\t"<<endl;

  for(int i=0;i<10;i++)
    {
      int fId = TMath::BinarySearch(10,svz,svz[i]);
      cout<<"NewVz : " <<svz[i]<<"\t"<<"Closest Vz : "<<vz[id[fId]]<<"\t"<<fId<<endl;
    } 

}