Fix Coverity 10974-82
[u/mrichter/AliRoot.git] / PWG2 / FEMTOSCOPY / AliFemtoUser / AliFemtoYlm.h
1 ////////////////////////////////////////////////////////////////////////////////
2 //                                                                            //
3 // AliFemtoYlm - the class to calculate varous components of spherical        //
4 //  harmonics                                                                 //
5 //                                                                            //
6 // Authors: Adam Kisiel kisiel@mps.ohio-state.edu                             //
7 //                                                                            //
8 ////////////////////////////////////////////////////////////////////////////////
9
10 #ifndef ALIFEMTOYLM_H
11 #define ALIFEMTOYLM_H
12 #include <cstdlib>
13 #include <cmath>
14 #include <complex>
15 #include <TMath.h>
16
17 class AliFemtoYlm {
18  public:
19   AliFemtoYlm();
20   ~AliFemtoYlm();
21
22   AliFemtoYlm(const AliFemtoYlm& aYlm);
23   AliFemtoYlm& operator=(const AliFemtoYlm& aYlm);
24
25   static double Legendre(int ell, int emm, double ctheta);
26   static void   LegendreUpToYlm(int lmax, double ctheta, double *lbuf);
27
28   static std::complex<double> Ylm(int ell,int m,double theta,double phi);
29   static std::complex<double> Ylm(int ell, int m, double x, double y, double z);
30
31   static void YlmUpToL(int lmax, double x, double y, double z, std::complex<double> *ylms);
32   static void YlmUpToL(int lmax, double ctheta, double phi, std::complex<double> *ylms);
33
34   static double ReYlm(int ell, int m, double theta, double phi);
35   static double ReYlm(int ell, int m, double x, double y, double z);
36   static double ImYlm(int ell, int m, double theta, double phi);
37   static double ImYlm(int ell, int m, double x, double y, double z);
38
39   static void InitializeYlms();
40   
41  private:
42   static std::complex<double> Ceiphi(double phi);
43
44   static double *fgPrefactors;
45   static int    *fgPrefshift;
46   static int    *fgPlmshift;
47 };
48
49 #endif