*/
}
-void AliL3Hough::Init(Char_t *path,Bool_t binary,Int_t n_eta_segments,Bool_t bit8=kFALSE)
+void AliL3Hough::Init(Char_t *path,Bool_t binary,Int_t n_eta_segments,Bool_t bit8=kFALSE,Int_t tv=0)
{
fBinary = binary;
strcpy(fPath,path);
fDoIterative = kFALSE;
fWriteDigits = kFALSE;
fUse8bits = bit8;
-
+ fversion = tv;
AliL3Transform::Init(fPath);
fNPatches = AliL3Transform::GetNPatches();
ProcessPatchIter(i);
fMerger->FillTracks(fTracks[i],i); //Copy tracks to merger
}
-
}
void AliL3Hough::ProcessPatchIter(Int_t patch)
<<AliL3Log::kDec<<"Found "<<tracks->GetNTracks()<<" tracks in patch "<<patch<<ENDLOG;
}
-
void AliL3Hough::AddAllHistograms()
{
//Add the histograms within one etaslice.
Double_t initTime,cpuTime;
initTime = GetCpuTime();
+
for(Int_t i=0; i<n_patches; i++)
{
AliL3HoughBaseTransformer *tr = fHoughTransformer[i];
Double_t eta_slice = (tr->GetEtaMax()-tr->GetEtaMin())/tr->GetNEtaSegments();
fTracks[i]->Reset();
+
for(Int_t j=0; j<fNEtaSegments; j++)
{
AliL3Histogram *hist = tr->GetHistogram(j);
fPeakFinder->Reset();
fPeakFinder->SetHistogram(hist);
fPeakFinder->FindMaxima(0,0); //Simple maxima finder
-
//fPeakFinder->FindAbsMaxima();
+
for(Int_t k=0; k<fPeakFinder->GetEntries(); k++)
{
if(fPeakFinder->GetWeight(k) == 0) continue;
return tv.tv_sec+(((Double_t)tv.tv_usec)/1000000.);
//return (Double_t)(clock()) / CLOCKS_PER_SEC;
}
+
AliL3Hough(Char_t *path,Bool_t binary,Int_t n_eta_segments=100,Int_t tversion=0);
virtual ~AliL3Hough();
- void Init(Char_t *path,Bool_t binary,Int_t n_eta_segments,Bool_t bit8=kFALSE);
+ void Init(Char_t *path,Bool_t binary,Int_t n_eta_segments,Bool_t bit8=kFALSE,Int_t tv=0);
void Process(Int_t minslice,Int_t maxslice);
void ReadData(Int_t slice,Int_t eventnr=0);
void Transform(Int_t row_range = -1);