coverity 15108 fixed
[u/mrichter/AliRoot.git] / FMD / flow / TestFlowSimple.C
1 Double_t GeneratePhi(Double_t psiR, Double_t v1, Double_t v2)
2 {
3   Double_t phi  =  gRandom->Uniform(0, 2 * TMath::Pi());
4   Double_t rel  =  phi - psiR;
5   Double_t dphi =  -2 * TMath::Sin(rel) * v1;
6   dphi          -= TMath::Sin(2 * rel) * v2;
7   phi           += dphi;
8   return phi;
9 }
10
11 void TestFlowSimple()
12 {
13   gSystem->Load("libFMDflow.so");
14   AliFMDFlowAxis      axis(10, -5, 5);
15   AliFMDFlowBinned1D  flow("flow", "analysed", 2, 1, axis);
16   TArrayD             phis(20000);
17   TArrayD             xs(20000);
18
19   for (Int_t i = 0; i < 100; i++) { 
20     Double_t psiR = gRandom->Uniform(0, 2*TMath::Pi());
21     Int_t    nObs = gRandom->Integer(20000);
22     
23     for (Int_t j = 0; j < nObs; j++) { 
24       xs[j]   = gRandom->Uniform(-5, 5);
25       phis[j] = GeneratePhi(psiR, 0, 0.05);
26     }
27     flow.Event(nObs, phis.fArray, xs.fArray);
28   }
29   flow.Draw("ht");
30 }
31