[u/mrichter/AliRoot.git] / HBTAN / AliHBTCrab.h

/* $Id$ */

// This class introduces the weight's calculation
// according to the Lednicky's algorithm.
// The detailed description of the algorithm can be found
// in comments to fortran code:
// fsiw.f, fsiini.f

#ifndef ALIHBTCrab_H
#define ALIHBTCrab_H

#include "AliHBTWeights.h"

class Complex;
typedef Complex double_complex;

class AliHBTPair;

class AliHBTCrab: public AliHBTWeights
 {
   public:

     virtual ~AliHBTCrab(){fgCrab =0x0;}
     static AliHBTCrab* Instance();
     void Set();

     Double_t GetWeight(const AliHBTPair* partpair);
     void Init(Int_t pid1,Int_t pid2); //put the initial values in fortran commons fsiini, led_bldata
   private:
     AliHBTCrab();
     AliHBTCrab(const AliHBTCrab &/*source*/);
     AliHBTCrab & operator=(const AliHBTCrab& /*source*/);

     void get_com_quantities(const AliHBTPair* pair, double *qred,double *r,double *qdotr,double *mom, int *test);
     double  corrcalc(double trueqred,double trueqdotr,double truer);
     
     Bool_t fBreitWigner;
     Bool_t fReducedMom;
     Float_t fMaxMomentum;
     
     Bool_t  SetConfig(const AliHBTPair* pair);
     
     Int_t fPid1;
     Int_t fPid2;
     
     Double_t MASS1;
     Double_t MASS2;
 
     Float_t INTERACTION_WSYM;/* fractions of symmetric and antisym weights of the various spin channels */
     Float_t INTERACTION_WANTI;
     Float_t INTERACTION_WNOSYM;
     
     Float_t INTERACTION_DELK;
     Int_t INTERACTION_NKMAX;/* number of momentum points in mesh for strong/coul. interaction */
     
     static const double_complex ci;
     static const Double_t fgkROOT2;//! some const
     static const Double_t fgkWcons; //constant for fm->GeV conversion 1/0.1973
     
     double_complex cgamma(double_complex c);
     
     static AliHBTCrab* fgCrab;
     ClassDef(AliHBTCrab,1)
 };
 
#endif
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
	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