]> git.uio.no Git - u/mrichter/AliRoot.git/blame - TRD/AliTRDReconstructor.cxx
Change in AliESDEvent fwd declaration (Ch. Klein Boesing)
[u/mrichter/AliRoot.git] / TRD / AliTRDReconstructor.cxx
CommitLineData
121a60bd 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// //
4fad09c9 20// Class for TRD reconstruction //
121a60bd 21// //
22///////////////////////////////////////////////////////////////////////////////
23
e3b2b5e5 24#include <TFile.h>
121a60bd 25
121a60bd 26#include "AliRunLoader.h"
3bc9d03e 27#include "AliRawReader.h"
28#include "AliLog.h"
3bc9d03e 29#include "AliESDTrdTrack.h"
af885e0f 30#include "AliESDEvent.h"
3bc9d03e 31
32#include "AliTRDReconstructor.h"
121a60bd 33#include "AliTRDclusterizerV1.h"
34#include "AliTRDtracker.h"
b0f03c34 35#include "AliTRDpidESD.h"
e3b2b5e5 36#include "AliTRDgtuTrack.h"
a8fbb161 37#include "AliTRDrawData.h"
38#include "AliTRDdigitsManager.h"
121a60bd 39
40ClassImp(AliTRDReconstructor)
41
3bc9d03e 42Bool_t AliTRDReconstructor::fgkSeedingOn = kFALSE;
43Int_t AliTRDReconstructor::fgStreamLevel = 0; // Stream (debug) level
0397ce53 44
121a60bd 45//_____________________________________________________________________________
a8fbb161 46void AliTRDReconstructor::ConvertDigits(AliRawReader *rawReader
47 , TTree *digitsTree) const
121a60bd 48{
3bc9d03e 49 //
25ca55ce 50 // Convert raw data digits into digit objects in a root tree
3bc9d03e 51 //
121a60bd 52
a8fbb161 53 AliTRDrawData rawData;
3a0f6479 54 rawReader->Reset();
55 rawReader->Select("TRD");
a8fbb161 56 AliTRDdigitsManager *manager = rawData.Raw2Digits(rawReader);
57 manager->MakeBranch(digitsTree);
58 manager->WriteDigits();
0ee00e25 59
121a60bd 60}
61
62//_____________________________________________________________________________
3bc9d03e 63void AliTRDReconstructor::Reconstruct(AliRunLoader *runLoader
64 , AliRawReader *rawReader) const
928e9fae 65{
3bc9d03e 66 //
67 // Reconstruct clusters
68 //
928e9fae 69
70 AliInfo("Reconstruct TRD clusters from RAW data");
71
3bc9d03e 72 AliLoader *loader = runLoader->GetLoader("TRDLoader");
928e9fae 73 loader->LoadRecPoints("recreate");
74
928e9fae 75 runLoader->CdGAFile();
928e9fae 76 Int_t nEvents = runLoader->GetNumberOfEvents();
77
3a0f6479 78 rawReader->Reset();
79 rawReader->Select("TRD");
80
928e9fae 81 for (Int_t iEvent = 0; iEvent < nEvents; iEvent++) {
82 if (!rawReader->NextEvent()) break;
3bc9d03e 83 AliTRDclusterizerV1 clusterer("clusterer","TRD clusterizer");
84 clusterer.Open(runLoader->GetFileName(),iEvent);
928e9fae 85 clusterer.ReadDigits(rawReader);
86 clusterer.MakeClusters();
87 clusterer.WriteClusters(-1);
88 }
89
90 loader->UnloadRecPoints();
0ee00e25 91
25ca55ce 92}
a8fbb161 93
25ca55ce 94//_____________________________________________________________________________
95void AliTRDReconstructor::Reconstruct(AliRawReader *rawReader
96 , TTree *clusterTree) const
97{
98 //
99 // Reconstruct clusters
100 //
101
102 AliInfo("Reconstruct TRD clusters from RAW data");
103
3a0f6479 104 rawReader->Reset();
105 rawReader->Select("TRD");
106
25ca55ce 107 AliTRDclusterizerV1 clusterer("clusterer","TRD clusterizer");
108 clusterer.OpenOutput(clusterTree);
109 clusterer.ReadDigits(rawReader);
110 clusterer.MakeClusters();
25ca55ce 111
25ca55ce 112}
113
114//_____________________________________________________________________________
a8fbb161 115void AliTRDReconstructor::Reconstruct(TTree *digitsTree
116 , TTree *clusterTree) const
25ca55ce 117{
118 //
119 // Reconstruct clusters
120 //
121
122 AliTRDclusterizerV1 clusterer("clusterer","TRD clusterizer");
123 clusterer.OpenOutput(clusterTree);
124 clusterer.ReadDigits(digitsTree);
125 clusterer.MakeClusters();
25ca55ce 126
25ca55ce 127}
128
129//_____________________________________________________________________________
130void AliTRDReconstructor::Reconstruct(AliRunLoader *runLoader) const
131{
132 //
133 // Reconstruct clusters
134 //
135
136 AliLoader *loader = runLoader->GetLoader("TRDLoader");
137 loader->LoadRecPoints("recreate");
138
139 runLoader->CdGAFile();
140 Int_t nEvents = runLoader->GetNumberOfEvents();
141
142 for (Int_t iEvent = 0; iEvent < nEvents; iEvent++) {
143 AliTRDclusterizerV1 clusterer("clusterer","TRD clusterizer");
144 clusterer.Open(runLoader->GetFileName(),iEvent);
145 clusterer.ReadDigits();
146 clusterer.MakeClusters();
147 clusterer.WriteClusters(-1);
148 }
149
150 loader->UnloadRecPoints();
151
928e9fae 152}
153
154//_____________________________________________________________________________
3bc9d03e 155AliTracker *AliTRDReconstructor::CreateTracker(AliRunLoader *runLoader) const
121a60bd 156{
3bc9d03e 157 //
158 // Create a TRD tracker
159 //
121a60bd 160
161 runLoader->CdGAFile();
3bc9d03e 162
121a60bd 163 return new AliTRDtracker(gFile);
3bc9d03e 164
121a60bd 165}
166
25ca55ce 167//_____________________________________________________________________________
168void AliTRDReconstructor::FillESD(AliRunLoader* /*runLoader*/
169 , AliRawReader* /*rawReader*/
eee6e868 170 , AliESDEvent* /*esd*/) const
25ca55ce 171{
172 //
173 // Make PID
174 //
175
eee6e868 176 //AliTRDpidESD trdPID;
177 //trdPID.MakePID(esd);
25ca55ce 178
25ca55ce 179}
180
181//_____________________________________________________________________________
182void AliTRDReconstructor::FillESD(AliRawReader* /*rawReader*/
183 , TTree* /*clusterTree*/
eee6e868 184 , AliESDEvent* /*esd*/) const
25ca55ce 185{
186 //
187 // Make PID
188 //
189
eee6e868 190 //AliTRDpidESD trdPID;
191 //trdPID.MakePID(esd);
25ca55ce 192
25ca55ce 193}
194
195//_____________________________________________________________________________
196void AliTRDReconstructor::FillESD(TTree* /*digitsTree*/
197 , TTree* /*clusterTree*/
eee6e868 198 , AliESDEvent* /*esd*/) const
25ca55ce 199{
200 //
201 // Make PID
202 //
203
eee6e868 204 //AliTRDpidESD trdPID;
205 //trdPID.MakePID(esd);
25ca55ce 206
25ca55ce 207}
208
121a60bd 209//_____________________________________________________________________________
36698bd7 210void AliTRDReconstructor::FillESD(AliRunLoader* /*runLoader*/
eee6e868 211 , AliESDEvent* /*esd*/) const
121a60bd 212{
3bc9d03e 213 //
214 // Make PID
215 //
b0f03c34 216
eee6e868 217 //AliTRDpidESD trdPID;
218 //trdPID.MakePID(esd);
0ee00e25 219
121a60bd 220}