[u/mrichter/AliRoot.git] / PWGCF / 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, int nBins=5000);
  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);}
  virtual void Fill(const AliFemtoParticleCollection* aCollection1,const AliFemtoParticleCollection* aCollection2) {AliFemtoCutMonitor::Fill(aCollection1, aCollection2);}

  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

  TH1D *fPsiVZERO;    // psi from vzero

};

#endif
Commit	Line	Data
76ce4b5b	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
	24	public:
	25	AliFemtoCutMonitorEventMult();
4eac0b05	26	AliFemtoCutMonitorEventMult(const char *aName, int nBins=5000);
76ce4b5b	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);}
973a91f8	39	virtual void Fill(const AliFemtoEvent* aEvent,const AliFemtoParticleCollection* aCollection) {AliFemtoCutMonitor::Fill(aEvent, aCollection);}
973a91f8	40	virtual void Fill(const AliFemtoParticleCollection* aCollection1,const AliFemtoParticleCollection* aCollection2) {AliFemtoCutMonitor::Fill(aCollection1, aCollection2);}
76ce4b5b	41
	42	void SetReadMC(Bool_t mc);
	43	void AdditionalMultHistsOn(Bool_t addhists);
	44	void Write();
	45
	46	virtual TList *GetOutputList();
	47
	48	private:
	49	TH1D *fEvMult; // Multiplicity distribution
	50	TH1D *fNormEvMult; // Normalized event multiplicity distribution
	51	TH1D *fSPDMult; // SPD tracklet multiplicity
	52	TH2D *fMultSumPt; // Event total pT vs. multiplicity
	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
	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
5e2038a9	68	TH1D *fPsiVZERO; // psi from vzero
5e2038a9	69
76ce4b5b	70	};
	71
	72	#endif