[u/mrichter/AliRoot.git] / RAW / AliRawReader.cxx

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

///////////////////////////////////////////////////////////////////////////////
///
/// This is the base class for reading raw data.
///
/// The derived classes, which operate on concrete raw data formats,
/// should implement
/// - ReadHeader to read the next (data/equipment) header
/// - ReadNextData to read the next raw data block (=1 DDL)
/// - ReadNext to read a given number of bytes
/// - several getters like GetType
///
/// Sequential access to the raw data is provided by the methods
/// ReadHeader, ReadNextData, ReadNextInt, ReadNextShort, ReadNextChar
///
/// If only data from a specific detector (and a given range of DDL numbers)
/// should be read, this can be achieved by the Select method.
/// Several getters provide information about the current event and the
/// current type of raw data.
///
///////////////////////////////////////////////////////////////////////////////

#include <Riostream.h>
#include "AliRawReader.h"
#include "AliDAQ.h"

ClassImp(AliRawReader)


AliRawReader::AliRawReader() :
  fEquipmentIdsIn(NULL),
  fEquipmentIdsOut(NULL),
  fRequireHeader(kTRUE),
  fHeader(NULL),
  fCount(0),
  fSelectEquipmentType(-1),
  fSelectMinEquipmentId(-1),
  fSelectMaxEquipmentId(-1),
  fSkipInvalid(kFALSE),
  fSelectEventType(-1),
  fErrorCode(0)
{
// default constructor: initialize data members
}

Bool_t AliRawReader::LoadEquipmentIdsMap(const char *fileName)
{
  // Open the mapping file
  // and load the mapping data
  ifstream input(fileName);
  if (input.is_open()) {
    Warning("AliRawReader","Equipment ID mapping file is found !");
    const Int_t kMaxDDL = 256;
    fEquipmentIdsIn = new TArrayI(kMaxDDL);
    fEquipmentIdsOut = new TArrayI(kMaxDDL);
    Int_t equipIn, equipOut;
    Int_t nIds = 0;
    while (input >> equipIn >> equipOut) {
      if (nIds >= kMaxDDL) {
	Error("AliRawReader","Too many equipment Id mappings found ! Truncating the list !");
	break;
      }
      fEquipmentIdsIn->AddAt(equipIn,nIds); 
      fEquipmentIdsOut->AddAt(equipOut,nIds);
      nIds++;
    }
    fEquipmentIdsIn->Set(nIds);
    fEquipmentIdsOut->Set(nIds);
    input.close();
    return kTRUE;
  }
  else {
    Error("AliRawReader","equipment id map file is not found ! Skipping the mapping !");
    return kFALSE;
  }
}

AliRawReader::AliRawReader(const AliRawReader& rawReader) :
  TObject(rawReader),
  fEquipmentIdsIn(rawReader.fEquipmentIdsIn),
  fEquipmentIdsOut(rawReader.fEquipmentIdsOut),
  fRequireHeader(rawReader.fRequireHeader),
  fHeader(rawReader.fHeader),
  fCount(rawReader.fCount),
  fSelectEquipmentType(rawReader.fSelectEquipmentType),
  fSelectMinEquipmentId(rawReader.fSelectMinEquipmentId),
  fSelectMaxEquipmentId(rawReader.fSelectMaxEquipmentId),
  fSkipInvalid(rawReader.fSkipInvalid),
  fSelectEventType(rawReader.fSelectEventType),
  fErrorCode(0)
{
// copy constructor
}

AliRawReader& AliRawReader::operator = (const AliRawReader& rawReader)
{
// assignment operator
  fEquipmentIdsIn = rawReader.fEquipmentIdsIn;
  fEquipmentIdsOut = rawReader.fEquipmentIdsOut;

  fHeader = rawReader.fHeader;
  fCount = rawReader.fCount;

  fSelectEquipmentType = rawReader.fSelectEquipmentType;
  fSelectMinEquipmentId = rawReader.fSelectMinEquipmentId;
  fSelectMaxEquipmentId = rawReader.fSelectMaxEquipmentId;
  fSkipInvalid = rawReader.fSkipInvalid;
  fSelectEventType = rawReader.fSelectEventType;

  fErrorCode = rawReader.fErrorCode;

  return *this;
}

AliRawReader::~AliRawReader()
{
  // destructor
  // delete the mapping arrays if
  // initialized
  if (fEquipmentIdsIn) delete fEquipmentIdsIn;
  if (fEquipmentIdsOut) delete fEquipmentIdsOut;
}

Int_t AliRawReader::GetMappedEquipmentId() const
{
  if (!fEquipmentIdsIn || !fEquipmentIdsOut) {
    Error("AliRawReader","equipment Ids mapping is not initialized !");
    return GetEquipmentId();
  }
  Int_t equipmentId = GetEquipmentId();
  for(Int_t iId = 0; iId < fEquipmentIdsIn->GetSize(); iId++) {
    if (equipmentId == fEquipmentIdsIn->At(iId)) {
      equipmentId = fEquipmentIdsOut->At(iId);
      break;
    }
  }
  return equipmentId;
}

Int_t AliRawReader::GetDetectorID() const
{
  // Get the detector ID
  // The list of detector IDs
  // can be found in AliDAQ.h
  Int_t equipmentId;
  if (fEquipmentIdsIn && fEquipmentIdsIn)
    equipmentId = GetMappedEquipmentId();
  else
    equipmentId = GetEquipmentId();

  if (equipmentId >= 0) {
    Int_t ddlIndex;
    return AliDAQ::DetectorIDFromDdlID(equipmentId,ddlIndex);
  }
  else
    return -1;
}

Int_t AliRawReader::GetDDLID() const
{
  // Get the DDL ID (within one sub-detector)
  // The list of detector IDs
  // can be found in AliDAQ.h
  Int_t equipmentId;
  if (fEquipmentIdsIn && fEquipmentIdsIn)
    equipmentId = GetMappedEquipmentId();
  else
    equipmentId = GetEquipmentId();

  if (equipmentId >= 0) {
    Int_t ddlIndex;
    AliDAQ::DetectorIDFromDdlID(equipmentId,ddlIndex);
    return ddlIndex;
  }
  else
    return -1;
}

void AliRawReader::Select(const char *detectorName, Int_t minDDLID, Int_t maxDDLID)
{
// read only data of the detector with the given name and in the given
// range of DDLs (minDDLID <= DDLID <= maxDDLID).
// no selection is applied if a value < 0 is used.
  Int_t detectorID = AliDAQ::DetectorID(detectorName);
  if(detectorID >= 0)
    Select(detectorID,minDDLID,maxDDLID);
}

void AliRawReader::Select(Int_t detectorID, Int_t minDDLID, Int_t maxDDLID)
{
// read only data of the detector with the given ID and in the given
// range of DDLs (minDDLID <= DDLID <= maxDDLID).
// no selection is applied if a value < 0 is used.

  fSelectEquipmentType = -1;

  if (minDDLID < 0)
    fSelectMinEquipmentId = AliDAQ::DdlIDOffset(detectorID);
  else
    fSelectMinEquipmentId = AliDAQ::DdlID(detectorID,minDDLID);

  if (maxDDLID < 0)
    fSelectMaxEquipmentId = AliDAQ::DdlID(detectorID,AliDAQ::NumberOfDdls(detectorID)-1);
  else
    fSelectMaxEquipmentId = AliDAQ::DdlID(detectorID,maxDDLID);
}

void AliRawReader::SelectEquipment(Int_t equipmentType, 
				   Int_t minEquipmentId, Int_t maxEquipmentId)
{
// read only data of the equipment with the given type and in the given
// range of IDs (minEquipmentId <= EquipmentId <= maxEquipmentId).
// no selection is applied if a value < 0 is used.

  fSelectEquipmentType = equipmentType;
  fSelectMinEquipmentId = minEquipmentId;
  fSelectMaxEquipmentId = maxEquipmentId;
}

void AliRawReader::SelectEvents(Int_t type)
{
// read only events with the given type.
// no selection is applied if a value < 0 is used.

  fSelectEventType = type;
}

Bool_t AliRawReader::IsSelected() const
{
// apply the selection (if any)

  if (fSkipInvalid && !IsValid()) return kFALSE;

  if (fSelectEquipmentType >= 0)
    if (GetEquipmentType() != fSelectEquipmentType) return kFALSE;

  Int_t equipmentId;
  if (fEquipmentIdsIn && fEquipmentIdsIn)
    equipmentId = GetMappedEquipmentId();
  else
    equipmentId = GetEquipmentId();

  if ((fSelectMinEquipmentId >= 0) && 
      (equipmentId < fSelectMinEquipmentId))
    return kFALSE;
  if ((fSelectMaxEquipmentId >= 0) && 
      (equipmentId > fSelectMaxEquipmentId))
    return kFALSE;

  return kTRUE;
}

Bool_t AliRawReader::IsEventSelected() const
{
// apply the event selection (if any)

  if (fSelectEventType >= 0) {
    if (GetType() != (UInt_t) fSelectEventType) return kFALSE;
  }

  return kTRUE;
}


Bool_t AliRawReader::ReadNextInt(UInt_t& data)
{
// reads the next 4 bytes at the current position
// returns kFALSE if the data could not be read

  while (fCount == 0) {
    if (!ReadHeader()) return kFALSE;
  }
  if (fCount < (Int_t) sizeof(data)) {
    Error("ReadNextInt", 
	  "too few data left (%d bytes) to read an UInt_t!", fCount);
    return kFALSE;
  }
  if (!ReadNext((UChar_t*) &data, sizeof(data))) {
    Error("ReadNextInt", "could not read data!");
    return kFALSE;
  }
  return kTRUE;
}

Bool_t AliRawReader::ReadNextShort(UShort_t& data)
{
// reads the next 2 bytes at the current position
// returns kFALSE if the data could not be read

  while (fCount == 0) {
    if (!ReadHeader()) return kFALSE;
  }
  if (fCount < (Int_t) sizeof(data)) {
    Error("ReadNextShort", 
	  "too few data left (%d bytes) to read an UShort_t!", fCount);
    return kFALSE;
  }
  if (!ReadNext((UChar_t*) &data, sizeof(data))) {
    Error("ReadNextShort", "could not read data!");
    return kFALSE;
  }
  return kTRUE;
}

Bool_t AliRawReader::ReadNextChar(UChar_t& data)
{
// reads the next 1 byte at the current stream position
// returns kFALSE if the data could not be read

  while (fCount == 0) {
    if (!ReadHeader()) return kFALSE;
  }
  if (!ReadNext((UChar_t*) &data, sizeof(data))) {
    Error("ReadNextChar", "could not read data!");
    return kFALSE;
  }
  return kTRUE;
}


Int_t AliRawReader::CheckData() const
{
// check the consistency of the data
// derived classes should overwrite the default method which returns 0 (no err)

  return 0;
}


void AliRawReader::DumpData(Int_t limit)
{
// print the raw data
// if limit is not negative, only the first and last "limit" lines of raw data
// are printed

  Reset();
  if (!ReadHeader()) {
    Error("DumpData", "no header");
    return;
  }
  printf("header:\n"
	 " type = %d  run = %d  ", GetType(), GetRunNumber());
  if (GetEventId()) {
    printf("event = %8.8x %8.8x\n", GetEventId()[1], GetEventId()[0]);
  } else {
    printf("event = -------- --------\n");
  }
  if (GetTriggerPattern()) {
    printf(" trigger = %8.8x %8.8x  ",
	   GetTriggerPattern()[1], GetTriggerPattern()[0]);
  } else {
    printf(" trigger = -------- --------  ");
  }
  if (GetDetectorPattern()) {
    printf("detector = %8.8x\n", GetDetectorPattern()[0]);
  } else {
    printf("detector = --------\n");
  }
  if (GetAttributes()) {
    printf(" attributes = %8.8x %8.8x %8.8x  ",
	   GetAttributes()[2], GetAttributes()[1], GetAttributes()[0]);
  } else {
    printf(" attributes = -------- -------- --------  ");
  }
  printf("GDC = %d\n", GetGDCId());
  printf("\n");

  do {
    printf("-------------------------------------------------------------------------------\n");
    printf("LDC = %d\n", GetLDCId());

    printf("equipment:\n"
	   " size = %d  type = %d  id = %d\n",
	   GetEquipmentSize(), GetEquipmentType(), GetEquipmentId());
    if (GetEquipmentAttributes()) {
      printf(" attributes = %8.8x %8.8x %8.8x  ", GetEquipmentAttributes()[2],
	     GetEquipmentAttributes()[1], GetEquipmentAttributes()[0]);
    } else {
      printf(" attributes = -------- -------- --------  ");
    }
    printf("element size = %d\n", GetEquipmentElementSize());

    printf("data header:\n"
	   " size = %d  version = %d  valid = %d  compression = %d\n",
	   GetDataSize(), GetVersion(), IsValid(), IsCompressed());

    printf("\n");
    if (limit == 0) continue;

    Int_t size = GetDataSize();
    char line[70];
    for (Int_t i = 0; i < 70; i++) line[i] = ' ';
    line[69] = '\0';
    Int_t pos = 0;
    Int_t max = 16;
    UChar_t byte;

    for (Int_t n = 0; n < size; n++) {
      if (!ReadNextChar(byte)) {
	Error("DumpData", "couldn't read byte number %d\n", n);
	break;
      }
      if (pos >= max) {
	printf("%8.8x  %s\n", n-pos, line);
	for (Int_t i = 0; i < 70; i++) line[i] = ' ';
	line[69] = '\0';
	pos = 0;
	if ((limit > 0) && (n/max == limit)) {
	  Int_t nContinue = ((size-1)/max+1-limit) * max;
	  if (nContinue > n) {
	    printf(" [skipping %d bytes]\n", nContinue-n);
	    n = nContinue-1;
	    continue;
	  }
	}
      }
      Int_t offset = pos/4;
      if ((byte > 0x20) && (byte < 0x7f)) {
	line[pos+offset] = byte;
      } else {
	line[pos+offset] = '.';
      }
      char hex[3];
      sprintf(hex, "%2.2x", byte);
      line[max+max/4+3+2*pos+offset] = hex[0];
      line[max+max/4+4+2*pos+offset] = hex[1];
      pos++;
    }

    if (pos > 0) printf("%8.8x  %s\n", size-pos, line);
    printf("\n");
	   
  } while (ReadHeader());
}
Commit	Line	Data
a6e7b125	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
bea6b2a4	16	/* $Id$ */
bea6b2a4	17
a6e7b125	18	///////////////////////////////////////////////////////////////////////////////
bea6b2a4	19	///
	20	/// This is the base class for reading raw data.
	21	///
	22	/// The derived classes, which operate on concrete raw data formats,
	23	/// should implement
	24	/// - ReadHeader to read the next (data/equipment) header
	25	/// - ReadNextData to read the next raw data block (=1 DDL)
	26	/// - ReadNext to read a given number of bytes
	27	/// - several getters like GetType
	28	///
	29	/// Sequential access to the raw data is provided by the methods
	30	/// ReadHeader, ReadNextData, ReadNextInt, ReadNextShort, ReadNextChar
	31	///
	32	/// If only data from a specific detector (and a given range of DDL numbers)
	33	/// should be read, this can be achieved by the Select method.
	34	/// Several getters provide information about the current event and the
	35	/// current type of raw data.
	36	///
a6e7b125	37	///////////////////////////////////////////////////////////////////////////////
a6e7b125	38
0cfe2438	39	#include <Riostream.h>
a6e7b125	40	#include "AliRawReader.h"
362c9d61	41	#include "AliDAQ.h"
8de97894	42
a6e7b125	43	ClassImp(AliRawReader)
	44
	45
c946ab02	46	AliRawReader::AliRawReader() :
0cfe2438	47	fEquipmentIdsIn(NULL),
0cfe2438	48	fEquipmentIdsOut(NULL),
18882b3b	49	fRequireHeader(kTRUE),
39f9963f	50	fHeader(NULL),
c946ab02	51	fCount(0),
f224de6c	52	fSelectEquipmentType(-1),
	53	fSelectMinEquipmentId(-1),
	54	fSelectMaxEquipmentId(-1),
39f9963f	55	fSkipInvalid(kFALSE),
dd9a70fe	56	fSelectEventType(-1),
c946ab02	57	fErrorCode(0)
a6e7b125	58	{
42d20574	59	// default constructor: initialize data members
0cfe2438	60	}
42d20574	61
0cfe2438	62	Bool_t AliRawReader::LoadEquipmentIdsMap(const char *fileName)
	63	{
	64	// Open the mapping file
	65	// and load the mapping data
	66	ifstream input(fileName);
	67	if (input.is_open()) {
	68	Warning("AliRawReader","Equipment ID mapping file is found !");
	69	const Int_t kMaxDDL = 256;
	70	fEquipmentIdsIn = new TArrayI(kMaxDDL);
	71	fEquipmentIdsOut = new TArrayI(kMaxDDL);
	72	Int_t equipIn, equipOut;
	73	Int_t nIds = 0;
	74	while (input >> equipIn >> equipOut) {
	75	if (nIds >= kMaxDDL) {
	76	Error("AliRawReader","Too many equipment Id mappings found ! Truncating the list !");
	77	break;
	78	}
	79	fEquipmentIdsIn->AddAt(equipIn,nIds);
	80	fEquipmentIdsOut->AddAt(equipOut,nIds);
	81	nIds++;
	82	}
	83	fEquipmentIdsIn->Set(nIds);
	84	fEquipmentIdsOut->Set(nIds);
	85	input.close();
	86	return kTRUE;
	87	}
	88	else {
	89	Error("AliRawReader","equipment id map file is not found ! Skipping the mapping !");
	90	return kFALSE;
	91	}
a6e7b125	92	}
a6e7b125	93
42d20574	94	AliRawReader::AliRawReader(const AliRawReader& rawReader) :
c946ab02	95	TObject(rawReader),
0cfe2438	96	fEquipmentIdsIn(rawReader.fEquipmentIdsIn),
0cfe2438	97	fEquipmentIdsOut(rawReader.fEquipmentIdsOut),
18882b3b	98	fRequireHeader(rawReader.fRequireHeader),
39f9963f	99	fHeader(rawReader.fHeader),
c946ab02	100	fCount(rawReader.fCount),
f224de6c	101	fSelectEquipmentType(rawReader.fSelectEquipmentType),
	102	fSelectMinEquipmentId(rawReader.fSelectMinEquipmentId),
	103	fSelectMaxEquipmentId(rawReader.fSelectMaxEquipmentId),
39f9963f	104	fSkipInvalid(rawReader.fSkipInvalid),
dd9a70fe	105	fSelectEventType(rawReader.fSelectEventType),
c946ab02	106	fErrorCode(0)
42d20574	107	{
42d20574	108	// copy constructor
42d20574	109	}
	110
	111	AliRawReader& AliRawReader::operator = (const AliRawReader& rawReader)
	112	{
	113	// assignment operator
0cfe2438	114	fEquipmentIdsIn = rawReader.fEquipmentIdsIn;
0cfe2438	115	fEquipmentIdsOut = rawReader.fEquipmentIdsOut;
42d20574	116
39f9963f	117	fHeader = rawReader.fHeader;
42d20574	118	fCount = rawReader.fCount;
42d20574	119
39f9963f	120	fSelectEquipmentType = rawReader.fSelectEquipmentType;
	121	fSelectMinEquipmentId = rawReader.fSelectMinEquipmentId;
	122	fSelectMaxEquipmentId = rawReader.fSelectMaxEquipmentId;
	123	fSkipInvalid = rawReader.fSkipInvalid;
dd9a70fe	124	fSelectEventType = rawReader.fSelectEventType;
42d20574	125
b4857df7	126	fErrorCode = rawReader.fErrorCode;
b4857df7	127
42d20574	128	return *this;
	129	}
	130
0cfe2438	131	AliRawReader::~AliRawReader()
	132	{
	133	// destructor
	134	// delete the mapping arrays if
	135	// initialized
	136	if (fEquipmentIdsIn) delete fEquipmentIdsIn;
	137	if (fEquipmentIdsOut) delete fEquipmentIdsOut;
	138	}
	139
	140	Int_t AliRawReader::GetMappedEquipmentId() const
	141	{
	142	if (!fEquipmentIdsIn \|\| !fEquipmentIdsOut) {
	143	Error("AliRawReader","equipment Ids mapping is not initialized !");
	144	return GetEquipmentId();
	145	}
	146	Int_t equipmentId = GetEquipmentId();
	147	for(Int_t iId = 0; iId < fEquipmentIdsIn->GetSize(); iId++) {
	148	if (equipmentId == fEquipmentIdsIn->At(iId)) {
	149	equipmentId = fEquipmentIdsOut->At(iId);
	150	break;
	151	}
	152	}
	153	return equipmentId;
	154	}
	155
	156	Int_t AliRawReader::GetDetectorID() const
	157	{
	158	// Get the detector ID
	159	// The list of detector IDs
362c9d61	160	// can be found in AliDAQ.h
0cfe2438	161	Int_t equipmentId;
	162	if (fEquipmentIdsIn && fEquipmentIdsIn)
	163	equipmentId = GetMappedEquipmentId();
	164	else
	165	equipmentId = GetEquipmentId();
	166
362c9d61	167	if (equipmentId >= 0) {
	168	Int_t ddlIndex;
	169	return AliDAQ::DetectorIDFromDdlID(equipmentId,ddlIndex);
	170	}
0cfe2438	171	else
	172	return -1;
	173	}
	174
	175	Int_t AliRawReader::GetDDLID() const
	176	{
	177	// Get the DDL ID (within one sub-detector)
	178	// The list of detector IDs
362c9d61	179	// can be found in AliDAQ.h
0cfe2438	180	Int_t equipmentId;
	181	if (fEquipmentIdsIn && fEquipmentIdsIn)
	182	equipmentId = GetMappedEquipmentId();
	183	else
	184	equipmentId = GetEquipmentId();
	185
362c9d61	186	if (equipmentId >= 0) {
	187	Int_t ddlIndex;
	188	AliDAQ::DetectorIDFromDdlID(equipmentId,ddlIndex);
	189	return ddlIndex;
	190	}
0cfe2438	191	else
	192	return -1;
	193	}
8de97894	194
362c9d61	195	void AliRawReader::Select(const char *detectorName, Int_t minDDLID, Int_t maxDDLID)
	196	{
	197	// read only data of the detector with the given name and in the given
	198	// range of DDLs (minDDLID <= DDLID <= maxDDLID).
	199	// no selection is applied if a value < 0 is used.
	200	Int_t detectorID = AliDAQ::DetectorID(detectorName);
	201	if(detectorID >= 0)
	202	Select(detectorID,minDDLID,maxDDLID);
	203	}
	204
8de97894	205	void AliRawReader::Select(Int_t detectorID, Int_t minDDLID, Int_t maxDDLID)
a6e7b125	206	{
8de97894	207	// read only data of the detector with the given ID and in the given
f224de6c	208	// range of DDLs (minDDLID <= DDLID <= maxDDLID).
8de97894	209	// no selection is applied if a value < 0 is used.
a6e7b125	210
0cfe2438	211	fSelectEquipmentType = -1;
362c9d61	212
	213	if (minDDLID < 0)
	214	fSelectMinEquipmentId = AliDAQ::DdlIDOffset(detectorID);
	215	else
	216	fSelectMinEquipmentId = AliDAQ::DdlID(detectorID,minDDLID);
	217
	218	if (maxDDLID < 0)
	219	fSelectMaxEquipmentId = AliDAQ::DdlID(detectorID,AliDAQ::NumberOfDdls(detectorID)-1);
	220	else
	221	fSelectMaxEquipmentId = AliDAQ::DdlID(detectorID,maxDDLID);
8de97894	222	}
a6e7b125	223
f224de6c	224	void AliRawReader::SelectEquipment(Int_t equipmentType,
	225	Int_t minEquipmentId, Int_t maxEquipmentId)
	226	{
	227	// read only data of the equipment with the given type and in the given
	228	// range of IDs (minEquipmentId <= EquipmentId <= maxEquipmentId).
	229	// no selection is applied if a value < 0 is used.
	230
	231	fSelectEquipmentType = equipmentType;
	232	fSelectMinEquipmentId = minEquipmentId;
	233	fSelectMaxEquipmentId = maxEquipmentId;
	234	}
	235
dd9a70fe	236	void AliRawReader::SelectEvents(Int_t type)
	237	{
	238	// read only events with the given type.
	239	// no selection is applied if a value < 0 is used.
	240
	241	fSelectEventType = type;
	242	}
	243
42d20574	244	Bool_t AliRawReader::IsSelected() const
a6e7b125	245	{
8de97894	246	// apply the selection (if any)
8de97894	247
39f9963f	248	if (fSkipInvalid && !IsValid()) return kFALSE;
f224de6c	249
0cfe2438	250	if (fSelectEquipmentType >= 0)
f224de6c	251	if (GetEquipmentType() != fSelectEquipmentType) return kFALSE;
0cfe2438	252
	253	Int_t equipmentId;
	254	if (fEquipmentIdsIn && fEquipmentIdsIn)
	255	equipmentId = GetMappedEquipmentId();
	256	else
	257	equipmentId = GetEquipmentId();
	258
	259	if ((fSelectMinEquipmentId >= 0) &&
	260	(equipmentId < fSelectMinEquipmentId))
	261	return kFALSE;
	262	if ((fSelectMaxEquipmentId >= 0) &&
	263	(equipmentId > fSelectMaxEquipmentId))
	264	return kFALSE;
f224de6c	265
8de97894	266	return kTRUE;
a6e7b125	267	}
a6e7b125	268
dd9a70fe	269	Bool_t AliRawReader::IsEventSelected() const
	270	{
	271	// apply the event selection (if any)
	272
	273	if (fSelectEventType >= 0) {
	274	if (GetType() != (UInt_t) fSelectEventType) return kFALSE;
	275	}
	276
	277	return kTRUE;
	278	}
	279
a6e7b125	280
a6e7b125	281	Bool_t AliRawReader::ReadNextInt(UInt_t& data)
a6e7b125	282	{
8de97894	283	// reads the next 4 bytes at the current position
8de97894	284	// returns kFALSE if the data could not be read
a6e7b125	285
a6e7b125	286	while (fCount == 0) {
c946ab02	287	if (!ReadHeader()) return kFALSE;
a6e7b125	288	}
8de97894	289	if (fCount < (Int_t) sizeof(data)) {
	290	Error("ReadNextInt",
	291	"too few data left (%d bytes) to read an UInt_t!", fCount);
	292	return kFALSE;
	293	}
	294	if (!ReadNext((UChar_t*) &data, sizeof(data))) {
a6e7b125	295	Error("ReadNextInt", "could not read data!");
	296	return kFALSE;
	297	}
a6e7b125	298	return kTRUE;
	299	}
	300
	301	Bool_t AliRawReader::ReadNextShort(UShort_t& data)
	302	{
8de97894	303	// reads the next 2 bytes at the current position
8de97894	304	// returns kFALSE if the data could not be read
a6e7b125	305
a6e7b125	306	while (fCount == 0) {
c946ab02	307	if (!ReadHeader()) return kFALSE;
a6e7b125	308	}
8de97894	309	if (fCount < (Int_t) sizeof(data)) {
	310	Error("ReadNextShort",
	311	"too few data left (%d bytes) to read an UShort_t!", fCount);
	312	return kFALSE;
	313	}
	314	if (!ReadNext((UChar_t*) &data, sizeof(data))) {
a6e7b125	315	Error("ReadNextShort", "could not read data!");
	316	return kFALSE;
	317	}
a6e7b125	318	return kTRUE;
	319	}
	320
	321	Bool_t AliRawReader::ReadNextChar(UChar_t& data)
	322	{
	323	// reads the next 1 byte at the current stream position
8de97894	324	// returns kFALSE if the data could not be read
a6e7b125	325
a6e7b125	326	while (fCount == 0) {
c946ab02	327	if (!ReadHeader()) return kFALSE;
a6e7b125	328	}
8de97894	329	if (!ReadNext((UChar_t*) &data, sizeof(data))) {
a6e7b125	330	Error("ReadNextChar", "could not read data!");
	331	return kFALSE;
	332	}
a6e7b125	333	return kTRUE;
	334	}
	335
b4857df7	336
	337	Int_t AliRawReader::CheckData() const
	338	{
	339	// check the consistency of the data
	340	// derived classes should overwrite the default method which returns 0 (no err)
	341
	342	return 0;
	343	}
	344
509a7696	345
	346	void AliRawReader::DumpData(Int_t limit)
	347	{
	348	// print the raw data
	349	// if limit is not negative, only the first and last "limit" lines of raw data
	350	// are printed
	351
	352	Reset();
	353	if (!ReadHeader()) {
	354	Error("DumpData", "no header");
	355	return;
	356	}
	357	printf("header:\n"
	358	" type = %d run = %d ", GetType(), GetRunNumber());
	359	if (GetEventId()) {
	360	printf("event = %8.8x %8.8x\n", GetEventId()[1], GetEventId()[0]);
	361	} else {
	362	printf("event = -------- --------\n");
	363	}
	364	if (GetTriggerPattern()) {
	365	printf(" trigger = %8.8x %8.8x ",
	366	GetTriggerPattern()[1], GetTriggerPattern()[0]);
	367	} else {
	368	printf(" trigger = -------- -------- ");
	369	}
	370	if (GetDetectorPattern()) {
	371	printf("detector = %8.8x\n", GetDetectorPattern()[0]);
	372	} else {
	373	printf("detector = --------\n");
	374	}
	375	if (GetAttributes()) {
	376	printf(" attributes = %8.8x %8.8x %8.8x ",
	377	GetAttributes()[2], GetAttributes()[1], GetAttributes()[0]);
	378	} else {
	379	printf(" attributes = -------- -------- -------- ");
	380	}
	381	printf("GDC = %d\n", GetGDCId());
	382	printf("\n");
	383
	384	do {
	385	printf("-------------------------------------------------------------------------------\n");
	386	printf("LDC = %d\n", GetLDCId());
	387
	388	printf("equipment:\n"
	389	" size = %d type = %d id = %d\n",
	390	GetEquipmentSize(), GetEquipmentType(), GetEquipmentId());
	391	if (GetEquipmentAttributes()) {
	392	printf(" attributes = %8.8x %8.8x %8.8x ", GetEquipmentAttributes()[2],
	393	GetEquipmentAttributes()[1], GetEquipmentAttributes()[0]);
	394	} else {
	395	printf(" attributes = -------- -------- -------- ");
	396	}
	397	printf("element size = %d\n", GetEquipmentElementSize());
	398
39f9963f	399	printf("data header:\n"
	400	" size = %d version = %d valid = %d compression = %d\n",
	401	GetDataSize(), GetVersion(), IsValid(), IsCompressed());
509a7696	402
	403	printf("\n");
	404	if (limit == 0) continue;
	405
	406	Int_t size = GetDataSize();
	407	char line[70];
	408	for (Int_t i = 0; i < 70; i++) line[i] = ' ';
	409	line[69] = '\0';
	410	Int_t pos = 0;
	411	Int_t max = 16;
	412	UChar_t byte;
	413
	414	for (Int_t n = 0; n < size; n++) {
	415	if (!ReadNextChar(byte)) {
	416	Error("DumpData", "couldn't read byte number %d\n", n);
	417	break;
	418	}
	419	if (pos >= max) {
	420	printf("%8.8x %s\n", n-pos, line);
	421	for (Int_t i = 0; i < 70; i++) line[i] = ' ';
	422	line[69] = '\0';
	423	pos = 0;
	424	if ((limit > 0) && (n/max == limit)) {
	425	Int_t nContinue = ((size-1)/max+1-limit) * max;
	426	if (nContinue > n) {
	427	printf(" [skipping %d bytes]\n", nContinue-n);
	428	n = nContinue-1;
	429	continue;
	430	}
	431	}
	432	}
	433	Int_t offset = pos/4;
	434	if ((byte > 0x20) && (byte < 0x7f)) {
	435	line[pos+offset] = byte;
	436	} else {
	437	line[pos+offset] = '.';
	438	}
	439	char hex[3];
	440	sprintf(hex, "%2.2x", byte);
	441	line[max+max/4+3+2*pos+offset] = hex[0];
	442	line[max+max/4+4+2*pos+offset] = hex[1];
	443	pos++;
	444	}
	445
	446	if (pos > 0) printf("%8.8x %s\n", size-pos, line);
	447	printf("\n");
	448
	449	} while (ReadHeader());
	450	}