[u/mrichter/AliRoot.git] / TPC / AliTPCcalibGainMult.h

#ifndef ALITPCCALIBGAINMULT_H
#define ALITPCCALIBGAINMULT_H

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

#include "AliTPCcalibBase.h"
#include "AliTPCCalPad.h"
#include "TH3F.h"
#include "TF1.h"
#include "THnSparse.h"
#include "TMatrixD.h"
#include "TVectorD.h"
class TH1F;
class TList;
class AliESDEvent;
class AliESDtrack;
class AliTPCseed;

#include "TTreeStream.h"


class AliTPCcalibGainMult:public AliTPCcalibBase {
public:
  AliTPCcalibGainMult(); 
  AliTPCcalibGainMult(const Text_t *name, const Text_t *title);
  virtual ~AliTPCcalibGainMult();
  void SetBBParam(TVectorD * param) {fBBParam=param;}
  //  virtual void Terminate();  
  //
  virtual void           Process(AliESDEvent *event);
  virtual void           ProcessV0s(AliESDEvent *event);
  virtual void           ProcessCosmic(const AliESDEvent *event);
  virtual void           ProcessKinks(const AliESDEvent *event);
  virtual void           ProcessTOF(const AliESDEvent *event);  
  virtual void           DumpHPT(const AliESDEvent *event);
  virtual Long64_t       Merge(TCollection *li);
  virtual void           Analyze();
  void                   DumpTrack(AliESDtrack * track, AliESDfriendTrack *ftrack, AliTPCseed * seed, Int_t index);
  //
  TH1F   *          GetHistNTracks() const {return fHistNTracks;};
  TH1F   *          GetHistClusterShape() const {return fHistClusterShape;};
  TH3F   *          GetHistQA() const {return fHistQA;};
  //
  THnSparseF *      GetHistGainSector() const {return fHistGainSector;};
  THnSparseF *      GetHistPadEqual() const {return fHistPadEqual;};
  THnSparseF *      GetHistGainMult() const {return fHistGainMult;};  
  THnSparseF * GetHistdEdxMap() const { return fHistdEdxMap;}      // 4D dedx histogram
  THnSparseF * GetHistdEdxMax() const { return fHistdEdxMax;}      // 4D dedx histogram
  THnSparseF * GetHistdEdxTot() const { return fHistdEdxTot;}      // 4D dedx histogram
  TTree *      GetdEdxTree() const {return fdEdxTree;}         // tree for the later minimization

  //
  void SetMIPvalue(Float_t mip){fMIP = mip;};
  void SetLowerTrunc(Float_t lowerTrunc){fLowerTrunc = lowerTrunc;};
  void SetUpperTrunc(Float_t upperTrunc){fUpperTrunc = upperTrunc;};
  void SetUseMax(Bool_t useMax){fUseMax = useMax;};
  //
  //
  void     Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
  void     Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
  //
  void     MakeLookup(THnSparse * hist, Char_t * outputFile);
  //
  void     UpdateGainMap();
  void     UpdateClusterParam();


private:
  //
  // parameter specifications
  //
  Float_t fMIP;                  // MIP position to be in fMIP
  Float_t fLowerTrunc;           // lower truncation for dEdx
  Float_t fUpperTrunc;           // upper truncation for dEdx
  //
  Bool_t fUseMax;                 // flag if Qmax or Qtot should be used
  //
  // histograms
  //
  TH1F  *fHistNTracks;            //  histogram showing number of ESD tracks per event
  TH1F  *fHistClusterShape;       //  histogram to check the cluster shape
  TH3F  *fHistQA;                 //  dE/dx histogram showing the final spectrum
  //
  THnSparseF * fHistGainSector;   //  histogram which shows MIP peak for each of the 3x36 sectors (pad region)
  THnSparseF * fHistPadEqual;     //  histogram for the equalization of the gain in the different pad regions -> pass0
  THnSparseF * fHistGainMult;     //  histogram which shows decrease of MIP signal as a function
  TMatrixD *fPIDMatrix;           //! custom PID matrix
  //
  THnSparseF * fHistdEdxMap;      // 4D dedx histogram - per sector/phi
  THnSparseF * fHistdEdxMax;      // 5D dedx histogram - per 1/dedx, tan(theta), tan(phi), mult, pad-type
  THnSparseF * fHistdEdxTot;      // 5D dedx histogram - per 1/dedx, tan(theta), tan(phi), mult, pad-type
  TTree *      fdEdxTree;         // tree for the later minimization
  TVectorD    *fBBParam;          // BetheBloch parameterization used for the dedx expected calculation
  //
  AliTPCcalibGainMult(const AliTPCcalibGainMult&); 
  AliTPCcalibGainMult& operator=(const AliTPCcalibGainMult&); 

  ClassDef(AliTPCcalibGainMult, 2); 
};

#endif
Commit	Line	Data
f72219cb	1	#ifndef ALITPCCALIBGAINMULT_H
	2	#define ALITPCCALIBGAINMULT_H
	3
	4	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	5	* See cxx source for full Copyright notice */
	6
	7	#include "AliTPCcalibBase.h"
	8	#include "AliTPCCalPad.h"
	9	#include "TH3F.h"
	10	#include "TF1.h"
	11	#include "THnSparse.h"
6feb400f	12	#include "TMatrixD.h"
6feb400f	13	#include "TVectorD.h"
f72219cb	14	class TH1F;
	15	class TList;
	16	class AliESDEvent;
	17	class AliESDtrack;
	18	class AliTPCseed;
	19
	20	#include "TTreeStream.h"
	21
	22
	23	class AliTPCcalibGainMult:public AliTPCcalibBase {
	24	public:
	25	AliTPCcalibGainMult();
	26	AliTPCcalibGainMult(const Text_t name, const Text_t title);
	27	virtual ~AliTPCcalibGainMult();
6feb400f	28	void SetBBParam(TVectorD * param) {fBBParam=param;}
6feb400f	29	// virtual void Terminate();
f72219cb	30	//
f72219cb	31	virtual void Process(AliESDEvent *event);
6feb400f	32	virtual void ProcessV0s(AliESDEvent *event);
	33	virtual void ProcessCosmic(const AliESDEvent *event);
	34	virtual void ProcessKinks(const AliESDEvent *event);
	35	virtual void ProcessTOF(const AliESDEvent *event);
	36	virtual void DumpHPT(const AliESDEvent *event);
f72219cb	37	virtual Long64_t Merge(TCollection *li);
f72219cb	38	virtual void Analyze();
6feb400f	39	void DumpTrack(AliESDtrack * track, AliESDfriendTrack ftrack, AliTPCseed seed, Int_t index);
f72219cb	40	//
	41	TH1F * GetHistNTracks() const {return fHistNTracks;};
	42	TH1F * GetHistClusterShape() const {return fHistClusterShape;};
	43	TH3F * GetHistQA() const {return fHistQA;};
	44	//
	45	THnSparseF * GetHistGainSector() const {return fHistGainSector;};
	46	THnSparseF * GetHistPadEqual() const {return fHistPadEqual;};
6feb400f	47	THnSparseF * GetHistGainMult() const {return fHistGainMult;};
	48	THnSparseF * GetHistdEdxMap() const { return fHistdEdxMap;} // 4D dedx histogram
	49	THnSparseF * GetHistdEdxMax() const { return fHistdEdxMax;} // 4D dedx histogram
	50	THnSparseF * GetHistdEdxTot() const { return fHistdEdxTot;} // 4D dedx histogram
	51	TTree * GetdEdxTree() const {return fdEdxTree;} // tree for the later minimization
	52
f72219cb	53	//
	54	void SetMIPvalue(Float_t mip){fMIP = mip;};
	55	void SetLowerTrunc(Float_t lowerTrunc){fLowerTrunc = lowerTrunc;};
	56	void SetUpperTrunc(Float_t upperTrunc){fUpperTrunc = upperTrunc;};
	57	void SetUseMax(Bool_t useMax){fUseMax = useMax;};
	58	//
	59	//
f72219cb	60	void Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
	61	void Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
	62	//
19f63078	63	void MakeLookup(THnSparse * hist, Char_t * outputFile);
19f63078	64	//
f72219cb	65	void UpdateGainMap();
	66	void UpdateClusterParam();
	67
	68
	69	private:
	70	//
	71	// parameter specifications
	72	//
	73	Float_t fMIP; // MIP position to be in fMIP
	74	Float_t fLowerTrunc; // lower truncation for dEdx
	75	Float_t fUpperTrunc; // upper truncation for dEdx
	76	//
	77	Bool_t fUseMax; // flag if Qmax or Qtot should be used
	78	//
	79	// histograms
	80	//
	81	TH1F *fHistNTracks; // histogram showing number of ESD tracks per event
	82	TH1F *fHistClusterShape; // histogram to check the cluster shape
	83	TH3F *fHistQA; // dE/dx histogram showing the final spectrum
	84	//
	85	THnSparseF * fHistGainSector; // histogram which shows MIP peak for each of the 3x36 sectors (pad region)
	86	THnSparseF * fHistPadEqual; // histogram for the equalization of the gain in the different pad regions -> pass0
	87	THnSparseF * fHistGainMult; // histogram which shows decrease of MIP signal as a function
6feb400f	88	TMatrixD *fPIDMatrix; //! custom PID matrix
	89	//
	90	THnSparseF * fHistdEdxMap; // 4D dedx histogram - per sector/phi
	91	THnSparseF * fHistdEdxMax; // 5D dedx histogram - per 1/dedx, tan(theta), tan(phi), mult, pad-type
	92	THnSparseF * fHistdEdxTot; // 5D dedx histogram - per 1/dedx, tan(theta), tan(phi), mult, pad-type
	93	TTree * fdEdxTree; // tree for the later minimization
	94	TVectorD *fBBParam; // BetheBloch parameterization used for the dedx expected calculation
f72219cb	95	//
	96	AliTPCcalibGainMult(const AliTPCcalibGainMult&);
	97	AliTPCcalibGainMult& operator=(const AliTPCcalibGainMult&);
	98
6feb400f	99	ClassDef(AliTPCcalibGainMult, 2);
f72219cb	100	};
	101
	102	#endif
	103
	104