ITS/AliITSCorrMapSDD.h

   1 #ifndef ALIITSCORRMAPSDD_H
   2 #define ALIITSCORRMAPSDD_H
   3 /* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 /* $Id$ */
   7
   8 ///////////////////////////////////////////////////////////////////
   9 //                                                               //
  10 // Mother class for SDD maps used to correct for                 //
  11 // voltage divider shape and doping fluctuations                 //
  12 // Origin: F.Prino, Torino, prino@to.infn.it                     //
  13 //                                                               //
  14 ///////////////////////////////////////////////////////////////////
  15
  16 #include<TNamed.h>
  17 #include "AliLog.h"
  18 class TH1F;
  19 class TH2F;
  20 class AliITSsegmentationSDD;
  21
  22 class AliITSCorrMapSDD : public TNamed {
  23
  24  public:
  25   // maps produced from residuals stores Xtrue-Xmeasured in drift lenght.
  26   // Since the map computes xtrue-xmeas in drift coordinate
  27   // and the difference is added to measured local coordinate, we have
  28   // Left:  Xmeas_loc_corr = Xmeas_loc + (Xtrue-Xmeas)_loc = Xmeas_loc - (Xtrue-Xmeas)_drift
  29   // Right: Xmeas_loc_corr = Xmeas_loc + (Xtrue-Xmeas)_loc = Xmeas_loc + (Xtrue-Xmeas)_drift
  30   // hence, for the left side the sign of the correction need to inverted
  31   enum {kLeftInvBit = BIT(14)};
  32   AliITSCorrMapSDD();
  33   AliITSCorrMapSDD(Char_t *mapname);
  34   virtual ~AliITSCorrMapSDD(){};
  35   //
  36   void   SetInversionBit(Bool_t v=kTRUE)           {SetBit(kLeftInvBit,v);}
  37   Bool_t GetInversionBit()               const     {return TestBit(kLeftInvBit);}
  38   //
  39   Int_t GetNBinsAnode() const {return fNAnodePts;}
  40   Int_t GetNBinsDrift() const {return fNDriftPts;}
  41   void SetNBinsAnode(Int_t nbins) {
  42     if(nbins<=kMaxNAnodePts) fNAnodePts=nbins;
  43     else AliError(Form("Max. number of anode bins = %d",kMaxNAnodePts));
  44   }
  45   void SetNBinsDrift(Int_t nbins) {
  46     if(nbins<=kMaxNDriftPts) fNDriftPts=nbins;
  47     else AliError(Form("Max. number of drift bins = %d",kMaxNDriftPts));
  48   }
  49
  50   Bool_t CheckAnodeBounds(Int_t iAn) const {
  51     if(iAn<0 || iAn>=fNAnodePts)return kFALSE;
  52     else return kTRUE;
  53   }
  54   Bool_t CheckDriftBounds(Int_t iTb) const {
  55     if(iTb<0 || iTb >= fNDriftPts)return kFALSE;
  56     else return kTRUE;
  57   }
  58
  59   virtual void Set1DMap(TH1F* /*hmap*/){
  60     AliError("Not implemented");
  61   }
  62   virtual void Set2DMap(TH2F* /*hmap*/){
  63     AliError("Not implemented");
  64   }
  65
  66   virtual void ResetMap(){
  67     AliError("Not implemented");
  68   }
  69   virtual void SetCellContent(Int_t /*iAn*/, Int_t /*iTb*/, Float_t /*devMicron*/){
  70     AliError("Not implemented");
  71   }
  72   virtual Float_t GetCellContent(Int_t /*iAn*/, Int_t /*iTb*/) const {
  73     AliError("Not implemented");
  74     return -99999.;
  75   }
  76
  77   void    ComputeGridPoints(Float_t z, Float_t x, AliITSsegmentationSDD *seg, Bool_t isReco=kTRUE);
  78   Float_t GetCorrection(Float_t z, Float_t x, AliITSsegmentationSDD *seg);
  79   Float_t GetShiftForSimulation(Float_t z, Float_t x, AliITSsegmentationSDD *seg);
  80   TH2F* GetMapHisto() const;
  81   TH1F* GetMapProfile() const;
  82   TH1F* GetResidualDistr(Float_t dmin=-300., Float_t dmax=300.) const;
  83
  84
  85  protected:
  86
  87   enum {kMaxNAnodePts=256};// max number of map points along anodes
  88   enum {kMaxNDriftPts=291};// max number of map points along drift
  89
  90   static const Int_t fgkNAnodePtsDefault; // default value for fNAnodePts
  91   static const Int_t fgkNDriftPtsDefault; // default value for fNDriftPts
  92   Int_t fNAnodePts; // number of map points along anodes
  93   Int_t fNDriftPts; // number of map points along anodes
  94
  95   Float_t fXt1;   // true coordinate in lower grid point
  96   Float_t fXt2;   // true coordinate in upper grid point
  97   Float_t fXm1;   // measured coordinate in lower grid point
  98   Float_t fXm2;   // measured coordinate in upper grid point
  99   Float_t fDrLen; // drift length
 100
 101   ClassDef(AliITSCorrMapSDD,2);
 102 };
 103 #endif
 104