1 void houghtest(Int_t slice,Char_t *path,Int_t n_eta_segments)
5 hough = new AliL3Hough();
6 hough->Init(path,kTRUE,n_eta_segments,kTRUE);
8 hough->ReadData(slice);
12 hough->SetPeakThreshold(1);
13 hough->AddAllHistograms();
14 hough->FindTrackCandidates();
16 hough->WriteTracks(slice);
19 tracks = (AliL3TrackArray*)hough->GetTracks(0);
20 for(int i=0; i<tracks->GetNTracks(); i++)
22 track = (AliL3HoughTrack*)tracks->GetCheckedTrack(i);
24 cout<<"pt "<<track->GetPt()<<" psi "<<track->GetPsi()<<" etaindex "<<track->GetEtaIndex()<<" weight "<<track->GetWeight()<<endl;
32 void display(AliL3Hough *hough,Int_t eta_index)
34 //Display the data/tracks in eta_index
36 hough->InitEvaluate();
37 digitd = new AliL3Histogram("Digits display","",250,0,250,250,-125,125);
38 trackd = new AliL3Histogram("Found tracks display","",250,0,250,250,-125,125);
39 for(int i=0; i<6; i++)
40 hough->GetEval(i)->DisplayEtaSlice(eta_index,digitd);
42 tracks = (AliL3TrackArray*)hough->GetTracks(0);
44 for(int i=0; i<tracks->GetNTracks(); i++)
46 AliL3HoughTrack *track = (AliL3HoughTrack*)tracks->GetCheckedTrack(i);
48 if(track->GetEtaIndex() != eta_index) continue;
50 for(int j=0; j<176; j++)
52 track->GetCrossingPoint(j,xyz);
53 trackd->Fill(xyz[0],xyz[1],1);
57 //Draw the parameter space
58 c1 = new TCanvas("c1","",2);
59 hough->GetTransformer(0)->GetHistogram(eta_index)->Draw("lego");
62 c2 = new TCanvas("c2","",2);
65 ((TH1F*)trackd->GetRootHisto())->SetMarkerColor(2);