[u/mrichter/AliRoot.git] / PWG / FLOW / Base / AliFlowEventSimple.h

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

/*****************************************************************
  AliFlowEventSimple: A simple event 
  for flow analysis                  
                                     
  origin: Naomi van der Kolk (kolk@nikhef.nl)           
          Ante Bilandzic     (anteb@nikhef.nl)         
          Raimond Snellings  (Raimond.Snellings@nikhef.nl)    
  mods:   Mikolaj Krzewicki  (mikolaj.krzewicki@cern.ch)
*****************************************************************/

#ifndef ALIFLOWEVENTSIMPLE_H
#define ALIFLOWEVENTSIMPLE_H

#include "TObject.h"
#include "TParameter.h"
#include "TMath.h"
class TTree;
class TF1;
class AliFlowVector;
class AliFlowTrackSimple;
class AliFlowTrackSimpleCuts;

class AliFlowEventSimple: public TObject {

 public:

  enum ConstructionMethod {kEmpty,kGenerate};

  AliFlowEventSimple();
  AliFlowEventSimple( Int_t nParticles,
                      ConstructionMethod m=kEmpty,
                      TF1* ptDist=NULL,
                      Double_t phiMin=0.0,
                      Double_t phiMax=TMath::TwoPi(),
                      Double_t etaMin=-1.0,
                      Double_t etaMax= 1.0 );
  AliFlowEventSimple(TTree* anInput, const AliFlowTrackSimpleCuts* rpCuts, const AliFlowTrackSimpleCuts* poiCuts);
  AliFlowEventSimple(const AliFlowEventSimple& anEvent);
  AliFlowEventSimple& operator=(const AliFlowEventSimple& anEvent);
  virtual  ~AliFlowEventSimple();

  Bool_t  IsFolder() const {return kTRUE;};
  void    Browse(TBrowser *b); 
  void    Print(Option_t* option = "") const;      //method to print stats
  
  Int_t    NumberOfTracks() const                   { return fNumberOfTracks; }
  Int_t    GetReferenceMultiplicity() const         { return fReferenceMultiplicity; }
  void     SetReferenceMultiplicity( Int_t m )      { fReferenceMultiplicity = m; }
  Int_t    GetEventNSelTracksRP() const             { return fNumberOfRPs; } 
  void     SetEventNSelTracksRP(Int_t nr)           { fNumberOfRPs = nr; }  
  Int_t    GetEventNSelTracksPOI() const            { return fNumberOfPOIs; } 
  void     SetEventNSelTracksPOI(Int_t np)          { fNumberOfPOIs = np; }  
  Int_t    GetNumberOfRPs() const                   { return fNumberOfRPs; }
  void     SetNumberOfRPs( Int_t nr )               { fNumberOfRPs=nr; }
  Int_t    GetNumberOfPOIs() const                  { return fNumberOfPOIs; }
  void     SetNumberOfPOIs( Int_t np )              { fNumberOfPOIs=np; }
  Double_t GetMCReactionPlaneAngle() const          { return fMCReactionPlaneAngle; }
  void     SetMCReactionPlaneAngle(Double_t fPhiRP) { fMCReactionPlaneAngle=fPhiRP; fMCReactionPlaneAngleIsSet=kTRUE; }
  Bool_t   IsSetMCReactionPlaneAngle() const        { return fMCReactionPlaneAngleIsSet; }
  void     SetAfterBurnerPrecision(Double_t p)      { fAfterBurnerPrecision=p; }
  Double_t GetAfterBurnerPrecision() const          { return fAfterBurnerPrecision; }
  void     SetUserModified(Bool_t s=kTRUE)          { fUserModified=s; }
  Bool_t   IsUserModified() const                   { return fUserModified; }
  void     SetShuffleTracks(Bool_t b)               {fShuffleTracks=b;}
  void     ShuffleTracks();

  void ResolutionPt(Double_t res);
  void TagSubeventsInEta(Double_t etaMinA, Double_t etaMaxA, Double_t etaMinB, Double_t etaMaxB );
  void TagSubeventsByCharge();
  void TagRP(const AliFlowTrackSimpleCuts* cuts );
  void TagPOI(const AliFlowTrackSimpleCuts* cuts );
  void CloneTracks(Int_t n);
  void AddV1( Double_t v1 );
  void AddV2( Double_t v2 );
  void AddV3( Double_t v3 );
  void AddV4( Double_t v4 );
  void AddV5( Double_t v5 );
  void AddFlow( Double_t v1, Double_t v2, Double_t v3, Double_t v4, Double_t v5 );
  void AddFlow(Double_t v1, Double_t v2, Double_t v3, Double_t v4, Double_t v5,
               Double_t rp1, Double_t rp2, Double_t rp3, Double_t rp4, Double_t rp5 );
  void AddV2( TF1* ptDepV2 );
  void DefineDeadZone( Double_t etaMin, Double_t etaMax, Double_t phiMin, Double_t phiMax );
  Int_t CleanUpDeadTracks();
  void ClearFast();
 
  static TF1* SimplePtSpectrum();
  static TF1* SimplePtDepV2();

  AliFlowTrackSimple* GetTrack(Int_t i);
  void AddTrack( AliFlowTrackSimple* track ); 
  AliFlowTrackSimple* MakeNewTrack();
 
  AliFlowVector GetQ(Int_t n=2, TList *weightsList=NULL, Bool_t usePhiWeights=kFALSE, Bool_t usePtWeights=kFALSE, Bool_t useEtaWeights=kFALSE);
  void Get2Qsub(AliFlowVector* Qarray, Int_t n=2, TList *weightsList=NULL, Bool_t usePhiWeights=kFALSE, Bool_t usePtWeights=kFALSE, Bool_t useEtaWeights=kFALSE);  

 protected:
  virtual void Generate( Int_t nParticles,
                         TF1* ptDist=NULL,
                         Double_t phiMin=0.0,
                         Double_t phiMax=TMath::TwoPi(),
                         Double_t etaMin=-1.0,
                         Double_t etaMax= 1.0 );

  //data members
  TObjArray*              fTrackCollection;           //-> collection of tracks
  Int_t                   fReferenceMultiplicity;     // reference multiplicity
  Int_t                   fNumberOfTracks;            // number of tracks
  Int_t                   fNumberOfRPs;               // number of tracks that have passed the RP selection
  Int_t                   fNumberOfPOIs;              // number of tracks that have passed the POI selection
  Double_t                fMCReactionPlaneAngle;      // the angle of the reaction plane from the MC truth
  Bool_t                  fMCReactionPlaneAngleIsSet; // did we set it from MC?
  Double_t                fAfterBurnerPrecision;      // iteration precision in afterburner
  Bool_t                  fUserModified;              // did we modify the event in any way (afterburner etc) ?
  TParameter<Int_t>*      fNumberOfTracksWrap;        //! number of tracks in TBrowser
  TParameter<Int_t>*      fNumberOfRPsWrap;           //! number of tracks that have passed the RP selection in TBrowser
  TParameter<Int_t>*      fNumberOfPOIsWrap;          //! number of tracks that have passed the POI selection in TBrowser
  TParameter<Double_t>*   fMCReactionPlaneAngleWrap;  //! the angle of the reaction plane from the MC truth in TBrowser
  Int_t*                  fShuffledIndexes;           //! placeholder for randomized indexes
  Bool_t                  fShuffleTracks;             // do we shuffle tracks on get?

  ClassDef(AliFlowEventSimple,1)
};

#endif
Commit	Line	Data
f1d945a1	1	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	2	* See cxx source for full Copyright notice */
	3	/* $Id$ */
	4
929098e4	5	/*****************************************************************
	6	AliFlowEventSimple: A simple event
	7	for flow analysis
	8
	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	*****************************************************************/
f1d945a1	14
929098e4	15	#ifndef ALIFLOWEVENTSIMPLE_H
	16	#define ALIFLOWEVENTSIMPLE_H
	17
	18	#include "TObject.h"
c076fda8	19	#include "TParameter.h"
bc231a12	20	#include "TMath.h"
929098e4	21	class TTree;
bc231a12	22	class TF1;
929098e4	23	class AliFlowVector;
f1d945a1	24	class AliFlowTrackSimple;
929098e4	25	class AliFlowTrackSimpleCuts;
f1d945a1	26
	27	class AliFlowEventSimple: public TObject {
	28
	29	public:
9fa64edc	30
	31	enum ConstructionMethod {kEmpty,kGenerate};
	32
46bec39c	33	AliFlowEventSimple();
bc231a12	34	AliFlowEventSimple( Int_t nParticles,
9fa64edc	35	ConstructionMethod m=kEmpty,
	36	TF1* ptDist=NULL,
	37	Double_t phiMin=0.0,
bc231a12	38	Double_t phiMax=TMath::TwoPi(),
	39	Double_t etaMin=-1.0,
	40	Double_t etaMax= 1.0 );
929098e4	41	AliFlowEventSimple(TTree* anInput, const AliFlowTrackSimpleCuts* rpCuts, const AliFlowTrackSimpleCuts* poiCuts);
e35ddff0	42	AliFlowEventSimple(const AliFlowEventSimple& anEvent);
e35ddff0	43	AliFlowEventSimple& operator=(const AliFlowEventSimple& anEvent);
f1d945a1	44	virtual ~AliFlowEventSimple();
c076fda8	45
c076fda8	46	Bool_t IsFolder() const {return kTRUE;};
29195b69	47	void Browse(TBrowser *b);
29195b69	48	void Print(Option_t* option = "") const; //method to print stats
f1d945a1	49
929098e4	50	Int_t NumberOfTracks() const { return fNumberOfTracks; }
df23c5ae	51	Int_t GetReferenceMultiplicity() const { return fReferenceMultiplicity; }
df23c5ae	52	void SetReferenceMultiplicity( Int_t m ) { fReferenceMultiplicity = m; }
bc231a12	53	Int_t GetEventNSelTracksRP() const { return fNumberOfRPs; }
bc231a12	54	void SetEventNSelTracksRP(Int_t nr) { fNumberOfRPs = nr; }
df23c5ae	55	Int_t GetEventNSelTracksPOI() const { return fNumberOfPOIs; }
df23c5ae	56	void SetEventNSelTracksPOI(Int_t np) { fNumberOfPOIs = np; }
59dab33a	57	Int_t GetNumberOfRPs() const { return fNumberOfRPs; }
59dab33a	58	void SetNumberOfRPs( Int_t nr ) { fNumberOfRPs=nr; }
df23c5ae	59	Int_t GetNumberOfPOIs() const { return fNumberOfPOIs; }
df23c5ae	60	void SetNumberOfPOIs( Int_t np ) { fNumberOfPOIs=np; }
929098e4	61	Double_t GetMCReactionPlaneAngle() const { return fMCReactionPlaneAngle; }
929098e4	62	void SetMCReactionPlaneAngle(Double_t fPhiRP) { fMCReactionPlaneAngle=fPhiRP; fMCReactionPlaneAngleIsSet=kTRUE; }
124fb262	63	Bool_t IsSetMCReactionPlaneAngle() const { return fMCReactionPlaneAngleIsSet; }
bc231a12	64	void SetAfterBurnerPrecision(Double_t p) { fAfterBurnerPrecision=p; }
bc231a12	65	Double_t GetAfterBurnerPrecision() const { return fAfterBurnerPrecision; }
9fa64edc	66	void SetUserModified(Bool_t s=kTRUE) { fUserModified=s; }
9fa64edc	67	Bool_t IsUserModified() const { return fUserModified; }
ad4d7e4e	68	void SetShuffleTracks(Bool_t b) {fShuffleTracks=b;}
ba4dd29f	69	void ShuffleTracks();
929098e4	70
7382279b	71	void ResolutionPt(Double_t res);
940a5ed1	72	void TagSubeventsInEta(Double_t etaMinA, Double_t etaMaxA, Double_t etaMinB, Double_t etaMaxB );
076df7bf	73	void TagSubeventsByCharge();
cc0afcfc	74	void TagRP(const AliFlowTrackSimpleCuts* cuts );
cc0afcfc	75	void TagPOI(const AliFlowTrackSimpleCuts* cuts );
929098e4	76	void CloneTracks(Int_t n);
bc231a12	77	void AddV1( Double_t v1 );
244c607a	78	void AddV2( Double_t v2 );
54089829	79	void AddV3( Double_t v3 );
bc231a12	80	void AddV4( Double_t v4 );
1a80f9f6	81	void AddV5( Double_t v5 );
	82	void AddFlow( Double_t v1, Double_t v2, Double_t v3, Double_t v4, Double_t v5 );
	83	void AddFlow(Double_t v1, Double_t v2, Double_t v3, Double_t v4, Double_t v5,
	84	Double_t rp1, Double_t rp2, Double_t rp3, Double_t rp4, Double_t rp5 );
99ff691b	85	void AddV2( TF1* ptDepV2 );
34b15925	86	void DefineDeadZone( Double_t etaMin, Double_t etaMax, Double_t phiMin, Double_t phiMax );
34b15925	87	Int_t CleanUpDeadTracks();
7d27a354	88	void ClearFast();
26c4cbb9	89
99ff691b	90	static TF1* SimplePtSpectrum();
	91	static TF1* SimplePtDepV2();
	92
f1d945a1	93	AliFlowTrackSimple* GetTrack(Int_t i);
7382279b	94	void AddTrack( AliFlowTrackSimple* track );
7b5556ef	95	AliFlowTrackSimple* MakeNewTrack();
03a02aca	96
03a02aca	97	AliFlowVector GetQ(Int_t n=2, TList *weightsList=NULL, Bool_t usePhiWeights=kFALSE, Bool_t usePtWeights=kFALSE, Bool_t useEtaWeights=kFALSE);
29195b69	98	void Get2Qsub(AliFlowVector* Qarray, Int_t n=2, TList *weightsList=NULL, Bool_t usePhiWeights=kFALSE, Bool_t usePtWeights=kFALSE, Bool_t useEtaWeights=kFALSE);
395fadba	99
7382279b	100	protected:
9fa64edc	101	virtual void Generate( Int_t nParticles,
	102	TF1* ptDist=NULL,
	103	Double_t phiMin=0.0,
	104	Double_t phiMax=TMath::TwoPi(),
	105	Double_t etaMin=-1.0,
	106	Double_t etaMax= 1.0 );
	107
7d27a354	108	//data members
bc231a12	109	TObjArray* fTrackCollection; //-> collection of tracks
7d27a354	110	Int_t fReferenceMultiplicity; // reference multiplicity
929098e4	111	Int_t fNumberOfTracks; // number of tracks
bc231a12	112	Int_t fNumberOfRPs; // number of tracks that have passed the RP selection
df23c5ae	113	Int_t fNumberOfPOIs; // number of tracks that have passed the POI selection
929098e4	114	Double_t fMCReactionPlaneAngle; // the angle of the reaction plane from the MC truth
929098e4	115	Bool_t fMCReactionPlaneAngleIsSet; // did we set it from MC?
bc231a12	116	Double_t fAfterBurnerPrecision; // iteration precision in afterburner
9fa64edc	117	Bool_t fUserModified; // did we modify the event in any way (afterburner etc) ?
929098e4	118	TParameter<Int_t>* fNumberOfTracksWrap; //! number of tracks in TBrowser
bc231a12	119	TParameter<Int_t>* fNumberOfRPsWrap; //! number of tracks that have passed the RP selection in TBrowser
df23c5ae	120	TParameter<Int_t>* fNumberOfPOIsWrap; //! number of tracks that have passed the POI selection in TBrowser
bc231a12	121	TParameter<Double_t>* fMCReactionPlaneAngleWrap; //! the angle of the reaction plane from the MC truth in TBrowser
ad4d7e4e	122	Int_t* fShuffledIndexes; //! placeholder for randomized indexes
ad4d7e4e	123	Bool_t fShuffleTracks; // do we shuffle tracks on get?
29195b69	124
bc231a12	125	ClassDef(AliFlowEventSimple,1)
f1d945a1	126	};
	127
	128	#endif
5fef318d	129
5fef318d	130