void Digits2RecParticles(){ gSystem->Setenv("CONFIG_SPLIT_FILE","1") ; TFile * fRootFile = TFile::Open("galice.root","update"); //========== Get AliRun object from file gAlice = (AliRun*) fRootFile->Get("gAlice") ; //=========== Get the PHOS object and associated geometry from the file AliPHOSv1 * fPHOS = (AliPHOSv1 *)gAlice->GetDetector("PHOS") ; AliPHOSGeometry * fGeom = AliPHOSGeometry::GetInstance( fPHOS->GetGeometry()->GetName(), fPHOS->GetGeometry()->GetTitle() ); AliPHOSIndexToObject * fObjGetter = AliPHOSIndexToObject::GetInstance(fPHOS) ; //========== Create the Clusterizer AliPHOSClusterizerv1 * fClu = new AliPHOSClusterizerv1() ; //========== Creates the track segment maker AliPHOSTrackSegmentMakerv1 * fTrs = new AliPHOSTrackSegmentMakerv1() ; // fTrs->UnsetUnfoldFlag() ; //========== Creates the particle identifier AliPHOSPIDv1* fPID = new AliPHOSPIDv1() ; fPID->SetShowerProfileCuts(0.3, 1.8, 0.3, 1.8 ) ; //========== Creates the Reconstructioner AliPHOSReconstructor * fRec = new AliPHOSReconstructor(fClu, fTrs, fPID) ; // fRec -> SetDebugReconstruction(kTRUE); gAlice->GetEvent(0) ; gAlice->TreeD()->GetEvent(0) ; if (gSystem->Getenv("CONFIG_SPLIT_FILE")) if(gAlice->TreeR() == 0) gAlice->MakeTree("R","Reco.PHOS.root"); else if(gAlice->TreeR() == 0) gAlice->MakeTree("R"); //=========== Do the reconstruction fPHOS->Reconstruction(fRec); }