void polarcolor(){

	TCanvas * c1 = new TCanvas("c1", "", 600, 600);
//	c1->SetFrameLineWidth(0);

	gStyle->SetOptStat(0);
//	gStyle->SetPalette(kDeepSea);

	TH2D* hist[2];

	hist[0] = new TH2D("h1","h1", 8, 0, TMath::Pi()*2, 4, 0, TMath::Pi()*2);
	hist[1] = new TH2D("h2","h2", 4, 0, TMath::Pi()*2, 4, 0, TMath::Pi()*2);

	double y1[8]={0.05,0.15,0.25,0.35,0.45,0.55,0.65,0.75};
	double y2[4]={0.05,     0.25,     0.45,     0.65};

	double x1[8]={0, TMath::Pi()/4, TMath::Pi()/2, TMath::Pi()*3/4, TMath::Pi(), TMath::Pi()*5/4, TMath::Pi()*6/4, TMath::Pi()*7/4};
	double x2[4]={0,                TMath::Pi()/2,                  TMath::Pi(),                  TMath::Pi()*6/4};	

	for(int i=0; i<8; i++){hist[0]->Fill(x1[i], TMath::Pi()/2, y1[i]);}
	for(int j=0; j<4; j++){hist[1]->Fill(x2[j], TMath::Pi(),   y2[j]);}

	hist[0]->Draw("colz");
	hist[1]->Draw("col same");

}