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

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

/********************************** 
 * functions and equations needed * 
 * for calculation of Q-cumulants *
 * and final flow estimates       *
 *                                *   
 * author:  Ante Bilandzic        * 
 *           (anteb@nikhef.nl)    *
 *********************************/ 

#ifndef ALIQCUMULANTSFUNCTIONS_H
#define ALIQCUMULANTSFUNCTIONS_H

#include "AliFlowCommonConstants.h"

class TH1;
class TProfile;
class TProfile2D;
class TProfile3D;

class TObjArray;
class TList;
class TFile;

class AliFlowCommonHistResults;

//================================================================================================================

class AliQCumulantsFunctions{
 public:
  AliQCumulantsFunctions();
  virtual ~AliQCumulantsFunctions();
  AliQCumulantsFunctions(TH1D *intRes, TH1D *diffRes2nd, TH1D *diffRes4th, TH1D *covar, TProfile *AvMult, TProfile *QVector, TProfile *QCorr, TProfile *QProd, TProfile *Direct, TProfile *binPt2p1n1nRP, TProfile *binPt2p2n2nRP, TProfile *binPt3p2n1n1nRP, TProfile *binPt3p1n1n2nRP, TProfile *binPt4p1n1n1n1nRP, TProfile *binEta2p1n1nRP, TProfile *binEta2p2n2nRP, TProfile *binEta3p2n1n1nRP, TProfile *binEta3p1n1n2nRP, TProfile *binEta4p1n1n1n1nRP, TProfile *binPt2p1n1nPOI, TProfile *binPt2p2n2nPOI, TProfile *binPt3p2n1n1nPOI, TProfile *binPt3p1n1n2nPOI, TProfile *binPt4p1n1n1n1nPOI, TProfile *binEta2p1n1nPOI, TProfile *binEta2p2n2nPOI, TProfile *binEta3p2n1n1nPOI, TProfile *binEta3p1n1n2nPOI, TProfile *binEta4p1n1n1n1nPOI, AliFlowCommonHist *ch2nd, AliFlowCommonHist *ch4th, AliFlowCommonHist *ch6th, AliFlowCommonHist *ch8th, AliFlowCommonHistResults *chr2nd, AliFlowCommonHistResults *chr4th, AliFlowCommonHistResults *chr6th, AliFlowCommonHistResults *chr8th);
 
  void Calculate();

 private:
  AliQCumulantsFunctions(const AliQCumulantsFunctions& Qfun);
  AliQCumulantsFunctions& operator=(const AliQCumulantsFunctions& Qfun);
  
  TH1D     *fIntRes;          //results for integrated flow
  TH1D     *fDiffRes2nd;      //results for differential flow (2nd order)
  TH1D     *fDiffRes4th;      //results for differential flow (4th order)
  TH1D     *fCovar;           //results for covariances (1st bin: <2*4>-<2>*<4>, 2nd bin: <2*6>-<2>*<6>, ...)  
          
  TProfile *fAvMult;          //avarage selected multiplicity for int. flow
  TProfile *fQVector;         //avarage values of Q-vector components
  TProfile *fQCorr;           //multi-particle correlations calculated from Q-vectors 
  TProfile *fQProd;           //average of products: 1st bin: <2*4>, 2nd bin: <2*6>, ...

  TProfile *fDirect;          //direct correlations (correlations calculated with nested loopps)
  
  TProfile *fbinPt2p1n1nRP;     //<<2'>>_{n|n} per pt-bin
  TProfile *fbinPt2p2n2nRP;     //<<2'>>_{2n|2n} per pt-bin
  TProfile *fbinPt3p2n1n1nRP;   //<<3'>>_{2n,n|n} per pt-bin
  TProfile *fbinPt3p1n1n2nRP;   //<<3'>>_{n,n|2n} per pt-bin
  TProfile *fbinPt4p1n1n1n1nRP; //<<4'>>_{n,n|n,n} per pt-bin
  
  TProfile *fbinEta2p1n1nRP;     //<<2'>>_{n|n} per eta-bin
  TProfile *fbinEta2p2n2nRP;     //<<2'>>_{2n|2n} per eta-bin
  TProfile *fbinEta3p2n1n1nRP;   //<<3'>>_{2n,n|n} per eta-bin
  TProfile *fbinEta3p1n1n2nRP;   //<<3'>>_{n,n|2n} per eta-bin
  TProfile *fbinEta4p1n1n1n1nRP; //<<4'>>_{n,n|n,n} per eta-bin  
  
  TProfile *fbinPt2p1n1nPOI;     //<<2'>>_{n|n} per pt-bin
  TProfile *fbinPt2p2n2nPOI;     //<<2'>>_{2n|2n} per pt-bin
  TProfile *fbinPt3p2n1n1nPOI;   //<<3'>>_{2n,n|n} per pt-bin
  TProfile *fbinPt3p1n1n2nPOI;   //<<3'>>_{n,n|2n} per pt-bin
  TProfile *fbinPt4p1n1n1n1nPOI; //<<4'>>_{n,n|n,n} per pt-bin
  
  TProfile *fbinEta2p1n1nPOI;     //<<2'>>_{n|n} per eta-bin
  TProfile *fbinEta2p2n2nPOI;     //<<2'>>_{2n|2n} per eta-bin
  TProfile *fbinEta3p2n1n1nPOI;   //<<3'>>_{2n,n|n} per eta-bin
  TProfile *fbinEta3p1n1n2nPOI;   //<<3'>>_{n,n|2n} per eta-bin
  TProfile *fbinEta4p1n1n1n1nPOI; //<<4'>>_{n,n|n,n} per eta-bin  
  
  AliFlowCommonHist *fch2nd; //common control histograms (taking into account only the events with 2 and more particles) 
  AliFlowCommonHist *fch4th; //common control histograms (taking into account only the events with 4 and more particles) 
  AliFlowCommonHist *fch6th; //common control histograms (taking into account only the events with 6 and more particles) 
  AliFlowCommonHist *fch8th; //common control histograms (taking into account only the events with 8 and more particles) 
  
  AliFlowCommonHistResults *fchr2nd; //final results for 2nd order int. and diff. flow stored in the common histograms
  AliFlowCommonHistResults *fchr4th; //final results for 4th order int. and diff. flow stored in the common histograms
  AliFlowCommonHistResults *fchr6th; //final results for 6th order int. and diff. flow stored in the common histograms
  AliFlowCommonHistResults *fchr8th; //final results for 8th order int. and diff. flow stored in the common histograms
              
  ClassDef(AliQCumulantsFunctions, 0);
};

//================================================================================================================

#endif
Commit	Line	Data
bc92c0cb	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	* functions and equations needed *
	9	* for calculation of Q-cumulants *
	10	* and final flow estimates *
	11	* *
	12	* author: Ante Bilandzic *
	13	* (anteb@nikhef.nl) *
	14	*********************************/
	15
52021ae2	16	#ifndef ALIQCUMULANTSFUNCTIONS_H
52021ae2	17	#define ALIQCUMULANTSFUNCTIONS_H
bc92c0cb	18
	19	#include "AliFlowCommonConstants.h"
	20
	21	class TH1;
	22	class TProfile;
	23	class TProfile2D;
	24	class TProfile3D;
	25
	26	class TObjArray;
	27	class TList;
	28	class TFile;
	29
1315fe58	30	class AliFlowCommonHistResults;
1315fe58	31
bc92c0cb	32	//================================================================================================================
	33
	34	class AliQCumulantsFunctions{
	35	public:
	36	AliQCumulantsFunctions();
	37	virtual ~AliQCumulantsFunctions();
7e58a232	38	AliQCumulantsFunctions(TH1D intRes, TH1D diffRes2nd, TH1D diffRes4th, TH1D covar, TProfile AvMult, TProfile QVector, TProfile QCorr, TProfile QProd, TProfile Direct, TProfile binPt2p1n1nRP, TProfile binPt2p2n2nRP, TProfile binPt3p2n1n1nRP, TProfile binPt3p1n1n2nRP, TProfile binPt4p1n1n1n1nRP, TProfile binEta2p1n1nRP, TProfile binEta2p2n2nRP, TProfile binEta3p2n1n1nRP, TProfile binEta3p1n1n2nRP, TProfile binEta4p1n1n1n1nRP, TProfile binPt2p1n1nPOI, TProfile binPt2p2n2nPOI, TProfile binPt3p2n1n1nPOI, TProfile binPt3p1n1n2nPOI, TProfile binPt4p1n1n1n1nPOI, TProfile binEta2p1n1nPOI, TProfile binEta2p2n2nPOI, TProfile binEta3p2n1n1nPOI, TProfile binEta3p1n1n2nPOI, TProfile binEta4p1n1n1n1nPOI, AliFlowCommonHist ch2nd, AliFlowCommonHist ch4th, AliFlowCommonHist ch6th, AliFlowCommonHist ch8th, AliFlowCommonHistResults chr2nd, AliFlowCommonHistResults chr4th, AliFlowCommonHistResults chr6th, AliFlowCommonHistResults *chr8th);
bc92c0cb	39
	40	void Calculate();
	41
	42	private:
	43	AliQCumulantsFunctions(const AliQCumulantsFunctions& Qfun);
	44	AliQCumulantsFunctions& operator=(const AliQCumulantsFunctions& Qfun);
	45
1dfa3c16	46	TH1D *fIntRes; //results for integrated flow
	47	TH1D *fDiffRes2nd; //results for differential flow (2nd order)
	48	TH1D *fDiffRes4th; //results for differential flow (4th order)
	49	TH1D fCovar; //results for covariances (1st bin: <24>-<2><4>, 2nd bin: <26>-<2>*<6>, ...)
bc92c0cb	50
1dfa3c16	51	TProfile *fAvMult; //avarage selected multiplicity for int. flow
	52	TProfile *fQVector; //avarage values of Q-vector components
	53	TProfile *fQCorr; //multi-particle correlations calculated from Q-vectors
	54	TProfile fQProd; //average of products: 1st bin: <24>, 2nd bin: <2*6>, ...
bc92c0cb	55
1dfa3c16	56	TProfile *fDirect; //direct correlations (correlations calculated with nested loopps)
bc92c0cb	57
1dfa3c16	58	TProfile *fbinPt2p1n1nRP; //<<2'>>_{n\|n} per pt-bin
	59	TProfile *fbinPt2p2n2nRP; //<<2'>>_{2n\|2n} per pt-bin
	60	TProfile *fbinPt3p2n1n1nRP; //<<3'>>_{2n,n\|n} per pt-bin
	61	TProfile *fbinPt3p1n1n2nRP; //<<3'>>_{n,n\|2n} per pt-bin
	62	TProfile *fbinPt4p1n1n1n1nRP; //<<4'>>_{n,n\|n,n} per pt-bin
	63
	64	TProfile *fbinEta2p1n1nRP; //<<2'>>_{n\|n} per eta-bin
	65	TProfile *fbinEta2p2n2nRP; //<<2'>>_{2n\|2n} per eta-bin
	66	TProfile *fbinEta3p2n1n1nRP; //<<3'>>_{2n,n\|n} per eta-bin
	67	TProfile *fbinEta3p1n1n2nRP; //<<3'>>_{n,n\|2n} per eta-bin
	68	TProfile *fbinEta4p1n1n1n1nRP; //<<4'>>_{n,n\|n,n} per eta-bin
	69
	70	TProfile *fbinPt2p1n1nPOI; //<<2'>>_{n\|n} per pt-bin
	71	TProfile *fbinPt2p2n2nPOI; //<<2'>>_{2n\|2n} per pt-bin
	72	TProfile *fbinPt3p2n1n1nPOI; //<<3'>>_{2n,n\|n} per pt-bin
	73	TProfile *fbinPt3p1n1n2nPOI; //<<3'>>_{n,n\|2n} per pt-bin
	74	TProfile *fbinPt4p1n1n1n1nPOI; //<<4'>>_{n,n\|n,n} per pt-bin
	75
	76	TProfile *fbinEta2p1n1nPOI; //<<2'>>_{n\|n} per eta-bin
	77	TProfile *fbinEta2p2n2nPOI; //<<2'>>_{2n\|2n} per eta-bin
	78	TProfile *fbinEta3p2n1n1nPOI; //<<3'>>_{2n,n\|n} per eta-bin
	79	TProfile *fbinEta3p1n1n2nPOI; //<<3'>>_{n,n\|2n} per eta-bin
	80	TProfile *fbinEta4p1n1n1n1nPOI; //<<4'>>_{n,n\|n,n} per eta-bin
1315fe58	81
7e58a232	82	AliFlowCommonHist *fch2nd; //common control histograms (taking into account only the events with 2 and more particles)
	83	AliFlowCommonHist *fch4th; //common control histograms (taking into account only the events with 4 and more particles)
	84	AliFlowCommonHist *fch6th; //common control histograms (taking into account only the events with 6 and more particles)
	85	AliFlowCommonHist *fch8th; //common control histograms (taking into account only the events with 8 and more particles)
	86
1315fe58	87	AliFlowCommonHistResults *fchr2nd; //final results for 2nd order int. and diff. flow stored in the common histograms
	88	AliFlowCommonHistResults *fchr4th; //final results for 4th order int. and diff. flow stored in the common histograms
	89	AliFlowCommonHistResults *fchr6th; //final results for 6th order int. and diff. flow stored in the common histograms
	90	AliFlowCommonHistResults *fchr8th; //final results for 8th order int. and diff. flow stored in the common histograms
1315fe58	91
bc92c0cb	92	ClassDef(AliQCumulantsFunctions, 0);
	93	};
	94
	95	//================================================================================================================
	96
	97	#endif
	98
	99
	100
	101
	102