[u/mrichter/AliRoot.git] / PWG2 / FEMTOSCOPY / AliFemto / AliFemtoCutMonitorEventMult.h

////////////////////////////////////////////////////////////////////////////////
///                                                                          ///
/// AliFemtoCutMonitorEventMult - the cut monitor for events to study        ///
/// the multiplicity distribution of events                                  ///
///                                                                          ///
////////////////////////////////////////////////////////////////////////////////
#ifndef AliFemtoCutMonitorEventMult_hh
#define AliFemtoCutMonitorEventMult_hh

class AliFemtoEvent;
class AliFemtoTrack;
class AliFemtoV0;
class AliFemtoKink;
class AliFemtoPair; 
class TH1D;
class TH2D;
class TList;
#include "AliFemtoString.h"
#include "AliFemtoParticleCollection.h"
#include "AliFemtoCutMonitor.h"

class AliFemtoCutMonitorEventMult : public AliFemtoCutMonitor{
  
 public:
  AliFemtoCutMonitorEventMult();
  AliFemtoCutMonitorEventMult(const char *aName);
  AliFemtoCutMonitorEventMult(const AliFemtoCutMonitorEventMult &aCut);
  virtual ~AliFemtoCutMonitorEventMult();

  AliFemtoCutMonitorEventMult& operator=(const AliFemtoCutMonitorEventMult& aCut);

  virtual AliFemtoString Report();
  virtual void Fill(const AliFemtoEvent* aEvent);
  virtual void Fill(const AliFemtoTrack* aTrack) {AliFemtoCutMonitor::Fill(aTrack);}
  virtual void Fill(const AliFemtoV0* aV0) {AliFemtoCutMonitor::Fill(aV0);}
  virtual void Fill(const AliFemtoKink* aKink) {AliFemtoCutMonitor::Fill(aKink);}
  virtual void Fill(const AliFemtoPair* aPair) {AliFemtoCutMonitor::Fill(aPair);}
  virtual void Fill(const AliFemtoParticleCollection* aCollection) {AliFemtoCutMonitor::Fill(aCollection);}
  virtual void Fill(const AliFemtoEvent* aEvent,const AliFemtoParticleCollection* aCollection)
  {AliFemtoCutMonitor::Fill(aEvent, aCollection);}

  void SetReadMC(Bool_t mc);
  void AdditionalMultHistsOn(Bool_t addhists);
  void Write();

  virtual TList *GetOutputList();

 private:
  TH1D *fEvMult;     // Multiplicity distribution
  TH1D *fNormEvMult; // Normalized event multiplicity distribution
  TH1D *fSPDMult;    // SPD tracklet multiplicity
  TH2D *fMultSumPt;  // Event total pT vs. multiplicity

  Bool_t freadMC;     // If true - add only one histogram to the output
  Bool_t faddhists;   // If true - add only additional multiplicity histograms
  
  TH1D *fEstimateITSTPC;     // Multiplicity estimate ITS+TPC
  TH1D *fEstimateTracklets;  // Multiplicity estimate Tracklets
  TH1D *fEstimateITSPure;    // Multiplicity estimate ITS Pure

  TH2D *fEst1Est2;           // ITS+TPC vs Tracklets
  TH2D *fEst1Est3;           // ITS+TPC vs ITS Pure
  TH2D *fEst2Est3;           // Tracklets vs ITS Pure
  TH2D *fEst1Norm;           // ITS+TPC vs Normalized
  TH2D *fEst2Norm;           // Tracklets vs Normalized
  TH2D *fEst3Norm;           // ITS Pure vs Normalized

};

#endif
Commit	Line	Data
00114ecb	1	////////////////////////////////////////////////////////////////////////////////
	2	/// ///
	3	/// AliFemtoCutMonitorEventMult - the cut monitor for events to study ///
	4	/// the multiplicity distribution of events ///
	5	/// ///
	6	////////////////////////////////////////////////////////////////////////////////
	7	#ifndef AliFemtoCutMonitorEventMult_hh
	8	#define AliFemtoCutMonitorEventMult_hh
	9
	10	class AliFemtoEvent;
	11	class AliFemtoTrack;
	12	class AliFemtoV0;
	13	class AliFemtoKink;
	14	class AliFemtoPair;
	15	class TH1D;
	16	class TH2D;
	17	class TList;
	18	#include "AliFemtoString.h"
	19	#include "AliFemtoParticleCollection.h"
	20	#include "AliFemtoCutMonitor.h"
	21
	22	class AliFemtoCutMonitorEventMult : public AliFemtoCutMonitor{
	23
73695088	24	public:
00114ecb	25	AliFemtoCutMonitorEventMult();
	26	AliFemtoCutMonitorEventMult(const char *aName);
	27	AliFemtoCutMonitorEventMult(const AliFemtoCutMonitorEventMult &aCut);
	28	virtual ~AliFemtoCutMonitorEventMult();
	29
	30	AliFemtoCutMonitorEventMult& operator=(const AliFemtoCutMonitorEventMult& aCut);
	31
	32	virtual AliFemtoString Report();
	33	virtual void Fill(const AliFemtoEvent* aEvent);
	34	virtual void Fill(const AliFemtoTrack* aTrack) {AliFemtoCutMonitor::Fill(aTrack);}
	35	virtual void Fill(const AliFemtoV0* aV0) {AliFemtoCutMonitor::Fill(aV0);}
	36	virtual void Fill(const AliFemtoKink* aKink) {AliFemtoCutMonitor::Fill(aKink);}
	37	virtual void Fill(const AliFemtoPair* aPair) {AliFemtoCutMonitor::Fill(aPair);}
	38	virtual void Fill(const AliFemtoParticleCollection* aCollection) {AliFemtoCutMonitor::Fill(aCollection);}
	39	virtual void Fill(const AliFemtoEvent* aEvent,const AliFemtoParticleCollection* aCollection)
	40	{AliFemtoCutMonitor::Fill(aEvent, aCollection);}
	41
859f1d8d	42	void SetReadMC(Bool_t mc);
859f1d8d	43	void AdditionalMultHistsOn(Bool_t addhists);
00114ecb	44	void Write();
	45
	46	virtual TList *GetOutputList();
	47
73695088	48	private:
	49	TH1D *fEvMult; // Multiplicity distribution
	50	TH1D *fNormEvMult; // Normalized event multiplicity distribution
e2b7a66d	51	TH1D *fSPDMult; // SPD tracklet multiplicity
782723e1	52	TH2D *fMultSumPt; // Event total pT vs. multiplicity
859f1d8d	53
	54	Bool_t freadMC; // If true - add only one histogram to the output
	55	Bool_t faddhists; // If true - add only additional multiplicity histograms
73c932f8	56
	57	TH1D *fEstimateITSTPC; // Multiplicity estimate ITS+TPC
	58	TH1D *fEstimateTracklets; // Multiplicity estimate Tracklets
	59	TH1D *fEstimateITSPure; // Multiplicity estimate ITS Pure
	60
	61	TH2D *fEst1Est2; // ITS+TPC vs Tracklets
	62	TH2D *fEst1Est3; // ITS+TPC vs ITS Pure
	63	TH2D *fEst2Est3; // Tracklets vs ITS Pure
	64	TH2D *fEst1Norm; // ITS+TPC vs Normalized
	65	TH2D *fEst2Norm; // Tracklets vs Normalized
	66	TH2D *fEst3Norm; // ITS Pure vs Normalized
	67
00114ecb	68	};
	69
	70	#endif