void MCsimang()
{
gROOT->Reset();
const auto PbR = 7;
const auto PbN = 1.3;
const auto kboxsize = 5;
const auto events = 10000;
double a[events];
double b[events];
double c[events];
double m[events];
double n[events];
double o[events];
double x[events];
TCanvas c11("c11", "c11", 800, 800);
for (auto i=0; i <events; i++)
{
	a[i] = (gRandom->Rndm() - 0.5) * 2 * PbR;
	b[i] = (gRandom->Rndm() - 0.5) * 2 * PbR;
	c[i] = (gRandom->Rndm() - 0.5) * 2 * PbR;
	x[i] = sqrt(a[i]*a[i] + b[i]*b[i] + c[i]*c[i]);
		if (x[i] < PbR)
		{
			m[i] = a[i];
			n[i] = b[i];
			o[i] = c[i];
			
		}

}
c11.Divide(1, 2);
TGraph gr(400, m, n);
gr.SetMarkerStyle(20);
gr.SetMarkerColor(30);
//gr.SetMarkerSize(3);
c11.cd(1);
gr.Draw("AP");
TGraph2D grr(400, m, n, o);
c11.cd(2);
grr.Draw();
c11.Update();
system("pause");
	}