void show_sim()
{
	TFile *fr = new TFile("sim.root");
	TTree *tr = (TTree*)fr->Get("cbmsim");

	TCanvas *c1 = new TCanvas();
	c1->Divide(3,2);

	c1->cd(1);
	tr->Draw("RTelescope1Si1Point.fY:RTelescope1Si1Point.fX", "", "box");
//	tr->Draw("RTelescope1Si1Point.fX", "", "box");

	c1->cd(3);
	tr->Draw("RTelescope1Si1Point.fELoss+RTelescope1Si2Point.fELoss+RTelescope1CsIPoint.fELoss",
			"RTelescope1Si1Point.fPID<3000 && RTelescope1Si2Point.fPID<3000 && RTelescope1CsIPoint.fPID<3000");

	c1->cd(2);
	tr->Draw("RTelescope1Si1Point.fPID", "RTelescope1Si1Point.fPID>800000000");

	c1->cd(3);
	tr->Draw("RTelescope1Si1Point.fPID", "RTelescope1Si1Point.fPID<800000000");

	c1->cd(4);
	tr->Draw("RTelescope1Si1Point.fELoss+RTelescope1Si2Point.fELoss+RTelescope1CsIPoint.fELoss", "");

	c1->cd(5);
		tr->Draw("RTelescope1Si1Point.fELoss+RTelescope1Si2Point.fELoss+RTelescope1CsIPoint.fELoss",
				"RTelescope1Si1Point.fPID<3000 && RTelescope1Si2Point.fPID<3000 && RTelescope1CsIPoint.fPID<3000");

	c1->cd(6);
	tr->Draw("RTelescope1Si1Point.fELoss+RTelescope1Si2Point.fELoss+RTelescope1CsIPoint.fELoss",
			"RTelescope1Si1Point.fPID>3000 && RTelescope1Si2Point.fPID>3000 && RTelescope1CsIPoint.fPID>3000");

}