TOF raw data: preliminary implementation and style changes

[u/mrichter/AliRoot.git] / TOF / AliTOFSDigit.cxx

diff --git a/TOF/AliTOFSDigit.cxx b/TOF/AliTOFSDigit.cxx
index 5d34990a4ec27a66034091f52d47a850b02c4225..ece39a8b5675eaa780e20bcec361b81c0b352271 100644 (file)
--- a/TOF/AliTOFSDigit.cxx
+++ b/TOF/AliTOFSDigit.cxx
@@ -13,6 +13,7 @@
  * provided "as is" without express or implied warranty.                  * 
  **************************************************************************/
 
+/* $Id$ */
 
 //_________________________________________________________________________
 //  TOF sdigit: member variables 
@@ -28,15 +29,14 @@
 //
 //*-- Authors: F. Pierella, A. Seganti, D. Vicinanza
 
-#include <iostream.h>
+#include <Riostream.h>
 #include "TArrayF.h"
 #include "TArrayI.h"
 
+#include "AliRun.h"
 #include "AliTOF.h"
+#include "AliTOFGeometry.h"
 #include "AliTOFSDigit.h"
-#include "AliTOFConstants.h"
-#include "AliRun.h"
-#include "AliMC.h"
 
 ClassImp(AliTOFSDigit)
 
@@ -54,10 +54,12 @@ ClassImp(AliTOFSDigit)
 
 ////////////////////////////////////////////////////////////////////////
 AliTOFSDigit::AliTOFSDigit(Int_t tracknum, Int_t *vol,Float_t *digit)
+:TObject()
 {
   //
   // Constructor of digit object
   //
+
   fSector = vol[0];
   fPlate  = vol[1];
   fStrip  = vol[2];
@@ -77,6 +79,7 @@ AliTOFSDigit::AliTOFSDigit(Int_t tracknum, Int_t *vol,Float_t *digit)
 
 ////////////////////////////////////////////////////////////////////////
 AliTOFSDigit::AliTOFSDigit(const AliTOFSDigit & digit)
+:TObject()
 {
   // 
   // copy ctor for AliTOFSDigit object
@@ -139,7 +142,7 @@ void AliTOFSDigit::Update(Float_t tdcbin, Int_t tdc, Int_t adc, Int_t track)
   //
   
   Int_t sameTime = -1;
-  Float_t tdcwindow=((Float_t)AliTOFConstants::fgkTimeDiff)/tdcbin;
+  Float_t tdcwindow=((Float_t)AliTOFGeometry::TimeDiff())/tdcbin;
   for (Int_t i = 0; i < fNDigits; i++) {
     if (TMath::Abs(tdc-fTdc->At(i)) < tdcwindow) {
       sameTime = i;
@@ -178,6 +181,67 @@ void AliTOFSDigit::Update(Float_t tdcbin, Int_t tdc, Int_t adc, Int_t track)
   }
   
 }
+
+////////////////////////////////////////////////////////////////////////
+void AliTOFSDigit::Update(AliTOFSDigit* sdig)
+{
+
+  //
+  // Perform the sum with sdig
+  //
+
+  // start loop on all sdig locations
+  Int_t nlocations=sdig->GetNDigits();
+
+  for (Int_t j = 0; j < nlocations; j++) {
+    Float_t tdcbin=50.; // [ps] hardwired for the time being
+    Int_t tdc=(Int_t)sdig->GetTdc(j);
+    Int_t adc=(Int_t)sdig->GetAdc(j);
+    // getting here only the first track number
+    Int_t track=GetTrack(j,0);
+    
+    
+    Int_t sameTime = -1;
+    Float_t tdcwindow=((Float_t)AliTOFGeometry::TimeDiff())/tdcbin;
+    for (Int_t i = 0; i < fNDigits; i++) {
+      if (TMath::Abs(tdc-fTdc->At(i)) < tdcwindow) {
+       sameTime = i;
+       break;
+      }
+    }
+    
+    if (sameTime >= 0) {
+      (*fAdc)[sameTime] += static_cast<Float_t>(adc);
+      // update track - find the first -1  value and replace it by the
+      // track number
+      for (Int_t iTrack=0; iTrack<kMAXDIGITS; iTrack++) {
+       if ((*fTracks)[sameTime*kMAXDIGITS+iTrack] == -1) {
+         (*fTracks)[sameTime*kMAXDIGITS+iTrack] = track;
+         break;
+       }
+       // write warning about many tracks going to this pad
+       if (iTrack == kMAXDIGITS) {
+         cerr<<"WARNING: AliTOFSDigit::Update  Many hits in the padhit"<<endl;
+         cerr<<"         ";
+         //    PrintPad();
+       }
+      }
+    } else {
+      // add new time slot
+      fNDigits++;
+      fTdc->Set(fNDigits);
+      (*fTdc)[fNDigits-1] = tdc;
+      fAdc->Set(fNDigits);
+      (*fAdc)[fNDigits-1] = adc;
+      fTracks->Set(fNDigits*kMAXDIGITS);
+      (*fTracks)[(fNDigits-1)*kMAXDIGITS] = track;
+      for (Int_t i = 1; i <kMAXDIGITS; i++) {
+       (*fTracks)[(fNDigits-1)*kMAXDIGITS+i] = -1;
+      } // for (Int_t i = 1; i <kMAXDIGITS; i++)
+    } // if (sameTime >= 0)
+  } // end loop on sdig locations
+}
+
 ////////////////////////////////////////////////////////////////////////
 AliTOFSDigit::~AliTOFSDigit()
 {
@@ -198,24 +262,30 @@ Int_t AliTOFSDigit::GetTotPad() const
   // starting from the digits data.
   //
   
-  Int_t pad = fPadx+AliTOFConstants::fgkNpadX*(fPadz-1);
+  Int_t pad = 2*fPadx + fPadz;
+  //Int_t pad = fPadx+AliTOFGeometry::NpadX()*fPadz;
   Int_t before=0;
   
   switch(fPlate){ 
-  case 1: before = 0;
+  case 0:
+    //before = 0;
     break;
-  case 2: before = AliTOFConstants::fgkNStripC;
+  case 1:
+    before = AliTOFGeometry::NStripC();
     break;
-  case 3: before = AliTOFConstants::fgkNStripB + AliTOFConstants::fgkNStripC;
+  case 2:
+    before = AliTOFGeometry::NStripB() +   AliTOFGeometry::NStripC();
     break;
-  case 4: before = AliTOFConstants::fgkNStripA + AliTOFConstants::fgkNStripB + AliTOFConstants::fgkNStripC;
+  case 3:
+    before = AliTOFGeometry::NStripA() +   AliTOFGeometry::NStripB() + AliTOFGeometry::NStripC();
     break;
-  case 5: before = AliTOFConstants::fgkNStripA + 2*AliTOFConstants::fgkNStripB + AliTOFConstants::fgkNStripC;
+  case 4:
+    before = AliTOFGeometry::NStripA() + 2*AliTOFGeometry::NStripB() + AliTOFGeometry::NStripC();
     break;
   }
   
-  Int_t strip = fStrip+before;
-  Int_t padTot = AliTOFConstants::fgkPadXStrip*(strip-1)+pad;
+  Int_t strip = fStrip + before;
+  Int_t padTot = AliTOFGeometry::NpadXStrip()*strip + pad;
   return padTot;
 }