-// //____________________________________________________________________________
-// void AliPHOSPID::SetSplitFile(const TString splitFileName) const
-// {
-// // Diverts the Digits in a file separate from the hits file
-
-
-// TDirectory * cwd = gDirectory ;
-// TFile * splitFile = gAlice->InitTreeFile("R",splitFileName.Data());
-// splitFile->cd() ;
-// gAlice->Write(0, TObject::kOverwrite);
-
-// TTree *treeE = gAlice->TreeE();
-// if (!treeE) {
-// cerr << "ERROR: AliPHOSPID::SetSplitFile -> No TreeE found "<<endl;
-// abort() ;
-// }
-
-// // copy TreeE
-// AliHeader *header = new AliHeader();
-// treeE->SetBranchAddress("Header", &header);
-// treeE->SetBranchStatus("*",1);
-// TTree *treeENew = treeE->CloneTree();
-// treeENew->Write(0, TObject::kOverwrite);
-
-// // copy AliceGeom
-// TGeometry *AliceGeom = static_cast<TGeometry*>(cwd->Get("AliceGeom"));
-// if (!AliceGeom) {
-// cerr << "ERROR: AliPHOSPID::SetSplitFile -> AliceGeom was not found in the input file "<<endl;
-// abort() ;
-// }
-// AliceGeom->Write(0, TObject::kOverwrite) ;
-
-// gAlice->MakeTree("R",splitFile);
-// cwd->cd() ;
-// cout << "INFO: AliPHOSPID::SetSPlitMode -> RecParticles will be stored in " << splitFileName.Data() << endl ;
-// }
+//____________________________________________________________________________
+void AliPHOSPID::SetInput(TTree *clustersTree, TClonesArray *trackSegments)
+{
+ // Read the clusters tree and creates the
+ // arrays with the EMC and CPV
+ // clusters.
+ // and set the corresponding branch addresses
+
+ fTrackSegments = trackSegments;
+
+ TBranch *emcbranch = clustersTree->GetBranch("PHOSEmcRP");
+ if (!emcbranch) {
+ AliError("can't get the branch with the PHOS EMC clusters !");
+ return;
+ }
+ fEMCRecPoints = new TObjArray(100) ;
+ emcbranch->SetAddress(&fEMCRecPoints);
+ emcbranch->GetEntry(0);
+
+ TBranch *cpvbranch = clustersTree->GetBranch("PHOSCpvRP");
+ if (!cpvbranch) {
+ AliError("can't get the branch with the PHOS CPV clusters !");
+ return;
+ }
+ fCPVRecPoints = new TObjArray(100) ;
+ cpvbranch->SetAddress(&fCPVRecPoints);
+ cpvbranch->GetEntry(0);
+}