TPC/AliTPCClusterParam.h

   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>
  16 #include <TVectorD.h>
  17
  18 class TTree;
  19 class TObjArray;
  20 //_____________________________________________________________________________
  21 class AliTPCClusterParam : public TObject {
  22  public:
  23   static AliTPCClusterParam* Instance();
  24   AliTPCClusterParam(){fRatio=0.01; fQNorm=0;}
  25   virtual           ~AliTPCClusterParam(){;}
  26   virtual void  Print(Option_t* option = "") const;
  27   void SetInstance(AliTPCClusterParam*param){fgInstance = param;}
  28   void FitData(TTree * tree);
  29   Float_t GetError0(Int_t dim, Int_t type, Float_t z, Float_t angle);
  30   Float_t GetError0Par(Int_t dim, Int_t type, Float_t z, Float_t angle);
  31   Float_t GetError1(Int_t dim, Int_t type, Float_t z, Float_t angle);
  32   Float_t GetErrorQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
  33   Float_t GetErrorQPar(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
  34   Float_t GetErrorQParScaled(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
  35   Float_t GetRMS0(Int_t dim, Int_t type, Float_t z, Float_t angle);
  36   Float_t GetRMS1(Int_t dim, Int_t type, Float_t z, Float_t angle);
  37   Float_t GetRMSQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
  38   Float_t GetRMSSigma(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean);
  39   Float_t GetShapeFactor(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean, Float_t rmsL, Float_t rmsM);
  40
  41   void Test(TTree * tree, const char *output="TestClusterParam.root");
  42   //
  43   // static methods equivalents  - use instance of param object - useful for tree draw and TF2 visualization
  44   static Float_t SGetError0(Int_t dim, Int_t type, Float_t z, Float_t angle){
  45     return fgInstance->GetError0(dim,type,z,angle);
  46   }
  47   static Float_t SGetError0Par(Int_t dim, Int_t type, Float_t z, Float_t angle){
  48     return fgInstance->GetError0Par(dim,type,z,angle);
  49   }
  50   static Float_t SGetError1(Int_t dim, Int_t type, Float_t z, Float_t angle){
  51     return fgInstance->GetError1(dim,type,z,angle);
  52   }
  53   static Float_t SGetErrorQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
  54     return fgInstance->GetErrorQ(dim,type,z,angle,Qmean);
  55   }
  56   static Float_t SGetErrorQPar(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
  57     return fgInstance->GetErrorQPar(dim,type,z,angle,Qmean);
  58   }
  59   static Float_t SGetErrorQParScaled(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
  60     return fgInstance->GetErrorQParScaled(dim,type,z,angle,Qmean);
  61   }
  62
  63   static Float_t SGetRMS0(Int_t dim, Int_t type, Float_t z, Float_t angle){
  64     return fgInstance->GetRMS0(dim,type,z,angle);
  65   }
  66   static Float_t SGetRMS1(Int_t dim, Int_t type, Float_t z, Float_t angle){
  67     return fgInstance->GetRMS1(dim,type,z,angle);
  68   }
  69   static Float_t SGetRMSQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
  70     return fgInstance->GetRMSQ(dim,type,z,angle,Qmean);
  71   }
  72   static Float_t SGetRMSSigma(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean){
  73     return fgInstance->GetRMSSigma(dim,type,z,angle,Qmean);
  74   }
  75   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){
  76     return fgInstance->GetShapeFactor(dim,type,z,angle,Qmean, rmsL, rmsM);
  77   }
  78   //
  79   //
  80   void FitResol(TTree * tree);
  81   void FitRMS(TTree * tree);
  82   Float_t Qnorm(Int_t ipad, Int_t itype, Float_t dr, Float_t ty, Float_t tz);
  83   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);}
  84
  85   void SetQnorm(Int_t ipad, Int_t itype,  TVectorD * norm);
  86
  87
  88  protected:
  89   void FitResol0(TTree * tree, Int_t dim, Int_t type, Float_t *param0, Float_t *error);
  90   void FitResol0Par(TTree * tree, Int_t dim, Int_t type, Float_t *param0, Float_t *error);
  91   void FitResol1(TTree * tree, Int_t dim, Float_t *param0, Float_t *error);
  92   void FitResolQ(TTree * tree, Int_t dim, Int_t type, Float_t *param0, Float_t *error);
  93   void FitResolQPar(TTree * tree, Int_t dim, Int_t type, Float_t *param0, Float_t *error);
  94   void FitRMS0(TTree * tree, Int_t dim, Int_t type, Float_t *param0, Float_t *error);
  95   void FitRMS1(TTree * tree, Int_t dim, Float_t *param0, Float_t *error);
  96   void FitRMSQ(TTree * tree, Int_t dim, Int_t type, Float_t *param0, Float_t *error);
  97   void FitRMSSigma(TTree * tree, Int_t dim, Int_t type, Float_t *param0, Float_t *error);
  98   //
  99
 100   Float_t fRatio;               //ratio of values constibution to error
 101   Float_t fParamS0[2][3][4];    //error parameterization coeficients
 102   Float_t fErrorS0[2][3][4];    //error parameterization coeficients
 103   Float_t fParamS0Par[2][3][7];    //error parameterization coeficients
 104   Float_t fErrorS0Par[2][3][7];    //error parameterization coeficients
 105   Float_t fParamSQ[2][3][6];    //error parameterization coeficients
 106   Float_t fErrorSQ[2][3][6];    //error parameterization coeficients
 107   Float_t fParamSQPar[2][3][9];    //error parameterization coeficients
 108   Float_t fErrorSQPar[2][3][9];    //error parameterization coeficients
 109   Float_t fParamS1[2][4];       //error parameterization coeficients
 110   Float_t fErrorS1[2][4];       //error parameterization coeficients
 111   //
 112   Float_t fParamRMS0[2][3][4];   //shape parameterization coeficients
 113   Float_t fErrorRMS0[2][3][4];   //shape parameterization coeficients
 114   Float_t fParamRMSQ[2][3][6];   //shape parameterization coeficients
 115   Float_t fErrorRMSQ[2][3][6];   //shape parameterization coeficients
 116   Float_t fParamRMS1[2][5];      //shape parameterization coeficients
 117   Float_t fErrorRMS1[2][5];      //shape parameterization coeficients
 118   Float_t fErrorRMSSys[2];        // systematic relative error of the parametererization
 119   Float_t fRMSSigmaRatio[2][2];   // mean value of the varation of RMS to RMS
 120   Float_t fRMSSigmaFit[2][3][2];   // mean value of the varation of RMS to RMS
 121   //
 122   TObjArray *fQNorm;              // q norm paramters
 123  protected:
 124   static AliTPCClusterParam*   fgInstance; //! Instance of this class (singleton implementation)
 125   ClassDef(AliTPCClusterParam,2)    //  TPC Cluster parameter class
 126 };
 127
 128 #endif