First version of the long range correlation studies in the EbyE directory (Petr Naumenko)

[u/mrichter/AliRoot.git] / PWG2 / EBYE / LRC / Test / TestPtPt.C

{\r
    gSystem->Load("libPWG2LRC.so");\r
Int_t errybins = 20000;\r
	Double_t errymin = 0.000025,\r
		  errymax = 1.000025;\r
\r
 	TH1D *gaus = new TH1D("source","source hist",1000,-5,5);\r
    gaus->FillRandom("gaus",10000);\r
    TH2D *source = new TH2D("PtN","PtN Test",100,0,5,100,0,5);\r
    TH2D *err = new TH2D("err","errl hist",100,0,5,errybins,errymin,errymax);\r
	Double_t x, y;\r
    for (Int_t i=0;i<10000;i++) \r
	{   \r
	    x =(gaus->GetRandom()+5)/2.0;\r
	    y = 0.678 * (x + gaus->GetRandom());\r
	    //Create source 2D histogram with correlation coefficient 0.678\r
		source->Fill(x,y);\r
		err->Fill(x,1/((gaus->GetRandom()+5)));\r
    }\r
    //Create TNN class encapsulated NN correlation algorithms\r
    //2D histogram pass into TNN constructor\r
	AliLRCPtPt final1("name", source, 0.35, err);\r
\r
   TCanvas *c1 = new TCanvas("c1","c1",800,1000);\r
	c1->Divide(3);\r
	c1->cd(1);\r
	source->Draw();\r
	c1->cd(2);\r
	final1.Draw_abs();\r
	c1->cd(3);\r
	final1.Draw_rel();\r
	c1->cd();\r
\r
	AliLRCPtPt final2;\r
        final2.MakeHistogramm("name", source, 0.35, err);\r
\r
   TCanvas *c2 = new TCanvas("c2","c2",800,1000);\r
	c2->Divide(3);\r
	c2->cd(1);\r
	source->Draw();\r
	c2->cd(2);\r
	final2.Draw_abs();\r
	c2->cd(3);\r
	final2.Draw_rel();\r
	c2->cd();\r
\r
}