2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved.
3 * See cxx source for full Copyright notice
7 /************************************
8 * Create an event and perform full *
9 * flow analysis 'on the fly'. *
11 * author: Ante Bilandzic *
12 * (abilandzic@gmail.com) *
13 ************************************/
15 #ifndef ALIFLOWEVENTSIMPLEMAKERONTHEFLY_H
16 #define ALIFLOWEVENTSIMPLEMAKERONTHEFLY_H
21 class AliFlowEventSimple;
22 class AliFlowTrackSimple;
23 class AliFlowTrackSimpleCuts;
25 class AliFlowEventSimpleMakerOnTheFly{
27 AliFlowEventSimpleMakerOnTheFly(UInt_t uiSeed = 0); // constructor
28 virtual ~AliFlowEventSimpleMakerOnTheFly(); // destructor
30 Bool_t AcceptOrNot(AliFlowTrackSimple *pTrack);
31 AliFlowEventSimple* CreateEventOnTheFly(AliFlowTrackSimpleCuts const *cutsRP, AliFlowTrackSimpleCuts const *cutsPOI);
32 // Setters and getters:
33 void SetMinMult(Int_t iMinMult) {this->fMinMult = iMinMult;}
34 Int_t GetMinMult() const {return this->fMinMult;}
35 void SetMaxMult(Int_t iMaxMult) {this->fMaxMult = iMaxMult;}
36 Int_t GetMaxMult() const {return this->fMaxMult;}
37 void SetMass(Double_t dMass) {this->fMass = dMass;}
38 Double_t GetMass() const {return this->fMass;}
39 void SetTemperature(Double_t dT) {this->fTemperature = dT;}
40 Double_t GetTemperature() const {return this->fTemperature;}
41 void SetV1(Double_t dV1) {this->fV1 = dV1;}
42 Double_t GetV1() const {return this->fV1;}
43 void SetV2(Double_t dV2) {this->fV2 = dV2;}
44 Double_t GetV2() const {return this->fV2;}
45 void SetV3(Double_t dV3) {this->fV3 = dV3;}
46 Double_t GetV3() const {return this->fV3;}
47 void SetV4(Double_t dV4) {this->fV4 = dV4;}
48 Double_t GetV4() const {return this->fV4;}
49 void SetV5(Double_t dV5) {this->fV5 = dV5;}
50 Double_t GetV5() const {return this->fV5;}
51 void SetV6(Double_t dV6) {this->fV6 = dV6;}
52 Double_t GetV6() const {return this->fV6;}
53 void SetUniformFluctuationsV2(Bool_t b) {this->fUniformFluctuationsV2 = b;}
54 Bool_t GetUniformFluctuationsV2() const {return this->fUniformFluctuationsV2;}
55 void SetMinV2(Double_t dMinV2) {this->fMinV2 = dMinV2;}
56 Double_t GetMinV2() const {return this->fMinV2;}
57 void SetMaxV2(Double_t dMaxV2) {this->fMaxV2 = dMaxV2;}
58 Double_t GetMaxV2() const {return this->fMaxV2;}
59 void SetPtDependentV2(Bool_t b) {this->fPtDependentV2 = b;}
60 Bool_t GetPtDependentV2() const {return this->fPtDependentV2;}
61 void SetV2vsPtCutOff(Double_t dV2vsPtCutOff) {this->fV2vsPtCutOff = dV2vsPtCutOff;}
62 Double_t GetV2vsPtCutOff() const {return this->fV2vsPtCutOff;}
63 void SetV2vsPtMax(Double_t dV2vsPtMax) {this->fV2vsPtMax = dV2vsPtMax;}
64 Double_t GetV2vsPtMax() const {return this->fV2vsPtMax;}
65 void SetSubeventEtaRange(Double_t minA, Double_t maxA, Double_t minB, Double_t maxB)
66 {this->fEtaMinA = minA;this->fEtaMaxA = maxA;this->fEtaMinB = minB;this->fEtaMaxB = maxB;};
67 void SetNTimes(Int_t nt) {this->fNTimes = nt;}
68 Int_t GetNTimes() const {return this->fNTimes;}
69 void SetUniformAcceptance(Bool_t ua) {this->fUniformAcceptance = ua;}
70 Bool_t GetUniformAcceptance() const {return this->fUniformAcceptance;}
71 void SetFirstSectorPhiMin(Double_t dPhiMin1) {this->fPhiMin1 = dPhiMin1;}
72 Double_t GetFirstSectorPhiMin() const {return this->fPhiMin1;}
73 void SetFirstSectorPhiMax(Double_t dPhiMax1) {this->fPhiMax1 = dPhiMax1;}
74 Double_t GetFirstSectorPhiMax() const {return this->fPhiMax1;}
75 void SetFirstSectorProbability(Double_t dProbability1) {this->fProbability1 = dProbability1;}
76 Double_t GetFirstSectorProbability() const {return this->fProbability1;}
77 void SetSecondSectorPhiMin(Double_t dPhiMin2) {this->fPhiMin2 = dPhiMin2;}
78 Double_t GetSecondSectorPhiMin() const {return this->fPhiMin2;}
79 void SetSecondSectorPhiMax(Double_t dPhiMax2) {this->fPhiMax2 = dPhiMax2;}
80 Double_t GetSecondSectorPhiMax() const {return this->fPhiMax2;}
81 void SetSecondSectorProbability(Double_t dProbability2) {this->fProbability2 = dProbability2;}
82 Double_t GetSecondSectorProbability() const {return this->fProbability2;}
84 AliFlowEventSimpleMakerOnTheFly(const AliFlowEventSimpleMakerOnTheFly& anAnalysis); // copy constructor
85 AliFlowEventSimpleMakerOnTheFly& operator=(const AliFlowEventSimpleMakerOnTheFly& anAnalysis); // assignment operator
86 Int_t fCount; // count number of events
87 Int_t fMinMult; // uniformly sampled multiplicity is >= iMinMult
88 Int_t fMaxMult; // uniformly sampled multiplicity is < iMaxMult
89 TF1 *fPtSpectra; // transverse momentum distribution (pt is sampled from hardwired Boltzmann distribution)
90 Double_t fMass; // mass in pt distribution (hardwired is Boltzmann pt distribution)
91 Double_t fTemperature; // "temperature" in pt distribution (hardwired is Boltzmann pt distribution)
92 TF1 *fPhiDistribution; // azimuthal distribution (phi is sampled from hardwired Fourier-like distribution)
93 Double_t fV1; // harmonic v1
94 Double_t fV2; // harmonic v2
95 Double_t fV3; // harmonic v3
96 Double_t fV4; // harmonic v4
97 Double_t fV5; // harmonic v5
98 Double_t fV6; // harmonic v6
99 Bool_t fUniformFluctuationsV2; // v2 is sampled uniformly for each event and for all particles from [fMinV2,fMaxV2]
100 Double_t fMinV2; // if v2 is sampled uniformly for each event, this is lower boundary on its value
101 Double_t fMaxV2; // if v2 is sampled uniformly for each event, this is upper boundary on its value
102 Bool_t fPtDependentV2; // v2 is pt-dependent
103 Double_t fV2vsPtCutOff; // if v2 is pt-dependent: for v2 < fV2vsPtCutOff v2 is growing linearly, otherwise v2 = fV2vsPtMax
104 Double_t fV2vsPtMax; // if v2 is pt-dependent: v2 = fV2vsPtMax for v2 >= fV2vsPtCutOff
105 Double_t fEtaMinA; // minimum eta of subevent A
106 Double_t fEtaMaxA; // maximum eta of subevent A
107 Double_t fEtaMinB; // minimum eta of subevent B
108 Double_t fEtaMaxB; // maximum eta of subevent B
109 Int_t fNTimes; // number of times to use the same particle in the analysis (simulating nonflow)
110 Bool_t fUniformAcceptance; // detector has uniform azimuthal acceptance or not
111 Double_t fPhiMin1; // first sector with non-uniform acceptance starts at azimuth fPhiMin1
112 Double_t fPhiMax1; // first sector with non-uniform acceptance ends at azimuth fPhiMax1
113 Double_t fProbability1; // particles emitted in fPhiMin1 < phi < fPhiMax1 are taken with probability fProbability1
114 Double_t fPhiMin2; // second sector with non-uniform acceptance starts at azimuth fPhiMin2
115 Double_t fPhiMax2; // second sector with non-uniform acceptance ends at azimuth fPhiMax2
116 Double_t fProbability2; // particles emitted in fPhiMin2 < phi < fPhiMax2 are taken with probability fProbability2
119 ClassDef(AliFlowEventSimpleMakerOnTheFly,1) // macro for rootcint