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

/* $Id$ */

//__________________________________________________________________________
////////////////////////////////////////////////////////////////////////////
//
// class AliHBTCrab
//
// 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
//
// Piotr.Skowronski@cern.ch
////////////////////////////////////////////////////////////////////////////

#ifndef ALIHBTCrab_H
#define ALIHBTCrab_H

#include "AliHBTWeights.h"

#ifdef __DECCXX
 #include <complex.h>
#else
 class Complex;
 typedef Complex doublecomplex;
#endif
 
//#include <math.h>
 
class AliHBTPair;

class AliHBTCrab: public AliHBTWeights
 {
   public:

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

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

     void GetComQuantities(const AliHBTPair* pair, double *qred,double *r,double *qdotr,double *mom, int *test);
     double  CorrCalc(double trueqred,double trueqdotr,double truer);
     
     Bool_t fBreitWigner;//switch if to calculated BW
     Bool_t fReducedMom;//switch if
     Float_t fMaxMomentum;//switch if
     
     Bool_t  SetConfig(const AliHBTPair* pair);
     
     Int_t fPid1;//PID of the first particle
     Int_t fPid2;//PID of the second particle
     
     Double_t fMass1;//mass of the first particle
     Double_t fMass2;//mass of the second particle
 
     Float_t fInteractionWsym;// fractions of symmetric and antisym weights of the various spin channels 
     Float_t fInteractionWanti;//comment
     Float_t fInteractionWnosym;//comment
     
     Float_t fInteractionDelk;//comment
     Int_t fInteractionNkmax;// number of momentum points in mesh for strong/coul. interaction 
     
#ifdef __DECCXX
     static const complex fgkCI;//complex (1,0)
#else
     static const doublecomplex fgkCI;//complex (1,0)
#endif
     static const Double_t fgkROOT2;//! some const
     static const Double_t fgkWcons; //constant for fm->GeV conversion 1/0.1973
     
#ifdef __DECCXX
     complex CGamma(complex c);
#else
     doublecomplex CGamma(doublecomplex c);
#endif
     
     static AliHBTCrab* fgCrab; //pointer to instance of this class - singleton
     ClassDef(AliHBTCrab,1)
 };
 
#endif
Commit	Line	Data
	1	/* $Id$ */
	2
	3	//__________________________________________________________________________
	4	////////////////////////////////////////////////////////////////////////////
	5	//
	6	// class AliHBTCrab
	7	//
	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:
	12	// fsiw.f, fsiini.f
	13	//
	14	// Piotr.Skowronski@cern.ch
	15	////////////////////////////////////////////////////////////////////////////
	16
	17	#ifndef ALIHBTCrab_H
	18	#define ALIHBTCrab_H
	19
	20	#include "AliHBTWeights.h"
	21
	22	#ifdef __DECCXX
	23	#include <complex.h>
	24	#else
	25	class Complex;
	26	typedef Complex doublecomplex;
	27	#endif
	28
	29	//#include <math.h>
	30
	31	class AliHBTPair;
	32
	33	class AliHBTCrab: public AliHBTWeights
	34	{
	35	public:
	36
	37	AliHBTCrab();
	38	virtual ~AliHBTCrab(){fgCrab =0x0;}
	39	static AliHBTCrab* Instance();
	40	void Set();
	41
	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
	44
	45	private:
	46	AliHBTCrab(const AliHBTCrab &/source/);
	47	AliHBTCrab & operator=(const AliHBTCrab& /source/);
	48
	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);
	51
	52	Bool_t fBreitWigner;//switch if to calculated BW
	53	Bool_t fReducedMom;//switch if
	54	Float_t fMaxMomentum;//switch if
	55
	56	Bool_t SetConfig(const AliHBTPair* pair);
	57
	58	Int_t fPid1;//PID of the first particle
	59	Int_t fPid2;//PID of the second particle
	60
	61	Double_t fMass1;//mass of the first particle
	62	Double_t fMass2;//mass of the second particle
	63
	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
	67
	68	Float_t fInteractionDelk;//comment
	69	Int_t fInteractionNkmax;// number of momentum points in mesh for strong/coul. interaction
	70
	71	#ifdef __DECCXX
	72	static const complex fgkCI;//complex (1,0)
	73	#else
	74	static const doublecomplex fgkCI;//complex (1,0)
	75	#endif
	76	static const Double_t fgkROOT2;//! some const
	77	static const Double_t fgkWcons; //constant for fm->GeV conversion 1/0.1973
	78
	79	#ifdef __DECCXX
	80	complex CGamma(complex c);
	81	#else
	82	doublecomplex CGamma(doublecomplex c);
	83	#endif
	84
	85	static AliHBTCrab* fgCrab; //pointer to instance of this class - singleton
	86	ClassDef(AliHBTCrab,1)
	87	};
	88
	89	#endif