1. Adding a possibility to read and reconstruct an old rcu formatted raw data. This...
[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 #include "AliPHOSTracker.h"
36 #include "AliRawReader.h"
37
38  
39 ClassImp(AliPHOSReconstructor)
40
41 Bool_t AliPHOSReconstructor::fgDebug = kFALSE ; 
42
43 //____________________________________________________________________________
44   AliPHOSReconstructor::AliPHOSReconstructor() 
45 {
46   // ctor
47
48
49
50 //____________________________________________________________________________
51   AliPHOSReconstructor::~AliPHOSReconstructor()
52 {
53   // dtor
54
55
56
57 //____________________________________________________________________________
58 void AliPHOSReconstructor::Reconstruct(AliRunLoader* runLoader) const
59 {
60   // method called by AliReconstruction; 
61   // Only the clusterization is performed,; the rest of the reconstruction is done in FillESD because the track
62   // segment maker needs access to the AliESD object to retrieve the tracks reconstructed by 
63   // the global tracking.
64  
65   TString headerFile(runLoader->GetFileName()) ; 
66   TString branchName(runLoader->GetEventFolder()->GetName()) ;  
67   
68   AliPHOSClusterizerv1 clu(headerFile, branchName);
69   clu.SetEventRange(0, -1) ; // do all the events
70   if ( Debug() ) 
71     clu.ExecuteTask("deb all") ; 
72   else 
73     clu.ExecuteTask("") ;  
74
75 }
76
77 //____________________________________________________________________________
78 void AliPHOSReconstructor::Reconstruct(AliRunLoader* runLoader, AliRawReader* rawreader) const
79 {
80   // method called by AliReconstruction; 
81   // Only the clusterization is performed,; the rest of the reconstruction is done in FillESD because the track
82   // segment maker needs access to the AliESD object to retrieve the tracks reconstructed by 
83   // the global tracking.
84   // Here we reconstruct from Raw Data
85
86   rawreader->Reset() ; 
87   TString headerFile(runLoader->GetFileName()) ; 
88   TString branchName(runLoader->GetEventFolder()->GetName()) ;  
89   
90   AliPHOSClusterizerv1 clu(headerFile, branchName);
91   clu.SetEventRange(0, -1) ; // do all the events
92   clu.SetRawReader(rawreader);
93
94   TString option = GetOption();
95   if (option.Contains("OldRCUFormat"))
96     clu.SetOldRCUFormat(kTRUE);
97
98   if ( Debug() ) 
99     clu.ExecuteTask("deb all") ; 
100   else 
101     clu.ExecuteTask("") ;
102
103 }
104
105 //____________________________________________________________________________
106 void AliPHOSReconstructor::FillESD(AliRunLoader* runLoader, AliESD* esd) const
107 {
108   // This function creates AliESDtracks from AliPHOSRecParticles
109   //         and
110   // writes them to the ESD
111
112   Int_t eventNumber = runLoader->GetEventNumber() ;
113
114   AliPHOSGetter::Instance()->Event(eventNumber, "P") ; 
115   TClonesArray *recParticles = AliPHOSGetter::Instance()->RecParticles();
116   Int_t nOfRecParticles = recParticles->GetEntries();
117   esd->SetNumberOfPHOSClusters(nOfRecParticles) ; 
118   esd->SetFirstPHOSCluster(esd->GetNumberOfCaloClusters()) ; 
119
120   for (Int_t recpart = 0 ; recpart < nOfRecParticles ; recpart++) {
121     AliPHOSRecParticle * rp = dynamic_cast<AliPHOSRecParticle*>(recParticles->At(recpart));
122     if (Debug()) 
123       rp->Print();
124     AliESDCaloCluster * ec = new AliESDCaloCluster() ; 
125 //     AliESDtrack * et = new AliESDtrack() ; 
126     // fills the ESDCaloCluster
127     Float_t xyz[3];
128     for (Int_t ixyz=0; ixyz<3; ixyz++) 
129       xyz[ixyz] = rp->GetPos()[ixyz];
130     ec->SetGlobalPosition(xyz);
131     ec->SetClusterEnergy(rp->Energy());
132     ec->SetPid          (rp->GetPID()) ;
133 //     ec->SetPrimaryIndex (rp->GetPrimaryIndex()); // for dry events only
134     // add the track to the esd object
135     esd->AddCaloCluster(ec);
136     delete ec;
137   }
138 }
139
140 AliTracker* AliPHOSReconstructor::CreateTracker(AliRunLoader* runLoader) const
141 {
142 // creates the PHOS tracker
143   if (!runLoader) return NULL; 
144   return new AliPHOSTracker(runLoader);
145 }
146