Added option for different binning of DCAxy axis in THnSparse. Same width for all...
[u/mrichter/AliRoot.git] / 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