{
	TFile* file_1 = TFile::Open("cutTest.root");
	TH1D* zVtx = (TH1D*) file_1->Get("ZVertex");
	TH1D* zVtx_MC = (TH1D*) file_1->Get("ZVertex_MC");
	
	zVtx->Rebin(4);
	zVtx_MC->Rebin(4);
	
	zVtx_MC->Scale(zVtx->Integral()/zVtx_MC->Integral());
	zVtx->Divide(zVtx_MC);
	
	TCanvas myCanvas;
	myCanvas.Print("cutTest.ps[");
	
	zVtx->Draw();

	float leftCutPointsX[5] = {-32.,-32.,-40.,-40.,-32.};
	float leftCutPointsY[5] = {-12.,12.,12.,-12.,-12.};
	TCutG leftAreaCut("leftCut",5,leftCutPointsX,leftCutPointsY);	
	
	leftAreaCut.SetFillStyle(3005);
	leftAreaCut.Draw("f");

	float rightCutPointsX[5] = {32.,32.,40.,40.,32.};
	float rightCutPointsY[5] = {-12.,12.,12.,-12.,-12.};
	TCutG rightAreaCut("rightCut",5,rightCutPointsX,rightCutPointsY);	

	rightAreaCut.SetFillStyle(3005);
	rightAreaCut.Draw("f");

	myCanvas.Update();
	myCanvas.Print("cutTest.ps");
	myCanvas.Print("cutTest.ps]");
}