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

// Analysis task for basic QA exploiting symmetries
//of global, TPC, and ITS tracks

#ifndef ALIANALYSISTASKQASYM_H
#define ALIANALYSISTASKQASYM_H
 

class TH1F;
class TH2F;
class TH3F;
class TList;

class AliESDEvent;
class AliESDtrack;
class AliESDtrackCuts;


#include "AliAnalysisTaskSE.h"
#include "TFile.h"

class AliAnalysisTaskQASym : public AliAnalysisTaskSE {
 public:
  AliAnalysisTaskQASym(const char *name = "AliAnalysisTaskQASym");
  virtual ~AliAnalysisTaskQASym() {}
  
  virtual void   UserCreateOutputObjects();
  virtual void   UserExec(Option_t *option);
  virtual void   Terminate(Option_t *);
  virtual void   SetTrackType(Int_t type) {fTrackType = type;}  
  virtual void   SetStandAloneTrack(Bool_t standAlone = kFALSE) {fStandAlone = standAlone;}  //needed for ITS tracks
  
  
  virtual void   SetCuts(AliESDtrackCuts* cuts)
     {fCuts = cuts;}

  virtual void   SetFieldOn(Bool_t b = kTRUE){fFieldOn = b;} 

  
 private:

  Int_t       fTrackType;       // track type
  Int_t       fStandAlone;      // needed for ITS tracks
  Bool_t      fFieldOn;         // field status

  TList       *fHists;          // List of histos

  //old
  TH1F        *fHistRECpt;      // pt 
  TH1F        *fEta;            // eta
  TH2F        *fEtaPhi;         // eta-phi
  TH1F        *fEtaPt;          // eta over pt 
  TH1F        *fQPt;            // charge over pt 
  TH1F        *fDca;            // distance of closest approach
  TH1F        *fqRec;           // reconstrcuted charge
  TH1F        *fSigmaPtHist;    // sigma_pT

   //positive und negative tracks
  TH1F        *fRecPtPos;      // pt of pos tracks
  TH1F        *fRecPtNeg;      // pt of neg tracks
  TH1F        *fRecPhiPos;     // phi of pos. tracks
  TH1F        *fRecPhiNeg;     // phi of neg. tracks
  TH1F        *fRecEtaPos;     // eta of neg. tracks
  TH1F        *fRecEtaNeg;     // eta of neg. tracks
  TH1F        *fRecEtaPtPos;   // eta over pt of neg. tracks
  TH1F        *fRecEtaPtNeg;   // eta over pt of neg. tracks
  TH1F        *fRecDcaPos;     // distance of closest approach of neg. tracks
  TH1F        *fRecDcaNeg;     // distance of closest of neg. tracks
  TH1F        *fRecDcaNegInv;  // invers dca of neg. tracks
  TH1F        *fRecDPos;       // impact parameter of neg. tracks
  TH1F        *fRecDNeg;       // impact parameter of neg. tracks

  // two sides of TPC -> Eta/Theta
  TH1F        *fRecQPtPosEta;   //charge/pT for pos. eta
  TH1F        *fRecQPtNegEta;   //charge/pT for neg. eta
  TH1F        *fRecPtPosEta;    //pT        for pos. eta
  TH1F        *fRecPtNegEta;    //pT        for neg. eta
  TH1F        *fRecPhiPosEta;   //phi       for pos. eta
  TH1F        *fRecPhiNegEta;   //phi       for neg. eta
  TH1F        *fRecDcaPosEta;   //dca       for pos. eta
  TH1F        *fRecDcaNegEta;   //dca       for neg. eta
  TH1F        *fRecDPosEta;     //d         for pos. eta
  TH1F        *fRecDNegEta;     //d         for neg. eta

  // 2D: all measures as function of z of first trackpoint 
  TH2F        *fRecPtPosVz;     //pt-zfirst of pos tracks
  TH2F        *fRecPtNegVz;     //pt-zfirst of neg tracks
  TH2F        *fRecEtaPosVz;    //eta-zfirst of pos tracks
  TH2F        *fRecEtaNegVz;    //eta-zfirst of neg tracks
  TH2F        *fRecPhiPosVz;    //phi-zfirst of pos tracks
  TH2F        *fRecPhiNegVz;    //phi-zfirst of neg tracks
  TH2F        *fSignedDcaPosVz; //dca-zfirst of pos tracks
  TH2F        *fSignedDcaNegVz; //dca-zfirst of neg tracks
  TH2F        *fRecQPtPosEtaVz; //charge/pT-zfirst of pos tracks
  TH2F        *fRecQPtNegEtaVz; //charge/pT-zfirst of neg tracks
  TH2F        *fRecEtaPtPosVz;  //eta/pT-zfirst of pos tracks
  TH2F        *fRecEtaPtNegVz;  //eta/pT-zfirst of neg tracks


  //high
  TH1F * fDeltaPhiAll;         // phiLeaingTracks-phiOthers
  TH2F * fDeltaPhiLeading;     // phiLeaingTracks-phiOthers vs. phiLeading
  TH1F * fDiffDcaD;            // d-dca

  //sim
  TH1F * fPhiRec;              //phi
  TH1F * fThetaRec;            //theta
  TH1F * fNumber;              //Number of tracks per event
  TH1F * fNumberAfterCut;      //Number of tracks per event after cuts
  TH1F * fVx;                  // x of first track point
  TH1F * fVy;                  // y of first track point
  TH1F * fVz;                  // z of first track point
  TH1F * fVertexX;             // x of vertex
  TH1F * fVertexY;             // y of vertex
  TH1F * fVertexZ;             // z of vertex
  TH1F * fNVertexSPD;          //number of vertices SPD
  TH1F * fNVertexTracks;       //number of vertices of Tracks

  //new
  TH2F        *fRecDcaPosPhi;     //dca-phi for pos.
  TH2F        *fRecDcaNegPhi;     //dca-phi for neg.
  TH2F        *fRecPtPosPhi;      //pt-phi for pos.
  TH2F        *fRecPtNegPhi;      //pt-phi for neg.
  TH2F        *fRecEtaPosPhi;     //eta-phi for pos.
  TH2F        *fRecEtaNegPhi;     //eta-phi for neg.
  TH2F        *fRecQPtPhi;        //charge/pt-phi
  TH2F        *fRecEtaPtPosPhi;   //eta/pt-phi for neg.
  TH2F        *fRecEtaPtNegPhi;   //eta/pt-phi for pos.

  TH1F        *fRecPtPosEtaPos;   //pt for pos tracks and pos eta
  TH1F        *fRecPtNegEtaPos;   //pt for neg tracks and pos eta
  TH1F        *fRecPtPosEtaNeg;   //pt for pos tracks and neg eta
  TH1F        *fRecPtNegEtaNeg;   //pt for neg tracks and neg eta

  TH1F        *fRec1PtPosEtaPos;   //1/pt for pos tracks and pos eta
  TH1F        *fRec1PtNegEtaPos;   //1/pt for neg tracks and pos eta
  TH1F        *fRec1PtPosEtaNeg;   //1/pt for pos tracks and neg eta
  TH1F        *fRec1PtNegEtaNeg;   //1/pt for neg tracks and neg eta

  TH1F        *fRecPhiPosEtaPos;   //phi for pos tracks and pos eta
  TH1F        *fRecPhiNegEtaPos;   //phi for neg tracks and pos eta
  TH1F        *fRecPhiPosEtaNeg;   //phi for pos tracks and neg eta
  TH1F        *fRecPhiNegEtaNeg;   //phi for neg tracks and neg eta

  TH2F        *fRecDcaPosPhiEtaPos;  //dca-phi for pos tracks and pos eta
  TH2F        *fRecDcaNegPhiEtaPos;  //dca-phi for neg tracks and pos eta
  TH2F        *fRecDcaPosPhiEtaNeg;  //dca-phi for pos tracks and neg eta
  TH2F        *fRecDcaNegPhiEtaNeg;  //dca-phi for neg tracks and neg eta

  TH2F        *fRecDcaPosPtEtaPos;  //dca-pt for pos tracks and pos eta
  TH2F        *fRecDcaNegPtEtaPos;  //dca-pt for neg tracks and pos eta
  TH2F        *fRecDcaPosPtEtaNeg;  //dca-pt for pos tracks and neg eta
  TH2F        *fRecDcaNegPtEtaNeg;  //dca-pt for neg tracks and neg eta

  TH2F        *fRecPtPosPhiEtaPos;  //pt-phi for pos tracks and pos eta
  TH2F        *fRecPtNegPhiEtaPos;  //pt-phi for neg tracks and pos eta 
  TH2F        *fRecPtPosPhiEtaNeg;  //pt-phi for pos tracks and neg eta
  TH2F        *fRecPtNegPhiEtaNeg;  //pt-phi for neg tracks and neg eta

  //  TH3F        *fRecDcaPhiPtPosEtaPos; //dca-pt-phi for pos tracks and pos eta
  //  TH3F        *fRecDcaPhiPtNegEtaPos; //dca-pt-phi for neg tracks and pos eta
  //  TH3F        *fRecDcaPhiPtPosEtaNeg; //dca-pt-phi for pos tracks and neg eta
  //  TH3F        *fRecDcaPhiPtNegEtaNeg; //dca-pt-phi for neg tracks and neg eta

  TH2F        *fEtavPt;                 // eta vs pt 
  TH2F        *fPhivPt;                 // phi vs pt 

  TH2F        *fCompareTPCparam;        // TPC param

  TH1F        *fITSlayer;               // ITS layer
  TH2F        *fITSlayerEta;            // ITS layer vs eta
  TH2F        *fITSlayerPhi;            // ITS layer vs phi

  AliESDtrackCuts* fCuts;               // List of cuts

  // sectors of TPC 
  TH1F        *fRecPtTpcSector[18];     //pt for TPC sectors
  TH1F        *fRecEtaTpcSector[18];    //eta for TPC sectors
  TH1F        *fRecQPtTpcSector[18];    //charge/pt for TPC sectors
  TH1F        *fRecEtaPtTpcSector[18];  //eta/pt for TPC sectors
  TH1F        *fSignedDcaTpcSector[18]; //dca for TPC sectors


  // 7 different case of hit in ITS ladders
  TH1F        *fRecPtPosLadder[7];  //pt for pos tracks
  TH1F        *fRecPtNegLadder[7];  //pt for neg tracks
  TH1F        *fRecPhiPosLadder[7]; //phi for pos tracks
  TH1F        *fRecPhiNegLadder[7]; //phi for neg tracks
  TH1F        *fRecEtaPosLadder[7]; //eta for pos tracks
  TH1F        *fRecEtaNegLadder[7]; //eta for neg tracks
  TH1F        *fSignDcaPos[7];      //dca for pos tracks
  TH1F        *fSignDcaNeg[7];      //dca for neg tracks
  TH1F        *fSignDcaNegInv[7];   //-dca for neg tracks
  TH1F        *fPtSigmaPos[7];      //sigma_pT for pos tracks
  TH1F        *fPtSigmaNeg[7];      //sigma_pT for neg tracks
  TH1F        *fqPtRec[7];          // charge/pt 
  TH2F        *fDcaSigmaPos[7];     // dca - sigma_pT for pos tracks
  TH2F        *fDcaSigmaNeg[7];     // dca - sigma_pT for neg tracks

  
  AliAnalysisTaskQASym(const AliAnalysisTaskQASym&); // not implemented
  AliAnalysisTaskQASym& operator=(const AliAnalysisTaskQASym&); // not implemented
  
  ClassDef(AliAnalysisTaskQASym, 1); // Basic QA exploiting symmetries
};

#endif
Commit	Line	Data
ddc7ef22	1	// Analysis task for basic QA exploiting symmetries
	2	//of global, TPC, and ITS tracks
	3
8d49bc89	4	#ifndef ALIANALYSISTASKQASYM_H
8d49bc89	5	#define ALIANALYSISTASKQASYM_H
9f9631f3	6
	7
	8	class TH1F;
	9	class TH2F;
1d55ca60	10	class TH3F;
9f9631f3	11	class TList;
9f9631f3	12
	13	class AliESDEvent;
	14	class AliESDtrack;
	15	class AliESDtrackCuts;
	16
	17
	18	#include "AliAnalysisTaskSE.h"
1d55ca60	19	#include "TFile.h"
9f9631f3	20
	21	class AliAnalysisTaskQASym : public AliAnalysisTaskSE {
	22	public:
1d55ca60	23	AliAnalysisTaskQASym(const char *name = "AliAnalysisTaskQASym");
	24	virtual ~AliAnalysisTaskQASym() {}
	25
	26	virtual void UserCreateOutputObjects();
	27	virtual void UserExec(Option_t *option);
	28	virtual void Terminate(Option_t *);
	29	virtual void SetTrackType(Int_t type) {fTrackType = type;}
b676e981	30	virtual void SetStandAloneTrack(Bool_t standAlone = kFALSE) {fStandAlone = standAlone;} //needed for ITS tracks
9f9631f3	31
9f9631f3	32
1d55ca60	33	virtual void SetCuts(AliESDtrackCuts* cuts)
	34	{fCuts = cuts;}
	35
	36	virtual void SetFieldOn(Bool_t b = kTRUE){fFieldOn = b;}
	37
9f9631f3	38
	39	private:
	40
8d49bc89	41	Int_t fTrackType; // track type
	42	Int_t fStandAlone; // needed for ITS tracks
	43	Bool_t fFieldOn; // field status
9f9631f3	44
	45	TList *fHists; // List of histos
	46
	47	//old
	48	TH1F *fHistRECpt; // pt
	49	TH1F *fEta; // eta
e4ab7356	50	TH2F *fEtaPhi; // eta-phi
9f9631f3	51	TH1F *fEtaPt; // eta over pt
	52	TH1F *fQPt; // charge over pt
	53	TH1F *fDca; // distance of closest approach
9f9631f3	54	TH1F *fqRec; // reconstrcuted charge
ddc7ef22	55	TH1F *fSigmaPtHist; // sigma_pT
1d55ca60	56
	57	//positive und negative tracks
	58	TH1F *fRecPtPos; // pt of pos tracks
	59	TH1F *fRecPtNeg; // pt of neg tracks
	60	TH1F *fRecPhiPos; // phi of pos. tracks
	61	TH1F *fRecPhiNeg; // phi of neg. tracks
	62	TH1F *fRecEtaPos; // eta of neg. tracks
	63	TH1F *fRecEtaNeg; // eta of neg. tracks
	64	TH1F *fRecEtaPtPos; // eta over pt of neg. tracks
	65	TH1F *fRecEtaPtNeg; // eta over pt of neg. tracks
	66	TH1F *fRecDcaPos; // distance of closest approach of neg. tracks
	67	TH1F *fRecDcaNeg; // distance of closest of neg. tracks
	68	TH1F *fRecDcaNegInv; // invers dca of neg. tracks
	69	TH1F *fRecDPos; // impact parameter of neg. tracks
	70	TH1F *fRecDNeg; // impact parameter of neg. tracks
9f9631f3	71
9f9631f3	72	// two sides of TPC -> Eta/Theta
1d55ca60	73	TH1F *fRecQPtPosEta; //charge/pT for pos. eta
	74	TH1F *fRecQPtNegEta; //charge/pT for neg. eta
	75	TH1F *fRecPtPosEta; //pT for pos. eta
	76	TH1F *fRecPtNegEta; //pT for neg. eta
	77	TH1F *fRecPhiPosEta; //phi for pos. eta
	78	TH1F *fRecPhiNegEta; //phi for neg. eta
	79	TH1F *fRecDcaPosEta; //dca for pos. eta
	80	TH1F *fRecDcaNegEta; //dca for neg. eta
	81	TH1F *fRecDPosEta; //d for pos. eta
	82	TH1F *fRecDNegEta; //d for neg. eta
	83
	84	// 2D: all measures as function of z of first trackpoint
	85	TH2F *fRecPtPosVz; //pt-zfirst of pos tracks
8d49bc89	86	TH2F *fRecPtNegVz; //pt-zfirst of neg tracks
1d55ca60	87	TH2F *fRecEtaPosVz; //eta-zfirst of pos tracks
	88	TH2F *fRecEtaNegVz; //eta-zfirst of neg tracks
	89	TH2F *fRecPhiPosVz; //phi-zfirst of pos tracks
	90	TH2F *fRecPhiNegVz; //phi-zfirst of neg tracks
	91	TH2F *fSignedDcaPosVz; //dca-zfirst of pos tracks
	92	TH2F *fSignedDcaNegVz; //dca-zfirst of neg tracks
	93	TH2F *fRecQPtPosEtaVz; //charge/pT-zfirst of pos tracks
	94	TH2F *fRecQPtNegEtaVz; //charge/pT-zfirst of neg tracks
	95	TH2F *fRecEtaPtPosVz; //eta/pT-zfirst of pos tracks
	96	TH2F *fRecEtaPtNegVz; //eta/pT-zfirst of neg tracks
9f9631f3	97
9f9631f3	98
9f9631f3	99	//high
1d55ca60	100	TH1F * fDeltaPhiAll; // phiLeaingTracks-phiOthers
	101	TH2F * fDeltaPhiLeading; // phiLeaingTracks-phiOthers vs. phiLeading
	102	TH1F * fDiffDcaD; // d-dca
9f9631f3	103
9f9631f3	104	//sim
1d55ca60	105	TH1F * fPhiRec; //phi
	106	TH1F * fThetaRec; //theta
	107	TH1F * fNumber; //Number of tracks per event
ddc7ef22	108	TH1F * fNumberAfterCut; //Number of tracks per event after cuts
1d55ca60	109	TH1F * fVx; // x of first track point
	110	TH1F * fVy; // y of first track point
	111	TH1F * fVz; // z of first track point
cbac53a7	112	TH1F * fVertexX; // x of vertex
	113	TH1F * fVertexY; // y of vertex
	114	TH1F * fVertexZ; // z of vertex
ddc7ef22	115	TH1F * fNVertexSPD; //number of vertices SPD
ddc7ef22	116	TH1F * fNVertexTracks; //number of vertices of Tracks
1d55ca60	117
	118	//new
	119	TH2F *fRecDcaPosPhi; //dca-phi for pos.
	120	TH2F *fRecDcaNegPhi; //dca-phi for neg.
	121	TH2F *fRecPtPosPhi; //pt-phi for pos.
	122	TH2F *fRecPtNegPhi; //pt-phi for neg.
	123	TH2F *fRecEtaPosPhi; //eta-phi for pos.
	124	TH2F *fRecEtaNegPhi; //eta-phi for neg.
	125	TH2F *fRecQPtPhi; //charge/pt-phi
	126	TH2F *fRecEtaPtPosPhi; //eta/pt-phi for neg.
	127	TH2F *fRecEtaPtNegPhi; //eta/pt-phi for pos.
	128
	129	TH1F *fRecPtPosEtaPos; //pt for pos tracks and pos eta
	130	TH1F *fRecPtNegEtaPos; //pt for neg tracks and pos eta
	131	TH1F *fRecPtPosEtaNeg; //pt for pos tracks and neg eta
	132	TH1F *fRecPtNegEtaNeg; //pt for neg tracks and neg eta
	133
	134	TH1F *fRec1PtPosEtaPos; //1/pt for pos tracks and pos eta
	135	TH1F *fRec1PtNegEtaPos; //1/pt for neg tracks and pos eta
	136	TH1F *fRec1PtPosEtaNeg; //1/pt for pos tracks and neg eta
	137	TH1F *fRec1PtNegEtaNeg; //1/pt for neg tracks and neg eta
	138
	139	TH1F *fRecPhiPosEtaPos; //phi for pos tracks and pos eta
	140	TH1F *fRecPhiNegEtaPos; //phi for neg tracks and pos eta
	141	TH1F *fRecPhiPosEtaNeg; //phi for pos tracks and neg eta
	142	TH1F *fRecPhiNegEtaNeg; //phi for neg tracks and neg eta
	143
	144	TH2F *fRecDcaPosPhiEtaPos; //dca-phi for pos tracks and pos eta
	145	TH2F *fRecDcaNegPhiEtaPos; //dca-phi for neg tracks and pos eta
	146	TH2F *fRecDcaPosPhiEtaNeg; //dca-phi for pos tracks and neg eta
	147	TH2F *fRecDcaNegPhiEtaNeg; //dca-phi for neg tracks and neg eta
	148
98a5f582	149	TH2F *fRecDcaPosPtEtaPos; //dca-pt for pos tracks and pos eta
	150	TH2F *fRecDcaNegPtEtaPos; //dca-pt for neg tracks and pos eta
	151	TH2F *fRecDcaPosPtEtaNeg; //dca-pt for pos tracks and neg eta
	152	TH2F *fRecDcaNegPtEtaNeg; //dca-pt for neg tracks and neg eta
	153
1d55ca60	154	TH2F *fRecPtPosPhiEtaPos; //pt-phi for pos tracks and pos eta
	155	TH2F *fRecPtNegPhiEtaPos; //pt-phi for neg tracks and pos eta
	156	TH2F *fRecPtPosPhiEtaNeg; //pt-phi for pos tracks and neg eta
	157	TH2F *fRecPtNegPhiEtaNeg; //pt-phi for neg tracks and neg eta
	158
8ae6f5d0	159	// TH3F *fRecDcaPhiPtPosEtaPos; //dca-pt-phi for pos tracks and pos eta
	160	// TH3F *fRecDcaPhiPtNegEtaPos; //dca-pt-phi for neg tracks and pos eta
	161	// TH3F *fRecDcaPhiPtPosEtaNeg; //dca-pt-phi for pos tracks and neg eta
	162	// TH3F *fRecDcaPhiPtNegEtaNeg; //dca-pt-phi for neg tracks and neg eta
1d55ca60	163
8d49bc89	164	TH2F *fEtavPt; // eta vs pt
ddc7ef22	165	TH2F *fPhivPt; // phi vs pt
ddc7ef22	166
8d49bc89	167	TH2F *fCompareTPCparam; // TPC param
ddc7ef22	168
8d49bc89	169	TH1F *fITSlayer; // ITS layer
ddc7ef22	170	TH2F *fITSlayerEta; // ITS layer vs eta
ddc7ef22	171	TH2F *fITSlayerPhi; // ITS layer vs phi
1d55ca60	172
8d49bc89	173	AliESDtrackCuts* fCuts; // List of cuts
1d55ca60	174
	175	// sectors of TPC
	176	TH1F *fRecPtTpcSector[18]; //pt for TPC sectors
	177	TH1F *fRecEtaTpcSector[18]; //eta for TPC sectors
	178	TH1F *fRecQPtTpcSector[18]; //charge/pt for TPC sectors
	179	TH1F *fRecEtaPtTpcSector[18]; //eta/pt for TPC sectors
	180	TH1F *fSignedDcaTpcSector[18]; //dca for TPC sectors
	181
	182
	183	// 7 different case of hit in ITS ladders
	184	TH1F *fRecPtPosLadder[7]; //pt for pos tracks
	185	TH1F *fRecPtNegLadder[7]; //pt for neg tracks
	186	TH1F *fRecPhiPosLadder[7]; //phi for pos tracks
	187	TH1F *fRecPhiNegLadder[7]; //phi for neg tracks
	188	TH1F *fRecEtaPosLadder[7]; //eta for pos tracks
	189	TH1F *fRecEtaNegLadder[7]; //eta for neg tracks
	190	TH1F *fSignDcaPos[7]; //dca for pos tracks
	191	TH1F *fSignDcaNeg[7]; //dca for neg tracks
	192	TH1F *fSignDcaNegInv[7]; //-dca for neg tracks
	193	TH1F *fPtSigmaPos[7]; //sigma_pT for pos tracks
	194	TH1F *fPtSigmaNeg[7]; //sigma_pT for neg tracks
	195	TH1F *fqPtRec[7]; // charge/pt
	196	TH2F *fDcaSigmaPos[7]; // dca - sigma_pT for pos tracks
	197	TH2F *fDcaSigmaNeg[7]; // dca - sigma_pT for neg tracks
	198
	199
	200
	201
	202
9f9631f3	203	AliAnalysisTaskQASym(const AliAnalysisTaskQASym&); // not implemented
	204	AliAnalysisTaskQASym& operator=(const AliAnalysisTaskQASym&); // not implemented
	205
	206	ClassDef(AliAnalysisTaskQASym, 1); // Basic QA exploiting symmetries
	207	};
	208
	209	#endif