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