]> git.uio.no Git - u/mrichter/AliRoot.git/blame - TRD/AliTRDReconstructor.cxx
Coding convention + var inizialization corrected
[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"
3fe61b77 33#include "AliTRDclusterizer.h"
121a60bd 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
ca21baaa 53 AliInfo("Convert raw data digits into digit objects [RawReader -> Digit TTree]");
54
a8fbb161 55 AliTRDrawData rawData;
3a0f6479 56 rawReader->Reset();
57 rawReader->Select("TRD");
a8fbb161 58 AliTRDdigitsManager *manager = rawData.Raw2Digits(rawReader);
59 manager->MakeBranch(digitsTree);
60 manager->WriteDigits();
c0ab62ff 61 delete manager;
0ee00e25 62
121a60bd 63}
64
65//_____________________________________________________________________________
3bc9d03e 66void AliTRDReconstructor::Reconstruct(AliRunLoader *runLoader
67 , AliRawReader *rawReader) const
928e9fae 68{
3bc9d03e 69 //
70 // Reconstruct clusters
71 //
928e9fae 72
ca21baaa 73 AliInfo("Reconstruct TRD clusters from RAW data [RunLoader, RawReader]");
928e9fae 74
3bc9d03e 75 AliLoader *loader = runLoader->GetLoader("TRDLoader");
928e9fae 76 loader->LoadRecPoints("recreate");
77
928e9fae 78 runLoader->CdGAFile();
928e9fae 79 Int_t nEvents = runLoader->GetNumberOfEvents();
80
3a0f6479 81 rawReader->Reset();
82 rawReader->Select("TRD");
83
928e9fae 84 for (Int_t iEvent = 0; iEvent < nEvents; iEvent++) {
ca21baaa 85
928e9fae 86 if (!rawReader->NextEvent()) break;
ca21baaa 87
ca21baaa 88 // New (fast) cluster finder
3fe61b77 89 AliTRDclusterizer clusterer("clusterer","TRD clusterizer");
3bc9d03e 90 clusterer.Open(runLoader->GetFileName(),iEvent);
ca21baaa 91 clusterer.Raw2ClustersChamber(rawReader);
92
928e9fae 93 clusterer.WriteClusters(-1);
ca21baaa 94
928e9fae 95 }
96
97 loader->UnloadRecPoints();
0ee00e25 98
25ca55ce 99}
a8fbb161 100
25ca55ce 101//_____________________________________________________________________________
102void AliTRDReconstructor::Reconstruct(AliRawReader *rawReader
103 , TTree *clusterTree) const
104{
105 //
106 // Reconstruct clusters
107 //
108
ca21baaa 109 AliInfo("Reconstruct TRD clusters from RAW data [RawReader -> Cluster TTree]");
25ca55ce 110
3a0f6479 111 rawReader->Reset();
112 rawReader->Select("TRD");
113
ca21baaa 114 // New (fast) cluster finder
3fe61b77 115 AliTRDclusterizer clusterer("clusterer","TRD clusterizer");
25ca55ce 116 clusterer.OpenOutput(clusterTree);
ca21baaa 117 clusterer.SetAddLabels(kFALSE);
118 clusterer.Raw2ClustersChamber(rawReader);
d6e36038 119
25ca55ce 120}
121
122//_____________________________________________________________________________
a8fbb161 123void AliTRDReconstructor::Reconstruct(TTree *digitsTree
124 , TTree *clusterTree) const
25ca55ce 125{
126 //
127 // Reconstruct clusters
128 //
ca21baaa 129 AliInfo("Reconstruct TRD clusters from Digits [Digit TTree -> Cluster TTree]");
25ca55ce 130
3fe61b77 131 AliTRDclusterizer clusterer("clusterer","TRD clusterizer");
25ca55ce 132 clusterer.OpenOutput(clusterTree);
133 clusterer.ReadDigits(digitsTree);
134 clusterer.MakeClusters();
d6e36038 135
25ca55ce 136}
137
138//_____________________________________________________________________________
139void AliTRDReconstructor::Reconstruct(AliRunLoader *runLoader) const
140{
141 //
142 // Reconstruct clusters
143 //
144
ca21baaa 145 AliInfo("Reconstruct TRD clusters [AliRunLoader]");
25ca55ce 146 AliLoader *loader = runLoader->GetLoader("TRDLoader");
147 loader->LoadRecPoints("recreate");
148
149 runLoader->CdGAFile();
150 Int_t nEvents = runLoader->GetNumberOfEvents();
151
152 for (Int_t iEvent = 0; iEvent < nEvents; iEvent++) {
3fe61b77 153 AliTRDclusterizer clusterer("clusterer","TRD clusterizer");
25ca55ce 154 clusterer.Open(runLoader->GetFileName(),iEvent);
155 clusterer.ReadDigits();
156 clusterer.MakeClusters();
157 clusterer.WriteClusters(-1);
158 }
159
160 loader->UnloadRecPoints();
161
928e9fae 162}
163
164//_____________________________________________________________________________
3bc9d03e 165AliTracker *AliTRDReconstructor::CreateTracker(AliRunLoader *runLoader) const
121a60bd 166{
3bc9d03e 167 //
168 // Create a TRD tracker
169 //
121a60bd 170
171 runLoader->CdGAFile();
3bc9d03e 172
121a60bd 173 return new AliTRDtracker(gFile);
3bc9d03e 174
121a60bd 175}
176
25ca55ce 177//_____________________________________________________________________________
178void AliTRDReconstructor::FillESD(AliRunLoader* /*runLoader*/
179 , AliRawReader* /*rawReader*/
eee6e868 180 , AliESDEvent* /*esd*/) const
25ca55ce 181{
182 //
af26ce80 183 // Fill ESD
25ca55ce 184 //
185
25ca55ce 186}
187
188//_____________________________________________________________________________
189void AliTRDReconstructor::FillESD(AliRawReader* /*rawReader*/
190 , TTree* /*clusterTree*/
eee6e868 191 , AliESDEvent* /*esd*/) const
25ca55ce 192{
193 //
af26ce80 194 // Fill ESD
25ca55ce 195 //
196
25ca55ce 197}
198
199//_____________________________________________________________________________
200void AliTRDReconstructor::FillESD(TTree* /*digitsTree*/
201 , TTree* /*clusterTree*/
ca21baaa 202 , AliESDEvent* /*esd*/) const
25ca55ce 203{
204 //
af26ce80 205 // Fill ESD
25ca55ce 206 //
207
25ca55ce 208}
209
121a60bd 210//_____________________________________________________________________________
36698bd7 211void AliTRDReconstructor::FillESD(AliRunLoader* /*runLoader*/
eee6e868 212 , AliESDEvent* /*esd*/) const
121a60bd 213{
3bc9d03e 214 //
af26ce80 215 // Fill ESD
3bc9d03e 216 //
b0f03c34 217
121a60bd 218}