[u/mrichter/AliRoot.git] / PWG2 / FLOW / AliFlowCommon / AliFlowEventSimpleMakerOnTheFly.h

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

/********************************** 
 * create an event and perform    *
 * flow analysis 'on the fly'     * 
 *                                * 
 * authors: Raimond Snellings     *
 *           (snelling@nikhef.nl) * 
 *          Ante Bilandzic        * 
 *           (anteb@nikhef.nl)    *
 *********************************/ 

#ifndef ALIFLOWEVENTSIMPLEMAKERONTHEFLY_H
#define ALIFLOWEVENTSIMPLEMAKERONTHEFLY_H

class TF1;
class TRandom3;

#include "AliFlowEventSimple.h"  //needed as include
    
class AliFlowEventSimpleMakerOnTheFly {

 public:

  AliFlowEventSimpleMakerOnTheFly(UInt_t);    // constructor
  virtual ~AliFlowEventSimpleMakerOnTheFly(); // destructor

  virtual void Init(); 
  
  AliFlowEventSimple* CreateEventOnTheFly();  // create an event on the fly
 
    
  //                        *****************************
  //                        **** SETTERS AND GETTERS ****
  //                        *****************************
  //................................................................................................
  // setters and getters for global parameters:
  void SetMultDistrOfRPsIsGauss(Bool_t const mdorig) {this->fMultDistrOfRPsIsGauss = mdorig;};
  Bool_t GetMultDistrOfRPsIsGauss() const {return this->fMultDistrOfRPsIsGauss;};
  
  void SetMultiplicityOfRP(Int_t multRP) {this->fMultiplicityOfRP = multRP;}
  Int_t GetMultiplicityOfRP() const {return this->fMultiplicityOfRP;} 
  
  void SetMultiplicitySpreadOfRP(Double_t multSpreadRP) {this->fMultiplicitySpreadOfRP = multSpreadRP;}
  Double_t GetMultiplicitySpreadOfRP() const {return this->fMultiplicitySpreadOfRP;} 
  
  void SetMinMultOfRP(Int_t minmr) {this->fMinMultOfRP = minmr;}
  Int_t GetMinMultOfRP() const {return this->fMinMultOfRP;} 
  
  void SetMaxMultOfRP(Int_t maxmr) {this->fMaxMultOfRP = maxmr;}
  Int_t GetMaxMultOfRP() const {return this->fMaxMultOfRP;} 
  
  void SetTemperatureOfRP(Double_t temperatureRP) {this->fTemperatureOfRP = temperatureRP;}
  Double_t GetTemperatureOfRP() const {return this->fTemperatureOfRP;} 
  
  void SetUseConstantHarmonics(Bool_t const uch) {this->fUseConstantHarmonics = uch;};
  Bool_t GetUseConstantHarmonics() const {return this->fUseConstantHarmonics;};
  
  // constant harmonics:  
  void SetV1RP(Double_t dV1RP) {this->fV1RP = dV1RP;}
  Double_t GetV1RP() const {return this->fV1RP;} 
  
  void SetV1SpreadRP(Double_t dV1SpreadRP) {this->fV1SpreadRP = dV1SpreadRP;}
  Double_t GetV1SpreadRP() const {return this->fV1SpreadRP;} 
  
  void SetV2DistrOfRPsIsGauss(Bool_t const v2dorig) {this->fV2DistrOfRPsIsGauss = v2dorig;};
  Bool_t GetV2DistrOfRPsIsGauss() const {return this->fV2DistrOfRPsIsGauss;};
  
  void SetV2RP(Double_t dV2RP) {this->fV2RP = dV2RP;}
  Double_t GetV2RP() const {return this->fV2RP;} 
  
  void SetV2SpreadRP(Double_t dV2SpreadRP) {this->fV2SpreadRP = dV2SpreadRP;}
  Double_t GetV2SpreadRP() const {return this->fV2SpreadRP;} 
  
  void SetMinV2RP(Double_t dMinV2RP) {this->fMinV2RP = dMinV2RP;}
  Double_t GetMinV2RP() const {return this->fMinV2RP;} 
  
  void SetMaxV2RP(Double_t dMaxV2RP) {this->fMaxV2RP = dMaxV2RP;}
  Double_t GetMaxV2RP() const {return this->fMaxV2RP;} 
  
  void SetV4RP(Double_t dV4RP) {this->fV4RP = dV4RP;}
  Double_t GetV4RP() const {return this->fV4RP;} 
  
  void SetV4SpreadRP(Double_t dV4SpreadRP) {this->fV4SpreadRP = dV4SpreadRP;}
  Double_t GetV4SpreadRP() const {return this->fV4SpreadRP;} 
  
  // (pt,eta) dependent harmonics:
  void SetV2RPMax(Double_t dV2RPMax) {this->fV2RPMax = dV2RPMax;}
  Double_t GetV2RPMax() const {return this->fV2RPMax;} 
  
  void SetPtCutOff(Double_t dPtCutOff) {this->fPtCutOff = dPtCutOff;}
  Double_t GetPtCutOff() const {return this->fPtCutOff;} 
  
  void SetFirstSectorPhiMin(Double_t dPhiMin1) {this->fPhiMin1 = dPhiMin1;}
  Double_t GetFirstSectorPhiMin() const {return this->fPhiMin1;} 
  
  void SetFirstSectorPhiMax(Double_t dPhiMax1) {this->fPhiMax1 = dPhiMax1;}
  Double_t GetFirstSectorPhiMax() const {return this->fPhiMax1;}
  
  void SetFirstSectorProbability(Double_t dProbability1) {this->fProbability1 = dProbability1;}
  Double_t GetFirstProbability() const {return this->fProbability1;}  
  
  void SetSecondSectorPhiMin(Double_t dPhiMin2) {this->fPhiMin2 = dPhiMin2;}
  Double_t GetSecondSectorPhiMin() const {return this->fPhiMin2;} 
  
  void SetSecondSectorPhiMax(Double_t dPhiMax2) {this->fPhiMax2 = dPhiMax2;}
  Double_t GetSecondSectorPhiMax() const {return this->fPhiMax2;}
  
  void SetSecondSectorProbability(Double_t dProbability2) {this->fProbability2 = dProbability2;}
  Double_t GetSecondProbability() const {return this->fProbability2;}  
  //................................................................................................
  
  void SetNoOfLoops(Int_t noofl) {this->fNoOfLoops = noofl;}
  Int_t GetNoOfLoops() const {return this->fNoOfLoops;} 

  void SetSubeventEtaRange(Double_t minA,Double_t maxA,Double_t minB,Double_t maxB) 
  {this->fEtaMinA = minA; this->fEtaMaxA = maxA;this->fEtaMinB = minB; this->fEtaMaxB = maxB;};

 private:
 
  AliFlowEventSimpleMakerOnTheFly(const AliFlowEventSimpleMakerOnTheFly& anAnalysis);            // copy constructor
  AliFlowEventSimpleMakerOnTheFly& operator=(const AliFlowEventSimpleMakerOnTheFly& anAnalysis); // assignment operator
  
  //................................................................................................
  // global parameters:
  Bool_t    fMultDistrOfRPsIsGauss;  // 1.) if kTRUE  = multiplicitiy of RPs is sampled e-b-e from Gaussian distribution with
                                     //                 mean = fMultiplicityOfRP and spread = fMultiplicitySpreadOfRP
                                     // 2.) if kFALSE = multiplicitiy of RPs is sampled e-b-e uniformly from 
                                     //                 interval [fMinMultOfRP,fMaxMultOfRP]
  Int_t     fMultiplicityOfRP;       // mean multiplicity of RPs (if sampled from Gaussian)
  Double_t  fMultiplicitySpreadOfRP; // multiplicity spread of RPs (if sampled from Gaussian)
  Int_t     fMinMultOfRP;            // minimal multiplicity of RPs (if sampled uniformly)
  Int_t     fMaxMultOfRP;            // maximum multiplicity of RPs (if sampled uniformly)
  Double_t  fTemperatureOfRP;        // "temperature" of RPs in GeV/c (increase this parameter to get more high pt RPs) 
  Bool_t    fUseConstantHarmonics;   // harmonics V1, V2, V4... are constant (kTRUE) or functions of pt and eta (kFALSE)     
  // constant harmonics: 
  Double_t  fV1RP;                   // directed flow of RPs
  Double_t  fV1SpreadRP;             // directed flow spread of RPs
  
  Bool_t    fV2DistrOfRPsIsGauss;    // 1.) if kTRUE  = elliptic flow of RPs is sampled e-b-e from Gaussian distribution with
                                     //                 mean = fV2RP and spread = fV2SpreadRP
                                     // 2.) if kFALSE = elliptic flow of RPs is sampled e-b-e uniformly from 
                                     //                 interval [fMinV2RP,fMaxV2RP]
  Double_t  fV2RP;                   // mean elliptic flow of RPs (if sampled from Gaussian)
  Double_t  fV2SpreadRP;             // elliptic flow spread of RPs (if sampled from Gaussian)
  Double_t  fMinV2RP;                // minimal elliptic flow of RPs (if sampled uniformly)
  Double_t  fMaxV2RP;                // minimal elliptic flow of RPs (if sampled uniformly)
  
  Double_t  fV4RP;                   // harmonic V4 of RPs
  Double_t  fV4SpreadRP;             // harmonic V4's spread of RPs
  // (pt,eta) dependent harmonics:
  Double_t  fV2RPMax;                // elliptic flow of RPs
  Double_t  fPtCutOff;               // elliptic flow spread of RPs
  // non-uniform acceptance:
  Double_t  fPhiMin1;                // first non-uniform sector starts at azimuth fPhiMin1
  Double_t  fPhiMax1;                // first non-uniform sector ends at azimuth fPhiMax1
  Double_t  fProbability1;           // particles emitted in fPhiMin1 < phi < fPhiMax1 are taken with probability fProbability1 
  Double_t  fPhiMin2;                // second non-uniform sector starts at azimuth fPhiMin2
  Double_t  fPhiMax2;                // second non-uniform sector starts at azimuth fPhiMax2
  Double_t  fProbability2;           // particles emitted in fPhiMin2 < phi < fPhiMax2 are taken with probability fProbability2
  //................................................................................................
  
  //................................................................................................
  // equations for distributions: 
  TF1*      fPtSpectra;  // transverse momentum distribution
  TF1*      fPhiDistribution; // azimuthal distribution
  //................................................................................................
  
  TRandom3* fMyTRandom3; // our TRandom3 generator
  Int_t     fCount;      // count number of events 
  Int_t     fNoOfLoops;  // number of times to use the same particle (nonflow)
  Double_t  fEtaMinA;              // minimum eta of subevent A eta range
  Double_t  fEtaMaxA;              // maximum eta of subevent A eta range
  Double_t  fEtaMinB;              // minimum eta of subevent B eta range
  Double_t  fEtaMaxB;              // maximum eta of subevent B eta range  

  ClassDef(AliFlowEventSimpleMakerOnTheFly,0) // macro for rootcint
};
 
#endif
Commit	Line	Data
93ff27bd	1	/*
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved.
	3	* See cxx source for full Copyright notice
	4	* $Id$
	5	*/
	6
	7	/**********************************
	8	* create an event and perform *
	9	* flow analysis 'on the fly' *
	10	* *
	11	* authors: Raimond Snellings *
	12	* (snelling@nikhef.nl) *
	13	* Ante Bilandzic *
	14	* (anteb@nikhef.nl) *
	15	*********************************/
	16
	17	#ifndef ALIFLOWEVENTSIMPLEMAKERONTHEFLY_H
	18	#define ALIFLOWEVENTSIMPLEMAKERONTHEFLY_H
	19
	20	class TF1;
5b8dccde	21	class TRandom3;
93ff27bd	22
	23	#include "AliFlowEventSimple.h" //needed as include
	24
	25	class AliFlowEventSimpleMakerOnTheFly {
	26
	27	public:
	28
5b8dccde	29	AliFlowEventSimpleMakerOnTheFly(UInt_t); // constructor
93ff27bd	30	virtual ~AliFlowEventSimpleMakerOnTheFly(); // destructor
93ff27bd	31
4d603348	32	virtual void Init();
4d603348	33
93ff27bd	34	AliFlowEventSimple* CreateEventOnTheFly(); // create an event on the fly
	35
	36
	37	// *****************************
	38	// ** SETTERS AND GETTERS **
	39	// *****************************
	40	//................................................................................................
	41	// setters and getters for global parameters:
678f8252	42	void SetMultDistrOfRPsIsGauss(Bool_t const mdorig) {this->fMultDistrOfRPsIsGauss = mdorig;};
	43	Bool_t GetMultDistrOfRPsIsGauss() const {return this->fMultDistrOfRPsIsGauss;};
	44
93ff27bd	45	void SetMultiplicityOfRP(Int_t multRP) {this->fMultiplicityOfRP = multRP;}
93ff27bd	46	Int_t GetMultiplicityOfRP() const {return this->fMultiplicityOfRP;}
93ff27bd	47
93ff27bd	48	void SetMultiplicitySpreadOfRP(Double_t multSpreadRP) {this->fMultiplicitySpreadOfRP = multSpreadRP;}
93ff27bd	49	Double_t GetMultiplicitySpreadOfRP() const {return this->fMultiplicitySpreadOfRP;}
4d603348	50
678f8252	51	void SetMinMultOfRP(Int_t minmr) {this->fMinMultOfRP = minmr;}
	52	Int_t GetMinMultOfRP() const {return this->fMinMultOfRP;}
	53
	54	void SetMaxMultOfRP(Int_t maxmr) {this->fMaxMultOfRP = maxmr;}
	55	Int_t GetMaxMultOfRP() const {return this->fMaxMultOfRP;}
	56
cf90787f	57	void SetTemperatureOfRP(Double_t temperatureRP) {this->fTemperatureOfRP = temperatureRP;}
	58	Double_t GetTemperatureOfRP() const {return this->fTemperatureOfRP;}
	59
4d603348	60	void SetUseConstantHarmonics(Bool_t const uch) {this->fUseConstantHarmonics = uch;};
	61	Bool_t GetUseConstantHarmonics() const {return this->fUseConstantHarmonics;};
	62
	63	// constant harmonics:
e8c7e66c	64	void SetV1RP(Double_t dV1RP) {this->fV1RP = dV1RP;}
	65	Double_t GetV1RP() const {return this->fV1RP;}
	66
	67	void SetV1SpreadRP(Double_t dV1SpreadRP) {this->fV1SpreadRP = dV1SpreadRP;}
	68	Double_t GetV1SpreadRP() const {return this->fV1SpreadRP;}
	69
55381065	70	void SetV2DistrOfRPsIsGauss(Bool_t const v2dorig) {this->fV2DistrOfRPsIsGauss = v2dorig;};
	71	Bool_t GetV2DistrOfRPsIsGauss() const {return this->fV2DistrOfRPsIsGauss;};
	72
e8c7e66c	73	void SetV2RP(Double_t dV2RP) {this->fV2RP = dV2RP;}
	74	Double_t GetV2RP() const {return this->fV2RP;}
	75
	76	void SetV2SpreadRP(Double_t dV2SpreadRP) {this->fV2SpreadRP = dV2SpreadRP;}
	77	Double_t GetV2SpreadRP() const {return this->fV2SpreadRP;}
4d603348	78
55381065	79	void SetMinV2RP(Double_t dMinV2RP) {this->fMinV2RP = dMinV2RP;}
	80	Double_t GetMinV2RP() const {return this->fMinV2RP;}
	81
	82	void SetMaxV2RP(Double_t dMaxV2RP) {this->fMaxV2RP = dMaxV2RP;}
	83	Double_t GetMaxV2RP() const {return this->fMaxV2RP;}
	84
4d603348	85	void SetV4RP(Double_t dV4RP) {this->fV4RP = dV4RP;}
	86	Double_t GetV4RP() const {return this->fV4RP;}
	87
	88	void SetV4SpreadRP(Double_t dV4SpreadRP) {this->fV4SpreadRP = dV4SpreadRP;}
	89	Double_t GetV4SpreadRP() const {return this->fV4SpreadRP;}
	90
	91	// (pt,eta) dependent harmonics:
	92	void SetV2RPMax(Double_t dV2RPMax) {this->fV2RPMax = dV2RPMax;}
	93	Double_t GetV2RPMax() const {return this->fV2RPMax;}
	94
	95	void SetPtCutOff(Double_t dPtCutOff) {this->fPtCutOff = dPtCutOff;}
	96	Double_t GetPtCutOff() const {return this->fPtCutOff;}
e23cbd39	97
	98	void SetFirstSectorPhiMin(Double_t dPhiMin1) {this->fPhiMin1 = dPhiMin1;}
	99	Double_t GetFirstSectorPhiMin() const {return this->fPhiMin1;}
	100
	101	void SetFirstSectorPhiMax(Double_t dPhiMax1) {this->fPhiMax1 = dPhiMax1;}
	102	Double_t GetFirstSectorPhiMax() const {return this->fPhiMax1;}
	103
	104	void SetFirstSectorProbability(Double_t dProbability1) {this->fProbability1 = dProbability1;}
	105	Double_t GetFirstProbability() const {return this->fProbability1;}
	106
	107	void SetSecondSectorPhiMin(Double_t dPhiMin2) {this->fPhiMin2 = dPhiMin2;}
	108	Double_t GetSecondSectorPhiMin() const {return this->fPhiMin2;}
	109
	110	void SetSecondSectorPhiMax(Double_t dPhiMax2) {this->fPhiMax2 = dPhiMax2;}
	111	Double_t GetSecondSectorPhiMax() const {return this->fPhiMax2;}
	112
	113	void SetSecondSectorProbability(Double_t dProbability2) {this->fProbability2 = dProbability2;}
	114	Double_t GetSecondProbability() const {return this->fProbability2;}
e8c7e66c	115	//................................................................................................
e8c7e66c	116
aaa62cf0	117	void SetNoOfLoops(Int_t noofl) {this->fNoOfLoops = noofl;}
	118	Int_t GetNoOfLoops() const {return this->fNoOfLoops;}
	119
809f4e8b	120	void SetSubeventEtaRange(Double_t minA,Double_t maxA,Double_t minB,Double_t maxB)
	121	{this->fEtaMinA = minA; this->fEtaMaxA = maxA;this->fEtaMinB = minB; this->fEtaMaxB = maxB;};
	122
93ff27bd	123	private:
	124
	125	AliFlowEventSimpleMakerOnTheFly(const AliFlowEventSimpleMakerOnTheFly& anAnalysis); // copy constructor
	126	AliFlowEventSimpleMakerOnTheFly& operator=(const AliFlowEventSimpleMakerOnTheFly& anAnalysis); // assignment operator
	127
	128	//................................................................................................
	129	// global parameters:
678f8252	130	Bool_t fMultDistrOfRPsIsGauss; // 1.) if kTRUE = multiplicitiy of RPs is sampled e-b-e from Gaussian distribution with
	131	// mean = fMultiplicityOfRP and spread = fMultiplicitySpreadOfRP
	132	// 2.) if kFALSE = multiplicitiy of RPs is sampled e-b-e uniformly from
	133	// interval [fMinMultOfRP,fMaxMultOfRP]
	134	Int_t fMultiplicityOfRP; // mean multiplicity of RPs (if sampled from Gaussian)
	135	Double_t fMultiplicitySpreadOfRP; // multiplicity spread of RPs (if sampled from Gaussian)
	136	Int_t fMinMultOfRP; // minimal multiplicity of RPs (if sampled uniformly)
	137	Int_t fMaxMultOfRP; // maximum multiplicity of RPs (if sampled uniformly)
cf90787f	138	Double_t fTemperatureOfRP; // "temperature" of RPs in GeV/c (increase this parameter to get more high pt RPs)
cf90787f	139	Bool_t fUseConstantHarmonics; // harmonics V1, V2, V4... are constant (kTRUE) or functions of pt and eta (kFALSE)
4d603348	140	// constant harmonics:
aaa62cf0	141	Double_t fV1RP; // directed flow of RPs
aaa62cf0	142	Double_t fV1SpreadRP; // directed flow spread of RPs
55381065	143
	144	Bool_t fV2DistrOfRPsIsGauss; // 1.) if kTRUE = elliptic flow of RPs is sampled e-b-e from Gaussian distribution with
	145	// mean = fV2RP and spread = fV2SpreadRP
	146	// 2.) if kFALSE = elliptic flow of RPs is sampled e-b-e uniformly from
	147	// interval [fMinV2RP,fMaxV2RP]
	148	Double_t fV2RP; // mean elliptic flow of RPs (if sampled from Gaussian)
	149	Double_t fV2SpreadRP; // elliptic flow spread of RPs (if sampled from Gaussian)
	150	Double_t fMinV2RP; // minimal elliptic flow of RPs (if sampled uniformly)
	151	Double_t fMaxV2RP; // minimal elliptic flow of RPs (if sampled uniformly)
	152
4d603348	153	Double_t fV4RP; // harmonic V4 of RPs
	154	Double_t fV4SpreadRP; // harmonic V4's spread of RPs
	155	// (pt,eta) dependent harmonics:
	156	Double_t fV2RPMax; // elliptic flow of RPs
	157	Double_t fPtCutOff; // elliptic flow spread of RPs
e23cbd39	158	// non-uniform acceptance:
	159	Double_t fPhiMin1; // first non-uniform sector starts at azimuth fPhiMin1
	160	Double_t fPhiMax1; // first non-uniform sector ends at azimuth fPhiMax1
	161	Double_t fProbability1; // particles emitted in fPhiMin1 < phi < fPhiMax1 are taken with probability fProbability1
	162	Double_t fPhiMin2; // second non-uniform sector starts at azimuth fPhiMin2
	163	Double_t fPhiMax2; // second non-uniform sector starts at azimuth fPhiMax2
	164	Double_t fProbability2; // particles emitted in fPhiMin2 < phi < fPhiMax2 are taken with probability fProbability2
93ff27bd	165	//................................................................................................
	166
	167	//................................................................................................
	168	// equations for distributions:
aaa62cf0	169	TF1* fPtSpectra; // transverse momentum distribution
aaa62cf0	170	TF1* fPhiDistribution; // azimuthal distribution
93ff27bd	171	//................................................................................................
93ff27bd	172
4d603348	173	TRandom3* fMyTRandom3; // our TRandom3 generator
4d603348	174	Int_t fCount; // count number of events
4d603348	175	Int_t fNoOfLoops; // number of times to use the same particle (nonflow)
809f4e8b	176	Double_t fEtaMinA; // minimum eta of subevent A eta range
	177	Double_t fEtaMaxA; // maximum eta of subevent A eta range
	178	Double_t fEtaMinB; // minimum eta of subevent B eta range
	179	Double_t fEtaMaxB; // maximum eta of subevent B eta range
	180
93ff27bd	181	ClassDef(AliFlowEventSimpleMakerOnTheFly,0) // macro for rootcint
	182	};
	183
	184	#endif
	185
	186
	187