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

#ifndef ALITPCCALIBTIMEGAIN_H
#define ALITPCCALIBTIMEGAIN_H

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

#include "AliTPCcalibBase.h"
#include "TH2F.h"
#include "TF1.h"
#include "TArrayD.h"
#include "TObjArray.h"
#include "AliSplineFit.h"

class TH1F;
class TH3F;
class TH2F;
class THnSparse;
class TList;
class TGraphErrors;
class AliESDEvent;
class AliESDtrack;
class AliTPCcalibLaser;
class AliTPCseed;

#include "TTreeStream.h"


class AliTPCcalibTimeGain:public AliTPCcalibBase {
public:
  AliTPCcalibTimeGain(); 
  AliTPCcalibTimeGain(const Text_t *name, const Text_t *title, UInt_t StartTime, UInt_t EndTime, Int_t deltaIntegrationTimeGain);
  virtual ~AliTPCcalibTimeGain();
  //
  virtual void           Process(AliESDEvent *event);
  virtual Long64_t       Merge(TCollection *li);
  virtual void           AnalyzeRun(Int_t minEntries);
  //
  void                   ProcessCosmicEvent(AliESDEvent *event);
  void                   ProcessBeamEvent(AliESDEvent *event);
  //
  void                   CalculateBetheAlephParams(TH2F *hist, Double_t * ini);
  static void            BinLogX(THnSparse *h, Int_t axisDim);
  static void            BinLogX(TH1 *h);
  //
  THnSparse *            GetHistGainTime() const {return (THnSparse*) fHistGainTime;};
  TH2F      *            GetHistDeDxTotal() const {return (TH2F*) fHistDeDxTotal;};
  //
  TGraphErrors *         GetGraphGainVsTime(Int_t runNumber = 0, Int_t minEntries = 2000);
  static AliSplineFit *  MakeSplineFit(TGraphErrors * graph);
  TGraphErrors *         GetGraphAttachment(Int_t minEntries, Int_t nmaxBin, Float_t fracLow=0.1, Float_t fracUp=0.9);
  //
  void SetMIP(Float_t MIP){fMIP = MIP;};
  void SetUseMax(Bool_t UseMax){fUseMax = UseMax;};
  void SetLowerTrunc(Float_t LowerTrunc){fLowerTrunc = LowerTrunc;};
  void SetUpperTrunc(Float_t UpperTrunc){fUpperTrunc = UpperTrunc;};
  void SetUseShapeNorm(Bool_t UseShapeNorm){fUseShapeNorm = UseShapeNorm;};
  void SetUsePosNorm(Bool_t UsePosNorm){fUsePosNorm = UsePosNorm;};
  void SetUsePadNorm(Int_t UsePadNorm){fUsePadNorm = UsePadNorm;};
  void SetIsCosmic(Bool_t IsCosmic){fIsCosmic = IsCosmic;};
  void SetLowMemoryConsumption(Bool_t LowMemoryConsumption){fLowMemoryConsumption = LowMemoryConsumption;};
  void SetUseCookAnalytical(Bool_t UseCookAnalytical){fUseCookAnalytical = UseCookAnalytical;};

private:
  //
  Float_t GetTPCdEdx(AliTPCseed * seed);   // wrapper for CookdEdxNorm or analytical
  //
  THnSparse    * fHistGainTime;            // dEdx vs. time, type, Driftlength, momentum P
  TGraphErrors * fGainVsTime;              // multiplication factor vs. time
  TH2F         * fHistDeDxTotal;           // dEdx vs. momentum for quality assurance
  //
  Float_t fIntegrationTimeDeDx;         // required statistics for each dEdx time bin
  //
  Float_t fMIP;                         // rough MIP position in order to have scaleable histograms
  //
  Bool_t  fUseMax;                      // true: use max charge for dE/dx calculation, false: use total charge for dE/dx calculation
  Float_t fLowerTrunc;                  // lower truncation of dE/dx ; at most 5%
  Float_t fUpperTrunc;                  // upper truncation of dE/dx ; ca. 70%
  Bool_t  fUseShapeNorm;                // use empirical correction of dependencies
  Bool_t  fUsePosNorm;                  // charge correction (analytical?)
  Int_t   fUsePadNorm;                  // normalization of pad geometries
  Bool_t  fUseCookAnalytical;           // true if CookdEdxAnalytical should be used
  //
  Bool_t  fIsCosmic;                    // kTRUE if the analyzed runs contain cosmic events
  Bool_t  fLowMemoryConsumption;        // set this option kTRUE if the momenta information should not be stored in order to save memory
  //
  AliTPCcalibTimeGain(const AliTPCcalibTimeGain&); 
  AliTPCcalibTimeGain& operator=(const AliTPCcalibTimeGain&); 
  void     Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
  void     Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}

  ClassDef(AliTPCcalibTimeGain, 1); 
};

#endif
Commit	Line	Data
	1	#ifndef ALITPCCALIBTIMEGAIN_H
	2	#define ALITPCCALIBTIMEGAIN_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 "TH2F.h"
	9	#include "TF1.h"
	10	#include "TArrayD.h"
	11	#include "TObjArray.h"
	12	#include "AliSplineFit.h"
	13
	14	class TH1F;
	15	class TH3F;
	16	class TH2F;
	17	class THnSparse;
	18	class TList;
	19	class TGraphErrors;
	20	class AliESDEvent;
	21	class AliESDtrack;
	22	class AliTPCcalibLaser;
	23	class AliTPCseed;
	24
	25	#include "TTreeStream.h"
	26
	27
	28	class AliTPCcalibTimeGain:public AliTPCcalibBase {
	29	public:
	30	AliTPCcalibTimeGain();
	31	AliTPCcalibTimeGain(const Text_t name, const Text_t title, UInt_t StartTime, UInt_t EndTime, Int_t deltaIntegrationTimeGain);
	32	virtual ~AliTPCcalibTimeGain();
	33	//
	34	virtual void Process(AliESDEvent *event);
	35	virtual Long64_t Merge(TCollection *li);
	36	virtual void AnalyzeRun(Int_t minEntries);
	37	//
	38	void ProcessCosmicEvent(AliESDEvent *event);
	39	void ProcessBeamEvent(AliESDEvent *event);
	40	//
	41	void CalculateBetheAlephParams(TH2F hist, Double_t ini);
	42	static void BinLogX(THnSparse *h, Int_t axisDim);
	43	static void BinLogX(TH1 *h);
	44	//
	45	THnSparse * GetHistGainTime() const {return (THnSparse*) fHistGainTime;};
	46	TH2F * GetHistDeDxTotal() const {return (TH2F*) fHistDeDxTotal;};
	47	//
	48	TGraphErrors * GetGraphGainVsTime(Int_t runNumber = 0, Int_t minEntries = 2000);
	49	static AliSplineFit * MakeSplineFit(TGraphErrors * graph);
	50	TGraphErrors * GetGraphAttachment(Int_t minEntries, Int_t nmaxBin, Float_t fracLow=0.1, Float_t fracUp=0.9);
	51	//
	52	void SetMIP(Float_t MIP){fMIP = MIP;};
	53	void SetUseMax(Bool_t UseMax){fUseMax = UseMax;};
	54	void SetLowerTrunc(Float_t LowerTrunc){fLowerTrunc = LowerTrunc;};
	55	void SetUpperTrunc(Float_t UpperTrunc){fUpperTrunc = UpperTrunc;};
	56	void SetUseShapeNorm(Bool_t UseShapeNorm){fUseShapeNorm = UseShapeNorm;};
	57	void SetUsePosNorm(Bool_t UsePosNorm){fUsePosNorm = UsePosNorm;};
	58	void SetUsePadNorm(Int_t UsePadNorm){fUsePadNorm = UsePadNorm;};
	59	void SetIsCosmic(Bool_t IsCosmic){fIsCosmic = IsCosmic;};
	60	void SetLowMemoryConsumption(Bool_t LowMemoryConsumption){fLowMemoryConsumption = LowMemoryConsumption;};
	61	void SetUseCookAnalytical(Bool_t UseCookAnalytical){fUseCookAnalytical = UseCookAnalytical;};
	62
	63	private:
	64	//
	65	Float_t GetTPCdEdx(AliTPCseed * seed); // wrapper for CookdEdxNorm or analytical
	66	//
	67	THnSparse * fHistGainTime; // dEdx vs. time, type, Driftlength, momentum P
	68	TGraphErrors * fGainVsTime; // multiplication factor vs. time
	69	TH2F * fHistDeDxTotal; // dEdx vs. momentum for quality assurance
	70	//
	71	Float_t fIntegrationTimeDeDx; // required statistics for each dEdx time bin
	72	//
	73	Float_t fMIP; // rough MIP position in order to have scaleable histograms
	74	//
	75	Bool_t fUseMax; // true: use max charge for dE/dx calculation, false: use total charge for dE/dx calculation
	76	Float_t fLowerTrunc; // lower truncation of dE/dx ; at most 5%
	77	Float_t fUpperTrunc; // upper truncation of dE/dx ; ca. 70%
	78	Bool_t fUseShapeNorm; // use empirical correction of dependencies
	79	Bool_t fUsePosNorm; // charge correction (analytical?)
	80	Int_t fUsePadNorm; // normalization of pad geometries
	81	Bool_t fUseCookAnalytical; // true if CookdEdxAnalytical should be used
	82	//
	83	Bool_t fIsCosmic; // kTRUE if the analyzed runs contain cosmic events
	84	Bool_t fLowMemoryConsumption; // set this option kTRUE if the momenta information should not be stored in order to save memory
	85	//
	86	AliTPCcalibTimeGain(const AliTPCcalibTimeGain&);
	87	AliTPCcalibTimeGain& operator=(const AliTPCcalibTimeGain&);
	88	void Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
	89	void Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
	90
	91	ClassDef(AliTPCcalibTimeGain, 1);
	92	};
	93
	94	#endif
	95
	96