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