void liubianxing1()
{
double asyinzi=1;
	const int bin=40;
const double pitch=5e-3*asyinzi;
const double sid=2.88675e-3*asyinzi;
 double xstart=-0.1*asyinzi;
 double ystart=(-0.1-pitch/sqrt(3)/6*5)*asyinzi;

TH2Poly *hc = new TH2Poly();
   hc->SetTitle("Honeycomb example");
//   hc->Honeycomb(xstart,ystart,sid,bin,bin);
double xc,yc;double ii,jj;
   for(int i=0;i<10;i++){
xstart=-0.1+i*pitch;
ystart=-0.1;
	hc->Honeycomb(xstart,ystart,sid,bin,bin);	
cout<<xstart<<"\t"<<ystart<<endl;
	for(int j=0;j<1;j++){
	hc->Fill(0.,0.);
}
jj= hc->GetNumberOfBins();
                                for(int k=1;k<jj+1;k++){ii= hc->GetBinContent(k);
                                        if(ii>0){
                                                yc=0+(k/bin+1);
                                                if((k/bin)%2==0){xc=0+((k/bin)/2+k)%bin;}
                                                else{xc=0+((k/bin)/2+k)%bin;}
                                          cout<<xc<<"\t"<<yc<<endl;
                                        }
                                }
hc->Clear();
}
//hc->Draw("COLZ0");
}