3 // Author: Constantin Loizides <mailto:loizides@fi.uib.no>
4 //*-- Copyright © ALICE HLT Group
6 #include "AliL3StandardIncludes.h"
8 #include "AliL3RootTypes.h"
9 #include "AliL3Logging.h"
10 #include "AliL3MemHandler.h"
11 #include "AliL3Transform.h"
12 #include "AliL3DigitData.h"
13 #include "AliL3HoughTransformerVhdl.h"
14 #include "AliL3FFloat.h"
20 /** \class AliL3HoughTransformerVhdl
22 //_____________________________________________________________
23 // AliL3HoughTransformerVhdl
25 // Hough transformation class for VHDL comparism.
30 ClassImp(AliL3HoughTransformerVhdl)
33 AliL3HoughTransformerVhdl::AliL3HoughTransformerVhdl()
34 : AliL3HoughTransformerLUT()
43 AliL3HoughTransformerVhdl::AliL3HoughTransformerVhdl(Int_t slice,Int_t patch,Int_t n_eta_segments,Int_t n_its)
44 : AliL3HoughTransformerLUT(slice,patch,n_eta_segments)
52 AliL3HoughTransformerVhdl::~AliL3HoughTransformerVhdl()
57 void AliL3HoughTransformerVhdl::CreateHistograms(Int_t nxbin,Float_t pt_min,Int_t nybin,Float_t phimin,Float_t phimax)
59 AliL3HoughTransformerLUT::CreateHistograms(nxbin,pt_min,nybin,phimin,phimax);
62 void AliL3HoughTransformerVhdl::CreateHistograms(Int_t nxbin,Float_t xmin,Float_t xmax,
63 Int_t nybin,Float_t ymin,Float_t ymax)
65 AliL3HoughTransformerLUT::CreateHistograms(nxbin,xmin,xmax,nybin,ymin,ymax);
67 fEpsilon=(ymax-ymin)/nybin;
68 fSinEpsilon=sin(fEpsilon);
69 fCosEpsilon=cos(fEpsilon);
77 //cout << fEpsilon << " - " << (xmax-xmin)/nxbin << endl;
80 void AliL3HoughTransformerVhdl::Init(Int_t slice,Int_t patch,Int_t n_eta_segments,Int_t n_its)
82 AliL3HoughTransformerLUT::Init(slice,patch,n_eta_segments);
85 void AliL3HoughTransformerVhdl::TransformCircle()
87 //Transform the input data with a circle HT.
88 //The function loops over all the data, and transforms each pixel with the equation:
90 //kappa = lastkappa +- epsilon * lastkappaprime
92 //kappaprime = lastkappaprime -+ epsilon * lastkappa
94 //Each pixel then transforms into a curve in the (kappa,phi0)-space. In order to find
95 //which histogram in which the pixel should be transformed, the eta-value is calcluated
96 //and the proper histogram index is found by GetEtaIndex(eta).
99 AliL3DigitRowData *tempPt = GetDataPointer();
102 LOG(AliL3Log::kError,"AliL3HoughTransformerVhdl::TransformCircle","Data")
103 <<"No input data "<<ENDLOG;
114 //Loop over the padrows:
115 for(Int_t i=fMinRow, row=0; i<=fMaxRow; i++, row++)
117 //Get the data on this padrow:
118 AliL3DigitData *digPt = tempPt->fDigitData;
119 if(i != (Int_t)tempPt->fRow)
121 LOG(AliL3Log::kError,"AliL3HoughTransformerVhdl::TransformCircle","Data")
122 <<"AliL3HoughTransformerLUT::TransformCircle : Mismatching padrow numbering "<<i<<" != "<<(Int_t)tempPt->fRow<<ENDLOG;
126 Float_t x = CalcX(row);
134 //Loop over the data on this padrow:
135 for(UInt_t j=0; j<tempPt->fNDigit; j++)
137 UShort_t charge = digPt[j].fCharge;
140 if((Int_t)charge <= GetLowerThreshold() || (Int_t)charge > GetUpperThreshold())
143 UChar_t pad = digPt[j].fPad;
144 UShort_t time = digPt[j].fTime;
146 if(fLastPad!=pad){ //only update if necessary
152 for(Int_t b=0; b<fNPhi0; b++)
153 fLUTKappa[b]=R2*(y*fLUT2cosphi0[b]-x*fLUT2sinphi0[b]);
155 //Fill the histogram along the phirange
157 Float_t kappaprime=0;
159 Float_t lastkappaprime=0;
161 Float_t R2=1/(x*x+y*y);
166 Float_t phi = fLUTphi0[fNPhi0/2];
167 kappa=A*cos(phi)+B*sin(phi);
168 kappaprime=B*cos(phi)-A*sin(phi);
171 lastkappaprime=kappaprime;
172 // hist->Fill(kappa,phi,charge);
174 for(Int_t b=fNPhi0/2+1; b<fNPhi0; b++){
175 //Float_t exact_kappa=R2*(y*fLUT2cosphi0[b]-x*fLUT2sinphi0[b]);
179 if(its==0) { //initialize or re-adjust values
180 kappa=A*cos(phi)+B*sin(phi); //equals exact_kappa!
181 kappaprime=B*cos(phi)-A*sin(phi);
184 //kappa=fCosEpsilon*lastkappa+fSinEpsilon*lastkappaprime;
185 //kappaprime=fCosEpsilon*lastkappaprime-fSinEpsilon*lastkappa;
186 kappa=lastkappa+fEpsilon*lastkappaprime;
187 kappaprime=lastkappaprime-fEpsilon*lastkappa;
192 lastkappaprime=kappaprime;
193 //hist->Fill(kappa,phi,charge);
194 } // end positive running values
196 phi = fLUTphi0[fNPhi0/2];
197 kappa=A*cos(phi)+B*sin(phi);
198 kappaprime=B*cos(phi)-A*sin(phi);
201 lastkappaprime=kappaprime;
202 //hist->Fill(kappa,fLUTphi0[b],charge);
204 for(Int_t b=fNPhi0/2-1; b>=0; b--){
205 //Float_t exact_kappa=R2*(y*fLUT2cosphi0[b]-x*fLUT2sinphi0[b]);
207 Float_t phi = fLUTphi0[b];
209 if(its==0) { //initialize or re-adjust values
210 kappa=A*cos(phi)+B*sin(phi); //equals exact_kappa!
211 kappaprime=B*cos(phi)-A*sin(phi);
214 //kappa=fCosEpsilon*lastkappa-fSinEpsilon*lastkappaprime;
215 //kappaprime=fCosEpsilon*lastkappaprime+fSinEpsilon*lastkappa;
216 kappa=lastkappa-fEpsilon*lastkappaprime;
217 kappaprime=lastkappaprime+fEpsilon*lastkappa;
222 lastkappaprime=kappaprime;
223 //hist->Fill(kappa,phi,charge);
229 Float_t z = CalcZ(time);
232 Float_t rz2 = 1 + (x*x+y*y)/(z*z);
233 Int_t eta_index = FindIndex(rz2);
235 if(eta_index < 0 || eta_index >= fNEtas){
236 //LOG(AliL3Log::kWarning,"AliL3HoughTransformerVhdl::TransformCircle","Histograms")<<"No histograms corresponding to eta index value of "<<eta_index<<"."<<ENDLOG;
240 //Get the correct histogrampointer:
241 AliL3Histogram *hist = fParamSpace[eta_index];
243 //LOG(AliL3Log::kWarning,"AliL3HoughTransformerVhdl::TransformCircle","Histograms")<<"Error getting histogram in index "<<eta_index<<"."<<ENDLOG;
247 for(Int_t b=0; b<fNPhi0; b++){
250 Float_t err=fabs((exact_kappa-kappa)/exact_kappa);
253 //cout << max_error << " - " << err << " " << kappa << " " << exact_kappa << " " << kappa/exact_kappa<< endl;
260 //hist->Fill(kappa,fLUTphi0[b],charge);
262 //cout << kappa << " " << fLUTphi0[b] << " " << charge << endl;
265 //Move the data pointer to the next padrow:
266 AliL3MemHandler::UpdateRowPointer(tempPt);
270 cout <<"Max Error: " << max_error << endl;
271 cout <<"Rel Error average: " << rel_error/counter << endl;
276 void AliL3HoughTransformerVhdl::Print()
278 AliL3HoughTransformerLUT::Print();
280 cout << "fEpsilon: " << fEpsilon << endl;
281 cout << "fIts: " << fIts << endl;
284 void AliL3HoughTransformerVhdl::PrintVhdl()
286 cout << "fSlice := " << GetSlice() << ";" << endl;
287 cout << "fPatch := " << GetPatch() << ";" << endl;
288 //cout << "fSector := " << fSector << ";" << endl;
289 //cout << "fSectorRow := " << fSectorRow << ";" << endl;
290 //cout << "fMinRow := " << fMinRow << ";" << endl;
291 //cout << "fMaxRow := " << fMaxRow << ";" << endl;
292 //cout << "fNRows := " << fNRows << ";" << endl;
293 //cout << "fNEtas := " << fNEtas << ";" << endl;
294 //cout << "fNPhi0 := " << fNPhi0 << ";" << endl;
295 cout << "fZSign := " << fZSign << ";" << endl;
296 cout << "fZLengthPlusOff := " << fZLengthPlusOff << ";" << endl;
297 cout << "fPadPitch := " << fPadPitch << ";" << endl;
298 cout << "fTimeWidth := " << fTimeWidth << ";" << endl;
301 cout << "fNLUTX :=" << fNRows << ";" << endl;
302 cout << "fLUTX := ( ";
303 for(Int_t i=0;i<fNRows-1;i++) cout << fLUTX[i] << ", ";
304 cout << fLUTX[fNRows-1] << " );\n" << endl;
305 cout << "fNLUTY := " << fNRows << ";" << endl;
306 cout << "fLUTY := ( ";
307 for(Int_t i=0;i<fNRows-1;i++) cout << fLUTY[i] << ", ";
308 cout << fLUTY[fNRows-1] << " );\n" << endl;
311 cout << "fNLUTEta := " << fNEtas << ";" << endl;
312 cout << "fLUTEta := (";
313 for(Int_t i=0;i<fNEtas-1;i++) cout << fLUTEta[i] << ", ";
314 cout << fLUTEta[fNEtas-1] << " );\n" << endl;
317 cout << "-- Kappa -- " << endl;
318 cout << "fNxbin := " << fNxbin << endl;
319 cout << "fxmin := " << fXmin << endl;
320 cout << "fxmax := " << fXmax << endl;
321 Float_t kdiff=(fXmax-fXmin)/fNxbin;
323 cout << "fKappa := ";
324 for(Int_t i=0;i<fNxbin;i++){
328 cout << k+kdiff << " );\n" << endl;
330 cout << "-- Phi --" << endl;
331 cout << "fNybin := " << fNybin << endl;
332 cout << "fymin := " << fYmin << endl;
333 cout << "fymax := " << fYmax << endl;
335 cout << "\nfNLUTphi0 := " << fNPhi0 << ";" << endl;
336 cout << "fLUTphi0 := (";
337 for(Int_t i=0;i<fNPhi0-1;i++) cout << fLUTphi0[i] << ", ";
338 cout << fLUTphi0[fNPhi0-1] << " );\n" << endl;
340 cout << "\nfNLUT2sinphi0 := " << fNPhi0 << ";" << endl;
341 cout << "fLUT2sinphi0 := (";
342 for(Int_t i=0;i<fNPhi0-1;i++) cout << fLUT2sinphi0[i] << ", ";
343 cout << fLUT2sinphi0[fNPhi0-1] << " );\n" << endl;
345 cout << "\nfNLUT2cosphi0 := " << fNPhi0 << ";" << endl;
346 cout << "fLUT2cosphi0 := (";
347 for(Int_t i=0;i<fNPhi0-1;i++) cout << fLUT2cosphi0[i] << ", ";
348 cout << fLUT2cosphi0[fNPhi0-1] << " );\n" << endl;
350 //cout << "\nfEpsilon := " << fEpsilon << endl;
351 //cout << "\nfSinEpsilon := " << fSinEpsilon << endl;
352 //cout << "\nfCosEpsilon := " << fCosEpsilon << endl;