{
  gStyle->SetPalette(1);
  gStyle->SetOptStat(0);

  Int_t nlon = 360;
  Int_t nlat = 180;
  Int_t nlat1 = 200;

  Float_t vlon[nlon+1];
  for (Int_t i=0; i<nlon+1; i++){
    vlon[i] = -180. + float(i)*360./float(nlon);
  }

  Float_t vlat[nlat+1];
  for (Int_t i=0; i<nlat+1; i++){
    vlat[i] =  -90. + float(i)*180./float(nlat); 
  }

  Float_t aaa;
  Float_t vlat1[nlat1+1];
  for (Int_t i=0; i<nlat1+1; i++){
    aaa = -1. + float(i)*2.0/float(nlat1);
    vlat1[i] =  TMath::ASin(aaa)*180./TMath::Pi(); 
  }

  TH2D *htest1 = new TH2D("htest1", "htest1", nlon, vlon, nlat, vlat);
  TH2D *htest2 = new TH2D("htest2", "htest2", nlon, vlon, nlat1, vlat1);

  TRandom *ran = new TRandom();
  ran.SetSeed(0);

  Float_t xlat, xlon;

  Int_t nevt = 100000;
  for (Int_t ieve=0; ieve<nevt; ieve++){
    xlon = -180. + ran.Rndm()*360.;
    xlat = -90. + ran.Rndm()*180.;
    htest1.Fill(xlon,xlat,1.);
    htest2.Fill(xlon,xlat,1.);
  }

  TCanvas *c1 = new TCanvas("c1", "Test sky maps", 0, 0, 900, 700);
  c1->SetFillStyle(4000);
  c1->SetFillColor(0);
  c1->SetBorderSize(2);
  c1->SetFrameBorderMode(0);
  c1->SetFrameBorderSize(2);
  c1->Divide(1,2);
  c1->cd(1);
  htest1->Draw("zaitoff");
  c1->cd(2);
  htest2->Draw("zaitoff");
  


}