[u/mrichter/AliRoot.git] / PWGCF / FEMTOSCOPY / AliFemtoUser / AliFemtoQinvCorrFctnEMCIC.h

/***************************************************************************
 *
 * $Id: AliFemtoQinvCorrFctnEMCIC.h  $
 *
 * Author: Nicolas Bock, Ohio State University, bock@mps.ohio-state.edu
 ***************************************************************************
 *
 * Description: Calculates of the Qinv Correlation Function, and also
 *              produces histograms to calculate Energy Momentum Conservation
 *              Induced Correlations  (EMCICs)
 *
 * This Class produces the following histograms as function of Qinv
 * (for both real and mixed pairs):
 *        1)   E1 + E2
 *        2)   E1 * E2
 *        3)   Pt1*Pt2
 *        4)   Pz1*Pz2
 *  
 * The class is derived from AliFemtoQinvCorrFctn, therefore it produces
 * also the histograms in that class. 
 * 
 * NOTE: The EMCIC histograms are not averaged in this class, to obtain 
 * the average, the user needs to divide the real pair histograms by 
 * the numerator, and the mixed pairs by denominator
 *
 ***************************************************************************
 *
 **************************************************************************/

#ifndef ALIFEMTOQINVCORRFCTNEMCIC_H
#define ALIFEMTOQINVCORRFCTNEMCIC_H

#include "TH1D.h"
#include "AliFemtoCorrFctn.h"
#include "AliFemtoQinvCorrFctn.h"


class AliFemtoQinvCorrFctnEMCIC : public AliFemtoQinvCorrFctn 
{
 public:
  AliFemtoQinvCorrFctnEMCIC(char* title, const int& nbins, 
			    const float& QinvLo, const float& QinvHi);
  AliFemtoQinvCorrFctnEMCIC(const AliFemtoQinvCorrFctnEMCIC& aCorrFctn);
  virtual ~AliFemtoQinvCorrFctnEMCIC();
  
  AliFemtoQinvCorrFctnEMCIC& operator=(const AliFemtoQinvCorrFctnEMCIC& aCorrFctn);

  virtual void AddRealPair(AliFemtoPair* aPair);
  virtual void AddMixedPair(AliFemtoPair* aPair);

  
  virtual TList* GetOutputList();
  void Write();
  
 private:
  //Emcic histograms:  
  /*TH1D* fESumReal;   //  <E1+E2>   from real Pairs 
  TH1D* fEMultReal;  //  <E1*E2>   from real Pairs 
  TH1D* fPtMultReal; //  <Pt1*Pt2> from real Pairs 
  TH1D* fPzMultReal; //  <Pz1*Pz2> from real Pairs */
  TH1D* fESumMix;    //  <E1+E2>   from mixed Pairs 
  TH1D* fEMultMix;   //  <E1*E2>   from mixed Pairs
  TH1D* fPtMultMix;  //  <PT1*Pt2> from mixed Pairs
  TH1D* fPzMultMix;  //  <Pz1*Pz2> from mixed Pairs


#ifdef __ROOT__
  ClassDef(AliFemtoQinvCorrFctnEMCIC, 1)
#endif
};


#endif
Commit	Line	Data
76ce4b5b	1	/***************************************************************************
	2	*
	3	* $Id: AliFemtoQinvCorrFctnEMCIC.h $
	4	*
	5	* Author: Nicolas Bock, Ohio State University, bock@mps.ohio-state.edu
	6	***************************************************************************
	7	*
	8	* Description: Calculates of the Qinv Correlation Function, and also
	9	* produces histograms to calculate Energy Momentum Conservation
	10	* Induced Correlations (EMCICs)
	11	*
	12	* This Class produces the following histograms as function of Qinv
	13	* (for both real and mixed pairs):
	14	* 1) E1 + E2
	15	* 2) E1 * E2
	16	* 3) Pt1*Pt2
	17	* 4) Pz1*Pz2
	18	*
	19	* The class is derived from AliFemtoQinvCorrFctn, therefore it produces
	20	* also the histograms in that class.
	21	*
	22	* NOTE: The EMCIC histograms are not averaged in this class, to obtain
	23	* the average, the user needs to divide the real pair histograms by
	24	* the numerator, and the mixed pairs by denominator
	25	*
	26	***************************************************************************
	27	*
	28	**************************************************************************/
	29
	30	#ifndef ALIFEMTOQINVCORRFCTNEMCIC_H
	31	#define ALIFEMTOQINVCORRFCTNEMCIC_H
	32
	33	#include "TH1D.h"
	34	#include "AliFemtoCorrFctn.h"
	35	#include "AliFemtoQinvCorrFctn.h"
	36
	37
	38	class AliFemtoQinvCorrFctnEMCIC : public AliFemtoQinvCorrFctn
	39	{
	40	public:
	41	AliFemtoQinvCorrFctnEMCIC(char* title, const int& nbins,
	42	const float& QinvLo, const float& QinvHi);
	43	AliFemtoQinvCorrFctnEMCIC(const AliFemtoQinvCorrFctnEMCIC& aCorrFctn);
	44	virtual ~AliFemtoQinvCorrFctnEMCIC();
	45
	46	AliFemtoQinvCorrFctnEMCIC& operator=(const AliFemtoQinvCorrFctnEMCIC& aCorrFctn);
	47
	48	virtual void AddRealPair(AliFemtoPair* aPair);
	49	virtual void AddMixedPair(AliFemtoPair* aPair);
	50
	51
	52
	53	virtual TList* GetOutputList();
	54	void Write();
	55
	56	private:
	57	//Emcic histograms:
	58	/TH1D fESumReal; // <E1+E2> from real Pairs
	59	TH1D* fEMultReal; // <E1*E2> from real Pairs
	60	TH1D* fPtMultReal; // <Pt1*Pt2> from real Pairs
	61	TH1D* fPzMultReal; // <Pz1Pz2> from real Pairs /
	62	TH1D* fESumMix; // <E1+E2> from mixed Pairs
	63	TH1D* fEMultMix; // <E1*E2> from mixed Pairs
	64	TH1D* fPtMultMix; // <PT1*Pt2> from mixed Pairs
65	TH1D* fPzMultMix; // <Pz1*Pz2> from mixed Pairs
66
67
68
69
70	#ifdef __ROOT__
71	ClassDef(AliFemtoQinvCorrFctnEMCIC, 1)
72	#endif
73	};
74
75
76
77	#endif
78