[u/mrichter/AliRoot.git] / TPC / CLFinderDemo.C

//Begin_Html
/*
<img src="../../../../pic/CLFinderDemo.gif">
*/
//End_Html
void CLFinderDemo()
{

  AliTPCDigitsArray *gDigitsArray= GetDigitsArray();
  AliTPCClustersArray *gExactClusters = GetExactClustersArray();
  
  
  AliH2F * his = gDigitsArray->LoadRow(0,0)->GenerHisto()->GetSubrange2d(440,470,60,90);
  
  AliClusterFinder cf;
  cf.SetThreshold(3);
  cf.SetNoise(1);
  Int_t index[3]={1,0,0};
  cf.SetDetectorIndex(index);
  cf.SetDetectorParam(gTPCParam);
  cf.GetHisto(his);
    
  cf.SetBFit(kFALSE);
  cf.GetMinuit()->SetPrintLevel(-1);
  cf.GetMinuit()->SetMaxIterations(100);
  cf.FindPeaks1()->GetEntriesFast();
    
    
  TCanvas *c0 = new TCanvas("signals", "Cluster finder drawing", 600,1100); 
  c0->Divide(1,3);
  c0->cd(1);
  cf.Draw("lego2");

  c0->cd(2);
  cf.Draw("cont1");
  cf.DrawCluster(4,5,5);  
  c0->cd(3);
  cf.DrawBorders("lego2");
    
}
Commit	Line	Data
cc80f89e	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