- prfinner->Read("prf_07504_Gati_056068_d02");
- prfouter->Read("prf_10006_Gati_047051_d03");
- f->Close();
- savedir->cd();
-
- param->SetInnerPRF(prfinner);
- param->SetOuterPRF(prfouter);
- param->SetTimeRF(rf);
- TPC->SetParam(param);
-
- cerr<<"Digitizing TPC...\n";
-
- //setup TPCDigitsArray
- AliTPCDigitsArray *arr = new AliTPCDigitsArray;
- arr->SetClass("AliSimDigits");
- arr->Setup(param);
- TPC->SetParam(param);
- arr->MakeTree();
-
- TPC->SetDigitsArray(arr);
- TPC->Hits2Digits();
- /* TPC->Hits2DigitsSector(1);
- TPC->Hits2DigitsSector(2);
- TPC->Hits2DigitsSector(3);
- TPC->Hits2DigitsSector(1+18);
- TPC->Hits2DigitsSector(2+18);
- TPC->Hits2DigitsSector(3+18);
-
- TPC->Hits2DigitsSector(36+1);
- TPC->Hits2DigitsSector(36+2);
- TPC->Hits2DigitsSector(36+3);
- TPC->Hits2DigitsSector(36+1+18);
- TPC->Hits2DigitsSector(36+2+18);
- TPC->Hits2DigitsSector(36+3+18); */
-
- //write results
-
- char treeName[100];
- sprintf(treeName,"TreeD_%s",param->GetTitle());
- TPC->GetDigitsArray()->GetTree()->Write(treeName);