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

/* $Id$ */

////////////////////////////////////////////////
//  Manager class for TPC parameters          //
////////////////////////////////////////////////
#include "AliTPCParam.h"

class AliTPCRF1D;
class AliTPCPRF2D;

class AliTPCParamSR : public AliTPCParam {
public:
  AliTPCParamSR();
  AliTPCParamSR(const AliTPCParamSR &param); // copy constructor
  AliTPCParamSR &operator = (const AliTPCParamSR & param); //assignment operator
  virtual ~AliTPCParamSR();

  Int_t  CalcResponse(Float_t* x, Int_t * index, Int_t row);
  Int_t CalcResponseFast(Float_t* x, Int_t * index, Int_t row,Float_t phase);
  //calculate bin response as function of the input position -x 
  //return number of valid response bin

  
  void XYZtoCRXYZ(Float_t *xyz, 
		  Int_t &sector, Int_t &padrow, Int_t option=3) const;
  //transform global position to the position relative to the sector padrow
  //if option=0  X calculate absolute            calculate sector
  //if option=1  X           absolute            use input sector
  //if option=2  X           relative to pad row calculate sector
  //if option=3  X           relative            use input sector

  void CRXYZtoXYZ(Float_t *xyz,
	    const Int_t &sector, const Int_t & padrow, Int_t option=3) const;  
  //transform relative position  to the gloabal position
  void TransformTo8(Float_t *xyz, Int_t *index) const;
  void TransformTo2(Float_t *xyz, Int_t *index) const;
  Bool_t Update();            //recalculate and check geometric parameters 
  void SetDefault();          //set default parameters
  void   SetInnerPRF(AliTPCPRF2D * prf) {fInnerPRF = prf;}
  void   SetOuter1PRF(AliTPCPRF2D * prf) {fOuter1PRF = prf;} //e.k
  void   SetOuter2PRF(AliTPCPRF2D * prf) {fOuter2PRF = prf;} //e.k
  void   SetTimeRF(AliTPCRF1D * timerf) {fTimeRF = timerf;}

  AliTPCPRF2D * GetInnerPRF() const {return fInnerPRF;}
  AliTPCPRF2D * GetOuter1PRF() const {return fOuter1PRF;} //e.k
  AliTPCPRF2D * GetOuter2PRF() const {return fOuter2PRF;} //e.k
  AliTPCRF1D  * GetTimeRF()   const {return fTimeRF;}
  void SetFacSigmaPadRow(Float_t fac=3.) {fFacSigmaPadRow=fac;}
  void SetFacSigmaPad(Float_t fac=3.) {fFacSigmaPad=fac;}
  void SetFacSigmaTime(Float_t fac=3.) {fFacSigmaTime=fac;}

  //  Float_t GetPadRowRadiiLow(Int_t irow) const;
  //  Float_t GetPadRowRadiiUp(Int_t irow) const;
  Float_t GetYInner(Int_t irow) const; //e,k
  Float_t GetYOuter(Int_t irow) const; //e.k

  virtual Float_t GetPrimaryLoss(Float_t *x, Int_t *index, Float_t *angle);
  virtual Float_t GetTotalLoss(Float_t *x, Int_t *index, Float_t *angle);

  virtual void GetClusterSize(Float_t *x, Int_t *index, Float_t *angle, Int_t mode, Float_t *sigma);
  virtual void GetSpaceResolution(Float_t *x, Int_t *index, Float_t *angle, Float_t amplitude, Int_t mode,Float_t *sigma);
  virtual Float_t  GetAmp(Float_t *x, Int_t *index, Float_t *angle);
  virtual Float_t * GetAnglesAccMomentum(Float_t *x, Int_t * index, Float_t* momentum, Float_t *angle); 
 
protected:
  AliTPCPRF2D * fInnerPRF;         //pad response function  for inner sector
  AliTPCPRF2D * fOuter1PRF;        //pad response function  for outer sector  
  AliTPCPRF2D * fOuter2PRF; 
  AliTPCRF1D  * fTimeRF;           //time response function object
  Float_t      fFacSigmaPadRow;    //factor-how many sigma of response I accept
  Float_t      fFacSigmaPad;       //factor-how many sigma of response I accept
  Float_t      fFacSigmaTime;      //factor-how many sigma of response I accept
  ClassDef(AliTPCParamSR,2)  //parameter  object for set:TPC
};

#endif