Do not include from subdirectories
[u/mrichter/AliRoot.git] / FASTSIM / PathLengthTest.C
1 void PathLengthTest() {
2   //
3   // It's a test macro for path length calculation in AliFastGlauber
4   //
5   // Andrea Dainese
6   //
7
8   AliFastGlauber g;
9
10   g.Init(2);
11
12   // centrality (fraction of geometrical cross section
13   g.SetCentralityClass(0.00,0.10); 
14
15   // plot b,ell, and ell1 vs ell2 (back-to-back partons) in that centrality 
16   // class
17   Int_t nEntries = 100;
18   g.PlotBDistr(nEntries);
19   g.PlotLengthDistr(nEntries);
20   g.PlotLengthB2BDistr(nEntries);
21
22
23   // examples on getting stuff out of it
24   Double_t b,ell,ell1,ell2;
25
26   g.GetRandomBHard(b);
27   printf(" Random b in cetrality class (according to hard cross section):\n %f fm\n",b);
28
29   printf(" Length:\n");
30   g.GetLength(ell,b);
31   printf("   for this b: %f fm\n",ell);
32   g.GetLength(ell);
33   printf("   for random b : %f fm\n",ell);
34
35   printf(" Lengths for two partons back-to-back:\n");
36   g.GetLengthsBackToBack(b,ell1,ell2);
37   printf("   for this b: %f fm and %f fm\n",ell1,ell2);
38   g.GetLengthsBackToBack(ell1,ell2);
39   printf("   for random b: %f fm and %f fm\n",ell1,ell2);
40
41
42   Double_t phis[3]={2.,4.,6.};
43   Double_t ells[3];
44   printf(" Lengths for N partons from PYTHIA:\n");
45   g.GetLengthsForPythia(3,phis,ells,b);
46   printf("   for this b: %f fm, %f fm and %f fm\n",ells[0],ells[1],ells[2]);
47   g.GetLengthsForPythia(3,phis,ells);
48   printf("   for random b: %f fm, %f fm and %f fm\n",ells[0],ells[1],ells[2]);
49
50   return;
51 }