const Double_t AliHBTCrab::fgkWcons = 1./0.1973;
const Double_t AliHBTCrab::fgkROOT2=1.41421356237309504880;
+#ifdef __DECCXX
+const complex AliHBTCrab::ci(0.0,1.0);
+#else
const double_complex AliHBTCrab::ci(0.0,1.0);
+#endif
/************************************************************/
double wsym_leftover,wanti_leftover,wnosym_leftover;
double qred,qdotr,r;
// const double rmass=MASS1*MASS2/(MASS1+MASS2);
+#ifdef __DECCXX
+ complex cphi1,cphi2,cphis,cphia;
+#else
double_complex cphi1,cphi2,cphis,cphia;
+#endif
arg=trueqdotr/197.323-2.0*TMath::Pi()*TMath::Floor(trueqdotr/(197.323*2.0*TMath::Pi()));
cphi1=exp(ci*arg);
return corr0;
}
+#ifdef __DECCXX
+complex AliHBTCrab::cgamma(complex c){
+#else
double_complex AliHBTCrab::cgamma(double_complex c){
+#endif
/* This calc.s gamma functions which are in the form gamma(n+i*y)
where n is an int and y is real. */
// This code is written by Scott Pratt
// taken from http://www.nscl.msu.edu/~pratt/freecodes/crab/home.html
+#ifdef __DECCXX
+ complex cg,cphase;
+#else
double_complex cg,cphase;
+#endif
int mm,j;
double x,y,phase,delp,cgmag;
x=real(c);
#include "AliHBTWeights.h"
-class Complex;
-typedef Complex double_complex;
+//class Complex;
+//typedef Complex double_complex;
+#ifdef __DECCXX
+#include <complex.h>
+#endif
+#include <math.h>
class AliHBTPair;
Float_t INTERACTION_DELK;
Int_t INTERACTION_NKMAX;/* number of momentum points in mesh for strong/coul. interaction */
+#ifdef __DECCXX
+ static const complex ci;
+#else
static const double_complex ci;
+#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);
+#endif
static AliHBTCrab* fgCrab;
ClassDef(AliHBTCrab,1)