[u/mrichter/AliRoot.git] / PWG2 / FEMTOSCOPY / AliFemto / AliFemtoModelWeightGeneratorLednicky.h

///////////////////////////////////////////////////////////////////////////
//                                                                       //
// AliFemtoModelWeightGeneratorLednicky : the most advanced weight       //
// generator available. Supports a large number of different pair types  //
// and interaction types. Can calculate pair weights coming from         //
// quantum statistics, coulomb interation and strong interaction ot any  //
// combination of the three, as applicable.                              //
//                                                                       //
///////////////////////////////////////////////////////////////////////////

#ifndef ALIFEMTOMODELWEIGHTGENERATORLEDNICKY_H
#define ALIFEMTOMODELWEIGHTGENERATORLEDNICKY_H

#include "AliFemtoTypes.h"
#include "AliFemtoModelWeightGenerator.h"

class AliFemtoModelWeightGeneratorLednicky : public  AliFemtoModelWeightGenerator {
 public: 
// --- Constructor
  AliFemtoModelWeightGeneratorLednicky(); // call SetDefaultCalcPar
  AliFemtoModelWeightGeneratorLednicky(const AliFemtoModelWeightGeneratorLednicky &aWeight); // call SetDefaultCalcPar
// --- Destructor : nothing to explicitly delete
  AliFemtoModelWeightGeneratorLednicky& operator=(const AliFemtoModelWeightGeneratorLednicky& aWeight);
  ~AliFemtoModelWeightGeneratorLednicky();

  virtual Double_t GenerateWeight(AliFemtoPair *aPair);

  virtual void     SetPairType(Int_t aPairType);
  virtual void     SetPairTypeFromPair(AliFemtoPair *aPair);
  virtual Int_t    GetPairType() const;

  virtual Double_t GetKStar() const;
  virtual Double_t GetKStarOut() const;
  virtual Double_t GetKStarSide() const;
  virtual Double_t GetKStarLong() const;
  virtual Double_t GetRStar() const;
  virtual Double_t GetRStarOut() const;
  virtual Double_t GetRStarSide() const;
  virtual Double_t GetRStarLong() const;

  virtual AliFemtoModelWeightGenerator* Clone() const;

// --- Setting

// >>> Calculation mode
  void SetDefaultCalcPar(); // Default is CoulOn, QuantumOn, StrongOn, 3BodyOff, Square, T0ApproxOff
  void SetCoulOn();
  void SetCoulOff();

  void SetQuantumOn();
  void SetQuantumOff();
  void SetStrongOn();
  void SetStrongOff();
  void Set3BodyOn();
  void Set3BodyOff();
  void SetSphere(); // use Spherical wave approximation
  void SetSquare(); // use use Square potential (only for p-p and pi+Pi-) otherwise, use spherical wave approx
  void SetT0ApproxOff();//only with  Spherical wave Approximation - this is default mode
  void SetT0ApproxOn(); 
 
// Test Lambda parameters
  void PrintLambdas();
  
  void SetNuclCharge(const double aNuclCharge); // for 3-body calculation
  void SetNuclMass(const double aNuclMass);

  virtual AliFemtoString Report();

protected:
  // Fsi weight output
  double  fWei;  // normal weight
  double  fWein; // weight with nuclear influence
  double  fWeif; // weight
  double  fWeightDen; // weight for the denominator

  // Setting parameters
  int fItest;    // if set to 1 default parameters will be used

  //int mNs;
  int    fIch;        // switch coulomb interaction on/off
  int    fIqs;        // switch quantum statistics on/off
  int    fIsi;        // switch strong interaction on/off
  int    fI3c;        // switch 3rd body influence on/off
  double fNuclMass;   // mass of the third body
  double fNuclCharge; // charge of the third body

  bool   fSphereApp;       // use spherical approximation
  bool   fT0App;           // use square well approximation

  //Pair identification
  int       fLL;             // internal pair type code
  short     fNuclChargeSign; // sign of the 3rd body charge
  bool      fSwap;           // are particle in right order ? 
  int const fLLMax;          // number of supported pairs
  char**    fLLName;         // name of the system
  int *     fNumProcessPair; // number of process pairs of each type
  int       fNumbNonId;      // Number of unidentified pairs

  // Interface to the fortran functions
  void FsiInit();
  void FsiSetLL();
  void FsiNucl();
  bool SetPid(const int aPid1,const int aPid2);

#ifdef __ROOT__
  ClassDef(AliFemtoModelWeightGeneratorLednicky,1)
#endif
};

#endif
Commit	Line	Data
37e4484e	1	///////////////////////////////////////////////////////////////////////////
	2	// //
	3	// AliFemtoModelWeightGeneratorLednicky : the most advanced weight //
	4	// generator available. Supports a large number of different pair types //
	5	// and interaction types. Can calculate pair weights coming from //
	6	// quantum statistics, coulomb interation and strong interaction ot any //
	7	// combination of the three, as applicable. //
	8	// //
	9	///////////////////////////////////////////////////////////////////////////
	10
	11	#ifndef ALIFEMTOMODELWEIGHTGENERATORLEDNICKY_H
	12	#define ALIFEMTOMODELWEIGHTGENERATORLEDNICKY_H
	13
	14	#include "AliFemtoTypes.h"
	15	#include "AliFemtoModelWeightGenerator.h"
	16
	17	class AliFemtoModelWeightGeneratorLednicky : public AliFemtoModelWeightGenerator {
	18	public:
	19	// --- Constructor
	20	AliFemtoModelWeightGeneratorLednicky(); // call SetDefaultCalcPar
	21	AliFemtoModelWeightGeneratorLednicky(const AliFemtoModelWeightGeneratorLednicky &aWeight); // call SetDefaultCalcPar
	22	// --- Destructor : nothing to explicitly delete
	23	AliFemtoModelWeightGeneratorLednicky& operator=(const AliFemtoModelWeightGeneratorLednicky& aWeight);
	24	~AliFemtoModelWeightGeneratorLednicky();
	25
	26	virtual Double_t GenerateWeight(AliFemtoPair *aPair);
	27
	28	virtual void SetPairType(Int_t aPairType);
	29	virtual void SetPairTypeFromPair(AliFemtoPair *aPair);
	30	virtual Int_t GetPairType() const;
	31
	32	virtual Double_t GetKStar() const;
	33	virtual Double_t GetKStarOut() const;
	34	virtual Double_t GetKStarSide() const;
	35	virtual Double_t GetKStarLong() const;
	36	virtual Double_t GetRStar() const;
	37	virtual Double_t GetRStarOut() const;
	38	virtual Double_t GetRStarSide() const;
	39	virtual Double_t GetRStarLong() const;
	40
	41	virtual AliFemtoModelWeightGenerator* Clone() const;
	42
	43	// --- Setting
	44
	45	// >>> Calculation mode
	46	void SetDefaultCalcPar(); // Default is CoulOn, QuantumOn, StrongOn, 3BodyOff, Square, T0ApproxOff
	47	void SetCoulOn();
	48	void SetCoulOff();
	49
	50	void SetQuantumOn();
	51	void SetQuantumOff();
	52	void SetStrongOn();
	53	void SetStrongOff();
	54	void Set3BodyOn();
	55	void Set3BodyOff();
	56	void SetSphere(); // use Spherical wave approximation
	57	void SetSquare(); // use use Square potential (only for p-p and pi+Pi-) otherwise, use spherical wave approx
	58	void SetT0ApproxOff();//only with Spherical wave Approximation - this is default mode
	59	void SetT0ApproxOn();
	60
	61	// Test Lambda parameters
	62	void PrintLambdas();
	63
	64	void SetNuclCharge(const double aNuclCharge); // for 3-body calculation
65	void SetNuclMass(const double aNuclMass);
66
67	virtual AliFemtoString Report();
68
69	protected:
70	// Fsi weight output
71	double fWei; // normal weight
72	double fWein; // weight with nuclear influence
73	double fWeif; // weight
74	double fWeightDen; // weight for the denominator
75
76	// Setting parameters
77	int fItest; // if set to 1 default parameters will be used
78
79	//int mNs;
80	int fIch; // switch coulomb interaction on/off
81	int fIqs; // switch quantum statistics on/off
82	int fIsi; // switch strong interaction on/off
83	int fI3c; // switch 3rd body influence on/off
84	double fNuclMass; // mass of the third body
85	double fNuclCharge; // charge of the third body
86
87	bool fSphereApp; // use spherical approximation
88	bool fT0App; // use square well approximation
89
90	//Pair identification
91	int fLL; // internal pair type code
92	short fNuclChargeSign; // sign of the 3rd body charge
93	bool fSwap; // are particle in right order ?
94	int const fLLMax; // number of supported pairs
95	char** fLLName; // name of the system
96	int * fNumProcessPair; // number of process pairs of each type
97	int fNumbNonId; // Number of unidentified pairs
98
99	// Interface to the fortran functions
100	void FsiInit();
101	void FsiSetLL();
102	void FsiNucl();
103	bool SetPid(const int aPid1,const int aPid2);
104
105	#ifdef __ROOT__
106	ClassDef(AliFemtoModelWeightGeneratorLednicky,1)
107	#endif
108	};
109
110	#endif