MUONRawStreamTrigger.C

Go to the documentation of this file.

/**************************************************************************
 * 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.                  *
 **************************************************************************/

// $Id$



#if !defined(__CINT__) || defined(__MAKECINT__)

// RAW includes
#include "AliRawReader.h"
#include "AliRawReaderFile.h"
#include "AliRawReaderRoot.h"

// MUON includes
#include "AliMUONRawStreamTrigger.h"
#include "AliMUONRawStreamTriggerHP.h"
#include "AliMUONDarcHeader.h"
#include "AliMUONRegHeader.h"
#include "AliMUONLocalStruct.h"
#include "AliMUONDDLTrigger.h"
#include "AliMpTriggerCrate.h"
#include "AliMpDDLStore.h"
#include "AliMpCDB.h"

#include "TStopwatch.h"


#endif

void MUONRawStreamTrigger(Int_t maxEvent = 1, Int_t minDDL = 0, Int_t maxDDL = 1, TString fileName = "./")
{
   
   TStopwatch timer;
   timer.Start(kTRUE);

   AliRawReader* rawReader = AliRawReader::Create(fileName.Data());

   // Load mapping
     if ( ! AliMpCDB::LoadDDLStore() ) {
       printf("Could not access mapping from OCDB !\n");
     }

   // raw stream
   AliMUONRawStreamTrigger* rawStream   = new AliMUONRawStreamTrigger(rawReader);

   // set the number ofreg & local that are PRESENT in the rawdata file
   // it's NOT the number to be read.
   // default wise set to 8, 16 respectively.
   //    rawStream->SetMaxReg(2);
   //    rawStream->SetMaxLoc(xx);

   // containers
   AliMUONDDLTrigger*       ddlTrigger  = 0x0;
   AliMUONDarcHeader*       darcHeader  = 0x0;
   AliMUONRegHeader*        regHeader   = 0x0;
   AliMUONLocalStruct*      localStruct = 0x0;


   // Loop over events  
   Int_t iEvent = 0;

   while (rawReader->NextEvent()) {

     if (iEvent == maxEvent)
       break;

     printf("Event %d\n",iEvent++);

     // read DDL while < 2 DDL
     while(rawStream->NextDDL()) {

      if (rawStream->GetDDL() < minDDL || rawStream->GetDDL() > maxDDL)
   continue;

       printf("\niDDL %d\n", rawStream->GetDDL());

       ddlTrigger = rawStream->GetDDLTrigger();
       darcHeader = ddlTrigger->GetDarcHeader();

       printf("Global output %x\n", (Int_t)darcHeader->GetGlobalOutput());

       // loop over regional structures
       Int_t nReg = darcHeader->GetRegHeaderEntries();
       for(Int_t iReg = 0; iReg < nReg ;iReg++){   //REG loop

//   printf("RegionalId %d\n", iReg);

   regHeader =  darcHeader->GetRegHeaderEntry(iReg);
   //  printf("Reg length %d\n",regHeader->GetHeaderLength());

   // crate info  
   AliMpTriggerCrate* crate = AliMpDDLStore::Instance()->
                              GetTriggerCrate(rawStream->GetDDL(), iReg);

   // loop over local structures
   Int_t nLocal = regHeader->GetLocalEntries();
   for(Int_t iLocal = 0; iLocal < nLocal; iLocal++) {  

     localStruct = regHeader->GetLocalEntry(iLocal);

     Int_t iLocCard = crate->GetLocalBoardId(localStruct->GetId());

     if ( !iLocCard ) continue; // empty slot

     // check if trigger 
     if (localStruct->GetTriggerX() 
         || localStruct->GetTriggerY()) { // no empty data

       printf("LocalId %d\n", localStruct->GetId());

       Int_t loStripX  = (Int_t)localStruct->GetXPos();
       Int_t loStripY  = (Int_t)localStruct->GetYPos();
       Int_t loDev     = (Int_t)localStruct->GetXDev();
         
       printf("iLocCard: %d, XPos: %d, YPos: %d Dev: %d\n", iLocCard, loStripX, loStripY, loDev);

     }
   } // iLocal
       } // iReg
     } // NextDDL
   }// NextEvent

   delete rawReader;
   delete rawStream;

   timer.Print();

}


void MUONRawStreamTriggerHPExpert(Int_t maxEvent = 1, Int_t minDDL = 0, Int_t maxDDL = 1, TString fileName = "./")
{
   
   TStopwatch timer;
   timer.Start(kTRUE);

   AliRawReader* rawReader = AliRawReader::Create(fileName.Data());

   // Load mapping
   if ( ! AliMpCDB::LoadDDLStore() ) {
     printf("Could not access mapping from OCDB !\n");
   }

   // raw stream
   AliMUONRawStreamTriggerHP* rawStream = new AliMUONRawStreamTriggerHP(rawReader);

   // Light weight interfaces to the raw data.
   const AliMUONRawStreamTriggerHP::AliHeader*          darcHeader  = 0x0;
   const AliMUONRawStreamTriggerHP::AliRegionalHeader*  regHeader   = 0x0;
   const AliMUONRawStreamTriggerHP::AliLocalStruct*     localStruct = 0x0;

   // Loop over events.
   Int_t iEvent = 0;

   while (rawReader->NextEvent()) {

     if (iEvent == maxEvent)
       break;

     printf("Event %d\n",iEvent++);

     // read DDL while < 2 DDL
     while(rawStream->NextDDL()) {

      if (rawStream->GetDDL() < minDDL || rawStream->GetDDL() > maxDDL)
   continue;

       printf("\niDDL %d\n", rawStream->GetDDL());

       darcHeader = rawStream->GetHeaders();

       printf("Global output %x\n", (Int_t)darcHeader->GetGlobalOutput());

       // loop over regional structures
       Int_t nReg = (Int_t) rawStream->GetRegionalHeaderCount();
       for(Int_t iReg = 0; iReg < nReg ;iReg++){   //REG loop

//   printf("RegionalId %d\n", iReg);

   regHeader =  rawStream->GetRegionalHeader(iReg);

   // crate info
   AliMpTriggerCrate* crate = AliMpDDLStore::Instance()->
                              GetTriggerCrate(rawStream->GetDDL(), iReg);

   // loop over local structures
   Int_t nLocal = rawStream->GetLocalStructCount(iReg);
   for(Int_t iLocal = 0; iLocal < nLocal; iLocal++) {

     localStruct = rawStream->GetLocalStruct(iReg, iLocal);

     Int_t iLocCard = crate->GetLocalBoardId(localStruct->GetId());

     if ( !iLocCard ) continue; // empty slot
    
     // check if trigger 
     if (localStruct->GetTriggerX() 
         || localStruct->GetTriggerY()) { // no empty data

       printf("LocalId %d\n", localStruct->GetId());

       Int_t loStripX  = (Int_t)localStruct->GetXPos();
       Int_t loStripY  = (Int_t)localStruct->GetYPos();
       Int_t loDev     = (Int_t)localStruct->GetXDev();

       printf("iLocCard: %d, XPos: %d, YPos: %d Dev: %d\n", iLocCard, loStripX, loStripY, loDev);

     }
   } // iLocal
       } // iReg
     } // NextDDL
   }// NextEvent

   delete rawReader;
   delete rawStream;

   timer.Print();

}


void MUONRawStreamTriggerHPExpert2(Int_t maxEvent = 1, Int_t minDDL = 0, Int_t maxDDL = 1, TString fileName = "./")
{
   
   TStopwatch timer;
   timer.Start(kTRUE);

   AliRawReader* rawReader = AliRawReader::Create(fileName.Data());

   // Load mapping
   if ( ! AliMpCDB::LoadDDLStore() ) {
     printf("Could not access mapping from OCDB !\n");
   }

   // raw stream
   AliMUONRawStreamTriggerHP* rawStream = new AliMUONRawStreamTriggerHP(rawReader);

   // Light weight interfaces to the raw data.
   const AliMUONRawStreamTriggerHP::AliHeader*          darcHeader  = 0x0;
   const AliMUONRawStreamTriggerHP::AliRegionalHeader*  regHeader   = 0x0;
   const AliMUONRawStreamTriggerHP::AliLocalStruct*     localStruct = 0x0;

   // Loop over events.
   Int_t iEvent = 0;

   while (rawReader->NextEvent()) {

     if (iEvent == maxEvent)
       break;

     printf("Event %d\n",iEvent++);

     // read DDL while < 2 DDL
     while(rawStream->NextDDL()) {

      if (rawStream->GetDDL() < minDDL || rawStream->GetDDL() > maxDDL)
   continue;

       printf("\niDDL %d\n", rawStream->GetDDL());

       darcHeader = rawStream->GetHeaders();

       printf("Global output %x\n", (Int_t)darcHeader->GetGlobalOutput());

       // loop over regional structures
       Int_t nReg = (Int_t) rawStream->GetRegionalHeaderCount();
       for(Int_t iReg = 0; iReg < nReg ;iReg++){   //REG loop

//   printf("RegionalId %d\n", iReg);

   regHeader =  rawStream->GetRegionalHeader(iReg);

   // crate info
   AliMpTriggerCrate* crate = AliMpDDLStore::Instance()->
                              GetTriggerCrate(rawStream->GetDDL(), iReg);

   // loop over local structures
   Int_t nLocal = regHeader->GetLocalStructCount();
   for(Int_t iLocal = 0; iLocal < nLocal; iLocal++) {

     localStruct = regHeader->GetLocalStruct(iLocal);

     Int_t iLocCard = crate->GetLocalBoardId(localStruct->GetId());

     if ( !iLocCard ) continue; // empty slot
    
     // check if trigger 
     if (localStruct->GetTriggerX() 
         || localStruct->GetTriggerY()) { // no empty data

       printf("LocalId %d\n", localStruct->GetId());

       Int_t loStripX  = (Int_t)localStruct->GetXPos();
       Int_t loStripY  = (Int_t)localStruct->GetYPos();
       Int_t loDev     = (Int_t)localStruct->GetXDev();

       printf("iLocCard: %d, XPos: %d, YPos: %d Dev: %d\n", iLocCard, loStripX, loStripY, loDev);

     }
   } // iLocal
       } // iReg
     } // NextDDL
   }// NextEvent

   delete rawReader;
   delete rawStream;

   timer.Print();

}


void MUONRawStreamTriggerHPExpert3(Int_t maxEvent = 1, Int_t minDDL = 0, Int_t maxDDL = 1, TString fileName = "./")
{
   
   TStopwatch timer;
   timer.Start(kTRUE);

   AliRawReader* rawReader = AliRawReader::Create(fileName.Data());

   // Load mapping
   if ( ! AliMpCDB::LoadDDLStore() ) {
     printf("Could not access mapping from OCDB !\n");
   }

   // raw stream
   AliMUONRawStreamTriggerHP* rawStream = new AliMUONRawStreamTriggerHP(rawReader);

   // Light weight interfaces to the raw data.
   const AliMUONRawStreamTriggerHP::AliHeader*          darcHeader  = 0x0;
   const AliMUONRawStreamTriggerHP::AliRegionalHeader*  regHeader   = 0x0;
   const AliMUONRawStreamTriggerHP::AliLocalStruct*     localStruct = 0x0;

   // Loop over events.
   Int_t iEvent = 0;

   while (rawReader->NextEvent()) {

     if (iEvent == maxEvent)
       break;

     printf("Event %d\n",iEvent++);

     // read DDL while < 2 DDL
     while(rawStream->NextDDL()) {

      if (rawStream->GetDDL() < minDDL || rawStream->GetDDL() > maxDDL)
   continue;

       printf("\niDDL %d\n", rawStream->GetDDL());

       darcHeader = rawStream->GetHeaders();

       printf("Global output %x\n", (Int_t)darcHeader->GetGlobalOutput());

       // loop over regional structures
       Int_t iReg = 0;
       regHeader = rawStream->GetFirstRegionalHeader();
       while (regHeader != NULL)
       {
//   printf("RegionalId %d\n", iReg);

   regHeader =  rawStream->GetRegionalHeader(iReg);

   // crate info
   AliMpTriggerCrate* crate = AliMpDDLStore::Instance()->
                              GetTriggerCrate(rawStream->GetDDL(), iReg);

   // loop over local structures
   Int_t iLocal = 0;
   localStruct = regHeader->GetFirstLocalStruct();
   while (localStruct != NULL)
   {
     Int_t iLocCard = crate->GetLocalBoardId(localStruct->GetId());
     if ( !iLocCard ) continue; // empty slot
    
     // check if trigger 
     if (localStruct->GetTriggerX() 
         || localStruct->GetTriggerY()) { // no empty data

       printf("LocalId %d\n", localStruct->GetId());

       Int_t loStripX  = (Int_t)localStruct->GetXPos();
       Int_t loStripY  = (Int_t)localStruct->GetYPos();
       Int_t loDev     = (Int_t)localStruct->GetXDev();

       printf("iLocCard: %d, XPos: %d, YPos: %d Dev: %d\n", iLocCard, loStripX, loStripY, loDev);

     }
     iLocal++;
     localStruct = localStruct->Next();
   } // iLocal
   iReg++;
   regHeader = regHeader->Next();
       } // iReg
     } // NextDDL
   }// NextEvent

   delete rawReader;
   delete rawStream;

   timer.Print();

}


void MUONRawStreamTriggerSimple(Int_t maxEvent = 1, TString fileName = "./")
{

  TStopwatch timer;
  timer.Start(kTRUE);

   AliRawReader* rawReader = 0x0;

   if (fileName.EndsWith("/")) {
     rawReader = new AliRawReaderFile(fileName);// DDL files
   } else if (fileName.EndsWith(".root")) {
     rawReader = new AliRawReaderRoot(fileName);
   } 

   // raw stream
   AliMUONRawStreamTrigger* rawStream   = new AliMUONRawStreamTrigger(rawReader);

   // set the number of reg & local that are PRESENT in the rawdata file
   // it's NOT the number to be read.
   // default wise set to 8, 16 respectively.
   //    rawStream->SetMaxReg(2);
   //    rawStream->SetMaxLoc(xx);

   UChar_t id;   
   UChar_t dec;
   Bool_t trigY; 
   UChar_t yPos; 
   UChar_t sXDev; 
   UChar_t xDev;
   UChar_t xPos;

   Bool_t triggerX; 
   Bool_t triggerY; 

 
   TArrayS xPattern; 
   TArrayS yPattern;

   // Loop over events  
   Int_t iEvent = 0;

   while (rawReader->NextEvent()) {

     if (iEvent == maxEvent)
   break;

     printf("Event %d\n",iEvent++);

     rawStream->First();

     // read while there are digits
     while( rawStream->Next(id, dec, trigY, yPos, sXDev, xDev, xPos,
          triggerX, triggerY, xPattern, yPattern) ) 
     {
       if ( triggerX || triggerY )  // no empty data
     printf("iLocCard: %d, XPos: %d, YPos: %d Dev: %d\n", id, xPos, yPos, xDev);

     }// Next
  
   }// NextEvent

   delete rawReader;
   delete rawStream;

   timer.Print();

}


void MUONRawStreamTriggerHPSimple(Int_t maxEvent = 1, TString fileName = "./")
{

  TStopwatch timer;
  timer.Start(kTRUE);

  AliRawReader* rawReader = AliRawReader::Create(fileName.Data());
  AliMUONRawStreamTriggerHP* rawStream   = new AliMUONRawStreamTriggerHP(rawReader);

  UChar_t id;
  UChar_t dec;
  Bool_t trigY;
  UChar_t yPos;
  UChar_t sXDev;
  UChar_t xDev;
  UChar_t xPos;
  Bool_t triggerX;
  Bool_t triggerY;
  TArrayS xPattern;
  TArrayS yPattern;

  // Loop over events
  Int_t iEvent = 0;

  while (rawReader->NextEvent())
  {

    if (iEvent == maxEvent)
       break;

     printf("Event %d\n", iEvent++);

     rawStream->First();

     // read while there are digits
     while ( rawStream->Next(id, dec, trigY, yPos, sXDev, xDev, xPos,
                            triggerX, triggerY, xPattern, yPattern) )
     {
       if ( triggerX || triggerY )  // no empty data
          printf("iLocCard: %d, XPos: %d, YPos: %d Dev: %d\n", id, xPos, yPos, xDev);
     }// Next
  }// NextEvent

  delete rawReader;
  delete rawStream;

  timer.Print();
}


void MUONRawStreamTriggerHPSimple2(Int_t maxEvent = 1, TString fileName = "./")
{

  TStopwatch timer;
  timer.Start(kTRUE);

  AliRawReader* rawReader = AliRawReader::Create(fileName.Data());
  AliMUONRawStreamTriggerHP* rawStream   = new AliMUONRawStreamTriggerHP(rawReader);

  const AliMUONRawStreamTriggerHP::AliLocalStruct* localStruct;

  // Loop over events
  Int_t iEvent = 0;

  while (rawReader->NextEvent())
  {

    if (iEvent == maxEvent)
       break;

     printf("Event %d\n", iEvent++);

     rawStream->First();

     // read while there are digits
     while ( (localStruct = rawStream->Next()) != NULL )
     {
       if ( localStruct->GetTriggerX() || localStruct->GetTriggerY() )  // no empty data
          printf("iLocCard: %d, XPos: %d, YPos: %d Dev: %d\n",
            localStruct->GetId(),
            localStruct->GetXPos(),
            localStruct->GetYPos(),
            localStruct->GetXDev()
          );
     }// Next
  }// NextEvent

  delete rawReader;
  delete rawStream;

  timer.Print();
}