[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"
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();
  //
  virtual void           Process(AliESDEvent *event);
  virtual Long64_t       Merge(TCollection *li);
  virtual void           Analyze();
  //
  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;};
  //
  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;};
  //
  //
  static void       BinLogX(const TH1 *h);   // method for correct histogram binning
  void     Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
  void     Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
  //
  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
  //
  AliTPCcalibGainMult(const AliTPCcalibGainMult&); 
  AliTPCcalibGainMult& operator=(const AliTPCcalibGainMult&); 

  ClassDef(AliTPCcalibGainMult, 1); 
};

#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"
	12	class TH1F;
	13	class TList;
	14	class AliESDEvent;
	15	class AliESDtrack;
	16	class AliTPCseed;
	17
	18	#include "TTreeStream.h"
	19
	20
	21	class AliTPCcalibGainMult:public AliTPCcalibBase {
	22	public:
	23	AliTPCcalibGainMult();
	24	AliTPCcalibGainMult(const Text_t name, const Text_t title);
	25	virtual ~AliTPCcalibGainMult();
	26	//
	27	virtual void Process(AliESDEvent *event);
	28	virtual Long64_t Merge(TCollection *li);
	29	virtual void Analyze();
	30	//
	31	TH1F * GetHistNTracks() const {return fHistNTracks;};
	32	TH1F * GetHistClusterShape() const {return fHistClusterShape;};
	33	TH3F * GetHistQA() const {return fHistQA;};
	34	//
	35	THnSparseF * GetHistGainSector() const {return fHistGainSector;};
	36	THnSparseF * GetHistPadEqual() const {return fHistPadEqual;};
	37	THnSparseF * GetHistGainMult() const {return fHistGainMult;};
	38	//
	39	void SetMIPvalue(Float_t mip){fMIP = mip;};
	40	void SetLowerTrunc(Float_t lowerTrunc){fLowerTrunc = lowerTrunc;};
	41	void SetUpperTrunc(Float_t upperTrunc){fUpperTrunc = upperTrunc;};
	42	void SetUseMax(Bool_t useMax){fUseMax = useMax;};
	43	//
	44	//
	45	static void BinLogX(const TH1 *h); // method for correct histogram binning
	46	void Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
	47	void Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
	48	//
	49	void UpdateGainMap();
	50	void UpdateClusterParam();
	51
	52
	53	private:
	54	//
	55	// parameter specifications
	56	//
	57	Float_t fMIP; // MIP position to be in fMIP
	58	Float_t fLowerTrunc; // lower truncation for dEdx
	59	Float_t fUpperTrunc; // upper truncation for dEdx
	60	//
	61	Bool_t fUseMax; // flag if Qmax or Qtot should be used
	62	//
	63	// histograms
	64	//
65	TH1F *fHistNTracks; // histogram showing number of ESD tracks per event
66	TH1F *fHistClusterShape; // histogram to check the cluster shape
67	TH3F *fHistQA; // dE/dx histogram showing the final spectrum
68	//
69	THnSparseF * fHistGainSector; // histogram which shows MIP peak for each of the 3x36 sectors (pad region)
70	THnSparseF * fHistPadEqual; // histogram for the equalization of the gain in the different pad regions -> pass0
71	THnSparseF * fHistGainMult; // histogram which shows decrease of MIP signal as a function
72	//
73	AliTPCcalibGainMult(const AliTPCcalibGainMult&);
74	AliTPCcalibGainMult& operator=(const AliTPCcalibGainMult&);
75
76	ClassDef(AliTPCcalibGainMult, 1);
77	};
78
79	#endif
80
81