TPC/CLFinderDemo.C

   1 //Begin_Html
   2 /*
   3 <img src="../../../../pic/CLFinderDemo.gif">
   4 */
   5 //End_Html
   6 void CLFinderDemo()
   7 {
   8
   9   AliTPCDigitsArray *gDigitsArray= GetDigitsArray();
  10   AliTPCClustersArray *gExactClusters = GetExactClustersArray();
  11
  12
  13   AliH2F * his = gDigitsArray->LoadRow(0,0)->GenerHisto()->GetSubrange2d(440,470,60,90);
  14
  15   AliClusterFinder cf;
  16   cf.SetThreshold(3);
  17   cf.SetNoise(1);
  18   Int_t index[3]={1,0,0};
  19   cf.SetDetectorIndex(index);
  20   cf.SetDetectorParam(gTPCParam);
  21   cf.GetHisto(his);
  22
  23   cf.SetBFit(kFALSE);
  24   cf.GetMinuit()->SetPrintLevel(-1);
  25   cf.GetMinuit()->SetMaxIterations(100);
  26   cf.FindPeaks1()->GetEntriesFast();
  27
  28
  29   TCanvas *c0 = new TCanvas("signals", "Cluster finder drawing", 600,1100);
  30   c0->Divide(1,3);
  31   c0->cd(1);
  32   cf.Draw("lego2");
  33
  34   c0->cd(2);
  35   cf.Draw("cont1");
  36   cf.DrawCluster(4,5,5);
  37   c0->cd(3);
  38   cf.DrawBorders("lego2");
  39
  40 }
  41