[u/mrichter/AliRoot.git] / EVGEN / AliGenGeVSim.h

#ifndef ALIGENGEVSIM_H
#define ALIGENGEVSIM_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

/* $Id$ */

////////////////////////////////////////////////////////////////////////////////
//
// AliGenGeVSim is a class implementing GeVSim event generator.
// 
// GeVSim is a simple Monte-Carlo event generator for testing detector and 
// algorythm performance especialy concerning flow and event-by-event studies
//
// In this event generator particles are generated from thermal distributions 
// without any dynamics and addicional constrains. Distribution parameters like
// multiplicity, particle type yields, inverse slope parameters, flow coeficients 
// and expansion velocities are expleicite defined by the user.
//
// GeVSim contains four thermal distributions the same as
// MevSim event generator developed for STAR experiment.
//
// In addition custom distributions can be used be the mean 
// either two dimensional formula (TF2), a two dimensional histogram or
// two one dimensional histograms.
//  
// Azimuthal distribution is deconvoluted from (Pt,Y) distribution
// and is described by two Fourier coefficients representing 
// Directed and Elliptic flow. 
// 
////////////////////////////////////////////////////////////////////////////////
//
// To apply flow to event ganerated by an arbitraly event generator
// refer to AliGenAfterBurnerFlow class.
//
////////////////////////////////////////////////////////////////////////////////
//
// For examples, parameters and testing macros refer to:
// http:/home.cern.ch/radomski
// 
// for more detailed description refer to ALICE NOTE
// "GeVSim Monte-Carlo Event Generator"
// S.Radosmki, P. Foka.
//  
// Author:
// Sylwester Radomski,
// GSI, March 2002
//  
// S.Radomski@gsi.de
//
////////////////////////////////////////////////////////////////////////////////
//
// Updated and revised: September 2002, S. Radomski, GSI
//
////////////////////////////////////////////////////////////////////////////////


class TFormula;
class TF1;
class TF2;
class TH1D;
class TH2D;
class TObjArray;
class AliGeVSimParticle;

#include "AliGenerator.h"


class AliGenGeVSim : public AliGenerator {

 public:
  
  AliGenGeVSim();
  AliGenGeVSim(Float_t psi, Bool_t isMultTotal = kTRUE);
  
  virtual ~AliGenGeVSim();
  
  /////////////////////////////////////////////////////////////////
  
  void AddParticleType(AliGeVSimParticle *part);
  void SetMultTotal(Bool_t isTotal = kTRUE);
  
  void Init();
  void Generate();
  
  /////////////////////////////////////////////////////////////////
  
 private:

  AliGenGeVSim(const AliGenGeVSim & ggs);
  AliGenGeVSim & operator=(const AliGenGeVSim & /*ggs*/);


  Int_t   fModel;            // Selected model (1-7)
  Float_t fPsi;              // Reaction Plane angle (0-2pi)
  Bool_t  fIsMultTotal;      // Mode od multiplicity: total, dN/dY

  TF1 *fPtFormula;           //! Pt formula for model (1)
  TF1 *fYFormula;            //! Y formula for model (1)
  TF2 *fPtYFormula[4];       //! Pt,Y formulae for model (2)-(4)
  TF1 *fPhiFormula;          //! phi formula 
  
  TFormula *fCurrentForm;    //! currently used formula
  TH1D *fHist[2];            //! two 1D histograms (fModel == 6)
  TH2D *fPtYHist;            //! two-dimensional histogram (fModel == 7)

  TObjArray *fPartTypes;     // Registered particles
  
  void InitFormula();
  void SetFormula(Int_t pdg);
  void AdjustFormula();
  void DetermineReactionPlane();
  void GetRandomPtY(Double_t &pt, Double_t &y);

  Float_t GetdNdYToTotal();
  
  Bool_t CheckPtYPhi(Float_t pt, Float_t y, Float_t phi) const;  // for histograms only
  Bool_t CheckAcceptance(Float_t p[3]);
  
  Float_t FindScaler(Int_t paramId, Int_t pdg);
  
  /////////////////////////////////////////////////////////////////


  ClassDef(AliGenGeVSim, 2)

};

#endif
Commit	Line	Data
	1	#ifndef ALIGENGEVSIM_H
	2	#define ALIGENGEVSIM_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	/* $Id$ */
	7
	8	////////////////////////////////////////////////////////////////////////////////
	9	//
	10	// AliGenGeVSim is a class implementing GeVSim event generator.
	11	//
	12	// GeVSim is a simple Monte-Carlo event generator for testing detector and
	13	// algorythm performance especialy concerning flow and event-by-event studies
	14	//
	15	// In this event generator particles are generated from thermal distributions
	16	// without any dynamics and addicional constrains. Distribution parameters like
	17	// multiplicity, particle type yields, inverse slope parameters, flow coeficients
	18	// and expansion velocities are expleicite defined by the user.
	19	//
	20	// GeVSim contains four thermal distributions the same as
	21	// MevSim event generator developed for STAR experiment.
	22	//
	23	// In addition custom distributions can be used be the mean
	24	// either two dimensional formula (TF2), a two dimensional histogram or
	25	// two one dimensional histograms.
	26	//
	27	// Azimuthal distribution is deconvoluted from (Pt,Y) distribution
	28	// and is described by two Fourier coefficients representing
	29	// Directed and Elliptic flow.
	30	//
	31	////////////////////////////////////////////////////////////////////////////////
	32	//
	33	// To apply flow to event ganerated by an arbitraly event generator
	34	// refer to AliGenAfterBurnerFlow class.
	35	//
	36	////////////////////////////////////////////////////////////////////////////////
	37	//
	38	// For examples, parameters and testing macros refer to:
	39	// http:/home.cern.ch/radomski
	40	//
	41	// for more detailed description refer to ALICE NOTE
	42	// "GeVSim Monte-Carlo Event Generator"
	43	// S.Radosmki, P. Foka.
	44	//
	45	// Author:
	46	// Sylwester Radomski,
	47	// GSI, March 2002
	48	//
	49	// S.Radomski@gsi.de
	50	//
	51	////////////////////////////////////////////////////////////////////////////////
	52	//
	53	// Updated and revised: September 2002, S. Radomski, GSI
	54	//
	55	////////////////////////////////////////////////////////////////////////////////
	56
	57
	58	class TFormula;
	59	class TF1;
	60	class TF2;
	61	class TH1D;
	62	class TH2D;
	63	class TObjArray;
	64	class AliGeVSimParticle;
	65
	66	#include "AliGenerator.h"
	67
	68
	69	class AliGenGeVSim : public AliGenerator {
	70
	71	public:
	72
	73	AliGenGeVSim();
	74	AliGenGeVSim(Float_t psi, Bool_t isMultTotal = kTRUE);
	75
	76	virtual ~AliGenGeVSim();
	77
	78	/////////////////////////////////////////////////////////////////
	79
	80	void AddParticleType(AliGeVSimParticle *part);
	81	void SetMultTotal(Bool_t isTotal = kTRUE);
	82
	83	void Init();
	84	void Generate();
	85
	86	/////////////////////////////////////////////////////////////////
	87
	88	private:
	89
	90	AliGenGeVSim(const AliGenGeVSim & ggs);
	91	AliGenGeVSim & operator=(const AliGenGeVSim & /ggs/);
	92
	93
	94	Int_t fModel; // Selected model (1-7)
	95	Float_t fPsi; // Reaction Plane angle (0-2pi)
	96	Bool_t fIsMultTotal; // Mode od multiplicity: total, dN/dY
	97
	98	TF1 *fPtFormula; //! Pt formula for model (1)
	99	TF1 *fYFormula; //! Y formula for model (1)
	100	TF2 *fPtYFormula[4]; //! Pt,Y formulae for model (2)-(4)
	101	TF1 *fPhiFormula; //! phi formula
	102
	103	TFormula *fCurrentForm; //! currently used formula
	104	TH1D *fHist[2]; //! two 1D histograms (fModel == 6)
	105	TH2D *fPtYHist; //! two-dimensional histogram (fModel == 7)
	106
	107	TObjArray *fPartTypes; // Registered particles
	108
	109	void InitFormula();
	110	void SetFormula(Int_t pdg);
	111	void AdjustFormula();
	112	void DetermineReactionPlane();
	113	void GetRandomPtY(Double_t &pt, Double_t &y);
	114
	115	Float_t GetdNdYToTotal();
	116
	117	Bool_t CheckPtYPhi(Float_t pt, Float_t y, Float_t phi) const; // for histograms only
	118	Bool_t CheckAcceptance(Float_t p[3]);
	119
	120	Float_t FindScaler(Int_t paramId, Int_t pdg);
	121
	122	/////////////////////////////////////////////////////////////////
	123
	124
	125	ClassDef(AliGenGeVSim, 2)
	126
	127	};
	128
	129	#endif