Transition to NewIO
[u/mrichter/AliRoot.git] / HBTAN / AliHBTLLWeights.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 ALIHBTLLWEIGHTS_H
10 #define ALIHBTLLWEIGHTS_H
11
12 #include <TObject.h>
13
14 class AliHBTPair;
15 class AliHBTLLWeights: public TObject
16  {
17    public:
18      virtual ~AliHBTLLWeights(){;}
19      static AliHBTLLWeights* Instance();
20      
21      Double_t GetWeight(const AliHBTPair* partpair); //get weight calculated by Lednicky's algorithm
22
23      void Init(); //put the initial values in fortran commons fsiini, led_bldata
24      void SetTest(Bool_t rtest = kTRUE);//Sets fTest member
25
26      void SetColoumb(Bool_t col = kTRUE);//: (ICH in fortran code) Coulomb interaction between the two particles ON (OFF)
27      void SetQuantumStatistics(Bool_t qss = kTRUE);//IQS: quantum statistics for the two particles ON (OFF) //if non-identical particles automatically off
28      void SetStrongInterSwitch(Bool_t sis = kTRUE);//ISI: strong interaction between the two particles ON (OFF)
29      void SetColWithResidNuclSwitch(Bool_t crn = kTRUE);//I3C: Coulomb interaction with residual nucleus ON (OFF)  
30      void SetApproxModel(Int_t ap);//sets  Model of Approximation (NS in Fortran code)
31      void SetRandomPosition(Bool_t rp = kTRUE); //ON=kTRUE(OFF=kFALSE)
32      void SetR1dw(Double_t R);   //spherical source model radii                                                                                                                                                                                                 
33      void SetParticlesTypes(Int_t pid1, Int_t pid2); //set AliRoot particles types   
34      void SetNucleusCharge(Double_t ch); // not used now  (see comments in fortran code)
35      void SetNucleusMass(Double_t mass); // (see comments in fortran code)
36
37    protected:
38      
39      Bool_t fTest;           //flag indicating if parameters listed below are to be calculated automatically (0)
40                              //or 
41      Bool_t fColoumbSwitch; 
42      Bool_t fQuantStatSwitch; 
43      Bool_t fStrongInterSwitch;//Switches strong interactions TRUE=ON
44      Bool_t fColWithResidNuclSwitch;//Switches couloumb interaction 
45                                     //with residual nucleus TRUE=ON          
46      Double_t fNuclMass;   //mass of nucleus
47      Double_t fNuclCharge; //charge of nucleus
48      
49      Bool_t   fRandomPosition;//flag indicating if positions of particles shuold be randomized according to parameters below
50      Double_t fRadius;        //radius used for randomizing particle vertex position
51
52      
53      Int_t fApproximationModel; //approximation used to calculate Bethe-Salpeter amplitude see SetApproxModel for more comments
54
55      Int_t fPID1;
56      Int_t fPID2;
57      Double_t fSigma; //constants for spherical source model in cm
58      static const Double_t fgkWcons; //constant for fm->GeV conversion 1/0.1973
59
60      static  AliHBTLLWeights *fgLLWeights;// pointer to wrapper of Fortran Lednicky code
61
62      static Int_t GetPairCode(Int_t pid1,Int_t pid2);
63      static Int_t GetPairCode(const AliHBTPair* partpair);//calculate automatically internal code 
64
65    private:
66      AliHBTLLWeights();
67      ClassDef(AliHBTLLWeights,1)
68  };
69
70 #endif