Changes to offer an example of the new simplified interface.
authorivana <ivana@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 27 Apr 2007 14:39:52 +0000 (14:39 +0000)
committerivana <ivana@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 27 Apr 2007 14:39:52 +0000 (14:39 +0000)
The old one is still there though (as it might be needed in some cases)
(Laurent)

MUON/MUONRawStreamTracker.C

index 12da17e..36e2ae7 100644 (file)
@@ -1,17 +1,17 @@
 /**************************************************************************
- * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- *                                                                        *
- * Author: The ALICE Off-line Project.                                    *
- * Contributors are mentioned in the code where appropriate.              *
- *                                                                        *
- * Permission to use, copy, modify and distribute this software and its   *
- * documentation strictly for non-commercial purposes is hereby granted   *
- * without fee, provided that the above copyright notice appears in all   *
- * copies and that both the copyright notice and this permission notice   *
- * appear in the supporting documentation. The authors make no claims     *
- * about the suitability of this software for any purpose. It is          *
- * provided "as is" without express or implied warranty.                  *
- **************************************************************************/
+* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+*                                                                        *
+* Author: The ALICE Off-line Project.                                    *
+* Contributors are mentioned in the code where appropriate.              *
+*                                                                        *
+* Permission to use, copy, modify and distribute this software and its   *
+* documentation strictly for non-commercial purposes is hereby granted   *
+* without fee, provided that the above copyright notice appears in all   *
+* copies and that both the copyright notice and this permission notice   *
+* appear in the supporting documentation. The authors make no claims     *
+* about the suitability of this software for any purpose. It is          *
+* provided "as is" without express or implied warranty.                  *
+**************************************************************************/
 
 //
 // Macro for reading tracker raw data
 #include "AliMUONBusStruct.h"
 #include "AliMUONDDLTracker.h"
 
-#endif
+#include "TStopwatch.h"
 
+#endif
 
 // Macro to read rawdata for tracker
 // Ch. Finck, Subatech, April. 2006
 //
-// This macro is interface with AliRawReader for RAW
-// The different stucture of the patload are readout and stored in TClonesArray
+// This macro is interfaced with AliRawReader for RAW
+//
+// There are 2 ways of reading the data : one where each intermediate structure
+// (block, dsp, buspatch) is looped over, and one, using an iterator, where
+// we're directly accessing the pad informations (charge).
+//
+// The different stucture of the payload are readout and stored in TClonesArray
 // with AliMUONRawStreamTracker classe.
 // The macro just simpy read again the TClonesArray contents.
 // The parameter of each structure could be seen in the container classes
 // AliMUONBlockHeader, AliMUONBlockHeader, AliMUONBusStruct.
 // The class AliMUONDDLTracker manages the structure containers.
 // The number of structures in the rawdata file could be set.
+//
+// 
 
-void MUONRawStreamTracker(Int_t maxEvent = 1000, Int_t minDDL = 0, Int_t maxDDL = 19, TString fileName = "./")
+void MUONRawStreamTrackerExpert(TString fileName = "./", Int_t maxEvent = 1000,  
+                                Int_t minDDL = 0, Int_t maxDDL = 19)
 {
-   AliRawReader* rawReader = 0x0;
-
-   // check extention to choose the rawdata file format
-   if (fileName.EndsWith("/")) {
-     rawReader = new AliRawReaderFile(fileName); // DDL files
-   } else if (fileName.EndsWith(".root")) {
-     rawReader = new AliRawReaderRoot(fileName);
-   } else if (!fileName.IsNull()) {
-     rawReader = new AliRawReaderDate(fileName); // DATE file
-   }
-
-   // raw stream
-   AliMUONRawStreamTracker* rawStream  = new AliMUONRawStreamTracker(rawReader);    
-
-   // set the number of DDL block Dsp & buspatch structures that are PRESENT in the rawdata file
-   // it's NOT the number to be read.
-   // default wise set to 20, 2, 5 ans 5 respectively.
-   //    rawStream->SetMaxDDL(xx);
-   //    rawStream->SetMaxBlock(xx);
-   //    rawStream->SetMaxDsp(xx);
-   //    rawStream->SetMaxBus(xx);
-
-   // containers
-   AliMUONDDLTracker*       ddlTracker = 0x0;
-   AliMUONBlockHeader*      blkHeader  = 0x0;
-   AliMUONDspHeader*        dspHeader  = 0x0;
-   AliMUONBusStruct*        busStruct  = 0x0;
-
-   //   Loop over events  
-   Int_t iEvent = 0;
-   Int_t dataSize;
-
-   while (rawReader->NextEvent()) {
-
-     if (iEvent == maxEvent)
-       break;
-
-     printf("Event %d\n",iEvent++);
-
-     // read DDL while < 20 DDL
-     while(rawStream->NextDDL()) {
-
-       if (rawStream->GetDDL() < minDDL || rawStream->GetDDL() > maxDDL)
-        continue;
-
-       printf("\niDDL %d\n", rawStream->GetDDL());
-
-       ddlTracker =  rawStream->GetDDLTracker();
-
-       // loop over block structure
-       Int_t nBlock = ddlTracker->GetBlkHeaderEntries();
-       for(Int_t iBlock = 0; iBlock < nBlock ;iBlock++){
-
-        blkHeader = ddlTracker->GetBlkHeaderEntry(iBlock);
-        printf("Block Total length %d\n",blkHeader->GetTotalLength());
-
-        // loop over DSP structure
-        Int_t nDsp = blkHeader->GetDspHeaderEntries();
-        for(Int_t iDsp = 0; iDsp < nDsp ;iDsp++){   //DSP loop
-
-          dspHeader =  blkHeader->GetDspHeaderEntry(iDsp);
-          //   printf("Dsp length %d even word %d\n",dspHeader->GetTotalLength(), dspHeader->GetEventWord());
-
-          // loop over BusPatch structure
-          Int_t nBusPatch = dspHeader->GetBusPatchEntries();
-          for(Int_t iBusPatch = 0; iBusPatch < nBusPatch; iBusPatch++) {  
-
-            busStruct = dspHeader->GetBusPatchEntry(iBusPatch);
-
-            printf("busPatchId %d", busStruct->GetBusPatchId());
-            printf(" BlockId %d", busStruct->GetBlockId());
-            printf(" DspId %d\n", busStruct->GetDspId());
-
-          // loop over data
-            dataSize = busStruct->GetLength();
-            for (Int_t iData = 0; iData < dataSize; iData++) {
-
-              Int_t  manuId    = busStruct->GetManuId(iData);
-              Int_t  channelId = busStruct->GetChannelId(iData);
-              Int_t  charge    = busStruct->GetCharge(iData);
-              printf("manuId: %d, channelId: %d charge: %d\n", manuId, 
-                     channelId, charge);
-
-            } // iData
-          } // iBusPatch
-        } // iDsp
-       } // iBlock
-     } // NextDDL
-   }// NextEvent
+  /// Reads the data from fileName, using an "expert" mode where all substructures
+  /// are looped upon.
+  
+  TStopwatch timer;
+  timer.Start(kTRUE);
+  
+  AliRawReader* rawReader = 0x0;
+  
+  // check extention to choose the rawdata file format
+  if (fileName.EndsWith("/")) {
+    rawReader = new AliRawReaderFile(fileName); // DDL files
+  } else if (fileName.EndsWith(".root")) {
+    rawReader = new AliRawReaderRoot(fileName);
+  } else if (!fileName.IsNull()) {
+    rawReader = new AliRawReaderDate(fileName); // DATE file
+  }
+  
+  // raw stream
+  AliMUONRawStreamTracker* rawStream  = new AliMUONRawStreamTracker(rawReader);    
+  
+  // set the number of DDL block Dsp & buspatch structures that are PRESENT in the rawdata file
+  // it's NOT the number to be read.
+  // default wise set to 20, 2, 5 ans 5 respectively.
+  //    rawStream->SetMaxDDL(xx);
+  //    rawStream->SetMaxBlock(xx);
+  //    rawStream->SetMaxDsp(xx);
+  //    rawStream->SetMaxBus(xx);
+  
+  // containers
+  AliMUONDDLTracker*       ddlTracker = 0x0;
+  AliMUONBlockHeader*      blkHeader  = 0x0;
+  AliMUONDspHeader*        dspHeader  = 0x0;
+  AliMUONBusStruct*        busStruct  = 0x0;
+  
+  //   Loop over events  
+  Int_t iEvent = 0;
+  Int_t dataSize;
+  
+  while (rawReader->NextEvent()) {
+    
+    if (iEvent == maxEvent)
+      break;
+    
+    printf("Event %d\n",iEvent++);
+    
+    // read DDL while < 20 DDL
+    while(rawStream->NextDDL()) {
+      
+      if (rawStream->GetDDL() < minDDL || rawStream->GetDDL() > maxDDL)
+        continue;
+      
+      printf("\niDDL %d\n", rawStream->GetDDL());
+      
+      ddlTracker =  rawStream->GetDDLTracker();
+      
+      // loop over block structure
+      Int_t nBlock = ddlTracker->GetBlkHeaderEntries();
+      for(Int_t iBlock = 0; iBlock < nBlock ;iBlock++){
+        
+        blkHeader = ddlTracker->GetBlkHeaderEntry(iBlock);
+        printf("Block %d Total length %d\n",iBlock,blkHeader->GetTotalLength());
+       
+        // loop over DSP structure
+        Int_t nDsp = blkHeader->GetDspHeaderEntries();
+        for(Int_t iDsp = 0; iDsp < nDsp ;iDsp++){   //DSP loop
+          
+          dspHeader =  blkHeader->GetDspHeaderEntry(iDsp);
+          printf("Dsp %d length %d error word %d\n",iDsp,dspHeader->GetTotalLength(), dspHeader->GetErrorWord());
+          
+          // loop over BusPatch structure
+          Int_t nBusPatch = dspHeader->GetBusPatchEntries();
+          for(Int_t iBusPatch = 0; iBusPatch < nBusPatch; iBusPatch++) {  
+            
+            busStruct = dspHeader->GetBusPatchEntry(iBusPatch);
+            
+            //      printf("busPatchId %d", busStruct->GetBusPatchId());
+            //      printf(" BlockId %d", busStruct->GetBlockId());
+            //      printf(" DspId %d\n", busStruct->GetDspId());
+            
+            // loop over data
+            dataSize = busStruct->GetLength();
+            for (Int_t iData = 0; iData < dataSize; iData++) {
+              
+              Int_t  manuId    = busStruct->GetManuId(iData);
+              Int_t  channelId = busStruct->GetChannelId(iData);
+              Int_t  charge    = busStruct->GetCharge(iData);
+              printf("buspatch %5d manuI %4d channel %3d charge %4d\n", 
+                     busStruct->GetBusPatchId(),
+                     manuId, 
+                     channelId, charge);
+            } // iData
+          } // iBusPatch
+        } // iDsp
+      } // iBlock
+    } // NextDDL
+  }// NextEvent
+  
+  delete rawReader;
+  delete rawStream;
+  timer.Print();
+}
 
-   delete rawReader;
-   delete rawStream;
+void MUONRawStreamTrackerSimple(TString fileName, Int_t maxEvent)
+{
+  /// Reads the raw data in fileName, using a simplified interface (iterator
+  /// over pads).
+  TStopwatch timer;
+  timer.Start(kTRUE);
+  
+  AliRawReader* rawReader = 0x0;
+  
+  // check extention to choose the rawdata file format
+  if (fileName.EndsWith("/")) {
+    rawReader = new AliRawReaderFile(fileName); // DDL files
+  } else if (fileName.EndsWith(".root")) {
+    rawReader = new AliRawReaderRoot(fileName);
+  } else if (!fileName.IsNull()) {
+    rawReader = new AliRawReaderDate(fileName); // DATE file
+  }
+  
+  // raw stream
+  AliMUONRawStreamTracker* rawStream  = new AliMUONRawStreamTracker(rawReader);    
+  
+  //   Loop over events  
+  Int_t iEvent = 0;
+  
+  while (rawReader->NextEvent()) {
+    
+    if (iEvent == maxEvent)
+      break;
+    
+    printf("Event %d\n",iEvent++);
+    
+    Int_t busPatch;
+    UShort_t manuId, adc;
+    UChar_t manuChannel;
+    
+    rawStream->First();
+    
+    while ( rawStream->Next(busPatch,manuId,manuChannel,adc) )
+    {      
+      printf("buspatch %5d manuI %4d channel %3d charge %4d\n", 
+             busPatch,manuId,manuChannel, adc);
+    }
+  }
+  
+  delete rawReader;
+  delete rawStream;
+  timer.Print();
 }
+