+
+ const char * inFile_new = "galice.root";
+ AliRunLoader *rl = AliRunLoader::Open(inFile_new,"Event","read");
+
+ Int_t nevents=rl->GetNumberOfEvents();
+ cout<<"Number of Events:"<<nevents<<endl;
+ Int_t choice=0;
+ do{
+ cout<<"Insert the event number: ";
+ cin>>choice;
+ }while (choice<=0 || choice>nevents);
+ rl->GetEvent(choice-1);
+ AliLoader *tpcloader=rl->GetLoader("TPCLoader");
+ tpcloader->LoadDigits();
+ TTree *digitsTree=tpcloader->TreeD();
+
+ AliSimDigits digrows, *dummy=&digrows;
+ digitsTree->GetBranch("Segment")->SetAddress(&dummy);
+ Stat_t nrows = digitsTree->GetEntries();
+ cout<<"Number of entries (rows):"<<nrows<<endl;
+ // get the TPC parameters
+ rl->CdGAFile();
+ AliTPCParamSR* param = AliTPC::LoadTPCParam(gFile);
+ if (!param)
+ cout<<"No TPC parameter"<<endl;
+ AliTPCDigitsArray *digarr=new AliTPCDigitsArray;
+ digarr->Setup(param);
+ digarr->ConnectTree(digitsTree);