Rename classes to read raw data 2004
[u/mrichter/AliRoot.git] / PHOS / AliPHOSReconstructor.cxx
CommitLineData
d15a28e7 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
b2a60966 16/* $Id$ */
17
d15a28e7 18//_________________________________________________________________________
a3dfe79c 19//*--
dfe0be07 20//*-- Yves Schutz (SUBATECH)
21// Reconstruction class. Redesigned from the old AliReconstructionner class and
22// derived from STEER/AliReconstructor.
23//
d15a28e7 24// --- ROOT system ---
25
d15a28e7 26// --- Standard library ---
364de5c6 27
d15a28e7 28// --- AliRoot header files ---
35293055 29#include "AliESD.h"
f444a19f 30#include "AliPHOSReconstructor.h"
7acf6008 31#include "AliPHOSClusterizerv1.h"
7acf6008 32#include "AliPHOSTrackSegmentMakerv1.h"
33#include "AliPHOSPIDv1.h"
35293055 34#include "AliPHOSGetter.h"
23904d16 35#include "AliPHOSTracker.h"
d22dd3b4 36#include "AliRawReader.h"
e957fea8 37
2e60107f 38
f444a19f 39ClassImp(AliPHOSReconstructor)
d15a28e7 40
2e60107f 41Bool_t AliPHOSReconstructor::fgDebug = kFALSE ;
42
d15a28e7 43//____________________________________________________________________________
2e60107f 44 AliPHOSReconstructor::AliPHOSReconstructor()
d15a28e7 45{
b2a60966 46 // ctor
0379a13e 47
6ad0bfa0 48}
49
0379a13e 50//____________________________________________________________________________
51 AliPHOSReconstructor::~AliPHOSReconstructor()
52{
53 // dtor
54
55}
7acf6008 56
7acf6008 57//____________________________________________________________________________
dfe0be07 58void AliPHOSReconstructor::Reconstruct(AliRunLoader* runLoader) const
35293055 59{
dfe0be07 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()) ;
9ee4fde2 66 TString branchName(runLoader->GetEventFolder()->GetName()) ;
a68156e6 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//____________________________________________________________________________
3255d660 78void AliPHOSReconstructor::Reconstruct(AliRunLoader* runLoader, AliRawReader* rawreader) const
a68156e6 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()) ;
35293055 89
dfe0be07 90 AliPHOSClusterizerv1 clu(headerFile, branchName);
d22dd3b4 91 clu.SetEventRange(0, 0) ; // reconstruct one events
92 clu.SetRawReader(rawreader);
dfe0be07 93 if ( Debug() )
94 clu.ExecuteTask("deb all") ;
95 else
d22dd3b4 96 clu.ExecuteTask("") ;
0379a13e 97
35293055 98}
7acf6008 99
7acf6008 100//____________________________________________________________________________
dfe0be07 101void AliPHOSReconstructor::FillESD(AliRunLoader* runLoader, AliESD* esd) const
7acf6008 102{
23904d16 103 // This function creates AliESDtracks from AliPHOSRecParticles
104 // and
105 // writes them to the ESD
bf72996e 106
23904d16 107 Int_t eventNumber = runLoader->GetEventNumber() ;
dfe0be07 108
0379a13e 109 AliPHOSGetter::Instance()->Event(eventNumber, "P") ;
110 TClonesArray *recParticles = AliPHOSGetter::Instance()->RecParticles();
dfe0be07 111 Int_t nOfRecParticles = recParticles->GetEntries();
8013c27e 112 esd->SetNumberOfPHOSParticles(nOfRecParticles) ;
113 esd->SetFirstPHOSParticle(esd->GetNumberOfTracks()) ;
114
dfe0be07 115 for (Int_t recpart = 0 ; recpart < nOfRecParticles ; recpart++) {
116 AliPHOSRecParticle * rp = dynamic_cast<AliPHOSRecParticle*>(recParticles->At(recpart));
117 if (Debug())
118 rp->Print();
119 AliESDtrack * et = new AliESDtrack() ;
120 // fills the ESDtrack
121 Double_t xyz[3];
8013c27e 122 for (Int_t ixyz=0; ixyz<3; ixyz++)
123 xyz[ixyz] = rp->GetPos()[ixyz];
dfe0be07 124 et->SetPHOSposition(xyz) ;
125 et->SetPHOSsignal (rp->Energy()) ;
126 et->SetPHOSpid (rp->GetPID()) ;
fd7be713 127 et->SetLabel (rp->GetPrimaryIndex());
dfe0be07 128 // add the track to the esd object
129 esd->AddTrack(et);
130 delete et;
7acf6008 131 }
d7d3b67b 132}
23904d16 133
134AliTracker* AliPHOSReconstructor::CreateTracker(AliRunLoader* runLoader) const
135{
136// creates the PHOS tracker
137 if (!runLoader) return NULL;
138 return new AliPHOSTracker(runLoader);
139}
140