[u/mrichter/AliRoot.git] / PWG1 / AliComparisonDraw.h


#ifndef AliComparisonDraw_h
#define AliComparisonDraw_h

#include <iostream>
#include <fstream>
using namespace std;
#include <TSelector.h>


class TFile;
class AliMCInfo;
class AliESDRecInfo;
class AliESDEvent; 
class AliESD;
class AliESDfriend;
class TH1I;
class TH3F;
class TH3;
class TProfile;
class TProfile2D;
class TGraph2D;
class TGraph; 

class AliComparisonDraw : public TObject {
public :
  AliComparisonDraw(); 
  virtual Bool_t    IsFolder(){return kTRUE;}
  void      InitHisto();
  void      Process(AliMCInfo* infoMC, AliESDRecInfo *infoRC);
  
  //
  //
  void      ProcessEff(AliMCInfo* infoMC, AliESDRecInfo *infoRC);
  void      ProcessResolConstrained(AliMCInfo* infoMC, AliESDRecInfo *infoRC);
  void      ProcessTPCdedx(AliMCInfo* infoMC, AliESDRecInfo *infoRC);
  void      ProcessDCA(AliMCInfo* infoMC, AliESDRecInfo *infoRC);

  void MakePlots();


  //TH1F            GetPtResol(Float_t pt0, Float_t pt1);
  static TH1F*       MakeResol(TH2F * his, Int_t integ, Bool_t type); 
  static TGraph2D *  MakeStat2D(TH3 * his, Int_t delta0, Int_t delta1, Int_t type);
  static TGraph *  MakeStat1D(TH3 * his, Int_t delta1, Int_t type);


protected:
  //
  // efficiency 
  //
  static    Bool_t    fgBDraw;         //option draw temporary results
  TProfile* fEffTPCPt;      // TPC efficiency as function of Pt (tan+-1)
  TProfile* fEffTPCPtMC;    // MC -TPC efficiency as function of Pt (tan+-1)
  TProfile* fEffTPCPtF;     // efficiency for findable tracks
  //
  TProfile* fEffTPCTan;   // TPC efficiency as function of Tan (pt>0.15
  TProfile* fEffTPCTanMC; // MC -TPC efficiency as function of Tan (pt>0.15)
  TProfile* fEffTPCTanF;  // efficiency for findable tracks Tan (pt>0.15)
  //
  TProfile2D* fEffTPCPtTan;    // TPC efficiency as function of Pt and tan
  TProfile2D* fEffTPCPtTanMC;  // MC -TPC efficiency as function of Pt and tan
  TProfile2D* fEffTPCPtTanF;  // TPC efficiency as function of Pt and tan
  //
  // dEdx resolution
  //
  TH2F* fTPCSignalNormTan; // tpc signal normalized to the mean signal - MC
  TH2F* fTPCSignalNormSPhi;   // tpc signal normalized to the mean signal - MC
  TH2F* fTPCSignalNormTPhi;   // tpc signal normalized to the mean signal - MC
  //
  TH3F* fTPCSignalNormTanSPhi;   // tpc signal normalized to the mean signal - MC
  TH3F* fTPCSignalNormTanTPhi;   // tpc signal normalized to the mean signal - MC
  TH3F* fTPCSignalNormTanSPt;   // tpc signal normalized to the mean signal - MC


  //
  //
  TH2F* fPtResolLPT;        // pt resolution - low pt
  TH2F* fPtResolHPT;        // pt resolution - high pt 
  TH2F* fPtPullLPT;         // pt resolution - low pt
  TH2F* fPtPullHPT;         // pt resolution - high pt 
  //
  // Resolution constrained param
  //
  TH2F   *fCPhiResolTan;   // angular resolution -  constrained
  TH2F   *fCTanResolTan;   // angular resolution -  constrained
  TH2F   *fCPtResolTan;    // pt resolution      -  constrained
  TH2F   *fCPhiPullTan;   // angular resolution -  constrained
  TH2F   *fCTanPullTan;   // angular resolution -  constrained
  TH2F   *fCPtPullTan;    // pt resolution      -  constrained
  //
  // DCA resolution
  //
  TH3F  *fD0TanSPtB1;   // distance to vertex y  
  TH3F  *fD1TanSPtB1;   // distance to vertex z  
  TH3F  *fD0TanSPtL1;   // distance to vertex y  
  TH3F  *fD1TanSPtL1;   // distance to vertex z  

protected:
   ClassDef(AliComparisonDraw,1);
};


#endif
Commit	Line	Data
06f2c0e2	1
	2	#ifndef AliComparisonDraw_h
	3	#define AliComparisonDraw_h
	4
	5	#include <iostream>
	6	#include <fstream>
	7	using namespace std;
	8	#include <TSelector.h>
	9
06f2c0e2	10
21486ec3	11	class TFile;
	12	class AliMCInfo;
	13	class AliESDRecInfo;
06f2c0e2	14	class AliESDEvent;
	15	class AliESD;
	16	class AliESDfriend;
	17	class TH1I;
2f396646	18	class TH3F;
	19	class TH3;
	20	class TProfile;
	21	class TProfile2D;
	22	class TGraph2D;
	23	class TGraph;
06f2c0e2	24
06f2c0e2	25	class AliComparisonDraw : public TObject {
	26	public :
	27	AliComparisonDraw();
	28	virtual Bool_t IsFolder(){return kTRUE;}
2f396646	29	void InitHisto();
	30	void Process(AliMCInfo* infoMC, AliESDRecInfo *infoRC);
	31
	32	//
	33	//
	34	void ProcessEff(AliMCInfo* infoMC, AliESDRecInfo *infoRC);
	35	void ProcessResolConstrained(AliMCInfo* infoMC, AliESDRecInfo *infoRC);
	36	void ProcessTPCdedx(AliMCInfo* infoMC, AliESDRecInfo *infoRC);
	37	void ProcessDCA(AliMCInfo* infoMC, AliESDRecInfo *infoRC);
	38
	39	void MakePlots();
	40
	41
06f2c0e2	42	//TH1F GetPtResol(Float_t pt0, Float_t pt1);
2f396646	43	static TH1F* MakeResol(TH2F * his, Int_t integ, Bool_t type);
	44	static TGraph2D * MakeStat2D(TH3 * his, Int_t delta0, Int_t delta1, Int_t type);
	45	static TGraph * MakeStat1D(TH3 * his, Int_t delta1, Int_t type);
	46
	47
21486ec3	48	protected:
2f396646	49	//
	50	// efficiency
	51	//
21486ec3	52	static Bool_t fgBDraw; //option draw temporary results
2f396646	53	TProfile* fEffTPCPt; // TPC efficiency as function of Pt (tan+-1)
	54	TProfile* fEffTPCPtMC; // MC -TPC efficiency as function of Pt (tan+-1)
	55	TProfile* fEffTPCPtF; // efficiency for findable tracks
	56	//
	57	TProfile* fEffTPCTan; // TPC efficiency as function of Tan (pt>0.15
	58	TProfile* fEffTPCTanMC; // MC -TPC efficiency as function of Tan (pt>0.15)
	59	TProfile* fEffTPCTanF; // efficiency for findable tracks Tan (pt>0.15)
	60	//
	61	TProfile2D* fEffTPCPtTan; // TPC efficiency as function of Pt and tan
	62	TProfile2D* fEffTPCPtTanMC; // MC -TPC efficiency as function of Pt and tan
	63	TProfile2D* fEffTPCPtTanF; // TPC efficiency as function of Pt and tan
	64	//
	65	// dEdx resolution
	66	//
	67	TH2F* fTPCSignalNormTan; // tpc signal normalized to the mean signal - MC
	68	TH2F* fTPCSignalNormSPhi; // tpc signal normalized to the mean signal - MC
	69	TH2F* fTPCSignalNormTPhi; // tpc signal normalized to the mean signal - MC
	70	//
	71	TH3F* fTPCSignalNormTanSPhi; // tpc signal normalized to the mean signal - MC
	72	TH3F* fTPCSignalNormTanTPhi; // tpc signal normalized to the mean signal - MC
	73	TH3F* fTPCSignalNormTanSPt; // tpc signal normalized to the mean signal - MC
	74
	75
	76	//
	77	//
06f2c0e2	78	TH2F* fPtResolLPT; // pt resolution - low pt
06f2c0e2	79	TH2F* fPtResolHPT; // pt resolution - high pt
2f396646	80	TH2F* fPtPullLPT; // pt resolution - low pt
	81	TH2F* fPtPullHPT; // pt resolution - high pt
	82	//
	83	// Resolution constrained param
	84	//
	85	TH2F *fCPhiResolTan; // angular resolution - constrained
	86	TH2F *fCTanResolTan; // angular resolution - constrained
	87	TH2F *fCPtResolTan; // pt resolution - constrained
	88	TH2F *fCPhiPullTan; // angular resolution - constrained
	89	TH2F *fCTanPullTan; // angular resolution - constrained
	90	TH2F *fCPtPullTan; // pt resolution - constrained
	91	//
	92	// DCA resolution
	93	//
	94	TH3F *fD0TanSPtB1; // distance to vertex y
	95	TH3F *fD1TanSPtB1; // distance to vertex z
	96	TH3F *fD0TanSPtL1; // distance to vertex y
	97	TH3F *fD1TanSPtL1; // distance to vertex z
	98
06f2c0e2	99	protected:
	100	ClassDef(AliComparisonDraw,1);
	101	};
	102
	103
	104
	105
	106
	107
	108
	109
	110
	111
	112
	113
	114
	115
	116
	117	#endif