AliTPCRF1D.h

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




// include files and class forward declarations


#include "TObject.h"
#include "TMath.h"
class TF1;


class AliTPCRF1D : public TObject {
public :
  AliTPCRF1D(Bool_t direct=kFALSE,Int_t np=0,Float_t step=0 );
  AliTPCRF1D(const AliTPCRF1D &prf);
  AliTPCRF1D & operator = (const AliTPCRF1D &prf);
  ~AliTPCRF1D();
  Float_t GetRF(Float_t xin); //return RF in point xin
  Float_t GetGRF(Float_t xin); //return generic response function  in xin
  void SetGauss(Float_t sigma,Float_t padWidth, Float_t kNorm);
  //adjust RF with GAUSIAN as generic GRF
  //if  direct = kTRUE then it does't convolute distribution
  void SetCosh(Float_t sigma,Float_t padWidth, Float_t kNorm);
  void SetGati(Float_t K3, Float_t padDistance, Float_t padWidth,
               Float_t kNorm);
  //adjust RF with 1/Cosh  as generic GRF
  void SetParam(TF1 * GRF,Float_t padwidth,Float_t kNorm,
                Float_t sigma=0);
  //adjust RF with general function
  void SetOffset(Float_t xoff) {fOffset=xoff;}
  //set offset value
  Float_t GetOffset(){return fOffset;}
  Float_t GetPadWidth(){ return fpadWidth;};
  //return  pad width
  Float_t  GetSigma(){return fSigma;}
  //return estimated sigma of RF
  void DrawRF(Float_t x1=-3 ,Float_t x2 =3.,Int_t N = 200);
  //draw RF it don't delete histograms after drawing
  void Update();
  static Double_t Gamma4(Double_t x, Double_t p0, Double_t p1);
private:
  Double_t funParam[5];
  Int_t  fNRF;      
  Float_t fDSTEPM1;    
  Float_t* fcharge; //[fNPRF]
  Float_t  forigsigma;
  Float_t fpadWidth;  
  Float_t fkNorm;     
  Float_t fInteg;     
  TF1 *  fGRF;        
  Float_t fSigma;     

  Float_t fOffset;    
  //have for expample shifted gauss)
  //calculated during update

  Bool_t fDirect;     

  Float_t fPadDistance;   
  char  fType[5];     
  static Int_t   fgNRF;
  static Float_t fgRFDSTEP;
  ClassDef(AliTPCRF1D,2)
};




#endif /* ALITPCRF1D_H */