From f0e88078adff245637669df1c4a3e4839c71fe3d Mon Sep 17 00:00:00 2001
From: kowal2 <kowal2@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Thu, 4 Apr 2002 16:27:23 +0000
Subject: [PATCH] New test macro by J.Chudoba

---
 TPC/AliTPCHits2SDigitsJCH.C | 60 +++++++++++++++++++++++++++++++++++++
 1 file changed, 60 insertions(+)
 create mode 100644 TPC/AliTPCHits2SDigitsJCH.C

diff --git a/TPC/AliTPCHits2SDigitsJCH.C b/TPC/AliTPCHits2SDigitsJCH.C
new file mode 100644
index 00000000000..bea5afeebc7
--- /dev/null
+++ b/TPC/AliTPCHits2SDigitsJCH.C
@@ -0,0 +1,60 @@
+#if !defined(__CINT__) || defined(__MAKECINT__)
+#include "iostream.h"
+#include "TFile.h"
+#include "TTree.h"
+#include "AliRun.h"
+#include "AliTPC.h"
+#endif
+
+
+Int_t AliTPCHits2SDigitsJCH(TString fileNameSDigits="sdigits.root", TString fileNameHits="rfio:galice.root", Int_t nEvents = 1, Int_t firstEvent = 0)
+{
+
+  TString fileMode = "read";
+  if (fileNameSDigits == fileNameHits || fileNameSDigits == "") fileMode = "update";
+  TFile *fileHits =  TFile::Open(fileNameHits.Data(),fileMode.Data());
+  if (!fileHits->IsOpen()) {
+    cerr<<"Can't open "<<fileNameHits.Data()<<" !\n";
+    return 1;
+  }
+
+  // Get AliRun object from file or create it if not on file
+  if (gAlice) delete gAlice; gAlice = 0;
+  gAlice = (AliRun*)fileHits->Get("gAlice");
+  if (!gAlice) {
+    cerr<<"AliTPCHits2Digits.C : AliRun object not found on file\n";
+    return 2;
+  }
+
+  AliTPC *TPC = (AliTPC*)gAlice->GetDetector("TPC");      
+
+  TStopwatch timer;
+  timer.Start();
+
+  if (fileNameSDigits != fileNameHits && fileNameSDigits != "")  gAlice->InitTreeFile("S",fileNameSDigits.Data());
+
+  // uncomment lines below to set active sectors 
+  //Int_t sec[10]={4,5,6,7,8,4+36,5+36,6+36,7+36,8+36};
+  //TPC->SetActiveSectors(sec,10);
+
+  for(Int_t iEvent = firstEvent;iEvent<firstEvent+nEvents;iEvent++){
+    printf("Processing event %d \n",iEvent);
+    gAlice->GetEvent(iEvent);
+//    TPC->SetActiveSectors(1); // all sectors set active
+    printf("\nActive sectors\n");
+    for (Int_t i=0;i<72;i++) {
+      if (i%10 == 0) printf("\n");
+      if (TPC->IsSectorActive(i)) printf("%2d ",i);
+    }
+    printf("\n");
+    TPC->Hits2SDigits2(iEvent);
+  } 
+
+  timer.Stop();
+  timer.Print();
+  delete gAlice; gAlice=0;
+  fileHits->Close(); delete fileHits;
+
+  return 0;
+};
+
-- 
2.39.3