macros
authorbnandi <bnandi@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 17 Jun 2004 09:14:39 +0000 (09:14 +0000)
committerbnandi <bnandi@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 17 Jun 2004 09:14:39 +0000 (09:14 +0000)
PMD/AliPMDDigits2RawData.C
PMD/AliPMDDigits2Recpoints.C

index fedd809..42963f5 100644 (file)
@@ -4,59 +4,11 @@
 //                                                     //
 // ----------------------------------------------------//
 
-#include <Riostream.h>
-#include "TBranch.h"
-#include "TStopwatch.h"
 
-extern AliRun *gAlice;
-
-Int_t AliPMDDigits2RawData(Int_t nevent = 1)
+Int_t AliPMDDigits2RawData()
 {
-  if (gAlice)
-    { 
-      delete gAlice->GetRunLoader();
-      delete gAlice;//if everything was OK here it is already NULL
-      gAlice = 0x0;
-    }
-
-  AliRunLoader *fRunLoader = AliRunLoader::Open("galice.root");
-  if (!fRunLoader)
-    {
-      cerr<<"Can't load RunLoader"<<endl;
-      return 1;
-    }
-  fRunLoader->LoadgAlice();
-  gAlice = fRunLoader->GetAliRun();
-
-  AliLoader *pmdloader = fRunLoader->GetLoader("PMDLoader");
-  //  Int_t nevent = fRunLoader->GetNumberOfEvents();
-  cout << " ************ nevent = " << nevent << endl;
-  
-  if (pmdloader == 0x0)
-    {
-      cerr<<" ===> Can not find PMD or PMDLoader <===\n";
-      delete fRunLoader;
-      return 2;
-    }
-  pmdloader->LoadDigits("READ");
-
-
-  for (Int_t ievt = 0; ievt < nevent; ievt++)
-    {
-
-      fRunLoader->GetEvent(ievt);
-      TTree *treeD = pmdloader->TreeD();
-      if (treeD == 0x0)
-       {
-         cout << " Can not get TreeD" << endl;
-         return 3;
-       }
-
-      AliPMDDDLRawData rawdata;
-      rawdata.WritePMDRawData(treeD, ievt);
-
-
-    }
+  AliSimulation sim;
+  sim.WriteRawData("PMD");
 
   return 0;
 }
index d4938e3..67c917d 100644 (file)
@@ -37,11 +37,18 @@ void AliPMDDigits2Recpoints(Int_t nevt=1)
   clus->SetDebug(1);
   clus->Load();
 
+
+
   for (Int_t ievt = 0; ievt < nevt; ievt++)
     {
-      clus->Digits2RecPoints(ievt);
+      // from digits data
+      //      clus->Digits2RecPoints(ievt);
+
+      // from raw data
+      AliRawReaderFile reader(ievt);
+      clus->Digits2RecPoints(ievt, &reader);
     }
-  clus->UnLoad("R");
+  clus->UnLoad();
 
   timer.Stop();
   timer.Print();