/* $Id$ */
+//__________________________________________________________________________
+////////////////////////////////////////////////////////////////////////////
+//
+// class AliHBTCrab
+//
// This class introduces the weight's calculation
// according to the Lednicky's algorithm.
// The detailed description of the algorithm can be found
// in comments to fortran code:
// fsiw.f, fsiini.f
+//
+// Piotr.Skowronski@cern.ch
+////////////////////////////////////////////////////////////////////////////
#ifndef ALIHBTCrab_H
#define ALIHBTCrab_H
#include "AliHBTWeights.h"
-class Complex;
-typedef Complex double_complex;
-
+#ifdef __DECCXX
+ #include <complex.h>
+#else
+ class Complex;
+ typedef Complex doublecomplex;
+#endif
+
+//#include <math.h>
+
class AliHBTPair;
class AliHBTCrab: public AliHBTWeights
{
public:
+ AliHBTCrab();
virtual ~AliHBTCrab(){fgCrab =0x0;}
static AliHBTCrab* Instance();
void Set();
- Double_t GetWeight(const AliHBTPair* partpair);
+ Double_t GetWeight(AliHBTPair* partpair);
void Init(Int_t pid1,Int_t pid2); //put the initial values in fortran commons fsiini, led_bldata
private:
- AliHBTCrab();
AliHBTCrab(const AliHBTCrab &/*source*/);
AliHBTCrab & operator=(const AliHBTCrab& /*source*/);
#ifdef __DECCXX
static const complex fgkCI;//complex (1,0)
#else
- static const double_complex fgkCI;//complex (1,0)
+ static const doublecomplex fgkCI;//complex (1,0)
#endif
static const Double_t fgkROOT2;//! some const
static const Double_t fgkWcons; //constant for fm->GeV conversion 1/0.1973
#ifdef __DECCXX
complex CGamma(complex c);
#else
- double_complex CGamma(double_complex c);
+ doublecomplex CGamma(doublecomplex c);
#endif
- static AliHBTCrab* fgCrab;
+ static AliHBTCrab* fgCrab; //pointer to instance of this class - singleton
ClassDef(AliHBTCrab,1)
};