+ // peaks->Draw("same");
+ c1->cd(2);
+ image->Draw("");
+ //road->Draw("same");
+}
+
+void process(char *path="/prog/alice/data/Rawdata/6_patch/hg_42105_s1-3/",bool bin=true)
+{
+ double torad = 3.1415/180;
+ a = new AliL3Hough(path,bin,1);
+ a->TransformSlice(1);
+
+ hist = (AliL3Histogram*)a->AddHistograms();
+
+ //hist->SetThreshold(10000);
+
+ b = new AliL3HoughMaxFinder("KappaPhi");
+ b->SetHistogram(hist);
+
+ Int_t xbin,ybin;
+
+ Int_t n=10;
+ Float_t x[10];
+ Float_t y[10];
+ b->FindPeak1(x,y,n);
+ printf("peak at pt %f phi0 %f\n",0.2*0.003/x[0],y[0]/torad);
+
+ track = new AliL3HoughTrack();
+ track->SetTrackParameters(x[0],y[0],1);
+
+ image = new AliL3Histogram("image","",250,0,250,250,-125,125);
+ a->Evaluate(image);
+ TH2F *road = new TH2F("road","",250,0,250,250,-125,125);
+ road->SetMarkerStyle(5);
+ road->SetMarkerColor(2);
+
+ float xyz[3];
+ for(int i=0; i<176; i++)
+ {
+ if(i%10) continue;
+ track->GetCrossingPoint(i,xyz);
+ road->Fill(xyz[0],xyz[1],1);
+ }
+ c1 = new TCanvas("c1","",1000,500);
+ c1->Divide(2);
+ c1->cd(1);
+ hist->Draw("box");