PWGCF/EBYE/LRC/Test/TestPtPt.C

   1 {
   2     gSystem->Load("libPWG2LRC.so");
   3 Int_t errybins = 20000;
   4         Double_t errymin = 0.000025,
   5                   errymax = 1.000025;
   6
   7         TH1D *gaus = new TH1D("source","source hist",1000,-5,5);
   8     gaus->FillRandom("gaus",10000);
   9     TH2D *source = new TH2D("PtN","PtN Test",100,0,5,100,0,5);
  10     TH2D *err = new TH2D("err","errl hist",100,0,5,errybins,errymin,errymax);
  11         Double_t x, y;
  12     for (Int_t i=0;i<10000;i++)
  13         {
  14             x =(gaus->GetRandom()+5)/2.0;
  15             y = 0.678 * (x + gaus->GetRandom());
  16             //Create source 2D histogram with correlation coefficient 0.678
  17                 source->Fill(x,y);
  18                 err->Fill(x,1/((gaus->GetRandom()+5)));
  19     }
  20     //Create TNN class encapsulated NN correlation algorithms
  21     //2D histogram pass into TNN constructor
  22         AliLRCPtPt final1("name", source, 0.35, err);
  23
  24    TCanvas *c1 = new TCanvas("c1","c1",800,1000);
  25         c1->Divide(3);
  26         c1->cd(1);
  27         source->Draw();
  28         c1->cd(2);
  29         final1.Draw_abs();
  30         c1->cd(3);
  31         final1.Draw_rel();
  32         c1->cd();
  33
  34         AliLRCPtPt final2;
  35         final2.MakeHistogramm("name", source, 0.35, err);
  36
  37    TCanvas *c2 = new TCanvas("c2","c2",800,1000);
  38         c2->Divide(3);
  39         c2->cd(1);
  40         source->Draw();
  41         c2->cd(2);
  42         final2.Draw_abs();
  43         c2->cd(3);
  44         final2.Draw_rel();
  45         c2->cd();
  46
  47 }