[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); 
  //
  // Charge parameterization
  //
  Float_t Qnorm(Int_t ipad, Int_t itype, Float_t dr, Float_t ty, Float_t tz);   

  Float_t QnormPos(Int_t ipad, Bool_t isMax,  Float_t pad, Float_t time, Float_t z, Float_t sy2, Float_t sz2, Float_t qm, Float_t qt);

  //
  // 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);}


 public: 
  //
  //
  //
  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
  //
  // charge normalization parametrization
  //
  TObjArray *fQNorm;              // q norm paramters
  TVectorD  *fPosQTnorm[3];        // q position normalization
  TVectorD  *fPosQMnorm[3];        // q position normalization
  TVectorD  *fQpadTnorm;           // q pad normalization - Total charge
  TVectorD  *fQpadMnorm;           // q pad normalization - Max charge
  //
 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);
11750d28	37	//
	38	// Charge parameterization
	39	//
dcf3a564	40	Float_t Qnorm(Int_t ipad, Int_t itype, Float_t dr, Float_t ty, Float_t tz);
b17540e4	41
	42	Float_t QnormPos(Int_t ipad, Bool_t isMax, Float_t pad, Float_t time, Float_t z, Float_t sy2, Float_t sz2, Float_t qm, Float_t qt);
	43
dcf3a564	44	//
	45	// Error parameterization
	46	//
12ca5da1	47	Float_t GetError0(Int_t dim, Int_t type, Float_t z, Float_t angle);
	48	Float_t GetError0Par(Int_t dim, Int_t type, Float_t z, Float_t angle);
	49	Float_t GetError1(Int_t dim, Int_t type, Float_t z, Float_t angle);
	50	Float_t GetErrorQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
	51	Float_t GetErrorQPar(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
	52	Float_t GetErrorQParScaled(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
dcf3a564	53	//
	54	// Shape parameterization
	55	//
12ca5da1	56	Float_t GetRMS0(Int_t dim, Int_t type, Float_t z, Float_t angle);
	57	Float_t GetRMS1(Int_t dim, Int_t type, Float_t z, Float_t angle);
	58	Float_t GetRMSQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
	59	Float_t GetRMSSigma(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
	60	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	61	//
	62	//
	63	//
12ca5da1	64	void Test(TTree * tree, const char *output="TestClusterParam.root");
	65	//
	66	// static methods equivalents - use instance of param object - useful for tree draw and TF2 visualization
	67	static Float_t SGetError0(Int_t dim, Int_t type, Float_t z, Float_t angle){
	68	return fgInstance->GetError0(dim,type,z,angle);
	69	}
	70	static Float_t SGetError0Par(Int_t dim, Int_t type, Float_t z, Float_t angle){
	71	return fgInstance->GetError0Par(dim,type,z,angle);
	72	}
	73	static Float_t SGetError1(Int_t dim, Int_t type, Float_t z, Float_t angle){
	74	return fgInstance->GetError1(dim,type,z,angle);
	75	}
	76	static Float_t SGetErrorQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
	77	return fgInstance->GetErrorQ(dim,type,z,angle,Qmean);
	78	}
	79	static Float_t SGetErrorQPar(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
	80	return fgInstance->GetErrorQPar(dim,type,z,angle,Qmean);
	81	}
	82	static Float_t SGetErrorQParScaled(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
	83	return fgInstance->GetErrorQParScaled(dim,type,z,angle,Qmean);
	84	}
	85
	86	static Float_t SGetRMS0(Int_t dim, Int_t type, Float_t z, Float_t angle){
	87	return fgInstance->GetRMS0(dim,type,z,angle);
	88	}
	89	static Float_t SGetRMS1(Int_t dim, Int_t type, Float_t z, Float_t angle){
	90	return fgInstance->GetRMS1(dim,type,z,angle);
	91	}
	92	static Float_t SGetRMSQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
	93	return fgInstance->GetRMSQ(dim,type,z,angle,Qmean);
	94	}
	95	static Float_t SGetRMSSigma(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
	96	return fgInstance->GetRMSSigma(dim,type,z,angle,Qmean);
	97	}
	98	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){
	99	return fgInstance->GetShapeFactor(dim,type,z,angle,Qmean, rmsL, rmsM);
	100	}
	101	//
	102	//
0a65832b	103	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	104
0a65832b	105
236a0d03	106	public:
dcf3a564	107	//
	108	//
	109	//
12ca5da1	110	void FitResol0(TTree * tree, Int_t dim, Int_t type, Float_t param0, Float_t error);
	111	void FitResol0Par(TTree * tree, Int_t dim, Int_t type, Float_t param0, Float_t error);
	112	void FitResol1(TTree * tree, Int_t dim, Float_t param0, Float_t error);
	113	void FitResolQ(TTree * tree, Int_t dim, Int_t type, Float_t param0, Float_t error);
	114	void FitResolQPar(TTree * tree, Int_t dim, Int_t type, Float_t param0, Float_t error);
	115	void FitRMS0(TTree * tree, Int_t dim, Int_t type, Float_t param0, Float_t error);
	116	void FitRMS1(TTree * tree, Int_t dim, Float_t param0, Float_t error);
	117	void FitRMSQ(TTree * tree, Int_t dim, Int_t type, Float_t param0, Float_t error);
	118	void FitRMSSigma(TTree * tree, Int_t dim, Int_t type, Float_t param0, Float_t error);
	119	//
0a65832b	120
12ca5da1	121	Float_t fRatio; //ratio of values constibution to error
	122	Float_t fParamS0[2][3][4]; //error parameterization coeficients
	123	Float_t fErrorS0[2][3][4]; //error parameterization coeficients
	124	Float_t fParamS0Par[2][3][7]; //error parameterization coeficients
	125	Float_t fErrorS0Par[2][3][7]; //error parameterization coeficients
	126	Float_t fParamSQ[2][3][6]; //error parameterization coeficients
	127	Float_t fErrorSQ[2][3][6]; //error parameterization coeficients
	128	Float_t fParamSQPar[2][3][9]; //error parameterization coeficients
	129	Float_t fErrorSQPar[2][3][9]; //error parameterization coeficients
	130	Float_t fParamS1[2][4]; //error parameterization coeficients
	131	Float_t fErrorS1[2][4]; //error parameterization coeficients
	132	//
	133	Float_t fParamRMS0[2][3][4]; //shape parameterization coeficients
	134	Float_t fErrorRMS0[2][3][4]; //shape parameterization coeficients
	135	Float_t fParamRMSQ[2][3][6]; //shape parameterization coeficients
	136	Float_t fErrorRMSQ[2][3][6]; //shape parameterization coeficients
	137	Float_t fParamRMS1[2][5]; //shape parameterization coeficients
	138	Float_t fErrorRMS1[2][5]; //shape parameterization coeficients
	139	Float_t fErrorRMSSys[2]; // systematic relative error of the parametererization
	140	Float_t fRMSSigmaRatio[2][2]; // mean value of the varation of RMS to RMS
	141	Float_t fRMSSigmaFit[2][3][2]; // mean value of the varation of RMS to RMS
	142	//
b17540e4	143	// charge normalization parametrization
11750d28	144	//
b17540e4	145	TObjArray *fQNorm; // q norm paramters
11750d28	146	TVectorD *fPosQTnorm[3]; // q position normalization
11750d28	147	TVectorD *fPosQMnorm[3]; // q position normalization
b17540e4	148	TVectorD *fQpadTnorm; // q pad normalization - Total charge
	149	TVectorD *fQpadMnorm; // q pad normalization - Max charge
	150	//
12ca5da1	151	protected:
12ca5da1	152	static AliTPCClusterParam* fgInstance; //! Instance of this class (singleton implementation)
0a65832b	153	ClassDef(AliTPCClusterParam,2) // TPC Cluster parameter class
12ca5da1	154	};
	155
	156	#endif