1 Int_t AliTPCHits2ExactClusters(Int_t nevent=1)
4 // new version by J.Belikov
6 // Connect the Root Galice file containing Geometry, Kine and Hits
8 const char * inFile_old = "galice.root";
9 const char * inFile_new = "galice.root";
10 TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject(inFile_old);
11 if (file) {file->Close(); delete file;}
12 file = TFile::Open(inFile_new,"UPDATE");
13 if (!file->IsOpen()) {
14 cerr<<"Can't open "<<inFile_new<<" !\n";
18 // Get AliRun object from file or create it if not on file
19 if (gAlice) delete gAlice;
20 gAlice = (AliRun*)file->Get("gAlice");
22 cerr<<"AliTPCHits2Digits.C : AliRun object not found on file\n";
28 // gAlice->GetEvent(0);
29 AliTPC *TPC = (AliTPC*)gAlice->GetDetector("TPC");
30 AliTPCParam * param = TPC->GetParam();
34 for(Int_t eventn =0;eventn<nevent;eventn++){
35 printf("Processing event %d\n",eventn);
36 gAlice->GetEvent(eventn);
37 //setup AliTPCClustersArray
38 AliTPCClustersArray * arr=new AliTPCClustersArray;
39 arr->SetClusterType("AliComplexCluster");
42 TPC->SetClustersArray(arr);
43 for (Int_t i=0;i<72;i++)
44 TPC->Hits2ExactClustersSector(i);
46 sprintf(treeName,"TreeCExact_%d",eventn);
47 TPC->GetClustersArray()->GetTree()->Write(treeName);
50 delete gAlice; gAlice=0;
51 file->Close(); delete file;