Adding the array of the recosntruction parameters (Marian)
[u/mrichter/AliRoot.git] / PHOS / Digits2RecParticles.C
1 void Digits2RecParticles(){
2   gSystem->Setenv("CONFIG_SPLIT_FILE","1") ;   
3
4   TFile * fRootFile = TFile::Open("galice.root","update");
5
6   //========== Get AliRun object from file 
7   gAlice = (AliRun*) fRootFile->Get("gAlice") ;
8
9   //=========== Get the PHOS object and associated geometry from the file      
10   AliPHOSv1 * fPHOS  = (AliPHOSv1 *)gAlice->GetDetector("PHOS") ;
11   AliPHOSGeometry * fGeom  = AliPHOSGeometry::GetInstance( fPHOS->GetGeometry()->GetName(), fPHOS->GetGeometry()->GetTitle() );
12  
13   AliPHOSIndexToObject * fObjGetter = AliPHOSIndexToObject::GetInstance(fPHOS) ; 
14
15
16   
17   //========== Create the Clusterizer
18   AliPHOSClusterizerv1 * fClu = new AliPHOSClusterizerv1() ; 
19   
20   //========== Creates the track segment maker
21   AliPHOSTrackSegmentMakerv1 * fTrs = new AliPHOSTrackSegmentMakerv1()  ;
22   //      fTrs->UnsetUnfoldFlag() ; 
23   
24   //========== Creates the particle identifier
25   AliPHOSPIDv1* fPID = new AliPHOSPIDv1() ;
26   fPID->SetShowerProfileCuts(0.3, 1.8, 0.3, 1.8 ) ;       
27       
28   //========== Creates the Reconstructioner
29   AliPHOSReconstructor * fRec = new AliPHOSReconstructor(fClu, fTrs, fPID) ; 
30   // fRec -> SetDebugReconstruction(kTRUE);     
31      
32   gAlice->GetEvent(0) ;
33
34   gAlice->TreeD()->GetEvent(0) ;
35
36   if (gSystem->Getenv("CONFIG_SPLIT_FILE"))    
37    if(gAlice->TreeR() == 0) gAlice->MakeTree("R","Reco.PHOS.root");
38   else 
39    if(gAlice->TreeR() == 0) gAlice->MakeTree("R");    
40
41   //=========== Do the reconstruction
42   fPHOS->Reconstruction(fRec);    
43   
44
45 }