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

#ifndef ALITPCKALMANDEDX_H
#define ALITPCKALMANDEDX_H

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

#include "TNamed.h"
#include "TMatrixD.h"
#include "TVectorD.h"
class TTreeSRedirector;
class TObjArray;

class AliTPCkalmandEdx: public TNamed{
public:
  AliTPCkalmandEdx();
  AliTPCkalmandEdx(const char* name, const char* title, Int_t sampleSize=50);
  AliTPCkalmandEdx(const AliTPCkalmandEdx & kalman);
  void UpdatedEdxPair(Int_t ip0, Int_t ip1,
		      Double_t dedx0, Double_t dedx1, 
		      Double_t s0, Double_t s1, 
		      Double_t kY0, Double_t kY1,
		      Double_t kZ0, Double_t kZ1,
		      Double_t dR0, Double_t dR1, 
		      TTreeSRedirector *debug=0);
  void UpdatedEdx(Int_t ip0,
		  Double_t dedx0, 
		  Double_t dedxRef, 
		  Double_t s0,
		  Double_t kY0,
		  Double_t kZ0,
		  Double_t dR0,
		  TTreeSRedirector *debug=0);
  static Int_t GetIndex(Int_t i, Int_t j) { return 15*i+j;}
  //
public:
  void AdddEdx(Int_t ip0,Double_t dedx0, Double_t dedxRef);
  void Init();
  TMatrixD * fState;           // state vector
  TMatrixD * fCovariance;      // covariance
  TMatrixD * fMat1;            // helper unit matrix
  Int_t      fNpad;            // number of pad types
  Int_t      fNpar;            // number of parameters
  Int_t      fNelem;           // number of elements
  //
  // initial parameters estimate
  //
  Int_t      fSampleSize;       // size of  starting  sample
  Int_t      fInit;             // number of initialized estimators
  //
  TVectorD   fSample[3];        // !training sample for robust estimate of initial parameters
  TVectorD   fSampleStat[3];    // sample statistic
  Int_t      fCounter[3];       // counter of samples
  //
private:
  AliTPCkalmandEdx&  operator=(const AliTPCkalmandEdx&);// not implemented
  AliTPCkalmandEdx(const AliTPCkalmandEdx&); //not implemented
  ClassDef(AliTPCkalmandEdx,1);
};

#endif
Commit	Line	Data
e032f934	1	#ifndef ALITPCKALMANDEDX_H
	2	#define ALITPCKALMANDEDX_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 "TNamed.h"
	8	#include "TMatrixD.h"
	9	#include "TVectorD.h"
	10	class TTreeSRedirector;
	11	class TObjArray;
	12
	13	class AliTPCkalmandEdx: public TNamed{
	14	public:
	15	AliTPCkalmandEdx();
	16	AliTPCkalmandEdx(const char* name, const char* title, Int_t sampleSize=50);
	17	AliTPCkalmandEdx(const AliTPCkalmandEdx & kalman);
	18	void UpdatedEdxPair(Int_t ip0, Int_t ip1,
	19	Double_t dedx0, Double_t dedx1,
	20	Double_t s0, Double_t s1,
	21	Double_t kY0, Double_t kY1,
	22	Double_t kZ0, Double_t kZ1,
	23	Double_t dR0, Double_t dR1,
	24	TTreeSRedirector *debug=0);
	25	void UpdatedEdx(Int_t ip0,
	26	Double_t dedx0,
	27	Double_t dedxRef,
	28	Double_t s0,
	29	Double_t kY0,
	30	Double_t kZ0,
	31	Double_t dR0,
	32	TTreeSRedirector *debug=0);
	33	static Int_t GetIndex(Int_t i, Int_t j) { return 15*i+j;}
	34	//
	35	public:
	36	void AdddEdx(Int_t ip0,Double_t dedx0, Double_t dedxRef);
	37	void Init();
	38	TMatrixD * fState; // state vector
	39	TMatrixD * fCovariance; // covariance
	40	TMatrixD * fMat1; // helper unit matrix
	41	Int_t fNpad; // number of pad types
	42	Int_t fNpar; // number of parameters
	43	Int_t fNelem; // number of elements
	44	//
	45	// initial parameters estimate
	46	//
	47	Int_t fSampleSize; // size of starting sample
	48	Int_t fInit; // number of initialized estimators
	49	//
	50	TVectorD fSample[3]; // !training sample for robust estimate of initial parameters
	51	TVectorD fSampleStat[3]; // sample statistic
	52	Int_t fCounter[3]; // counter of samples
	53	//
d9382d87	54	private:
	55	AliTPCkalmandEdx& operator=(const AliTPCkalmandEdx&);// not implemented
	56	AliTPCkalmandEdx(const AliTPCkalmandEdx&); //not implemented
e032f934	57	ClassDef(AliTPCkalmandEdx,1);
	58	};
	59
	60	#endif
	61