[u/mrichter/AliRoot.git] / PWGGA / CaloTrackCorrelations / AliAnaChargedParticles.h

#ifndef ALIANACHARGEDPARTICLES_H
#define ALIANACHARGEDPARTICLES_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice     */

//_________________________________________________________________________
//
// Class for track selection and identification (not done now)
// Tracks from the CTS are kept in the AOD.
// Few histograms produced.
//
//-- Author: Gustavo Conesa (INFN-LNF)

// Root system
class TH2F; 

// Analysis system
#include "AliAnaCaloTrackCorrBaseClass.h"
 
class AliAnaChargedParticles : public AliAnaCaloTrackCorrBaseClass {
  
 public: 
  AliAnaChargedParticles() ; // default ctor
  virtual ~AliAnaChargedParticles() { ; } //virtual dtor
  
  TList * GetCreateOutputObjects();
    
  void    Init();
  
  void    InitParameters();
  
  void    Print(const Option_t * opt) const;
  
  void    MakeAnalysisFillAOD()  ;
  
  void    MakeAnalysisFillHistograms() ; 
  
  void    SwitchOnFillPileUpHistograms()     { fFillPileUpHistograms    = kTRUE  ; }
  void    SwitchOffFillPileUpHistograms()    { fFillPileUpHistograms    = kFALSE ; }
  
  void    SwitchOnFillTrackBCHistograms()  { fFillVertexBC0Histograms = kTRUE  ; }
  void    SwitchOffFillTrackBCHistograms() { fFillVertexBC0Histograms = kFALSE ; }

  void    SwitchOnFillVertexBC0Histograms()  { fFillVertexBC0Histograms = kTRUE  ; }
  void    SwitchOffFillVertexBC0Histograms() { fFillVertexBC0Histograms = kFALSE ; }

 private:
  
  Bool_t  fFillPileUpHistograms;    // Fill pile-up related histograms
  Bool_t  fFillTrackBCHistograms;   // Fill histograms for tracks with TOF BC=0 or not related histograms
  Bool_t  fFillVertexBC0Histograms; // Fill histograms for tracks with vertex BC=0 or not related histograms
  
  //Histograms
  TH1F * fhNtracks;     //! track multiplicity distribution
  TH1F * fhPt;          //! pT distribution
  TH1F * fhPtNoCut;     //! pT distribution, no cut
  TH1F * fhPtCutDCA;    //! pT distribution, Apply DCA cut
  TH1F * fhPtCutDCABCOK;//! pT distribution, Apply DCA cut, BC=0 or -100

  TH1F * fhPtPileUp[7]; //! pT distribution, pile-up defined events
  TH2F * fhPhiNeg;      //! phi distribution vs pT, negative
  TH2F * fhEtaNeg;      //! eta distribution vs pT, negative
  TH2F * fhPhiPos;      //! phi distribution vs pT, positive
  TH2F * fhEtaPos;      //! eta distribution vs pT, positive
  TH2F * fhEtaPhiPos;   //! eta vs phi distribution of positive charge  
  TH2F * fhEtaPhiNeg;   //! eta vs phi distribution of negative charge
  
  TH1F * fhPtVtxOutBC0;    //! pT distribution of tracks from a vertex with BC!=0
  TH2F * fhEtaPhiVtxOutBC0;//! eta/phi distribution of tracks from a vertex with BC!=0
  TH1F * fhPtVtxInBC0;     //! pT distribution of tracks from a vertex with BC=0
  TH2F * fhEtaPhiVtxInBC0; //! eta/phi distribution of tracks from a vertex with BC=0

  TH1F * fhPtSPDRefit;     //! pT distribution of tracks with SPD and ITS refit
  TH1F * fhPtNoSPDRefit;   //! pT distribution of constrained tracks no SPD and with ITSRefit
  TH1F * fhPtNoSPDNoRefit; //! pT distribution of constrained tracks with no SPD requierement and without ITSRefit

  TH2F * fhEtaPhiSPDRefitPt02;     //! eta-phi distribution of tracks with SPD and ITS refit, 0 < pT < 2 GeV
  TH2F * fhEtaPhiNoSPDRefitPt02;   //! eta-phi distribution of constrained tracks no SPD and with ITSRefit,  0 < pT < 2 GeV
  TH2F * fhEtaPhiNoSPDNoRefitPt02; //! eta-phi distribution of constrained tracks with no SPD requierement and without ITSRefit,  0 < pT < 2 GeV

  TH2F * fhEtaPhiSPDRefitPt3;     //! eta-phi distribution of tracks with SPD and ITS refit, pT > 3 GeV
  TH2F * fhEtaPhiNoSPDRefitPt3;   //! eta-phi distribution of constrained tracks no SPD and with ITSRefit,  pT > 3 GeV
  TH2F * fhEtaPhiNoSPDNoRefitPt3; //! eta-phi distribution of constrained tracks with no SPD requierement and without ITSRefit,  pT > 3 GeV

  //MC
  TH1F * fhPtPion;      //! pT distribution
  TH2F * fhPhiPion;     //! phi distribution vs pT
  TH2F * fhEtaPion;     //! eta distribution vs pT
  
  TH1F * fhPtProton;    //! pT distribution
  TH2F * fhPhiProton;   //! phi distribution vs pT
  TH2F * fhEtaProton;   //! eta distribution vs pT
  
  TH1F * fhPtElectron;  //! pT distribution
  TH2F * fhPhiElectron; //! phi distribution vs pT
  TH2F * fhEtaElectron; //! eta distribution vs pT
  
  TH1F * fhPtKaon;      //! pT distribution
  TH2F * fhPhiKaon;     //! phi distribution vs pT
  TH2F * fhEtaKaon;     //! eta distribution vs pT
  
  TH1F * fhPtUnknown;   //! pT distribution
  TH2F * fhPhiUnknown;  //! phi distribution vs pT
  TH2F * fhEtaUnknown;  //! eta distribution vs pT

  TH1F * fhMCPt; //! Number of generated charged hadrons vs pT coming from MC particle
  TH2F * fhMCPhi; //! Number of generated charged hadrons vs phi coming from MC particle
  TH2F * fhMCEta; //! Number of generated charged hadrons vs eta coming from MC particle
  TH1F * fhMCRecPt; //! Number of reconstructed charged hadrons vs pT coming from MC particle

  // TOF
  TH1F * fhTOFSignal;                    //! TOF signal
  TH1F * fhTOFSignalPtCut;               //! TOF signal pt and acceptance cut
  TH1F * fhTOFSignalBCOK;                //! TOF signal pt and acceptance cut
  TH2F * fhPtTOFSignal;                  //! TOF signal vs track pT, good status
  TH2F * fhPtTOFSignalDCACut;            //! TOF signal vs track pT, good status
  TH2F * fhPtTOFSignalPileUp[7];         //! TOF signal vs track pT, good status, pile-up
  TH2F * fhPtTOFSignalVtxOutBC0;         //! TOF signal vs track pT, good status
  TH2F * fhPtTOFSignalVtxOutBC0PileUp[7];//! TOF signal vs track pT, good status, pile-up
  TH2F * fhPtTOFSignalVtxInBC0;          //! TOF signal vs track pT, good status
  TH2F * fhPtTOFSignalVtxInBC0PileUp[7]; //! TOF signal vs track pT, good status, pile-up
  TH1F * fhPtTOFStatus0;                 //! pT of tracks not passing TOF status selection
  TH2F * fhEtaPhiTOFStatus0;             //! eta/phi of tracks not passing TOF status selection
  TH2F * fhEtaPhiTOFBC0;                 //! eta/phi of tracks passing TOF status selection, tracks in BC=0
  TH2F * fhEtaPhiTOFBCPlus;              //! eta/phi of tracks passing TOF status selection, tracks in BC>0
  TH2F * fhEtaPhiTOFBCMinus;             //! eta/phi of tracks passing TOF status selection, tracks in BC<0
  TH2F * fhEtaPhiTOFBC0PileUpSPD;        //! eta/phi of tracks passing TOF status selection, tracks in BC=0, pile-up spd
  TH2F * fhEtaPhiTOFBCPlusPileUpSPD;     //! eta/phi of tracks passing TOF status selection, tracks in BC>0, pile-up spd
  TH2F * fhEtaPhiTOFBCMinusPileUpSPD;    //! eta/phi of tracks passing TOF status selection, tracks in BC<0, pile-up spd
  
  TH1F * fhProductionVertexBC;           //!  check BC of production vertex
  TH1F * fhProductionVertexBCPileUp[7];  //!  check BC of production vertex, pile-up

  TH2F * fhPtDCA[3];                     //! DCA (xy,z,constrained) of all tracks
  
  TH2F * fhPtDCASPDRefit[3];             //! DCA (xy,z,constrained) of tracks with SPD and ITS refit
  TH2F * fhPtDCANoSPDRefit[3];           //! DCA (xy,z,constrained) of constrained tracks no SPD and with ITSRefit
  TH2F * fhPtDCANoSPDNoRefit[3];         //! DCA (xy,z,constrained) of constrained tracks with no SPD requierement and without ITSRefit

  TH2F * fhPtDCAVtxOutBC0[3];            //! DCA (xy,z,constrained) of all tracks, vertex BC!=0
  TH2F * fhPtDCAVtxInBC0[3];             //! DCA (xy,z,constrained) of all tracks, vertex BC==0
  TH2F * fhPtDCAPileUp[3];               //! DCA (xy,z,constrained) of all tracks, SPD pile-up
  TH2F * fhPtDCAVtxOutBC0PileUp[3];      //! DCA (xy,z,constrained) of all tracks, vertex BC!=0, SPD pile-up
  TH2F * fhPtDCAVtxInBC0PileUp[3];       //! DCA (xy,z,constrained) of all tracks, vertex BC==0, SPD pile-up

  TH2F * fhPtDCATOFBC0[3];               //! DCA (xy,z,constrained) of all tracks, hit in TOF and BC=0
  TH2F * fhPtDCAPileUpTOFBC0[3];         //! DCA (xy,z,constrained) of all tracks, hit in TOF and BC=0
  TH2F * fhPtDCATOFBCOut[3];             //! DCA (xy,z,constrained) of all tracks, hit in TOF and BC!=0

  TH2F * fhPtDCANoTOFHit[3];                //! DCA (xy,z,constrained) of all tracks, no hit in TOF
  TH2F * fhPtDCAVtxOutBC0NoTOFHit[3];       //! DCA (xy,z,constrained) of all tracks, vertex BC!=0, no hit in TOF
  TH2F * fhPtDCAVtxInBC0NoTOFHit[3];        //! DCA (xy,z,constrained) of all tracks, vertex BC=0, no hit in TOF
  TH2F * fhPtDCAPileUpNoTOFHit[3];          //! DCA (xy,z,constrained) of all tracks, SPD pile-up, no hit in TOF
  TH2F * fhPtDCAVtxOutBC0PileUpNoTOFHit[3]; //! DCA (xy,z,constrained) of all tracks, vertex BC!=0, SPD pile-up, no hit in TOF
  TH2F * fhPtDCAVtxInBC0PileUpNoTOFHit[3];  //! DCA (xy,z,constrained) of all tracks, vertex BC=0, SPD pile-up, no hit in TOF
  
  
  AliAnaChargedParticles(              const AliAnaChargedParticles & ch) ; // cpy ctor
  AliAnaChargedParticles & operator = (const AliAnaChargedParticles & ch) ; // cpy assignment
  
  ClassDef(AliAnaChargedParticles,8)

} ;


#endif //ALIANACHARGEDPARTICLES_H
Commit	Line	Data
477d6cee	1	#ifndef ALIANACHARGEDPARTICLES_H
	2	#define ALIANACHARGEDPARTICLES_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
477d6cee	5
477d6cee	6	//_________________________________________________________________________
34c16486	7	//
	8	// Class for track selection and identification (not done now)
	9	// Tracks from the CTS are kept in the AOD.
	10	// Few histograms produced.
477d6cee	11	//
477d6cee	12	//-- Author: Gustavo Conesa (INFN-LNF)
34c16486	13
7cd4e982	14	// Root system
477d6cee	15	class TH2F;
05d0d05d	16
7cd4e982	17	// Analysis system
745913ae	18	#include "AliAnaCaloTrackCorrBaseClass.h"
477d6cee	19
745913ae	20	class AliAnaChargedParticles : public AliAnaCaloTrackCorrBaseClass {
477d6cee	21
477d6cee	22	public:
477d6cee	23	AliAnaChargedParticles() ; // default ctor
05d0d05d	24	virtual ~AliAnaChargedParticles() { ; } //virtual dtor
52bd777f	25
477d6cee	26	TList * GetCreateOutputObjects();
cc944149	27
05d0d05d	28	void Init();
477d6cee	29
05d0d05d	30	void InitParameters();
477d6cee	31
05d0d05d	32	void Print(const Option_t * opt) const;
477d6cee	33
05d0d05d	34	void MakeAnalysisFillAOD() ;
477d6cee	35
05d0d05d	36	void MakeAnalysisFillHistograms() ;
477d6cee	37
52bd777f	38	void SwitchOnFillPileUpHistograms() { fFillPileUpHistograms = kTRUE ; }
52bd777f	39	void SwitchOffFillPileUpHistograms() { fFillPileUpHistograms = kFALSE ; }
fedea415	40
9b4a1ee2	41	void SwitchOnFillTrackBCHistograms() { fFillVertexBC0Histograms = kTRUE ; }
	42	void SwitchOffFillTrackBCHistograms() { fFillVertexBC0Histograms = kFALSE ; }
	43
52bd777f	44	void SwitchOnFillVertexBC0Histograms() { fFillVertexBC0Histograms = kTRUE ; }
	45	void SwitchOffFillVertexBC0Histograms() { fFillVertexBC0Histograms = kFALSE ; }
	46
fedea415	47	private:
05d0d05d	48
cbd6236f	49	Bool_t fFillPileUpHistograms; // Fill pile-up related histograms
9b4a1ee2	50	Bool_t fFillTrackBCHistograms; // Fill histograms for tracks with TOF BC=0 or not related histograms
cbd6236f	51	Bool_t fFillVertexBC0Histograms; // Fill histograms for tracks with vertex BC=0 or not related histograms
52bd777f	52
fedea415	53	//Histograms
05d0d05d	54	TH1F * fhNtracks; //! track multiplicity distribution
05d0d05d	55	TH1F * fhPt; //! pT distribution
52bd777f	56	TH1F * fhPtNoCut; //! pT distribution, no cut
	57	TH1F * fhPtCutDCA; //! pT distribution, Apply DCA cut
	58	TH1F * fhPtCutDCABCOK;//! pT distribution, Apply DCA cut, BC=0 or -100
	59
fedea415	60	TH1F * fhPtPileUp[7]; //! pT distribution, pile-up defined events
05d0d05d	61	TH2F * fhPhiNeg; //! phi distribution vs pT, negative
	62	TH2F * fhEtaNeg; //! eta distribution vs pT, negative
	63	TH2F * fhPhiPos; //! phi distribution vs pT, positive
	64	TH2F * fhEtaPos; //! eta distribution vs pT, positive
	65	TH2F * fhEtaPhiPos; //! eta vs phi distribution of positive charge
fe9615b9	66	TH2F * fhEtaPhiNeg; //! eta vs phi distribution of negative charge
fe9615b9	67
52bd777f	68	TH1F * fhPtVtxOutBC0; //! pT distribution of tracks from a vertex with BC!=0
	69	TH2F * fhEtaPhiVtxOutBC0;//! eta/phi distribution of tracks from a vertex with BC!=0
	70	TH1F * fhPtVtxInBC0; //! pT distribution of tracks from a vertex with BC=0
	71	TH2F * fhEtaPhiVtxInBC0; //! eta/phi distribution of tracks from a vertex with BC=0
	72
8f467228	73	TH1F * fhPtSPDRefit; //! pT distribution of tracks with SPD and ITS refit
	74	TH1F * fhPtNoSPDRefit; //! pT distribution of constrained tracks no SPD and with ITSRefit
	75	TH1F * fhPtNoSPDNoRefit; //! pT distribution of constrained tracks with no SPD requierement and without ITSRefit
	76
	77	TH2F * fhEtaPhiSPDRefitPt02; //! eta-phi distribution of tracks with SPD and ITS refit, 0 < pT < 2 GeV
	78	TH2F * fhEtaPhiNoSPDRefitPt02; //! eta-phi distribution of constrained tracks no SPD and with ITSRefit, 0 < pT < 2 GeV
	79	TH2F * fhEtaPhiNoSPDNoRefitPt02; //! eta-phi distribution of constrained tracks with no SPD requierement and without ITSRefit, 0 < pT < 2 GeV
	80
	81	TH2F * fhEtaPhiSPDRefitPt3; //! eta-phi distribution of tracks with SPD and ITS refit, pT > 3 GeV
	82	TH2F * fhEtaPhiNoSPDRefitPt3; //! eta-phi distribution of constrained tracks no SPD and with ITSRefit, pT > 3 GeV
	83	TH2F * fhEtaPhiNoSPDNoRefitPt3; //! eta-phi distribution of constrained tracks with no SPD requierement and without ITSRefit, pT > 3 GeV
	84
477d6cee	85	//MC
05d0d05d	86	TH1F * fhPtPion; //! pT distribution
	87	TH2F * fhPhiPion; //! phi distribution vs pT
	88	TH2F * fhEtaPion; //! eta distribution vs pT
477d6cee	89
05d0d05d	90	TH1F * fhPtProton; //! pT distribution
	91	TH2F * fhPhiProton; //! phi distribution vs pT
	92	TH2F * fhEtaProton; //! eta distribution vs pT
477d6cee	93
05d0d05d	94	TH1F * fhPtElectron; //! pT distribution
477d6cee	95	TH2F * fhPhiElectron; //! phi distribution vs pT
	96	TH2F * fhEtaElectron; //! eta distribution vs pT
	97
05d0d05d	98	TH1F * fhPtKaon; //! pT distribution
	99	TH2F * fhPhiKaon; //! phi distribution vs pT
	100	TH2F * fhEtaKaon; //! eta distribution vs pT
477d6cee	101
05d0d05d	102	TH1F * fhPtUnknown; //! pT distribution
	103	TH2F * fhPhiUnknown; //! phi distribution vs pT
	104	TH2F * fhEtaUnknown; //! eta distribution vs pT
727a309a	105
	106	TH1F * fhMCPt; //! Number of generated charged hadrons vs pT coming from MC particle
	107	TH2F * fhMCPhi; //! Number of generated charged hadrons vs phi coming from MC particle
	108	TH2F * fhMCEta; //! Number of generated charged hadrons vs eta coming from MC particle
	109	TH1F * fhMCRecPt; //! Number of reconstructed charged hadrons vs pT coming from MC particle
	110
fe9615b9	111	// TOF
fedea415	112	TH1F * fhTOFSignal; //! TOF signal
	113	TH1F * fhTOFSignalPtCut; //! TOF signal pt and acceptance cut
	114	TH1F * fhTOFSignalBCOK; //! TOF signal pt and acceptance cut
	115	TH2F * fhPtTOFSignal; //! TOF signal vs track pT, good status
cc944149	116	TH2F * fhPtTOFSignalDCACut; //! TOF signal vs track pT, good status
fedea415	117	TH2F * fhPtTOFSignalPileUp[7]; //! TOF signal vs track pT, good status, pile-up
52bd777f	118	TH2F * fhPtTOFSignalVtxOutBC0; //! TOF signal vs track pT, good status
	119	TH2F * fhPtTOFSignalVtxOutBC0PileUp[7];//! TOF signal vs track pT, good status, pile-up
	120	TH2F * fhPtTOFSignalVtxInBC0; //! TOF signal vs track pT, good status
	121	TH2F * fhPtTOFSignalVtxInBC0PileUp[7]; //! TOF signal vs track pT, good status, pile-up
fedea415	122	TH1F * fhPtTOFStatus0; //! pT of tracks not passing TOF status selection
	123	TH2F * fhEtaPhiTOFStatus0; //! eta/phi of tracks not passing TOF status selection
	124	TH2F * fhEtaPhiTOFBC0; //! eta/phi of tracks passing TOF status selection, tracks in BC=0
	125	TH2F * fhEtaPhiTOFBCPlus; //! eta/phi of tracks passing TOF status selection, tracks in BC>0
	126	TH2F * fhEtaPhiTOFBCMinus; //! eta/phi of tracks passing TOF status selection, tracks in BC<0
	127	TH2F * fhEtaPhiTOFBC0PileUpSPD; //! eta/phi of tracks passing TOF status selection, tracks in BC=0, pile-up spd
	128	TH2F * fhEtaPhiTOFBCPlusPileUpSPD; //! eta/phi of tracks passing TOF status selection, tracks in BC>0, pile-up spd
	129	TH2F * fhEtaPhiTOFBCMinusPileUpSPD; //! eta/phi of tracks passing TOF status selection, tracks in BC<0, pile-up spd
52bd777f	130
	131	TH1F * fhProductionVertexBC; //! check BC of production vertex
	132	TH1F * fhProductionVertexBCPileUp[7]; //! check BC of production vertex, pile-up
fe9615b9	133
fedea415	134	TH2F * fhPtDCA[3]; //! DCA (xy,z,constrained) of all tracks
8f467228	135
	136	TH2F * fhPtDCASPDRefit[3]; //! DCA (xy,z,constrained) of tracks with SPD and ITS refit
	137	TH2F * fhPtDCANoSPDRefit[3]; //! DCA (xy,z,constrained) of constrained tracks no SPD and with ITSRefit
	138	TH2F * fhPtDCANoSPDNoRefit[3]; //! DCA (xy,z,constrained) of constrained tracks with no SPD requierement and without ITSRefit
	139
52bd777f	140	TH2F * fhPtDCAVtxOutBC0[3]; //! DCA (xy,z,constrained) of all tracks, vertex BC!=0
52bd777f	141	TH2F * fhPtDCAVtxInBC0[3]; //! DCA (xy,z,constrained) of all tracks, vertex BC==0
fedea415	142	TH2F * fhPtDCAPileUp[3]; //! DCA (xy,z,constrained) of all tracks, SPD pile-up
52bd777f	143	TH2F * fhPtDCAVtxOutBC0PileUp[3]; //! DCA (xy,z,constrained) of all tracks, vertex BC!=0, SPD pile-up
52bd777f	144	TH2F * fhPtDCAVtxInBC0PileUp[3]; //! DCA (xy,z,constrained) of all tracks, vertex BC==0, SPD pile-up
fedea415	145
	146	TH2F * fhPtDCATOFBC0[3]; //! DCA (xy,z,constrained) of all tracks, hit in TOF and BC=0
	147	TH2F * fhPtDCAPileUpTOFBC0[3]; //! DCA (xy,z,constrained) of all tracks, hit in TOF and BC=0
52bd777f	148	TH2F * fhPtDCATOFBCOut[3]; //! DCA (xy,z,constrained) of all tracks, hit in TOF and BC!=0
fedea415	149
fedea415	150	TH2F * fhPtDCANoTOFHit[3]; //! DCA (xy,z,constrained) of all tracks, no hit in TOF
52bd777f	151	TH2F * fhPtDCAVtxOutBC0NoTOFHit[3]; //! DCA (xy,z,constrained) of all tracks, vertex BC!=0, no hit in TOF
52bd777f	152	TH2F * fhPtDCAVtxInBC0NoTOFHit[3]; //! DCA (xy,z,constrained) of all tracks, vertex BC=0, no hit in TOF
fedea415	153	TH2F * fhPtDCAPileUpNoTOFHit[3]; //! DCA (xy,z,constrained) of all tracks, SPD pile-up, no hit in TOF
52bd777f	154	TH2F * fhPtDCAVtxOutBC0PileUpNoTOFHit[3]; //! DCA (xy,z,constrained) of all tracks, vertex BC!=0, SPD pile-up, no hit in TOF
52bd777f	155	TH2F * fhPtDCAVtxInBC0PileUpNoTOFHit[3]; //! DCA (xy,z,constrained) of all tracks, vertex BC=0, SPD pile-up, no hit in TOF
fedea415	156
fedea415	157
34c16486	158	AliAnaChargedParticles( const AliAnaChargedParticles & ch) ; // cpy ctor
34c16486	159	AliAnaChargedParticles & operator = (const AliAnaChargedParticles & ch) ; // cpy assignment
c5693f62	160
9b4a1ee2	161	ClassDef(AliAnaChargedParticles,8)
05d0d05d	162
05d0d05d	163	} ;
477d6cee	164
	165
	166	#endif //ALIANACHARGEDPARTICLES_H
	167
	168
	169