HBTAN/AliHBTCrab.h

   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
  14 class Complex;
  15 typedef Complex double_complex;
  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
  29    private:
  30      AliHBTCrab();
  31      AliHBTCrab(const AliHBTCrab &/*source*/);
  32      AliHBTCrab & operator=(const AliHBTCrab& /*source*/);
  33
  34      void get_com_quantities(const AliHBTPair* pair, double *qred,double *r,double *qdotr,double *mom, int *test);
  35      double  corrcalc(double trueqred,double trueqdotr,double truer);
  36
  37      Bool_t fBreitWigner;
  38      Bool_t fReducedMom;
  39      Float_t fMaxMomentum;
  40
  41      Bool_t  SetConfig(const AliHBTPair* pair);
  42
  43      Int_t fPid1;
  44      Int_t fPid2;
  45
  46      Double_t MASS1;
  47      Double_t MASS2;
  48
  49      Float_t INTERACTION_WSYM;/* fractions of symmetric and antisym weights of the various spin channels */
  50      Float_t INTERACTION_WANTI;
  51      Float_t INTERACTION_WNOSYM;
  52
  53      Float_t INTERACTION_DELK;
  54      Int_t INTERACTION_NKMAX;/* number of momentum points in mesh for strong/coul. interaction */
  55
  56      static const double_complex ci;
  57      static const Double_t fgkROOT2;//! some const
  58      static const Double_t fgkWcons; //constant for fm->GeV conversion 1/0.1973
  59
  60      double_complex cgamma(double_complex c);
  61
  62      static AliHBTCrab* fgCrab;
  63      ClassDef(AliHBTCrab,1)
  64  };
  65
  66 #endif