Program to run hough transform standalone.

diff --git a/HLT/programs/runhough.cxx b/HLT/programs/runhough.cxx
new file mode 100644 (file)
index 0000000..fbce379
--- /dev/null
+++ b/HLT/programs/runhough.cxx
@@ -0,0 +1,206 @@
+// $Id$
+
+// Author: Constantin Loizides <mailto: loizides@ikf.physik.uni-frankfurt.de
+
+#include <stream.h>
+#include <libgen.h>
+
+#include <AliL3RootTypes.h>
+#include <AliL3Transform.h>
+#include <AliL3Track.h>
+#include <AliL3TrackArray.h>
+#include <AliL3HoughTrack.h>
+#include <AliL3ClustFinderNew.h>
+#include <AliL3MemHandler.h>
+#include <AliL3SpacePointData.h>
+#include <AliL3Logging.h>
+#include <AliL3Logger.h>
+#include <AliL3HoughBaseTransformer.h>
+#include <AliL3HoughTransformerVhdl.h>
+#include <AliL3HoughTransformer.h>
+#include <AliL3Hough.h>
+//#include <AliL3FFloat.h>
+#include <AliL3HoughMaxFinder.h>
+#ifndef no_root
+#include <TROOT.h>
+#include <TApplication.h>
+#include <TCanvas.h>
+#include <TH2.h>
+#include <TGraph.h>
+#include <TGraphErrors.h>
+#endif
+
+
+int main(int argc,char **argv)
+{
+  Int_t sl=0;
+  Int_t sh=0;
+  Int_t segs=100;
+
+  AliL3Logger l;
+  l.Set(AliL3Logger::kAll);
+  l.UseStderr();
+  //l.UseStdout();
+  //l.UseStream();
+
+  Char_t path[1024];
+  if(argc<2){
+    cout<<"Usage: runhough path [slow] [shigh] [segs]"<<endl;
+    exit(1);
+  }
+  strcpy(path,argv[1]);
+  if (argc>2) {
+    sl=atoi(argv[2]);
+  }
+  if (argc>3) {
+    sh=atoi(argv[3]);
+  }
+  if (argc>4) {
+    segs=atoi(argv[4]);
+  }
+
+  //AliL3FFloat::SetParams(10000);
+  AliL3Transform::Init(path);
+
+#if 0
+  runhough(sl,sh,path,segs);
+#else
+  //VESTBO: Look here, init of six trafos
+  AliL3HoughBaseTransformer **fHoughTransformer = new AliL3HoughBaseTransformer*[6];
+
+  for(int i=0;i<6;i++){
+    cout << "----------------------------- " << i << "----------------------------- " << endl;
+
+    //VESTBO: init one of them
+    fHoughTransformer[i] = new AliL3HoughTransformerVhdl(0,i,segs);
+    //VESTBO: Print parameters
+    fHoughTransformer[i]->Print();
+
+    //VESTBO: Play around with combinations of creating histos or/and init values.
+    fHoughTransformer[i]->CreateHistograms(64,0.1,64,-30,30);
+    //fHoughTransformer[i]->Init(sl,i,segs);
+
+    //VESTBO: Print parameters. Error should happen here: crash or be in endless loop, that means 
+    //somehow data members get overwritten!!
+    fHoughTransformer[i]->Print();
+  }
+#endif
+
+  //AliL3FFloat::PrintStat();
+  exit(0);
+}
+
+
+
+
+//----------------------------------------------------------------------------------
+// dont look beyond...
+//----------------------------------------------------------------------------------
+
+#if 0
+void runhough(Int_t sl,Int_t sh,Char_t *path,Int_t n_eta_segments, Int_t show_seg=-1)
+{
+
+  Bool_t binary = kTRUE;
+  Bool_t bit8 = kTRUE;
+  Int_t tv=1;
+  Int_t th=14000;
+
+  AliL3Hough *hough = new AliL3Hough();
+  hough->Init(path,binary,n_eta_segments,bit8,tv);
+  hough->GetMaxFinder()->SetThreshold(th);
+  Int_t ntracks=0;
+
+  for(Int_t slice=sl;slice<=sh;slice++){
+    hough->ReadData(slice);
+    hough->Transform();
+    hough->AddAllHistograms();
+    hough->FindTrackCandidates();
+    //hough->Evaluate(5);
+
+    AliL3TrackArray *tracks = (AliL3TrackArray*)hough->GetTracks(0);
+    ntracks=tracks->GetNTracks();
+    for(int i=0; i<ntracks; i++)
+      {
+       AliL3HoughTrack *track = (AliL3HoughTrack*)tracks->GetCheckedTrack(i);
+       if(!track) continue;
+       if(sl==sh) cout<<"pt "<<track->GetPt()<<" psi "<<track->GetPsi()<<" eta "<<track->GetEta()<<" etaindex "<<track->GetEtaIndex()<<" weight "<<track->GetWeight()<<endl;
+       if(show_seg<0) show_seg=track->GetEtaIndex();
+      }
+
+    cout<<"Found total "<<tracks->GetNTracks()<<" tracks"<<endl;
+    hough->WriteTracks(slice);
+  }
+
+  //if((ntracks>0)&&(sl==sh)) display(hough,show_seg);
+}
+#endif
+
+#if 0
+void display(AliL3Hough *hough,Int_t eta_index)
+{
+  //Display the data/tracks in eta_index
+  
+  hough->InitEvaluate();
+  AliL3Histogram *digitd = new AliL3Histogram("Digits display","",250,0,250,250,-125,125);
+  AliL3Histogram *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);
+  
+  float xyz[3];
+  tracks = (AliL3TrackArray*)hough->GetTracks(0);
+  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
+  TCanvas *c1 = new TCanvas("c1","",2);
+  hough->GetTransformer(0)->GetHistogram(eta_index)->Draw("box");
+  
+  //Draw the tracks
+  TCanvas *c2 = new TCanvas("c2","",2);
+  digitd->Draw();
+  trackd->Draw("same");
+  ((TH1F*)trackd->GetRootHisto())->SetMarkerColor(2);
+}
+#endif
+
+#if 0
+struct GoodTrack
+{
+  Int_t event;
+  Int_t label;
+  Double_t eta;
+  Int_t code;
+  Double_t px,py,pz;
+  Double_t pt;
+  Int_t nhits;
+};
+
+void geteff(char *fname)
+{
+  GoodTrack gt[15000];
+  int counter=0;
+  ifstream in(fname);
+  if(!in)
+    {
+      cerr<<"Could not open "<<fname<<endl;
+      return;
+    }
+  while(in>>gt[counter].event>>gt[counter].label>>gt[counter].code
+       >>gt[counter].px>>gt[counter].py>>gt[counter].pz>>gt[counter].pt>>gt[counter].eta>>gt[counter].nhits)
+    counter++;
+  
+  char filename[100];
+  file = new AliL3MemHandler();
+}
+#endif