+ gime->PostSDigits( title, GetTitle()) ;
+}
+
+//__________________________________________________________________
+void AliPHOSSDigitizer::SetSplitFile(const TString splitFileName)
+{
+ // Diverts the SDigits in a file separate from the hits file
+
+ TDirectory * cwd = gDirectory ;
+
+ if ( !(gAlice->GetTreeSFileName() == splitFileName) ) {
+ if (gAlice->GetTreeSFile() )
+ gAlice->GetTreeSFile()->Close() ;
+ }
+
+ fSplitFile = gAlice->InitTreeFile("S",splitFileName.Data());
+ fSplitFile->cd() ;
+ gAlice->Write(0, TObject::kOverwrite);
+
+ TTree *treeE = gAlice->TreeE();
+ if (!treeE) {
+ cerr << "ERROR: AliPHOSSDigitizer::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: AliPHOSSDigitizer::SetSPlitFile -> AliceGeom was not found in the input file "<<endl;
+ abort() ;
+ }
+ AliceGeom->Write(0, TObject::kOverwrite);
+
+ gAlice->MakeTree("S",fSplitFile);
+ cwd->cd() ;
+ cout << "INFO: AliPHOSSDigitizer::SetSPlitMode -> SDigits will be stored in " << splitFileName.Data() << endl ;
+