PHOS/Digits2RecParticles.C

   1 void Digits2RecParticles(){
   2
   3   TFile * fRootFile = TFile::Open("galice.root","update");
   4
   5   //========== Get AliRun object from file
   6   gAlice = (AliRun*) fRootFile->Get("gAlice") ;
   7
   8   //=========== Get the PHOS object and associated geometry from the file
   9   AliPHOSv1 * fPHOS  = (AliPHOSv1 *)gAlice->GetDetector("PHOS") ;
  10   AliPHOSGeometry * fGeom  = AliPHOSGeometry::GetInstance( fPHOS->GetGeometry()->GetName(), fPHOS->GetGeometry()->GetTitle() );
  11
  12   AliPHOSIndexToObject * fObjGetter = AliPHOSIndexToObject::GetInstance(fPHOS) ;
  13
  14
  15
  16   //========== Create the Clusterizer
  17   AliPHOSClusterizerv1 * fClu = new AliPHOSClusterizerv1() ;
  18
  19   //========== Creates the track segment maker
  20   AliPHOSTrackSegmentMakerv1 * fTrs = new AliPHOSTrackSegmentMakerv1()  ;
  21   //      fTrs->UnsetUnfoldFlag() ;
  22
  23   //========== Creates the particle identifier
  24   AliPHOSPIDv1* fPID = new AliPHOSPIDv1() ;
  25   fPID->SetShowerProfileCuts(0.3, 1.8, 0.3, 1.8 ) ;
  26
  27   //========== Creates the Reconstructioner
  28   AliPHOSReconstructioner * fRec = new AliPHOSReconstructioner(fClu, fTrs, fPID) ;
  29   // fRec -> SetDebugReconstruction(kTRUE);
  30
  31   gAlice->GetEvent(0) ;
  32
  33   gAlice->TreeD()->GetEvent(0) ;
  34
  35   //=========== Do the reconstruction
  36   fPHOS->Reconstruction(fRec);
  37
  38
  39 }