]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - TOF/AliTOFRawStream.cxx
Adding the track fit residuals as a consequence of the ExB distortions (Marian)
[u/mrichter/AliRoot.git] / TOF / AliTOFRawStream.cxx
index d7b1364bc6989521e8d59483aa9a76310684c2e1..7648a59b017b9467894986570efbcd24d98dbcdd 100644 (file)
 
 /*
 $Log$
+Revision 1.19  2007/05/18 13:07:53  decaro
+Error messages stored in the global raw-reader error log (Cvetan, Chiara)
+
+Revision 1.18  2007/05/08 11:53:29  arcelli
+Improved class flexibility for further use (R.Preghenella)
+
 Revision 1.17  2007/05/03 08:53:50  decaro
 Coding convention: RS3 violation -> suppression
 
@@ -90,11 +96,11 @@ Revision 0.01  2005/07/22 A. De Caro
 #include "AliRawReader.h"
 
 #include "AliTOFGeometry.h"
-#include "AliTOFGeometryV5.h"
 #include "AliTOFrawData.h"
 #include "AliTOFRawMap.h"
 #include "AliTOFRawStream.h"
 
+#include "AliTOFHitData.h"
 
 ClassImp(AliTOFRawStream)
 
@@ -102,7 +108,6 @@ ClassImp(AliTOFRawStream)
 //_____________________________________________________________________________
 AliTOFRawStream::AliTOFRawStream(AliRawReader* rawReader):
   fRawReader(rawReader),
-  //  fTOFrawData(new TClonesArray("AliTOFrawData",1000)),
   fTOFrawData(0x0),
   fDecoder(new AliTOFDecoder()),
   fDDL(-1),
@@ -120,7 +125,6 @@ AliTOFRawStream::AliTOFRawStream(AliRawReader* rawReader):
   fStrip(-1),
   fPadX(-1),
   fPadZ(-1),
-  fTOFGeometry(new AliTOFGeometryV5()),
   fPackedDigits(0),
   fWordType(-1),
   fSlotID(-1),
@@ -152,7 +156,6 @@ AliTOFRawStream::AliTOFRawStream(AliRawReader* rawReader):
 //_____________________________________________________________________________
 AliTOFRawStream::AliTOFRawStream():
   fRawReader(0x0), 
-  //  fTOFrawData(new TClonesArray("AliTOFrawData",1000)),
   fTOFrawData(0x0),
   fDecoder(new AliTOFDecoder()),
   fDDL(-1),
@@ -170,7 +173,6 @@ AliTOFRawStream::AliTOFRawStream():
   fStrip(-1),
   fPadX(-1),
   fPadZ(-1),
-  fTOFGeometry(new AliTOFGeometryV5()),
   fPackedDigits(0),
   fWordType(-1),
   fSlotID(-1),
@@ -199,7 +201,6 @@ AliTOFRawStream::AliTOFRawStream():
 AliTOFRawStream::AliTOFRawStream(const AliTOFRawStream& stream) :
   TObject(stream),
   fRawReader(0x0),
-  //  fTOFrawData(new TClonesArray("AliTOFrawData",1000)),
   fTOFrawData(0x0),
   fDecoder(new AliTOFDecoder()),
   fDDL(-1),
@@ -217,7 +218,6 @@ AliTOFRawStream::AliTOFRawStream(const AliTOFRawStream& stream) :
   fStrip(-1),
   fPadX(-1),
   fPadZ(-1),
-  fTOFGeometry(new AliTOFGeometryV5()),
   fPackedDigits(0),
   fWordType(-1),
   fSlotID(-1),
@@ -256,8 +256,6 @@ AliTOFRawStream::AliTOFRawStream(const AliTOFRawStream& stream) :
   fPadX = stream.fPadX;
   fPadZ = stream.fPadZ;
 
-  fTOFGeometry = stream.fTOFGeometry;
-
   fPackedDigits = stream.fPackedDigits;
 
   fWordType = stream.fWordType;
@@ -308,8 +306,6 @@ AliTOFRawStream& AliTOFRawStream::operator = (const AliTOFRawStream& stream)
   fPadX = stream.fPadX;
   fPadZ = stream.fPadZ;
 
-  fTOFGeometry = stream.fTOFGeometry;
-
   fPackedDigits = stream.fPackedDigits;
 
   fWordType = stream.fWordType;
@@ -341,8 +337,6 @@ AliTOFRawStream::~AliTOFRawStream()
 
   fPackedDigits = 0;
 
-  delete fTOFGeometry;
   for (Int_t i=0;i<72;i++){ 
     delete fDataBuffer[i];
     delete fPackedDataBuffer[i];
@@ -359,15 +353,14 @@ AliTOFRawStream::~AliTOFRawStream()
 
 void AliTOFRawStream::LoadRawData(Int_t indexDDL)
 {
+  //
+  // To load raw data
+  //
 
-  /*
-    for (Int_t indexDDL = 0;
-    indexDDL < AliDAQ::NumberOfDdls("TOF");
-    indexDDL++) {
-  */
+  fTOFrawData->Clear();
+
+  TClonesArray &arrayTofRawData =  *fTOFrawData;
 
-  fTOFrawData = new TClonesArray("AliTOFrawData",1000); //potential memory leak
-    //fTOFrawData->Clear();
   fPackedDigits = 0;
 
   // create raw data map
@@ -381,13 +374,13 @@ void AliTOFRawStream::LoadRawData(Int_t indexDDL)
     
   Bool_t signal = kFALSE;
 
-  
+  AliTOFrawData *rawDigit = NULL;
 
   while(Next()) {
 
     signal = (fSector!=-1 && fPlate!=-1 && fStrip!=-1 && fPadZ!=-1 && fPadX!=-1);
     if (signal) {
-      //printf("  %2i  %1i  %2i  %1i  %2i  \n", fSector, fPlate, fStrip, fPadZ, fPadX);
+      AliDebug(2,Form("  %2i  %1i  %2i  %1i  %2i  \n", fSector, fPlate, fStrip, fPadZ, fPadX));
 
       slot[0] = fTRM;
       slot[1] = fTRMchain;
@@ -396,7 +389,7 @@ void AliTOFRawStream::LoadRawData(Int_t indexDDL)
 
       if (rawMap->TestHit(slot) != kEmpty) {
 
-       AliTOFrawData *rawDigit = static_cast<AliTOFrawData*>(rawMap->GetHit(slot));
+       rawDigit = static_cast<AliTOFrawData*>(rawMap->GetHit(slot));
 
        if (rawDigit->GetLeading()!=-1 && rawDigit->GetTrailing()==-1 &&
            fLeadingEdge==-1 && fTrailingEdge!=-1) {
@@ -406,24 +399,19 @@ void AliTOFRawStream::LoadRawData(Int_t indexDDL)
        else if ( ((rawDigit->GetTOF()!=-1 || rawDigit->GetLeading()!=-1 || rawDigit->GetTrailing()!=-1) &&
                   (fLeadingEdge!=-1 || fTrailingEdge!=-1 || fTime!=-1) )
 
-/*     else if ( ((rawDigit->GetTOF()!=-1 || rawDigit->GetLeading()!=-1 || rawDigit->GetTrailing()!=-1) &&
-                  (fTime!=-1 || fLeadingEdge!=-1 || fTrailingEdge!=-1))*/  /*||
-                 (rawDigit->GetLeading()==-1 && rawDigit->GetTrailing()!=-1 &&
-                 fLeadingEdge!=-1 && fTrailingEdge==-1) */)
+                 )
          {
 
-         TClonesArray &arrayTofRawData =  *fTOFrawData;
-         new (arrayTofRawData[fPackedDigits++]) AliTOFrawData(fTRM, fTRMchain, fTDC, fTDCchannel, fTime, fToT, fLeadingEdge, fTrailingEdge, fPSbit, fACQ, fErrorFlag);
+           new (arrayTofRawData[fPackedDigits++]) AliTOFrawData(fTRM, fTRMchain, fTDC, fTDCchannel, fTime, fToT, fLeadingEdge, fTrailingEdge, fPSbit, fACQ, fErrorFlag);
 
-       rawMap->SetHit(slot);
+           rawMap->SetHit(slot);
 
-       }
+         }
 
 
       }
       else {
 
-       TClonesArray &arrayTofRawData =  *fTOFrawData;
        new (arrayTofRawData[fPackedDigits++]) AliTOFrawData(fTRM, fTRMchain, fTDC, fTDCchannel, fTime, fToT, fLeadingEdge, fTrailingEdge, fPSbit, fACQ, fErrorFlag);
 
        rawMap->SetHit(slot);
@@ -434,34 +422,7 @@ void AliTOFRawStream::LoadRawData(Int_t indexDDL)
 
   } // closed -> while (Next())
 
-    /*
-      fDDL  = fRawReader->GetDDLID();
-
-      for (Int_t ii=0; ii<fTOFrawData->GetEntriesFast(); ii++) {
-
-      AliTOFrawData *dummy = (AliTOFrawData*)fTOFrawData->UncheckedAt(ii);
-
-      fTRM = dummy->GetTRM();
-      fTRMchain = dummy->GetTRMchain();
-      fTDC = dummy->GetTDC();
-      fTDCchannel = dummy->GetTDCchannel();
-
-      SetSector();
-      SetPlate();
-      SetStrip();
-      SetPadZ();
-      SetPadX();
-
-      printf("  %2i, %2i %1i, %2i, %1i  -->  %2i, %1i, %2i, %1i, %2i  \n",
-            fDDL, fTRM, fTRMchain, fTDC, fTDCchannel,
-            fSector, fPlate, fStrip, fPadZ, fPadX);
-
-            } // closed loop on TOF raw data TClonesArray
-    */
-
-    //} // closed loop on indexDDL
-
-
+  rawMap->Delete();
 
 }
 
@@ -1034,9 +995,10 @@ Int_t AliTOFRawStream::Equip2VolNpad(Int_t iDDL, Int_t iChain, Int_t nTDC,
   iPadAlongTheStrip = iTDClocal*AliTOFGeometry::NCh() + iCHlocal;
 
   if (((iDDL==1 || iDDL==2) && iPadAlongTheStrip< AliTOFGeometry::NpadX()) ||
-      ((iDDL==0 || iDDL==3) && iPadAlongTheStrip>=AliTOFGeometry::NpadX()))
-    AliError("Problems with the padX number!");
-
+      ((iDDL==0 || iDDL==3) && iPadAlongTheStrip>=AliTOFGeometry::NpadX())) {
+    fRawReader->AddMajorErrorLog(kPadXError);
+    AliWarning("Problems with the padX number!");
+  }
   return iPadAlongTheStrip;
 
 }
@@ -1099,14 +1061,22 @@ void AliTOFRawStream::EquipmentId2VolumeId(Int_t nDDL, Int_t nTRM, Int_t iChain,
   Int_t iSector = GetSectorNumber(nDDL);
 
   Int_t iPlate = Equip2VolNplate(iDDL, nTRM, nTDC);
-  if (iPlate==-1) AliError("Problems with the plate number!");
+  if (iPlate==-1) {
+    fRawReader->AddMajorErrorLog(kPlateError,"plate = -1");
+    AliWarning("Problems with the plate number!");
+  }
 
   Int_t iStrip = Equip2VolNstrip(iDDL, nTRM, nTDC);
-  if (iStrip==-1) AliError("Problems with the strip number!");
+  if (iStrip==-1) {
+    fRawReader->AddMajorErrorLog(kStripError,"strip = -1");
+    AliWarning("Problems with the strip number!");
+  }
 
   Int_t iPadAlongTheStrip  = Equip2VolNpad(iDDL, iChain, nTDC, iCH);
-  if (iPadAlongTheStrip==-1)
-    AliError("Problems with the pad number along the strip!");
+  if (iPadAlongTheStrip==-1){
+    fRawReader->AddMajorErrorLog(kPadAlongStripError,"pad = -1");
+    AliWarning("Problems with the pad number along the strip!");
+  }
 
   Int_t iPadX  = (Int_t)(iPadAlongTheStrip/(Float_t(AliTOFGeometry::NpadZ())));
   Int_t iPadZ  = iPadAlongTheStrip%AliTOFGeometry::NpadZ();
@@ -1118,80 +1088,38 @@ void AliTOFRawStream::EquipmentId2VolumeId(Int_t nDDL, Int_t nTRM, Int_t iChain,
   volume[4] = iPadZ;
 
 }
-//_____________________________________________________________________________
-
-Int_t AliTOFRawStream::GetIndex(Int_t *detId)
-{
-  //Retrieve calibration channel index
-  const Int_t kSectors = fTOFGeometry->NSectors();
-  const Int_t kPlates = fTOFGeometry->NPlates();
-  const Int_t kStripA = fTOFGeometry->NStripA();
-  const Int_t kStripB = fTOFGeometry->NStripB();
-  const Int_t kStripC = fTOFGeometry->NStripC();
-  const Int_t kPadX = fTOFGeometry->NpadX();
-  const Int_t kPadZ = fTOFGeometry->NpadZ();
-
-
-  Int_t isector = detId[0];
-  if (isector >= kSectors)
-    AliError(Form("Wrong sector number in TOF (%d) !",isector));
-  Int_t iplate = detId[1];
-  if (iplate >= kPlates)
-    AliError(Form("Wrong plate number in TOF (%d) !",iplate));
-  Int_t istrip = detId[2];
-  Int_t ipadz = detId[3];
-  Int_t ipadx = detId[4];
-  Int_t stripOffset = 0;
-  switch (iplate) {
-  case 0:
-    stripOffset = 0;
-    break;
-  case 1:
-    stripOffset = kStripC;
-    break;
-  case 2:
-    stripOffset = kStripC+kStripB;
-    break;
-  case 3:
-    stripOffset = kStripC+kStripB+kStripA;
-    break;
-  case 4:
-    stripOffset = kStripC+kStripB+kStripA+kStripB;
-    break;
-  default:
-    AliError(Form("Wrong plate number in TOF (%d) !",iplate));
-    break;
-  };
-
-  Int_t idet = ((2*(kStripC+kStripB)+kStripA)*kPadZ*kPadX)*isector +
-               (stripOffset*kPadZ*kPadX)+
-               (kPadZ*kPadX)*istrip+
-              (kPadX)*ipadz+
-               ipadx;
-  return idet;
-}
 //-----------------------------------------------------------------------------
-Bool_t AliTOFRawStream::DecodeDDL(Int_t DDLMin, Int_t DDLMax, Int_t verbose = 0){
+Bool_t AliTOFRawStream::DecodeDDL(Int_t nDDLMin, Int_t nDDLMax, Int_t verbose = 0) {
+  //
+  // To decode raw data for DDL number in [nDDLmin; nDDLmax]
+  //
+
   //check and fix valid DDL range
-  if (DDLMin < 0){
-    DDLMin = 0;
-    AliError("Wrong DDL range: setting first DDL ID to 0");
+  if (nDDLMin < 0){
+    nDDLMin = 0;
+    fRawReader->AddMinorErrorLog(kDDLMinError);
+    AliWarning("Wrong DDL range: setting first DDL ID to 0");
   }
-  if (DDLMax > 71){
-    DDLMax = 71;
-    AliError("Wrong DDL range: setting last DDL ID to 71");
+  if (nDDLMax > 71){
+    nDDLMax = 71;
+    fRawReader->AddMinorErrorLog(kDDLMaxError);
+    AliWarning("Wrong DDL range: setting last DDL ID to 71");
   }  
 
   //select required DDLs
-  fRawReader->Select("TOF", DDLMin, DDLMax);
+  fRawReader->Select("TOF", nDDLMin, nDDLMax);
 
   if (verbose)
-    AliInfo(Form("Selected TOF DDL range: %d-%d", DDLMin, DDLMax));
+    AliInfo(Form("Selected TOF DDL range: %d-%d", nDDLMin, nDDLMax));
 
   return(Decode(verbose));
 }
 //-----------------------------------------------------------------------------
-Bool_t AliTOFRawStream::Decode(Int_t verbose = 0){
+Bool_t AliTOFRawStream::Decode(Int_t verbose = 0) {
+  //
+  // New decoder method
+  //
+
   Int_t currentEquipment;
   Int_t currentDDL;
 
@@ -1221,8 +1149,9 @@ Bool_t AliTOFRawStream::Decode(Int_t verbose = 0){
     //read equipment payload
     if (!fRawReader->ReadNext(data, kDataSize))
       {
+       fRawReader->AddMajorErrorLog(kDDLdataReading);
        if (verbose)
-         AliError("Error while reading DDL data. Go to next equipment");
+         AliWarning("Error while reading DDL data. Go to next equipment");
        delete [] data;
        data = 0x0;
        continue;
@@ -1238,8 +1167,10 @@ Bool_t AliTOFRawStream::Decode(Int_t verbose = 0){
     fDecoder->SetPackedDataBuffer(fPackedDataBuffer[currentDDL]);
     
     //start decoding
-    if (fDecoder->Decode((UInt_t *)data, kDataWords) == kTRUE)
-      AliError(Form("Error while decoding DDL # %d: decoder returned with errors", currentDDL));
+    if (fDecoder->Decode((UInt_t *)data, kDataWords) == kTRUE) {
+      fRawReader->AddMajorErrorLog(kDDLDecoder,Form("DDL # = %d",currentDDL));
+      AliWarning(Form("Error while decoding DDL # %d: decoder returned with errors", currentDDL));
+    }
     
     delete [] data;
     data = 0x0;
@@ -1283,7 +1214,8 @@ AliTOFRawStream::LoadRawDataBuffers(Int_t indexDDL, Int_t verbose)
     AliInfo(Form("Decoding raw data for DDL # %d ...", indexDDL));
 
   if (DecodeDDL(indexDDL, indexDDL, verbose) != 0){ //decode DDL
-    AliError(Form("Error while decoding DDL # %d", indexDDL));
+    fRawReader->AddMajorErrorLog(kDDLDecoder,Form("DDL # = %d",indexDDL));
+    AliWarning(Form("Error while decoding DDL # %d", indexDDL));
     return kTRUE;
   }
   
@@ -1307,7 +1239,7 @@ AliTOFRawStream::LoadRawDataBuffers(Int_t indexDDL, Int_t verbose)
     Int_t   hitTimeBin = hitData->GetTimeBin();
     Int_t   hitTOTBin = hitData->GetTOTBin();
     
-    Int_t hitLeading = -1;
+    Int_t hitLeading = hitData->GetTimeBin();//-1; // adc
     Int_t hitTrailing = -1;
     Int_t hitError = -1;