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

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

/* $Id: AliTPCClusterParam.h,v */

////////////////////////////////////////////////////
//                                                //
//  TPC cluster error and shape parameterization  //
//                                                //
////////////////////////////////////////////////////


#include <TObject.h>
#include <TVectorD.h>

class TTree;
class TObjArray;
//_____________________________________________________________________________
class AliTPCClusterParam : public TObject {
 public:
  static AliTPCClusterParam* Instance();
  AliTPCClusterParam();
  AliTPCClusterParam(const AliTPCClusterParam& param);
  AliTPCClusterParam & operator=(const AliTPCClusterParam&);
  virtual           ~AliTPCClusterParam();
  virtual void	Print(Option_t* option = "") const;
  void SetInstance(AliTPCClusterParam*param){fgInstance = param;}
  //
  // Seting functions
  //
  void FitData(TTree * tree);
  void FitResol(TTree * tree);
  void FitRMS(TTree * tree);
  void SetQnorm(Int_t ipad, Int_t itype,  TVectorD * norm); 
  Float_t Qnorm(Int_t ipad, Int_t itype, Float_t dr, Float_t ty, Float_t tz);   
  //
  // Error parameterization
  //
  Float_t GetError0(Int_t dim, Int_t type, Float_t z, Float_t angle);
  Float_t GetError0Par(Int_t dim, Int_t type, Float_t z, Float_t angle);
  Float_t GetError1(Int_t dim, Int_t type, Float_t z, Float_t angle);
  Float_t GetErrorQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
  Float_t GetErrorQPar(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
  Float_t GetErrorQParScaled(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
  //
  // Shape parameterization
  //
  Float_t GetRMS0(Int_t dim, Int_t type, Float_t z, Float_t angle);
  Float_t GetRMS1(Int_t dim, Int_t type, Float_t z, Float_t angle);
  Float_t GetRMSQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
  Float_t GetRMSSigma(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
  Float_t GetShapeFactor(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean, Float_t rmsL, Float_t rmsM);
  //
  //
  //
  void Test(TTree * tree, const char *output="TestClusterParam.root");
  //
  // static methods equivalents  - use instance of param object - useful for tree draw and TF2 visualization 
  static Float_t SGetError0(Int_t dim, Int_t type, Float_t z, Float_t angle){
    return fgInstance->GetError0(dim,type,z,angle);
  }
  static Float_t SGetError0Par(Int_t dim, Int_t type, Float_t z, Float_t angle){
    return fgInstance->GetError0Par(dim,type,z,angle);
  }
  static Float_t SGetError1(Int_t dim, Int_t type, Float_t z, Float_t angle){
    return fgInstance->GetError1(dim,type,z,angle);
  }
  static Float_t SGetErrorQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
    return fgInstance->GetErrorQ(dim,type,z,angle,Qmean);
  }
  static Float_t SGetErrorQPar(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
    return fgInstance->GetErrorQPar(dim,type,z,angle,Qmean);
  }
  static Float_t SGetErrorQParScaled(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
    return fgInstance->GetErrorQParScaled(dim,type,z,angle,Qmean);
  }

  static Float_t SGetRMS0(Int_t dim, Int_t type, Float_t z, Float_t angle){
    return fgInstance->GetRMS0(dim,type,z,angle);
  }
  static Float_t SGetRMS1(Int_t dim, Int_t type, Float_t z, Float_t angle){
    return fgInstance->GetRMS1(dim,type,z,angle);
  }
  static Float_t SGetRMSQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
    return fgInstance->GetRMSQ(dim,type,z,angle,Qmean);
  }
  static Float_t SGetRMSSigma(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
    return fgInstance->GetRMSSigma(dim,type,z,angle,Qmean);
  }
  static Float_t SGetShapeFactor(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean, Float_t rmsL, Float_t rmsM){
    return fgInstance->GetShapeFactor(dim,type,z,angle,Qmean, rmsL, rmsM);
  }
  //
  //
  static Float_t SQnorm(Int_t ipad, Int_t itype,Float_t dr, Float_t ty, Float_t tz) {return fgInstance->Qnorm(ipad, itype, dr,ty,tz);}


 protected: 
  //
  //
  //
  void FitResol0(TTree * tree, Int_t dim, Int_t type, Float_t *param0, Float_t *error);
  void FitResol0Par(TTree * tree, Int_t dim, Int_t type, Float_t *param0, Float_t *error);
  void FitResol1(TTree * tree, Int_t dim, Float_t *param0, Float_t *error);
  void FitResolQ(TTree * tree, Int_t dim, Int_t type, Float_t *param0, Float_t *error);
  void FitResolQPar(TTree * tree, Int_t dim, Int_t type, Float_t *param0, Float_t *error);
  void FitRMS0(TTree * tree, Int_t dim, Int_t type, Float_t *param0, Float_t *error);
  void FitRMS1(TTree * tree, Int_t dim, Float_t *param0, Float_t *error);
  void FitRMSQ(TTree * tree, Int_t dim, Int_t type, Float_t *param0, Float_t *error);  
  void FitRMSSigma(TTree * tree, Int_t dim, Int_t type, Float_t *param0, Float_t *error);  
  //

  Float_t fRatio;               //ratio of values constibution to error
  Float_t fParamS0[2][3][4];    //error parameterization coeficients
  Float_t fErrorS0[2][3][4];    //error parameterization coeficients
  Float_t fParamS0Par[2][3][7];    //error parameterization coeficients
  Float_t fErrorS0Par[2][3][7];    //error parameterization coeficients  
  Float_t fParamSQ[2][3][6];    //error parameterization coeficients
  Float_t fErrorSQ[2][3][6];    //error parameterization coeficients
  Float_t fParamSQPar[2][3][9];    //error parameterization coeficients
  Float_t fErrorSQPar[2][3][9];    //error parameterization coeficients
  Float_t fParamS1[2][4];       //error parameterization coeficients
  Float_t fErrorS1[2][4];       //error parameterization coeficients
  //
  Float_t fParamRMS0[2][3][4];   //shape parameterization coeficients
  Float_t fErrorRMS0[2][3][4];   //shape parameterization coeficients
  Float_t fParamRMSQ[2][3][6];   //shape parameterization coeficients
  Float_t fErrorRMSQ[2][3][6];   //shape parameterization coeficients
  Float_t fParamRMS1[2][5];      //shape parameterization coeficients
  Float_t fErrorRMS1[2][5];      //shape parameterization coeficients
  Float_t fErrorRMSSys[2];        // systematic relative error of the parametererization
  Float_t fRMSSigmaRatio[2][2];   // mean value of the varation of RMS to RMS
  Float_t fRMSSigmaFit[2][3][2];   // mean value of the varation of RMS to RMS
  //
  TObjArray *fQNorm;              // q norm paramters
 protected:
  static AliTPCClusterParam*   fgInstance; //! Instance of this class (singleton implementation)
  ClassDef(AliTPCClusterParam,2)    //  TPC Cluster parameter class
};

#endif
Commit	Line	Data
12ca5da1	1	#ifndef ALITPCCLUSTERPARAM_H
	2	#define ALITPCCLUSTERPARAM_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	/* $Id: AliTPCClusterParam.h,v */
	7
	8	////////////////////////////////////////////////////
	9	// //
	10	// TPC cluster error and shape parameterization //
	11	// //
	12	////////////////////////////////////////////////////
	13
	14
	15	#include <TObject.h>
0a65832b	16	#include <TVectorD.h>
12ca5da1	17
0a65832b	18	class TTree;
0a65832b	19	class TObjArray;
12ca5da1	20	//_____________________________________________________________________________
	21	class AliTPCClusterParam : public TObject {
	22	public:
	23	static AliTPCClusterParam* Instance();
f1c2a4a3	24	AliTPCClusterParam();
38caa778	25	AliTPCClusterParam(const AliTPCClusterParam& param);
38caa778	26	AliTPCClusterParam & operator=(const AliTPCClusterParam&);
f1c2a4a3	27	virtual ~AliTPCClusterParam();
12ca5da1	28	virtual void Print(Option_t* option = "") const;
12ca5da1	29	void SetInstance(AliTPCClusterParam*param){fgInstance = param;}
dcf3a564	30	//
	31	// Seting functions
	32	//
12ca5da1	33	void FitData(TTree * tree);
dcf3a564	34	void FitResol(TTree * tree);
	35	void FitRMS(TTree * tree);
	36	void SetQnorm(Int_t ipad, Int_t itype, TVectorD * norm);
	37	Float_t Qnorm(Int_t ipad, Int_t itype, Float_t dr, Float_t ty, Float_t tz);
	38	//
	39	// Error parameterization
	40	//
12ca5da1	41	Float_t GetError0(Int_t dim, Int_t type, Float_t z, Float_t angle);
	42	Float_t GetError0Par(Int_t dim, Int_t type, Float_t z, Float_t angle);
	43	Float_t GetError1(Int_t dim, Int_t type, Float_t z, Float_t angle);
	44	Float_t GetErrorQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
	45	Float_t GetErrorQPar(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
	46	Float_t GetErrorQParScaled(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
dcf3a564	47	//
	48	// Shape parameterization
	49	//
12ca5da1	50	Float_t GetRMS0(Int_t dim, Int_t type, Float_t z, Float_t angle);
	51	Float_t GetRMS1(Int_t dim, Int_t type, Float_t z, Float_t angle);
	52	Float_t GetRMSQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
	53	Float_t GetRMSSigma(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
	54	Float_t GetShapeFactor(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean, Float_t rmsL, Float_t rmsM);
dcf3a564	55	//
	56	//
	57	//
12ca5da1	58	void Test(TTree * tree, const char *output="TestClusterParam.root");
	59	//
	60	// static methods equivalents - use instance of param object - useful for tree draw and TF2 visualization
	61	static Float_t SGetError0(Int_t dim, Int_t type, Float_t z, Float_t angle){
	62	return fgInstance->GetError0(dim,type,z,angle);
	63	}
	64	static Float_t SGetError0Par(Int_t dim, Int_t type, Float_t z, Float_t angle){
	65	return fgInstance->GetError0Par(dim,type,z,angle);
	66	}
	67	static Float_t SGetError1(Int_t dim, Int_t type, Float_t z, Float_t angle){
	68	return fgInstance->GetError1(dim,type,z,angle);
	69	}
	70	static Float_t SGetErrorQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
	71	return fgInstance->GetErrorQ(dim,type,z,angle,Qmean);
	72	}
	73	static Float_t SGetErrorQPar(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
	74	return fgInstance->GetErrorQPar(dim,type,z,angle,Qmean);
	75	}
	76	static Float_t SGetErrorQParScaled(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
	77	return fgInstance->GetErrorQParScaled(dim,type,z,angle,Qmean);
	78	}
	79
	80	static Float_t SGetRMS0(Int_t dim, Int_t type, Float_t z, Float_t angle){
	81	return fgInstance->GetRMS0(dim,type,z,angle);
	82	}
	83	static Float_t SGetRMS1(Int_t dim, Int_t type, Float_t z, Float_t angle){
	84	return fgInstance->GetRMS1(dim,type,z,angle);
	85	}
	86	static Float_t SGetRMSQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
	87	return fgInstance->GetRMSQ(dim,type,z,angle,Qmean);
	88	}
	89	static Float_t SGetRMSSigma(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
	90	return fgInstance->GetRMSSigma(dim,type,z,angle,Qmean);
	91	}
	92	static Float_t SGetShapeFactor(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean, Float_t rmsL, Float_t rmsM){
	93	return fgInstance->GetShapeFactor(dim,type,z,angle,Qmean, rmsL, rmsM);
	94	}
	95	//
	96	//
0a65832b	97	static Float_t SQnorm(Int_t ipad, Int_t itype,Float_t dr, Float_t ty, Float_t tz) {return fgInstance->Qnorm(ipad, itype, dr,ty,tz);}
0a65832b	98
0a65832b	99
12ca5da1	100	protected:
dcf3a564	101	//
	102	//
	103	//
12ca5da1	104	void FitResol0(TTree * tree, Int_t dim, Int_t type, Float_t param0, Float_t error);
	105	void FitResol0Par(TTree * tree, Int_t dim, Int_t type, Float_t param0, Float_t error);
	106	void FitResol1(TTree * tree, Int_t dim, Float_t param0, Float_t error);
	107	void FitResolQ(TTree * tree, Int_t dim, Int_t type, Float_t param0, Float_t error);
	108	void FitResolQPar(TTree * tree, Int_t dim, Int_t type, Float_t param0, Float_t error);
	109	void FitRMS0(TTree * tree, Int_t dim, Int_t type, Float_t param0, Float_t error);
	110	void FitRMS1(TTree * tree, Int_t dim, Float_t param0, Float_t error);
	111	void FitRMSQ(TTree * tree, Int_t dim, Int_t type, Float_t param0, Float_t error);
	112	void FitRMSSigma(TTree * tree, Int_t dim, Int_t type, Float_t param0, Float_t error);
	113	//
0a65832b	114
12ca5da1	115	Float_t fRatio; //ratio of values constibution to error
	116	Float_t fParamS0[2][3][4]; //error parameterization coeficients
	117	Float_t fErrorS0[2][3][4]; //error parameterization coeficients
	118	Float_t fParamS0Par[2][3][7]; //error parameterization coeficients
	119	Float_t fErrorS0Par[2][3][7]; //error parameterization coeficients
	120	Float_t fParamSQ[2][3][6]; //error parameterization coeficients
	121	Float_t fErrorSQ[2][3][6]; //error parameterization coeficients
	122	Float_t fParamSQPar[2][3][9]; //error parameterization coeficients
	123	Float_t fErrorSQPar[2][3][9]; //error parameterization coeficients
	124	Float_t fParamS1[2][4]; //error parameterization coeficients
	125	Float_t fErrorS1[2][4]; //error parameterization coeficients
	126	//
	127	Float_t fParamRMS0[2][3][4]; //shape parameterization coeficients
	128	Float_t fErrorRMS0[2][3][4]; //shape parameterization coeficients
	129	Float_t fParamRMSQ[2][3][6]; //shape parameterization coeficients
	130	Float_t fErrorRMSQ[2][3][6]; //shape parameterization coeficients
	131	Float_t fParamRMS1[2][5]; //shape parameterization coeficients
	132	Float_t fErrorRMS1[2][5]; //shape parameterization coeficients
	133	Float_t fErrorRMSSys[2]; // systematic relative error of the parametererization
	134	Float_t fRMSSigmaRatio[2][2]; // mean value of the varation of RMS to RMS
	135	Float_t fRMSSigmaFit[2][3][2]; // mean value of the varation of RMS to RMS
	136	//
0a65832b	137	TObjArray *fQNorm; // q norm paramters
12ca5da1	138	protected:
12ca5da1	139	static AliTPCClusterParam* fgInstance; //! Instance of this class (singleton implementation)
0a65832b	140	ClassDef(AliTPCClusterParam,2) // TPC Cluster parameter class
12ca5da1	141	};
	142
	143	#endif