]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TRD/AliTRDdigitsManager.cxx
Separated TOF libraries (base,rec,sim)
[u/mrichter/AliRoot.git] / TRD / AliTRDdigitsManager.cxx
index 4f86fc0a6acf97fdfc9da2c322eab82f3b1f570d..5d1e88fc69833c74a654c0b6015cc2962e2f2709 100644 (file)
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
 
-/*
-$Log$
-Revision 1.12  2001/05/16 14:57:28  alibrary
-New files for folders and Stack
-
-Revision 1.11  2001/03/13 09:30:35  cblume
-Update of digitization. Moved digit branch definition to AliTRD
-
-Revision 1.10  2001/01/26 19:56:57  hristov
-Major upgrade of AliRoot code
-
-Revision 1.9  2000/11/02 09:25:53  cblume
-Change also the dictionary to AliTRDdataArray
-
-Revision 1.8  2000/11/01 15:20:13  cblume
-Change AliTRDdataArrayI to AliTRDdataArray in MakeBranch()
-
-Revision 1.7  2000/11/01 14:53:20  cblume
-Merge with TRD-develop
-
-Revision 1.1.2.5  2000/10/17 02:27:34  cblume
-Get rid of global constants
-
-Revision 1.1.2.4  2000/10/15 23:40:01  cblume
-Remove AliTRDconst
-
-Revision 1.1.2.3  2000/10/06 16:49:46  cblume
-Made Getters const
-
-Revision 1.1.2.2  2000/10/04 16:34:58  cblume
-Replace include files by forward declarations
-
-Revision 1.5  2000/06/09 11:10:07  cblume
-Compiler warnings and coding conventions, next round
-
-Revision 1.4  2000/06/08 18:32:58  cblume
-Make code compliant to coding conventions
-
-Revision 1.3  2000/06/07 16:27:01  cblume
-Try to remove compiler warnings on Sun and HP
-
-Revision 1.2  2000/05/08 16:17:27  cblume
-Merge TRD-develop
-
-Revision 1.1.2.1  2000/05/08 14:44:01  cblume
-Add new class AliTRDdigitsManager
-
-*/
+/* $Id$ */
 
 ///////////////////////////////////////////////////////////////////////////////
 //                                                                           //
@@ -69,7 +22,7 @@ Add new class AliTRDdigitsManager
 //                                                                           //
 ///////////////////////////////////////////////////////////////////////////////
 
-#include <iostream.h>
+#include <Riostream.h>
  
 #include <TROOT.h>
 #include <TTree.h>                                                              
@@ -82,7 +35,6 @@ Add new class AliTRDdigitsManager
 #include "AliTRDdataArrayI.h"
 #include "AliTRDdigit.h"
 #include "AliTRDgeometry.h"
-#include "AliTRD.h"
 
 ClassImp(AliTRDdigitsManager)
 
@@ -98,19 +50,22 @@ AliTRDdigitsManager::AliTRDdigitsManager():TObject()
   // Default constructor
   //
 
-  fIsRaw = kFALSE;
-
-  fDigits = new AliTRDsegmentArray("AliTRDdataArrayI",AliTRDgeometry::Ndet());
+  fIsRaw   = kFALSE;
+  fEvent   = 0;
+  fDebug   = 0;
+  fSDigits = 0;
 
+  fTree    = NULL;
+  fDigits  = NULL;
   for (Int_t iDict = 0; iDict < kNDict; iDict++) {
-    fDictionary[iDict] = new AliTRDsegmentArray("AliTRDdataArrayI"
-                                               ,AliTRDgeometry::Ndet());
+    fDictionary[iDict] = NULL;
   }
 
 }
 
 //_____________________________________________________________________________
 AliTRDdigitsManager::AliTRDdigitsManager(const AliTRDdigitsManager &m)
+:TObject(m)
 {
   //
   // AliTRDdigitsManager copy constructor
@@ -130,11 +85,13 @@ AliTRDdigitsManager::~AliTRDdigitsManager()
   if (fDigits) {
     fDigits->Delete();
     delete fDigits;
+    fDigits = NULL;
   }
 
   for (Int_t iDict = 0; iDict < kNDict; iDict++) {
     fDictionary[iDict]->Delete();
     delete fDictionary[iDict];
+    fDictionary[iDict] = NULL;
   }
 
 }
@@ -146,15 +103,58 @@ void AliTRDdigitsManager::Copy(TObject &m)
   // Copy function
   //
 
-  ((AliTRDdigitsManager &) m).fIsRaw = fIsRaw;
+  ((AliTRDdigitsManager &) m).fIsRaw   = fIsRaw;
+  ((AliTRDdigitsManager &) m).fEvent   = fEvent;
+  ((AliTRDdigitsManager &) m).fDebug   = fDebug;
+  ((AliTRDdigitsManager &) m).fSDigits = fSDigits;
 
   TObject::Copy(m);
 
 }
 
+//_____________________________________________________________________________
+void AliTRDdigitsManager::CreateArrays()
+{
+  //
+  // Create the data arrays
+  //
+
+  fDigits = new AliTRDsegmentArray("AliTRDdataArrayI",AliTRDgeometry::Ndet());
+
+  for (Int_t iDict = 0; iDict < kNDict; iDict++) {
+    fDictionary[iDict] = new AliTRDsegmentArray("AliTRDdataArrayI"
+                                               ,AliTRDgeometry::Ndet());
+  }
+
+}
+//_____________________________________________________________________________
+void AliTRDdigitsManager::ResetArrays()
+{
+  //
+  // Reset the data arrays
+  //
+
+  if (fDigits) {
+    delete fDigits;
+  }
+  fDigits = new AliTRDsegmentArray("AliTRDdataArrayI",AliTRDgeometry::Ndet());
+
+  for (Int_t iDict = 0; iDict < kNDict; iDict++) {
+    if (fDictionary[iDict]) {  
+      delete fDictionary[iDict];
+    }
+    fDictionary[iDict] = new AliTRDsegmentArray("AliTRDdataArrayI"
+                                               ,AliTRDgeometry::Ndet());
+  }
+
+}
+
 //_____________________________________________________________________________
 void AliTRDdigitsManager::SetRaw()
 {
+  //
+  // Switch on the raw digits flag
+  //
 
   fIsRaw = kTRUE;
 
@@ -175,67 +175,65 @@ Short_t AliTRDdigitsManager::GetDigitAmp(Int_t row, Int_t col,Int_t time
 }
  
 //_____________________________________________________________________________
-Bool_t AliTRDdigitsManager::MakeBranch(char *file)
+Bool_t AliTRDdigitsManager::MakeBranch(TTree *tree)
 {
   //
-  // Creates the branches for the digits and the dictionary in the digits tree
+  // Creates the tree and branches for the digits and the dictionary
   //
 
   Int_t buffersize = 64000;
 
   Bool_t status = kTRUE;
 
-  //TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("TRD.Digits.root");
-
-  AliTRD * TRD = (AliTRD *) gAlice->GetDetector("TRD") ;
-
-  if (gAlice->TreeD()) {
+  if (tree) {
+    fTree = tree;
+  }
 
-    // Make the branch for the digits
-    if (fDigits) {
-      const AliTRDdataArray *kDigits = 
-           (AliTRDdataArray *) fDigits->At(0);
-      if (kDigits) {
-           TRD->MakeBranchInTree(gAlice->TreeD(), 
-                                "TRDdigits", kDigits->IsA()->GetName(),
-                                &kDigits,buffersize, 1,file);
-        printf("AliTRDdigitsManager::MakeBranch -- ");
+  // Make the branch for the digits
+  if (fDigits) {
+    const AliTRDdataArray *kDigits = (AliTRDdataArray *) fDigits->At(0);
+    if (kDigits) {
+      if (!fTree) return kFALSE;
+      TBranch* branch = fTree->GetBranch("TRDdigits");
+      if (!branch) fTree->Branch("TRDdigits",kDigits->IsA()->GetName(),
+                                 &kDigits,buffersize,99);
+      if (fDebug > 0) {
+        printf("<AliTRDdigitsManager::MakeBranch> ");
         printf("Making branch TRDdigits\n");
       }
-      else {
-        status = kFALSE;
-      }
     }
     else {
       status = kFALSE;
     }
+  }
+  else {
+    status = kFALSE;
+  }
 
-    // Make the branches for the dictionaries
-    for (Int_t iDict = 0; iDict < kNDict; iDict++) {
-      Char_t branchname[15];
-      sprintf(branchname,"TRDdictionary%d",iDict);
-      if (fDictionary[iDict]) {
-        const AliTRDdataArray *kDictionary = 
-             (AliTRDdataArray *) fDictionary[iDict]->At(0);
-        if (kDictionary) {
-            TRD->MakeBranchInTree(gAlice->TreeD(), 
-                             branchname,kDictionary->IsA()->GetName(),
-                             &kDictionary,buffersize, 1,file);
-          printf("AliTRDdigitsManager::MakeBranch -- ");
+  // Make the branches for the dictionaries
+  for (Int_t iDict = 0; iDict < kNDict; iDict++) {
+    Char_t branchname[15];
+    sprintf(branchname,"TRDdictionary%d",iDict);
+    if (fDictionary[iDict]) {
+      const AliTRDdataArray *kDictionary = 
+              (AliTRDdataArray *) fDictionary[iDict]->At(0);
+      if (kDictionary) {
+       if (!fTree) return kFALSE;
+       TBranch* branch = fTree->GetBranch(branchname);
+       if (!branch) fTree->Branch(branchname,kDictionary->IsA()->GetName(),
+                                  &kDictionary,buffersize,99);
+        if (fDebug > 0) {
+          printf("<AliTRDdigitsManager::MakeBranch> ");
           printf("Making branch %s\n",branchname);
        }
-        else {
-          status = kFALSE;
-       }
       }
       else {
         status = kFALSE;
       }
     }
-
-  }
-  else {
-    status = kFALSE;
+    else {
+      status = kFALSE;
+    }
   }
 
   return status;
@@ -243,7 +241,7 @@ Bool_t AliTRDdigitsManager::MakeBranch(char *file)
 }
 
 //_____________________________________________________________________________
-Bool_t AliTRDdigitsManager::ReadDigits()
+Bool_t AliTRDdigitsManager::ReadDigits(TTree *tree)
 {
   //
   // Reads the digit information from the input file
@@ -251,12 +249,26 @@ Bool_t AliTRDdigitsManager::ReadDigits()
 
   Bool_t status = kTRUE;
 
-  status = fDigits->LoadArray("TRDdigits");
+  if (tree) {
+
+    fTree = tree;
+
+  }
+
+  if (!fDigits) {
+    if (fDebug > 0) {
+      printf("<AliTRDdigitsManager::ReadDigits> ");
+      printf("Create the data arrays.\n");
+    }
+    CreateArrays();
+  }
+
+  status = fDigits->LoadArray("TRDdigits",fTree);
 
   for (Int_t iDict = 0; iDict < kNDict; iDict++) {
     Char_t branchname[15];
     sprintf(branchname,"TRDdictionary%d",iDict);
-    status = fDictionary[iDict]->LoadArray(branchname);
+    status = fDictionary[iDict]->LoadArray(branchname,fTree);
   }  
 
   if (fDigits->TestBit(AliTRDdigit::RawDigit())) {
@@ -277,29 +289,26 @@ Bool_t AliTRDdigitsManager::WriteDigits()
   // Writes out the TRD-digits and the dictionaries
   //
 
-  // Create the branches
-  if (!(gAlice->TreeD()->GetBranch("TRDdigits"))) { 
-    printf("AliTRDdigitsManager::WriteDigits -- ");
-    printf("Call MakeBranch\n");
-    if (!MakeBranch()) return kFALSE;
-  }
-
   // Store the contents of the segment array in the tree
-  if (!fDigits->StoreArray("TRDdigits")) {
-    printf("AliTRDdigitsManager::WriteDigits -- ");
+  if (!fDigits->StoreArray("TRDdigits",fTree)) {
+    printf("<AliTRDdigitsManager::WriteDigits> ");
     printf("Error while storing digits in branch TRDdigits\n");
     return kFALSE;
   }
   for (Int_t iDict = 0; iDict < kNDict; iDict++) {
     Char_t branchname[15];
     sprintf(branchname,"TRDdictionary%d",iDict);
-    if (!fDictionary[iDict]->StoreArray(branchname)) {
-      printf("AliTRDdigitsManager::WriteDigits -- ");
+    if (!fDictionary[iDict]->StoreArray(branchname,fTree)) {
+      printf("<AliTRDdigitsManager::WriteDigits> ");
       printf("Error while storing dictionary in branch %s\n",branchname);
       return kFALSE;
     }
   }
 
+  // Write the new tree to the output file
+  //fTree->Write();
+  fTree->AutoSave();  // Modification by Jiri
+
   return kTRUE;
 
 }