]>
Commit | Line | Data |
---|---|---|
dd82cadc | 1 | /* $Id$ */ |
2 | ||
3 | // This class introduces the weight's calculation | |
4 | // according to the Lednicky's algorithm. | |
5 | // The detailed description of the algorithm can be found | |
6 | // in comments to fortran code: | |
7 | // fsiw.f, fsiini.f | |
8 | ||
9 | #ifndef ALIHBTCrab_H | |
10 | #define ALIHBTCrab_H | |
11 | ||
12 | #include "AliHBTWeights.h" | |
13 | ||
227b510f | 14 | class Complex; |
15 | typedef Complex double_complex; | |
dd82cadc | 16 | |
17 | class AliHBTPair; | |
18 | ||
19 | class AliHBTCrab: public AliHBTWeights | |
20 | { | |
21 | public: | |
22 | ||
23 | virtual ~AliHBTCrab(){fgCrab =0x0;} | |
24 | static AliHBTCrab* Instance(); | |
25 | void Set(); | |
26 | ||
27 | Double_t GetWeight(const AliHBTPair* partpair); | |
28 | void Init(Int_t pid1,Int_t pid2); //put the initial values in fortran commons fsiini, led_bldata | |
c114f4f8 | 29 | |
dd82cadc | 30 | private: |
31 | AliHBTCrab(); | |
32 | AliHBTCrab(const AliHBTCrab &/*source*/); | |
33 | AliHBTCrab & operator=(const AliHBTCrab& /*source*/); | |
34 | ||
c114f4f8 | 35 | void GetComQuantities(const AliHBTPair* pair, double *qred,double *r,double *qdotr,double *mom, int *test); |
36 | double CorrCalc(double trueqred,double trueqdotr,double truer); | |
dd82cadc | 37 | |
c114f4f8 | 38 | Bool_t fBreitWigner;//switch if to calculated BW |
39 | Bool_t fReducedMom;//switch if | |
40 | Float_t fMaxMomentum;//switch if | |
dd82cadc | 41 | |
42 | Bool_t SetConfig(const AliHBTPair* pair); | |
43 | ||
c114f4f8 | 44 | Int_t fPid1;//PID of the first particle |
45 | Int_t fPid2;//PID of the second particle | |
dd82cadc | 46 | |
c114f4f8 | 47 | Double_t fMass1;//mass of the first particle |
48 | Double_t fMass2;//mass of the second particle | |
dd82cadc | 49 | |
c114f4f8 | 50 | Float_t fInteractionWsym;// fractions of symmetric and antisym weights of the various spin channels |
51 | Float_t fInteractionWanti;//comment | |
52 | Float_t fInteractionWnosym;//comment | |
dd82cadc | 53 | |
c114f4f8 | 54 | Float_t fInteractionDelk;//comment |
55 | Int_t fInteractionNkmax;// number of momentum points in mesh for strong/coul. interaction | |
dd82cadc | 56 | |
8882f543 | 57 | #ifdef __DECCXX |
c114f4f8 | 58 | static const complex fgkCI;//complex (1,0) |
8882f543 | 59 | #else |
c114f4f8 | 60 | static const double_complex fgkCI;//complex (1,0) |
8882f543 | 61 | #endif |
dd82cadc | 62 | static const Double_t fgkROOT2;//! some const |
63 | static const Double_t fgkWcons; //constant for fm->GeV conversion 1/0.1973 | |
64 | ||
8882f543 | 65 | #ifdef __DECCXX |
c114f4f8 | 66 | complex CGamma(complex c); |
8882f543 | 67 | #else |
c114f4f8 | 68 | double_complex CGamma(double_complex c); |
8882f543 | 69 | #endif |
dd82cadc | 70 | |
71 | static AliHBTCrab* fgCrab; | |
72 | ClassDef(AliHBTCrab,1) | |
73 | }; | |
74 | ||
75 | #endif |