1 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
2 * See cxx source for full Copyright notice */
5 /*****************************************************************
6 AliFlowEventSimple: A simple event
9 origin: Naomi van der Kolk (kolk@nikhef.nl)
10 Ante Bilandzic (anteb@nikhef.nl)
11 Raimond Snellings (Raimond.Snellings@nikhef.nl)
12 mods: Mikolaj Krzewicki (mikolaj.krzewicki@cern.ch)
13 *****************************************************************/
15 #ifndef ALIFLOWEVENTSIMPLE_H
16 #define ALIFLOWEVENTSIMPLE_H
19 #include "TParameter.h"
24 class AliFlowTrackSimple;
25 class AliFlowTrackSimpleCuts;
27 class AliFlowEventSimple: public TObject {
31 enum ConstructionMethod {kEmpty,kGenerate};
34 AliFlowEventSimple( Int_t nParticles,
35 ConstructionMethod m=kEmpty,
38 Double_t phiMax=TMath::TwoPi(),
40 Double_t etaMax= 1.0 );
41 AliFlowEventSimple(TTree* anInput, const AliFlowTrackSimpleCuts* rpCuts, const AliFlowTrackSimpleCuts* poiCuts);
42 AliFlowEventSimple(const AliFlowEventSimple& anEvent);
43 AliFlowEventSimple& operator=(const AliFlowEventSimple& anEvent);
44 virtual ~AliFlowEventSimple();
46 Bool_t IsFolder() const {return kTRUE;};
47 void Browse(TBrowser *b);
48 void Print(Option_t* option = "") const; //method to print stats
50 Int_t NumberOfTracks() const { return fNumberOfTracks; }
51 Int_t GetReferenceMultiplicity() const { return fReferenceMultiplicity; }
52 void SetReferenceMultiplicity( Int_t m ) { fReferenceMultiplicity = m; }
53 Int_t GetEventNSelTracksRP() const { return fNumberOfRPs; }
54 void SetEventNSelTracksRP(Int_t nr) { fNumberOfRPs = nr; }
55 Int_t GetEventNSelTracksPOI() const { return fNumberOfPOIs; }
56 void SetEventNSelTracksPOI(Int_t np) { fNumberOfPOIs = np; }
57 Int_t GetNumberOfRPs() const { return fNumberOfRPs; }
58 void SetNumberOfRPs( Int_t nr ) { fNumberOfRPs=nr; }
59 Int_t GetNumberOfPOIs() const { return fNumberOfPOIs; }
60 void SetNumberOfPOIs( Int_t np ) { fNumberOfPOIs=np; }
61 Double_t GetMCReactionPlaneAngle() const { return fMCReactionPlaneAngle; }
62 void SetMCReactionPlaneAngle(Double_t fPhiRP) { fMCReactionPlaneAngle=fPhiRP; fMCReactionPlaneAngleIsSet=kTRUE; }
63 Bool_t IsSetMCReactionPlaneAngle() const { return fMCReactionPlaneAngleIsSet; }
64 void SetAfterBurnerPrecision(Double_t p) { fAfterBurnerPrecision=p; }
65 Double_t GetAfterBurnerPrecision() const { return fAfterBurnerPrecision; }
66 void SetUserModified(Bool_t s=kTRUE) { fUserModified=s; }
67 Bool_t IsUserModified() const { return fUserModified; }
69 void ResolutionPt(Double_t res);
70 void TagSubeventsInEta(Double_t etaMinA, Double_t etaMaxA, Double_t etaMinB, Double_t etaMaxB );
71 void TagSubeventsByCharge();
72 void TagRP(const AliFlowTrackSimpleCuts* cuts );
73 void TagPOI(const AliFlowTrackSimpleCuts* cuts );
74 void CloneTracks(Int_t n);
75 void AddV1( Double_t v1 );
76 void AddV2( Double_t v2 );
77 void AddV3( Double_t v3 );
78 void AddV4( Double_t v4 );
79 void AddV5( Double_t v5 );
80 void AddFlow( Double_t v1, Double_t v2, Double_t v3, Double_t v4, Double_t v5 );
81 void AddFlow(Double_t v1, Double_t v2, Double_t v3, Double_t v4, Double_t v5,
82 Double_t rp1, Double_t rp2, Double_t rp3, Double_t rp4, Double_t rp5 );
83 void AddV2( TF1* ptDepV2 );
84 void DefineDeadZone( Double_t etaMin, Double_t etaMax, Double_t phiMin, Double_t phiMax );
85 Int_t CleanUpDeadTracks();
88 static TF1* SimplePtSpectrum();
89 static TF1* SimplePtDepV2();
91 AliFlowTrackSimple* GetTrack(Int_t i);
92 void AddTrack( AliFlowTrackSimple* track );
93 AliFlowTrackSimple* MakeNewTrack();
95 AliFlowVector GetQ(Int_t n=2, TList *weightsList=NULL, Bool_t usePhiWeights=kFALSE, Bool_t usePtWeights=kFALSE, Bool_t useEtaWeights=kFALSE);
96 void Get2Qsub(AliFlowVector* Qarray, Int_t n=2, TList *weightsList=NULL, Bool_t usePhiWeights=kFALSE, Bool_t usePtWeights=kFALSE, Bool_t useEtaWeights=kFALSE);
99 virtual void Generate( Int_t nParticles,
102 Double_t phiMax=TMath::TwoPi(),
103 Double_t etaMin=-1.0,
104 Double_t etaMax= 1.0 );
107 TObjArray* fTrackCollection; //-> collection of tracks
108 Int_t fReferenceMultiplicity; // reference multiplicity
109 Int_t fNumberOfTracks; // number of tracks
110 Int_t fNumberOfRPs; // number of tracks that have passed the RP selection
111 Int_t fNumberOfPOIs; // number of tracks that have passed the POI selection
112 Double_t fMCReactionPlaneAngle; // the angle of the reaction plane from the MC truth
113 Bool_t fMCReactionPlaneAngleIsSet; // did we set it from MC?
114 Double_t fAfterBurnerPrecision; // iteration precision in afterburner
115 Bool_t fUserModified; // did we modify the event in any way (afterburner etc) ?
116 TParameter<Int_t>* fNumberOfTracksWrap; //! number of tracks in TBrowser
117 TParameter<Int_t>* fNumberOfRPsWrap; //! number of tracks that have passed the RP selection in TBrowser
118 TParameter<Int_t>* fNumberOfPOIsWrap; //! number of tracks that have passed the POI selection in TBrowser
119 TParameter<Double_t>* fMCReactionPlaneAngleWrap; //! the angle of the reaction plane from the MC truth in TBrowser
121 ClassDef(AliFlowEventSimple,1)