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