Introduce extern for the common block references

[u/mrichter/AliRoot.git] / PMD / AliPMDDigits2Recpoints.C

diff --git a/PMD/AliPMDDigits2Recpoints.C b/PMD/AliPMDDigits2Recpoints.C
index d4938e35a26760858641847ee801723fdcb62a43..fe7d43a14346c03283a3d32b0521250c1a91def2 100644 (file)
--- a/PMD/AliPMDDigits2Recpoints.C
+++ b/PMD/AliPMDDigits2Recpoints.C
@@ -30,19 +30,47 @@ void AliPMDDigits2Recpoints(Int_t nevt=1)
   fRunLoader->LoadgAlice();
   gAlice = fRunLoader->GetAliRun();
 
+  printf(" Do you want reconstruction from Digits file or RAW data \n");
+  printf(" If RAW,    type 0 \n");
+  printf(" If Digits, type 1 \n");
+  Int_t itype;
+  cin >> itype;
 
   // Create the PMD Cluster Finder 
   AliPMDClusterFinder *clus = new AliPMDClusterFinder(fRunLoader);
-  
   clus->SetDebug(1);
-  clus->Load();
+  if (itype == 0)
+    {
+      clus->Load();
+    }
+  else if (itype == 1)
+    {
+      clus->LoadClusters();
+    }  
+
 
   for (Int_t ievt = 0; ievt < nevt; ievt++)
     {
-      clus->Digits2RecPoints(ievt);
+      if (itype == 0)
+       {
+         // from digits data
+         clus->Digits2RecPoints(ievt);
+       }
+      else if (itype == 1)
+       {
+         // from raw data
+         AliRawReaderFile reader(ievt);
+         clus->Digits2RecPoints(ievt, &reader);
+       }
+    }
+  if (itype == 0)
+    {
+      clus->UnLoad();
+    }
+  else if (itype == 1)
+    {
+      clus->UnLoadClusters();
     }
-  clus->UnLoad("R");
-
   timer.Stop();
   timer.Print();
 }