AliMiniHeader structure and executables addMiniHeader and printMiniHeader removed
[u/mrichter/AliRoot.git] / PHOS / AliPHOSReconstructor.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16 /* $Id$ */
17
18 //_________________________________________________________________________
19 //*--
20 //*-- Yves Schutz (SUBATECH) 
21 // Reconstruction class. Redesigned from the old AliReconstructionner class and 
22 // derived from STEER/AliReconstructor. 
23 // 
24 // --- ROOT system ---
25
26 // --- Standard library ---
27
28 // --- AliRoot header files ---
29 #include "AliESD.h"
30 #include "AliPHOSReconstructor.h"
31 #include "AliPHOSClusterizerv1.h"
32 #include "AliPHOSTrackSegmentMakerv1.h"
33 #include "AliPHOSPIDv1.h"
34 #include "AliPHOSGetter.h"
35
36
37 ClassImp(AliPHOSReconstructor)
38
39 //____________________________________________________________________________
40   AliPHOSReconstructor::AliPHOSReconstructor() : fDebug(kFALSE) 
41 {
42   // ctor
43
44
45
46 //____________________________________________________________________________
47 void AliPHOSReconstructor::Reconstruct(AliRunLoader* runLoader) const 
48 {
49   // method called by AliReconstruction; 
50   // Only the clusterization is performed,; the rest of the reconstruction is done in FillESD because the track
51   // segment maker needs access to the AliESD object to retrieve the tracks reconstructed by 
52   // the global tracking.
53  
54   TString headerFile(runLoader->GetFileName()) ; 
55   TString branchName("Default") ;  
56   
57   AliPHOSClusterizerv1 clu(headerFile, branchName);
58   clu.SetEventRange(0, -1) ; // do all the events
59   if ( Debug() ) 
60     clu.ExecuteTask("deb all") ; 
61   else 
62     clu.ExecuteTask("") ;  
63 }
64
65 //____________________________________________________________________________
66 void AliPHOSReconstructor::FillESD(AliRunLoader* runLoader, AliESD* esd) const
67 {
68   // Called by AliReconstruct after Reconstruct() and global tracking and vertxing 
69   //Creates the tracksegments and Recparticles
70   
71   TString headerFile(runLoader->GetFileName()) ; 
72   TString branchName("Default") ;  
73
74   AliPHOSTrackSegmentMakerv1 tsMaker(headerFile, branchName);
75   AliPHOSPIDv1 pID(headerFile, branchName);
76
77   AliPHOSGetter *gime = AliPHOSGetter::Instance( (runLoader->GetFileName()).Data() ) ;
78   Int_t eventNumber = gime->EventNumber();
79   // do current event; the loop over events is done by AliReconstruction::Run()
80   tsMaker.SetEventRange(eventNumber, eventNumber) ; 
81   pID.SetEventRange(eventNumber, eventNumber) ; 
82   if ( Debug() ) {
83    tsMaker.ExecuteTask("deb all") ;
84    pID.ExecuteTask("deb all") ;
85   }
86   else {
87     tsMaker.ExecuteTask("") ;
88     pID.ExecuteTask("") ;
89   }
90   
91   // Creates AliESDtrack from AliPHOSRecParticles 
92   gime->Event(eventNumber, "P") ; 
93   TClonesArray *recParticles = gime->RecParticles();
94   Int_t nOfRecParticles = recParticles->GetEntries();
95   for (Int_t recpart = 0 ; recpart < nOfRecParticles ; recpart++) {
96     AliPHOSRecParticle * rp = dynamic_cast<AliPHOSRecParticle*>(recParticles->At(recpart));
97     if (Debug()) 
98       rp->Print();
99     AliESDtrack * et = new AliESDtrack() ; 
100     // fills the ESDtrack
101     Double_t xyz[3];
102     for (Int_t ixyz=0; ixyz<3; ixyz++) xyz[ixyz] = rp->GetPos()[ixyz];
103     et->SetPHOSposition(xyz) ; 
104     et->SetPHOSsignal  (rp->Energy()) ; 
105     et->SetPHOSpid     (rp->GetPID()) ;
106     // add the track to the esd object
107     esd->AddTrack(et);
108     delete et;
109   }
110 }