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