]> git.uio.no Git - u/mrichter/AliRoot.git/blame - HLT/src/AliL3Transform.h
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Bug correction (energy recalculation when adding new incarnation of the particle).
[u/mrichter/AliRoot.git] / HLT / src / AliL3Transform.h
CommitLineData
3e87ef69 1// @(#) $Id$
2
108615fc 3#ifndef ALIL3TRANSFORM_H
4#define ALIL3TRANSFORM_H
5
108615fc 6#include "AliL3RootTypes.h"
108615fc 7
8class AliL3Transform {
494fad94 9
3e87ef69 10 public:
11 enum VersionType { fV_default=0, fV_deprecated=1, fV_aliroot=10, fV_cosmics=100};
12
108615fc 13 private:
af99f571 14 const static Double_t fBFACT;
3e87ef69 15 const static Double_t fPi;
16 const static Double_t fAnodeWireSpacing;
17
18 static Int_t fNPatches; //6 (dont change this)
19 static Int_t fRows[6][2];
20 static Int_t fNRows[6];
21
ef79795d 22 static Double_t fBField;
3e87ef69 23 static Double_t fBFieldFactor;
24 static Double_t fSolenoidBField;
494fad94 25 static Int_t fNTimeBins;
26 static Int_t fNRowLow;
27 static Int_t fNRowUp;
5923e698 28 static Int_t fNRowUp1;
29 static Int_t fNRowUp2;
494fad94 30 static Int_t fNSectorLow;
31 static Int_t fNSectorUp;
3e87ef69 32 static Int_t fSlice2Sector[36][2];
33 static Int_t fSector2Slice[72];
34 static Int_t fSectorLow[72];
494fad94 35 static Double_t fPadPitchWidthLow;
36 static Double_t fPadPitchWidthUp;
37 static Double_t fZWidth;
38 static Double_t fZSigma;
494fad94 39 static Double_t fZLength;
40 static Double_t fZOffset;
3e87ef69 41 static Int_t fNSector; //72 (dont change this)
42 static Int_t fNSlice; //36 (dont change this)
43 static Int_t fNRow; //159 (dont change this)
44 static Double_t fNRotShift; //Rotation shift (eg. 0.5 for 10 degrees)
45 static Int_t fNPads[159]; //fill this following Init and fVersion
46 static Double_t fX[159]; //X position in local coordinates
47 static Int_t fVersion; //flags the version
768861a0 48 static Double_t fDiffT; //Transversal diffusion constant
49 static Double_t fDiffL; //Longitudinal diffusion constant
3e87ef69 50 static Double_t fOmegaTau; //ExB effects
768861a0 51 static Double_t fInnerPadLength;
5923e698 52 static Double_t fOuter1PadLength;
53 static Double_t fOuter2PadLength;
768861a0 54 static Double_t fInnerPRFSigma;
5923e698 55 static Double_t fOuter1PRFSigma;
56 static Double_t fOuter2PRFSigma;
768861a0 57 static Double_t fTimeSigma; //Minimal longitudinal width
3e87ef69 58 static Int_t fADCSat; //ADC Saturation (1024 = 10 bit)
59 static Int_t fZeroSup; //Zero suppression threshold
60 static Double_t fCos[36]; //stores the cos value for local to global rotations
61 static Double_t fSin[36]; //stores the sin value for local to global rotations
af99f571 62
3e87ef69 63 public:
64 static Bool_t Init(Char_t* path,Bool_t UseAliTPCParam=kFALSE); //init transformer settings (versions)
02f030e3 65 static Bool_t MakeInitFile(Char_t *rootfilename,Char_t *filename); //create the init file from rootfile
66 static Bool_t ReadInit(Char_t *path); //read init (possibly from root file)
67 static Bool_t ReadInitFile(Char_t *path); //read init from text file
68 static Bool_t SaveInitFile(Char_t *filename); //save parameters in init file
2220326f 69
3e87ef69 70 //setters
71 static void SetNPatches(Int_t i){fNPatches = i;}
72 static void SetNRows(Int_t s[6]){for(Int_t i=0;i<fNPatches;i++) fNRows[i] = s[i];}
73 static void SetRows(Int_t s[6][2]){
74 for(Int_t i=0;i<fNPatches;i++){
75 fRows[i][0] = s[i][0];
76 fRows[i][1] = s[i][1];
77 }
78 }
79 static void SetBField(Double_t f) {fBField = f;} //careful, these 3 are not independent!
80 static void SetBFieldFactor(Double_t f) {fBFieldFactor = f;fBField=fBFieldFactor*fSolenoidBField*0.1;}
81 static void SetSolenoidBField(Double_t f){fSolenoidBField = f;fBField=fBFieldFactor*fSolenoidBField*0.1;}
82 static void SetNTimeBins(Int_t i){fNTimeBins = i;}
83 static void SetNRowLow(Int_t i){fNRowLow = i;}
84 static void SetNRowUp(Int_t i){fNRowUp = i;}
85 static void SetNRowUp1(Int_t i){fNRowUp1 = i;}
86 static void SetNRowUp2(Int_t i){fNRowUp2 = i;}
87 static void SetSlice2Sector(Int_t s[36][2]){
88 for(Int_t i=0;i<fNSlice;i++){
89 fSlice2Sector[i][0] = s[i][0];
90 fSlice2Sector[i][1] = s[i][1];
91 }
92 }
93 static void SetSector2Slice(Int_t s[72]){for(Int_t i=0;i<fNSector;i++) fSector2Slice[i] = s[i];}
94 static void SetSectorLow(Int_t s[72]){for(Int_t i=0;i<fNSector;i++) fSectorLow[i] = s[i];}
95 static void SetNSectorLow(Int_t i){fNSectorLow = i;}
96 static void SetNSectorUp(Int_t i){fNSectorUp = i;}
97 static void SetPadPitchWidthLow(Double_t f){fPadPitchWidthLow = f;}
98 static void SetPadPitchWidthUp(Double_t f){fPadPitchWidthUp = f;}
99 static void SetZWidth(Double_t f){fZWidth = f;}
100 static void SetZSigma(Double_t f){fZSigma = f;}
101 static void SetZLength(Double_t f){fZLength = f;}
102 static void SetZOffset(Double_t f){fZOffset = f;}
103 static void SetNSector(Int_t i){fNSector = i;}
104 static void SetNSlice(Int_t i){fNSlice = i;}
105 static void SetNRow(Int_t i){fNRow = i;}
106 static void SetNRotShift(Double_t f){fNRotShift = f;}
107 static void SetNPads(Int_t pads[159]){for(Int_t i=0;i<fNRow;i++) fNPads[i] = pads[i];}
108 static void SetX(Double_t xs[159]){for(Int_t i=0;i<fNRow;i++) fX[i] = xs[i];}
109 static void SetVersion(Int_t i){fVersion = i;}
110 static void SetDiffT(Double_t f){fDiffT = f;}
111 static void SetDiffL(Double_t f){fDiffL = f;}
112 static void SetOmegaTau(Double_t f){fOmegaTau = f;}
113 static void SetInnerPadLength(Double_t f){fInnerPadLength = f;}
114 static void SetOuter1PadLength(Double_t f){fOuter1PadLength = f;}
115 static void SetOuter2PadLength(Double_t f){fOuter2PadLength = f;}
116 static void SetInnerPRFSigma(Double_t f){fInnerPRFSigma = f;}
117 static void SetOuter1PRFSigma(Double_t f){fOuter1PRFSigma = f;}
118 static void SetOuter2PRFSigma(Double_t f){fOuter2PRFSigma = f;}
119 static void SetTimeSigma(Double_t f){fTimeSigma = f;}
120 static void SetADCSat(Int_t i) {fADCSat = i;}
121 static void SetZeroSup(Int_t i) {fZeroSup = i;}
af99f571 122
3e87ef69 123 //getters
5e0f9911 124 static const Char_t* GetParamName() {return "75x40_100x60_150x60";}
125
e376e129 126 static Int_t GetFirstRow(Int_t patch);
127 static Int_t GetLastRow(Int_t patch);
128 static Int_t GetNRows(Int_t patch);
3e87ef69 129 static Int_t GetPatch(Int_t padrow);
f500b93b 130 static Int_t GetNRows() {return fNRow;}
3e87ef69 131 static Int_t GetNRowLow() {return fNRowLow;}
132 static Int_t GetNRowUp1() {return fNRowUp1;}
133 static Int_t GetNRowUp2() {return fNRowUp2;}
134 static Int_t GetPadRow(Float_t x);
9ca67380 135 static Int_t GetNPatches() {return fNPatches;}
3e87ef69 136 static Int_t GetNPads(Int_t row);
137 static Int_t GetNTimeBins(){return fNTimeBins;}
ef79795d 138 static Double_t GetBField() {return fBField;}
3e87ef69 139 static Double_t GetSolenoidField() {return fSolenoidBField;}
af99f571 140 static Double_t GetBFact() {return fBFACT;}
3e87ef69 141 static Double_t GetBFactFactor() {return fBFieldFactor;}
af99f571 142 static Double_t GetBFieldValue() {return (fBField*fBFACT);}
3e44d77c 143 static Float_t Deg2Rad(Float_t angle) {return angle*fPi/180;}
9ca67380 144 static Double_t Pi() {return fPi;}
494fad94 145 static Int_t GetVersion(){return fVersion;}
146 static Double_t GetPadPitchWidthLow() {return fPadPitchWidthLow;}
147 static Double_t GetPadPitchWidthUp() {return fPadPitchWidthUp;}
3e87ef69 148 static Double_t GetPadPitchWidth(Int_t patch);
494fad94 149 static Double_t GetZWidth() {return fZWidth;}
150 static Double_t GetZLength() {return fZLength;}
151 static Double_t GetZOffset() {return fZOffset;}
768861a0 152 static Double_t GetDiffT() {return fDiffT;}
153 static Double_t GetDiffL() {return fDiffL;}
3e87ef69 154 static Double_t GetParSigmaY2(Int_t padrow,Float_t z,Float_t angle);
155 static Double_t GetParSigmaZ2(Int_t padrow,Float_t z,Float_t tgl);
156 static Double_t GetOmegaTau() {return fOmegaTau;}
768861a0 157 static Double_t GetAnodeWireSpacing() {return fAnodeWireSpacing;}
5923e698 158 static Double_t GetPadLength(Int_t padrow);
159 static Double_t GetPRFSigma(Int_t padrow);
768861a0 160 static Double_t GetTimeSigma() {return fTimeSigma;}
3e87ef69 161 static Double_t GetZSigma() {return fZSigma;}
6f3df1c8 162 static Int_t GetADCSat() {return fADCSat;}
3e87ef69 163 static Int_t GetZeroSup() {return fZeroSup;}
164 static Int_t GetNSlice() {return fNSlice;}
165 static Int_t GetNSector() {return fNSector;}
494fad94 166 static Int_t GetNSectorLow() {return fNSectorLow;}
167 static Int_t GetNSectorUp() {return fNSectorUp;}
95a00d93 168
494fad94 169 static Bool_t Slice2Sector(Int_t slice, Int_t slicerow, Int_t &sector, Int_t &row);
170 static Bool_t Sector2Slice(Int_t &slice, Int_t sector);
171 static Bool_t Sector2Slice(Int_t &slice, Int_t &slicerow, Int_t sector, Int_t row);
1727f1c9 172
494fad94 173 static Double_t Row2X(Int_t slicerow);
174 static Double_t GetMaxY(Int_t slicerow);
175 static Double_t GetEta(Float_t *xyz);
3ac82106 176 static Double_t GetEta(Int_t slice,Int_t padrow, Int_t pad, Int_t time);
494fad94 177 static Double_t GetPhi(Float_t *xyz);
1727f1c9 178
494fad94 179 static void XYZtoRPhiEta(Float_t *rpe, Float_t *xyz);
180 static void Local2Global(Float_t *xyz, Int_t slice);
181 static void Local2GlobalAngle(Float_t *angle, Int_t slice);
182 static void Global2LocalAngle(Float_t *angle, Int_t slice);
108615fc 183
494fad94 184 static void Raw2Local(Float_t *xyz, Int_t sector, Int_t row, Float_t pad, Float_t time);
185 static void Local2Global(Float_t *xyz, Int_t sector, Int_t row);
186 static void Global2Local(Float_t *xyz, Int_t sector, Bool_t isSlice=kFALSE);
187 static void Raw2Global(Float_t *xyz, Int_t sector, Int_t row, Float_t pad, Float_t time);
188 static void Local2Raw(Float_t *xyz, Int_t sector, Int_t row);
189 static void Global2Raw(Float_t *xyz, Int_t sector, Int_t row);
3e87ef69 190
191 static void PrintCompileOptions();
108615fc 192
4499ed26 193 ClassDef(AliL3Transform,1)
108615fc 194};
108615fc 195#endif