coverity defects fixed
[u/mrichter/AliRoot.git] / PMD / AliPMDDigits2Recpoints.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 AliPMDDigits2Recpoints(Int_t nevt=1) 
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   Int_t itype;
37   cin >> itype;
38
39   // Create the PMD Cluster Finder 
40   AliPMDClusterFinder *clus = new AliPMDClusterFinder(fRunLoader);
41
42   if (itype == 1)
43     {
44       clus->Load();
45     }
46   else if (itype == 0)
47     {
48       clus->LoadClusters();
49     }  
50
51
52   for (Int_t ievt = 0; ievt < nevt; ievt++)
53     {
54       if (itype == 1)
55         {
56           // from digits data
57           clus->Digits2RecPoints(ievt);
58         }
59       else if (itype == 0)
60         {
61           // from raw data
62           AliRawReaderFile reader(ievt);
63           clus->Digits2RecPoints(ievt, &reader);
64         }
65     }
66   if (itype == 1)
67     {
68       clus->UnLoad();
69     }
70   else if (itype == 0)
71     {
72       clus->UnLoadClusters();
73     }
74   timer.Stop();
75   timer.Print();
76 }
77