Recreating rec. points (Yu.Belikov)
[u/mrichter/AliRoot.git] / TRD / AliTRDReconstructor.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 // class for TRD reconstruction                                              //
21 //                                                                           //
22 ///////////////////////////////////////////////////////////////////////////////
23
24
25 #include "AliTRDReconstructor.h"
26 #include "AliRunLoader.h"
27 #include "AliTRDparameter.h"
28 #include "AliTRD.h"
29 #include "AliTRDclusterizerV1.h"
30 #include "AliTRDtracker.h"
31 #include <TFile.h>
32
33
34 ClassImp(AliTRDReconstructor)
35
36
37 //_____________________________________________________________________________
38 void AliTRDReconstructor::Reconstruct(AliRunLoader* runLoader) const
39 {
40 // reconstruct clusters
41
42   AliLoader *loader=runLoader->GetLoader("TRDLoader");
43   loader->LoadRecPoints("recreate");
44
45   AliTRDclusterizerV1 clusterer("clusterer", "TRD clusterizer");
46   runLoader->CdGAFile();
47   AliTRDparameter* trdParam = GetTRDparameter(runLoader); 
48   if (!trdParam) {
49     Error("Reconstruct", "no TRD parameters found");
50     return;
51   }
52   trdParam->ReInit();
53   clusterer.SetParameter(trdParam);
54   Int_t nEvents = runLoader->GetNumberOfEvents();
55
56   for (Int_t iEvent = 0; iEvent < nEvents; iEvent++) {
57     clusterer.Open(runLoader->GetFileName(), iEvent);
58     clusterer.ReadDigits();
59     clusterer.MakeClusters();
60     clusterer.WriteClusters(-1);
61   }
62 }
63
64 //_____________________________________________________________________________
65 AliTracker* AliTRDReconstructor::CreateTracker(AliRunLoader* runLoader) const
66 {
67 // create a TRD tracker
68
69   runLoader->CdGAFile();
70   return new AliTRDtracker(gFile);
71 }
72
73 //_____________________________________________________________________________
74 void AliTRDReconstructor::FillESD(AliRunLoader* /*runLoader*/, 
75                                   AliESD* /*esd*/) const
76 {
77 }
78
79
80 //_____________________________________________________________________________
81 AliTRDparameter* AliTRDReconstructor::GetTRDparameter(AliRunLoader* runLoader) const
82 {
83 // get the TRD parameters
84
85   runLoader->CdGAFile();
86   AliTRDparameter* trdParam = (AliTRDparameter*) gFile->Get("TRDparameter"); 
87   if (!trdParam) {
88     Error("GetTRDparameter", "no TRD parameters available");
89     return NULL;
90   }
91   return trdParam;
92 }
93
94