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();
  25   AliTPCClusterParam(const AliTPCClusterParam& param);
  26   AliTPCClusterParam & operator=(const AliTPCClusterParam&);
  27   virtual           ~AliTPCClusterParam();
  28   virtual void  Print(Option_t* option = "") const;
  29   void SetInstance(AliTPCClusterParam*param){fgInstance = param;}
  30   //
  31   // Seting functions
  32   //
  33   void FitData(TTree * tree);
  34   void FitResol(TTree * tree);
  35   void FitRMS(TTree * tree);
  36   void SetQnorm(Int_t ipad, Int_t itype,  TVectorD * norm);
  37   //
  38   // Charge parameterization
  39   //
  40   Float_t Qnorm(Int_t ipad, Int_t itype, Float_t dr, Float_t ty, Float_t tz);
  41
  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);
  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);}
  47   //
  48   // Error parameterization
  49   //
  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);
  56   //
  57   // Shape parameterization
  58   //
  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);
  64   //
  65   //
  66   //
  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   //
 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);}
 107
 108
 109  public:
 110   //
 111   //
 112   //
 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   //
 123
 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   //
 146   // charge normalization parametrization
 147   //
 148   TObjArray *fQNorm;              // q norm paramters
 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
 159   //
 160  protected:
 161   static AliTPCClusterParam*   fgInstance; //! Instance of this class (singleton implementation)
 162   ClassDef(AliTPCClusterParam,3)    //  TPC Cluster parameter class
 163 };
 164
 165 #endif