- changes due to deletion of files

[u/mrichter/AliRoot.git] / TPC / AliTPCSDigits2Digits.C

diff --git a/TPC/AliTPCSDigits2Digits.C b/TPC/AliTPCSDigits2Digits.C
index 4064cddaf135dd7abc90a662c5dbd15196f7566d..528f8c5b634a3bf84b6c69c25c27616dee26802d 100644 (file)
--- a/TPC/AliTPCSDigits2Digits.C
+++ b/TPC/AliTPCSDigits2Digits.C
@@ -9,44 +9,82 @@ Int_t AliTPCSDigits2Digits(Int_t nevent=1)
 
   // Connect the Root Galice file containing Geometry, Kine and Hits
 
-  const char * inFile_old = "galice.root"; 
+  const char * inFile_old = "galice.root";
   const char * inFile_new = "galice.root";
+
+   if (gAlice)
+    {
+      delete AliRunLoader::Instance();
+      delete gAlice;//if everything was OK here it is already NULL
+      gAlice = 0x0;
+    }
+
   TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject(inFile_old);
   if (file) {file->Close(); delete file;}
-  file =  TFile::Open(inFile_new,"UPDATE");
-  if (!file->IsOpen()) {
+
+  AliRunLoader *rl = AliRunLoader::Open(inFile_new,"Event","update");
+
+  if (!rl) {
     cerr<<"Can't open "<<inFile_new<<" !\n";
     return 1;
   }
 
   // Get AliRun object from file or create it if not on file
   //  if (gAlice) delete gAlice;
-  gAlice = (AliRun*)file->Get("gAlice");
+
+  rl->LoadgAlice();
+
+  gAlice = rl->GetAliRun();
   if (!gAlice) {
     cerr<<"AliTPCHits2Digits.C : AliRun object not found on file\n";
+    delete rl;
     return 2;
   }
 
-
-
   // gAlice->GetEvent(0);
   AliTPC *TPC = (AliTPC*)gAlice->GetDetector("TPC");      
 
+  AliLoader * tpcl = rl->GetLoader("TPCLoader");
+  
+  if ((TPC == 0x0) || (tpcl == 0x0))
+   {
+    cerr<<"AliTPCHits2Digits.C : Can not find TPC or TPCLoader\n";
+//    delete rl;
+    return 3;
+   }
+
+  tpcl->LoadSDigits("READ");
+  tpcl->LoadDigits("RECREATE");
+  rl->CdGAFile();
+  AliTPCParamSR *dig=(AliTPCParamSR *)gDirectory->Get("75x40_100x60");
+   if(dig){
+     cerr<<"2 pad-length geom hits with 3 pad-lengths geom digits\n";
+     delete dig;
+     dig = new AliTPCParamSR();
+   }
+   else
+   {
+     dig=(AliTPCParamSR *)gDirectory->Get("75x40_100x60_150x60");
+   }
+   if (!dig) {cerr<<"TPC parameters have not been found !\n"; return 4;}
+   TPC->SetParam(dig);
+       
+
   TStopwatch timer;
   timer.Start();
 
   for(Int_t eventn =0;eventn<nevent;eventn++){
     printf("Processing event %d\n",eventn);
-    gAlice->GetEvent(eventn);
-
+    cout<<"rl->GetEvent(eventn);\n";
+    rl->GetEvent(eventn);
+    cout<<"TPC->SDigits2Digits2(eventn);\n";
     TPC->SDigits2Digits2(eventn);
   } 
 
-  delete gAlice; gAlice=0;
-  file->Close(); delete file;
   timer.Stop();
   timer.Print();
 
+  delete rl;
   return 0;
 };