3 //__________________________________________________________________________
4 ////////////////////////////////////////////////////////////////////////////
8 // This class introduces the weight's calculation
9 // according to the Lednicky's algorithm.
10 // The detailed description of the algorithm can be found
11 // in comments to fortran code:
14 // Piotr.Skowronski@cern.ch
15 ////////////////////////////////////////////////////////////////////////////
20 #include "AliHBTWeights.h"
26 typedef Complex doublecomplex;
33 class AliHBTCrab: public AliHBTWeights
38 virtual ~AliHBTCrab(){fgCrab =0x0;}
39 static AliHBTCrab* Instance();
42 Double_t GetWeight(AliHBTPair* partpair);
43 void Init(Int_t pid1,Int_t pid2); //put the initial values in fortran commons fsiini, led_bldata
46 AliHBTCrab(const AliHBTCrab &/*source*/);
47 AliHBTCrab & operator=(const AliHBTCrab& /*source*/);
49 void GetComQuantities(const AliHBTPair* pair, double *qred,double *r,double *qdotr,double *mom, int *test);
50 double CorrCalc(double trueqred,double trueqdotr,double truer);
52 Bool_t fBreitWigner;//switch if to calculated BW
53 Bool_t fReducedMom;//switch if
54 Float_t fMaxMomentum;//switch if
56 Bool_t SetConfig(const AliHBTPair* pair);
58 Int_t fPid1;//PID of the first particle
59 Int_t fPid2;//PID of the second particle
61 Double_t fMass1;//mass of the first particle
62 Double_t fMass2;//mass of the second particle
64 Float_t fInteractionWsym;// fractions of symmetric and antisym weights of the various spin channels
65 Float_t fInteractionWanti;//comment
66 Float_t fInteractionWnosym;//comment
68 Float_t fInteractionDelk;//comment
69 Int_t fInteractionNkmax;// number of momentum points in mesh for strong/coul. interaction
72 static const complex fgkCI;//complex (1,0)
74 static const doublecomplex fgkCI;//complex (1,0)
76 static const Double_t fgkROOT2;//! some const
77 static const Double_t fgkWcons; //constant for fm->GeV conversion 1/0.1973
80 complex CGamma(complex c);
82 doublecomplex CGamma(doublecomplex c);
85 static AliHBTCrab* fgCrab; //pointer to instance of this class - singleton
86 ClassDef(AliHBTCrab,1)