Example macro how to run Hough.

diff --git a/HLT/exa/houghtest.C b/HLT/exa/houghtest.C
new file mode 100644 (file)
index 0000000..c633e9c
--- /dev/null
+++ b/HLT/exa/houghtest.C
@@ -0,0 +1,64 @@
+void houghtest(Int_t slice,Char_t *path,Int_t n_eta_segments)
+{
+  
+
+  hough = new AliL3Hough(path,kTRUE,n_eta_segments);
+  hough->ReadData(slice);
+  
+  hough->Transform();
+  
+  hough->SetPeakThreshold(1000);
+  hough->AddAllHistograms();
+  hough->FindTrackCandidates();
+  
+  hough->WriteTracks(slice);
+  
+  //hough->Evaluate();
+  tracks = (AliL3TrackArray*)hough->GetTracks(0);
+  for(int i=0; i<tracks->GetNTracks(); i++)
+    {
+      track = (AliL3HoughTrack*)tracks->GetCheckedTrack(i);
+      if(!track) continue;
+      cout<<"pt "<<track->GetPt()<<" psi "<<track->GetPsi()<<" etaindex "<<track->GetEtaIndex()<<" weight "<<track->GetWeight()<<endl;
+    }
+  
+
+  display(hough,0);
+  
+}
+
+void display(AliL3Hough *hough,Int_t eta_index)
+{
+  //Display the data/tracks in eta_index
+  
+  hough->InitEvaluate();
+  digitd = new AliL3Histogram("Digits display","",250,0,250,250,-125,125);
+  trackd = new AliL3Histogram("Found tracks display","",250,0,250,250,-125,125);
+  for(int i=0; i<6; i++)
+    hough->GetEval(i)->DisplayEtaSlice(eta_index,digitd);
+  
+  tracks = (AliL3TrackArray*)hough->GetTracks(0);
+  float xyz[3];
+  for(int i=0; i<tracks->GetNTracks(); i++)
+    {
+      AliL3HoughTrack *track = (AliL3HoughTrack*)tracks->GetCheckedTrack(i);
+      if(!track) continue;
+      if(track->GetEtaIndex() != eta_index) continue;
+
+      for(int j=0; j<176; j++)
+       {
+         track->GetCrossingPoint(j,xyz);
+         trackd->Fill(xyz[0],xyz[1],1);
+       }
+    }
+  
+  //Draw the parameter space
+  c1 = new TCanvas("c1","",2);
+  hough->GetTransformer(0)->GetHistogram(eta_index)->Draw("box");
+  
+  //Draw the tracks
+  c2 = new TCanvas("c2","",2);
+  digitd->Draw();
+  trackd->Draw("same");
+  ((TH1F*)trackd->GetRootHisto())->SetMarkerColor(2);
+}