PMD/AliPMDReconstruction.C

   1 // ----------------------------------------------------//
   2 //                                                     //
   3 //    This macro does Digits to Reconstructed Points   //
   4 //                                                     //
   5 // ----------------------------------------------------//
   6
   7 #include "Riostream.h"
   8 #include "TROOT.h"
   9 #include "TFile.h"
  10 #include "TNetFile.h"
  11 #include "TRandom.h"
  12 #include "TTree.h"
  13 #include "TBranch.h"
  14 #include "TClonesArray.h"
  15 #include "TStopwatch.h"
  16 #include <stdlib.h>
  17
  18 void AliPMDReconstruction()
  19 {
  20   TStopwatch timer;
  21   timer.Start();
  22
  23   // Open the AliRoot file
  24   AliRunLoader *fRunLoader = AliRunLoader::Open("galice.root");
  25   if (!fRunLoader)
  26     {
  27       cerr<<"Can't load RunLoader"<<endl;
  28       return 1;
  29     }
  30   fRunLoader->LoadgAlice();
  31   gAlice = fRunLoader->GetAliRun();
  32
  33   printf(" Do you want reconstruction from Digits file or RAW data \n");
  34   printf(" If RAW,    type 0 \n");
  35   printf(" If Digits, type 1 \n");
  36
  37   Int_t itype;
  38
  39   cin >> itype;
  40
  41   if (itype == 0)
  42     {
  43       AliRawReaderFile reader;
  44       AliPMDReconstructor pmdreco;
  45       pmdreco.Reconstruct(fRunLoader, &reader);
  46     }
  47   else if (itype == 1)
  48     {
  49       AliPMDReconstructor pmdreco;
  50       pmdreco.Reconstruct(fRunLoader);
  51     }
  52
  53   timer.Stop();
  54   timer.Print();
  55 }
  56