[u/mrichter/AliRoot.git] / HLT / RCU / AliHLTAltroChannelSelectorComponent.cxx

// $Id$

//**************************************************************************
//* This file is property of and copyright by the ALICE HLT Project        * 
//* ALICE Experiment at CERN, All rights reserved.                         *
//*                                                                        *
//* Primary Authors: Matthias Richter <Matthias.Richter@ift.uib.no>        *
//*                  for The ALICE HLT Project.                            *
//*                                                                        *
//* 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.                  *
//**************************************************************************

/// @file   AliHLTAltroChannelSelectorComponent.cxx
/// @author Matthias Richter edited by Jason Glyndwr Ulery
/// @date   
/// @brief  A filter/selective readout component for Altro data.
///

#include <cassert>
#include <memory>
#include "AliHLTAltroChannelSelectorComponent.h"
#include "AliHLTErrorGuard.h"
#include "AliHLTDAQ.h"
#include "AliRawReaderMemory.h"
#include "AliAltroRawStreamV3.h"
#include "TMath.h"

/** ROOT macro for the implementation of ROOT specific class methods */
ClassImp(AliHLTAltroChannelSelectorComponent)

AliHLTAltroChannelSelectorComponent::AliHLTAltroChannelSelectorComponent()
  :
  AliHLTProcessor(),
  fSkipCorrupted(true),
  fTalkative(false),
  fStartTimeBin(0),
  fEndTimeBin(1024),
  fSignalThreshold(0),
  fRMSThreshold(0)
{
  // see header file for class documentation
  // or
  // refer to README to build package
  // or
  // visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
}

AliHLTAltroChannelSelectorComponent::~AliHLTAltroChannelSelectorComponent()
{
  // see header file for class documentation
}

const char* AliHLTAltroChannelSelectorComponent::GetComponentID()
{
  // see header file for class documentation
  return "AltroChannelSelector";
}

void AliHLTAltroChannelSelectorComponent::GetInputDataTypes(AliHLTComponentDataTypeList& list)
{
  // see header file for class documentation
  list.clear();
  list.push_back(kAliHLTDataTypeDDLRaw|kAliHLTDataOriginTPC);
  list.push_back(kAliHLTDataTypeHwAddr16);
}

AliHLTComponentDataType AliHLTAltroChannelSelectorComponent::GetOutputDataType()
{
  // see header file for class documentation
  return kAliHLTDataTypeDDLRaw|kAliHLTDataOriginTPC;
}

void AliHLTAltroChannelSelectorComponent::GetOutputDataSize(unsigned long& constBase, double& inputMultiplier)
{
  // see header file for class documentation
  constBase=0;
  inputMultiplier=1.0;
}

AliHLTComponent* AliHLTAltroChannelSelectorComponent::Spawn()
{
  // see header file for class documentation
  return new AliHLTAltroChannelSelectorComponent;
}

int AliHLTAltroChannelSelectorComponent::DoInit(int argc, const char** argv)
{
  // see header file for class documentation
  int iResult=0;
  TString argument="";
  bool bMissingParam=0;
  char* cpErr=NULL;
  int i=0;
  for (; i<argc && iResult>=0; i++) {
    cpErr=NULL;
    argument=argv[i];
    if (argument.IsNull()) continue;

    // -skip-corrupted, just for backward compatibility, not announced
    if (argument.CompareTo("-skip-corrupted")==0) {
      fSkipCorrupted=true;

    // -keep-corrupted
    } else if (argument.CompareTo("-keep-corrupted")==0) {
      fSkipCorrupted=false;

    // -talkative
    } else if (argument.CompareTo("-talkative")==0) {
      fTalkative=true;

    // -start-timebin
    } else if (argument.CompareTo("-start-timebin")==0) {
      if ((bMissingParam=(++i>=argc))) break;
      fStartTimeBin = strtoul( argv[i], &cpErr ,0);
      if ( *cpErr ) break;

    // -end-timebin
    } else if (argument.CompareTo("-end-timebin")==0) {
      if ((bMissingParam=(++i>=argc))) break;
      fEndTimeBin = strtoul( argv[i], &cpErr ,0);
      if ( *cpErr ) break;

    // -signal-threshold
    } else if (argument.CompareTo("-signal-threshold")==0) {
      if ((bMissingParam=(++i>=argc))) break;
      fSignalThreshold = strtoul( argv[i], &cpErr ,0);
      if ( *cpErr ) break;

    // -rms-threshold
    } else if (argument.CompareTo("-rms-threshold")==0) {
      if ((bMissingParam=(++i>=argc))) break;
      fRMSThreshold = strtoul( argv[i], &cpErr ,0);
      if ( *cpErr ) break;

    } else {
      HLTError("unknown argument %s", argument.Data());
      iResult=-EINVAL;
    }
  }

  if (cpErr && *cpErr) {
    HLTError("Cannot convert specifier '%s' for argument '%s'", argv[i], argument.Data());
    iResult=-EINVAL;
  } else if (bMissingParam) {
    HLTError("missing parameter for argument %s", argument.Data());
    iResult=-EINVAL;
  }

  return iResult;
}

int AliHLTAltroChannelSelectorComponent::DoDeinit()
{
  // see header file for class documentation
  return 0;
}

int AliHLTAltroChannelSelectorComponent::DoEvent(const AliHLTComponentEventData& evtData,
						 const AliHLTComponentBlockData* blocks, 
						 AliHLTComponentTriggerData& /*trigData*/,
						 AliHLTUInt8_t* outputPtr, 
						 AliHLTUInt32_t& size,
						 AliHLTComponentBlockDataList& outputBlocks )
{
  // see header file for class documentation
  int iResult=0;
  AliHLTUInt32_t iCapacity=size;
  size=0;

  const int cdhSize=32;//8 32-bit words so 32 bytes
  const int rcuSize=36;//9 32-bit words

  if (!IsDataEvent()) {
    size=0;
    return 0;
  }

  // process the DLL input
  int blockno=0;
  const AliHLTComponentBlockData* pDesc=NULL;
  std::auto_ptr<AliRawReaderMemory> pRawReader(new AliRawReaderMemory);
  if (!pRawReader.get()) return -ENOMEM;

  for (pDesc=GetFirstInputBlock(kAliHLTDataTypeDDLRaw); pDesc!=NULL; pDesc=GetNextInputBlock(), blockno++) {
    iResult=0;
    if (pDesc->fSize<=32) {
      continue;
    }

    // search for the active pad information
    AliHLTUInt16_t* pActiveHwAddressArray=NULL;
    int iArraySize=0;
    if (fSignalThreshold==0 && fRMSThreshold==0) {
    for (int i=0; i<(int)evtData.fBlockCnt; i++ ) {
      // search for selection data of hw address type
      // which matches the data specification of the block
      if (blocks[i].fDataType==kAliHLTDataTypeHwAddr16 && blocks[i].fSpecification==pDesc->fSpecification) {
	pActiveHwAddressArray=reinterpret_cast<AliHLTUInt16_t*>(blocks[i].fPtr);
	iArraySize=blocks[i].fSize/sizeof(AliHLTUInt16_t);
	break;
      }
    }
    if (pActiveHwAddressArray==NULL) {
      HLTWarning("no block of type %s for specification 0x%08x available, data block skipped", 
		 DataType2Text(kAliHLTDataTypeHwAddr16).c_str(), 
		 pDesc->fSpecification);
      break;
    }
    }

    pRawReader->Reset();
    int ddlid=AliHLTDAQ::DdlIDFromHLTBlockData(pDesc->fDataType.fOrigin, pDesc->fSpecification);
    if (ddlid<0) {
      HLTError("unable to extract DDL Id for data block %s 0x%08x", DataType2Text(pDesc->fDataType).c_str(), pDesc->fSpecification);
      continue;
    }

    if (!pRawReader->AddBuffer((UChar_t*)pDesc->fPtr,pDesc->fSize, ddlid)) {
      ALIHLTERRORGUARD(1, "can not set up AltroDecoder for data block %s 0x%08x,"
		       " skipping data block and suppressing further messages",
		       DataType2Text(pDesc->fDataType).c_str(), pDesc->fSpecification);
      continue;
    }

    std::auto_ptr<AliAltroRawStreamV3> altroRawStream(new AliAltroRawStreamV3(pRawReader.get()));

    if (!altroRawStream.get()) {
      iResult=-ENOMEM;
      break;
    }

    altroRawStream->Reset();
    if (!altroRawStream->NextDDL()) {
      ALIHLTERRORGUARD(1, "internal error, can not read data from AliRawReaderMemory");
      continue;
    }

    unsigned int rcuTrailerLength=rcuSize;
    if(rcuTrailerLength>pDesc->fSize-cdhSize) HLTWarning("corrupted data block: RCU trailer length exceeds buffer size");
  
    if (iResult<0) {
      // TODO: here the trigger has to come into play. It is up to
      // policy if a corrupted data block should be kept (original
      // DDL) or discarded. In any case, the block is not going to
      // be part of HLTOUT
      HLTWarning("skipping corrupted data block for event %lu, data block %s 0x%08x", evtData.fEventID,
		 DataType2Text(pDesc->fDataType).c_str(), pDesc->fSpecification);
      iResult=0;
      continue;
    }

    int iSelected=0;
    int iTotal=0;
    int iCorrupted=0;
    AliHLTUInt32_t iOutputSize=0;

    //First add the common data header (CDH)
    memcpy(outputPtr+size,pDesc->fPtr,cdhSize);
    iOutputSize+=cdhSize;

    while (altroRawStream->NextChannel() && iResult>=0) {
      iTotal++;
      int hwAddress=altroRawStream->GetHWAddress();
  
      if (fSignalThreshold!=0) {
	// treshold by adc counts
	unsigned int sumSignals=0;
	unsigned int maxSignal=0;
	unsigned int nofSignals=0;
	while(altroRawStream->NextBunch()){
	  const UShort_t *bunchData=altroRawStream->GetSignals();
	  unsigned int time=altroRawStream->GetStartTimeBin();
	  for(Int_t i=0;i<altroRawStream->GetBunchLength();i++){
	    if(bunchData[i]>0){// disregarding 0 data.
	      if(time+i>=fStartTimeBin && time+i<=fEndTimeBin){
		sumSignals+=bunchData[i];
		if (maxSignal<bunchData[i]) maxSignal=bunchData[i];
		nofSignals++;
	      }
	    }
	  }
	}
	if (nofSignals==0 || maxSignal<=(sumSignals/nofSignals)+fSignalThreshold) {
	  continue;
	}

      } else if (fRMSThreshold!=0) {
	// treshold by adc counts
	unsigned int sumSignals=0;
	unsigned int maxSignal=0;
	unsigned int nofSignals=0;
	while(altroRawStream->NextBunch()){
	  const UShort_t *bunchData=altroRawStream->GetSignals();
	  unsigned int time=altroRawStream->GetStartTimeBin();
	  for(Int_t i=0;i<altroRawStream->GetBunchLength();i++){
	    if(bunchData[i]>0){// disregarding 0 data.
	      if(time+i>=fStartTimeBin && time+i<=fEndTimeBin){
		sumSignals+=bunchData[i]*bunchData[i];
		if (maxSignal<bunchData[i]) maxSignal=bunchData[i];
		nofSignals++;
	      }
	    }
	  }
	}
	if (nofSignals==0 || maxSignal<=TMath::Sqrt(sumSignals/nofSignals)*fRMSThreshold) {
	  continue;
	}
	
      } else {
      int active=0;
      for (active=0; active<iArraySize; active++) {
	if (pActiveHwAddressArray[active]==(AliHLTUInt16_t)hwAddress) {
	  break;
	}
      }
      if (active>=iArraySize) {
	HLTDebug("ALTRO block %#x (%d) discarded (inactive)", hwAddress, hwAddress);
	continue;
      }
      }

      // no of 10 bit words is without the fill words to fill complete 40 bit words
      // in addition, align to complete 40 bit words (the '+3')
      // also, the 4 bytes of the Altro trailer must be added to get the full size
      int channelSize=((altroRawStream->GetChannelPayloadSize()+2)/3)*4;
      if (channelSize==0) {
	if (fTalkative) HLTWarning("skipping zero length channel (hw address %d)", hwAddress);
	iCorrupted++;
	continue;
      }
      channelSize+=4;
      HLTDebug("ALTRO block hwAddress 0x%08x (%d) selected (active), size %d", hwAddress, hwAddress, channelSize);

      if (iCapacity>=(channelSize+iOutputSize)) {
	const UChar_t *ChannelData=altroRawStream->GetChannelPayload();
	iResult=channelSize;
	memcpy(outputPtr+size+iOutputSize,ChannelData,channelSize);
	if (channelSize == iResult){
	  if (channelSize%4 == 0) {
	    iOutputSize+=channelSize;
	  } else {
	    if (fTalkative) HLTWarning("corrupted ALTRO channel: incomplete 32 bit word (channel hw address %d)", hwAddress);
	    iCorrupted++;
	    continue;
	  }
	} else {
	  if (fTalkative) HLTWarning("internal error: failed to copy full channel: %d out of %d bytes (hw address %d)", iResult, channelSize, hwAddress);
	  iCorrupted++;
	  continue;
	}
      } else {
	if (fTalkative) HLTError("failed to write ALTRO channel of length %d for block %d  (hw address %d)", channelSize, blockno, hwAddress);
	// corrupted channel, but keep going
	iCorrupted++;
	iResult=0;
	continue;
      }  
      iSelected++;
    }
    
    if (iResult>=0) {
      //Write the RCU Trailer
      AliHLTUInt8_t* pSrc=reinterpret_cast<AliHLTUInt8_t*>(pDesc->fPtr);
      pSrc+=pDesc->fSize-rcuTrailerLength;
      memcpy(outputPtr+size+iOutputSize,pSrc,rcuTrailerLength);
      iOutputSize+=rcuTrailerLength;

      //Set new payload length of the new data bloock 
      AliHLTUInt32_t* pSize=reinterpret_cast<AliHLTUInt32_t*>(outputPtr+size+iOutputSize-rcuTrailerLength);
      (*pSize)&=~0x3FFFFFF;
      (*pSize)|=(iOutputSize-rcuTrailerLength-cdhSize+3)/4;//# of 32 bit words in payload
      
      //Write DDL length     
      AliHLTUInt32_t* pCdhSize=reinterpret_cast<AliHLTUInt32_t*>(outputPtr+size);
      *pCdhSize=iOutputSize;

      // insert block descriptor
      AliHLTComponentBlockData bd;
      FillBlockData(bd);
      bd.fOffset=size;
      bd.fSize=iOutputSize;
      bd.fDataType=pDesc->fDataType;
      bd.fSpecification=pDesc->fSpecification;
      outputBlocks.push_back(bd);
      size+=iOutputSize;
    }
    if (fTalkative) HLTImportant("data block %d (0x%08x): selected %d out of %d ALTRO channel(s), %d corrupted channels skipped", blockno, pDesc->fSpecification, iSelected, iTotal, iCorrupted);
  }

  if (iResult<0) {
    outputBlocks.clear();
  }
  return iResult;
}
Commit	Line	Data
3e6ec852	1	// $Id$
3e6ec852	2
13398559	3	//**************************************************************************
	4	//* This file is property of and copyright by the ALICE HLT Project *
	5	//* ALICE Experiment at CERN, All rights reserved. *
	6	//* *
	7	//* Primary Authors: Matthias Richter <Matthias.Richter@ift.uib.no> *
	8	//* for The ALICE HLT Project. *
	9	//* *
	10	//* Permission to use, copy, modify and distribute this software and its *
	11	//* documentation strictly for non-commercial purposes is hereby granted *
	12	//* without fee, provided that the above copyright notice appears in all *
	13	//* copies and that both the copyright notice and this permission notice *
	14	//* appear in the supporting documentation. The authors make no claims *
	15	//* about the suitability of this software for any purpose. It is *
	16	//* provided "as is" without express or implied warranty. *
	17	//**************************************************************************
3e6ec852	18
ff4a65b9	19	/// @file AliHLTAltroChannelSelectorComponent.cxx
1ea9f419	20	/// @author Matthias Richter edited by Jason Glyndwr Ulery
ff4a65b9	21	/// @date
	22	/// @brief A filter/selective readout component for Altro data.
	23	///
3e6ec852	24
e1440dab	25	#include <cassert>
ff4a65b9	26	#include <memory>
3e6ec852	27	#include "AliHLTAltroChannelSelectorComponent.h"
ff4a65b9	28	#include "AliHLTErrorGuard.h"
	29	#include "AliHLTDAQ.h"
	30	#include "AliRawReaderMemory.h"
	31	#include "AliAltroRawStreamV3.h"
2a32beb8	32	#include "TMath.h"
3e6ec852	33
	34	/** ROOT macro for the implementation of ROOT specific class methods */
	35	ClassImp(AliHLTAltroChannelSelectorComponent)
	36
	37	AliHLTAltroChannelSelectorComponent::AliHLTAltroChannelSelectorComponent()
	38	:
e1440dab	39	AliHLTProcessor(),
2a32beb8	40	fSkipCorrupted(true),
	41	fTalkative(false),
	42	fStartTimeBin(0),
f5bd9c68	43	fEndTimeBin(1024),
2a32beb8	44	fSignalThreshold(0),
2a32beb8	45	fRMSThreshold(0)
3e6ec852	46	{
	47	// see header file for class documentation
	48	// or
	49	// refer to README to build package
	50	// or
	51	// visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
	52	}
	53
	54	AliHLTAltroChannelSelectorComponent::~AliHLTAltroChannelSelectorComponent()
	55	{
	56	// see header file for class documentation
	57	}
	58
	59	const char* AliHLTAltroChannelSelectorComponent::GetComponentID()
	60	{
	61	// see header file for class documentation
	62	return "AltroChannelSelector";
	63	}
	64
	65	void AliHLTAltroChannelSelectorComponent::GetInputDataTypes(AliHLTComponentDataTypeList& list)
	66	{
	67	// see header file for class documentation
	68	list.clear();
	69	list.push_back(kAliHLTDataTypeDDLRaw\|kAliHLTDataOriginTPC);
64defa03	70	list.push_back(kAliHLTDataTypeHwAddr16);
3e6ec852	71	}
	72
	73	AliHLTComponentDataType AliHLTAltroChannelSelectorComponent::GetOutputDataType()
	74	{
	75	// see header file for class documentation
	76	return kAliHLTDataTypeDDLRaw\|kAliHLTDataOriginTPC;
	77	}
	78
	79	void AliHLTAltroChannelSelectorComponent::GetOutputDataSize(unsigned long& constBase, double& inputMultiplier)
	80	{
	81	// see header file for class documentation
	82	constBase=0;
	83	inputMultiplier=1.0;
	84	}
	85
	86	AliHLTComponent* AliHLTAltroChannelSelectorComponent::Spawn()
	87	{
	88	// see header file for class documentation
	89	return new AliHLTAltroChannelSelectorComponent;
	90	}
	91
	92	int AliHLTAltroChannelSelectorComponent::DoInit(int argc, const char** argv)
	93	{
	94	// see header file for class documentation
	95	int iResult=0;
	96	TString argument="";
e1440dab	97	bool bMissingParam=0;
2a32beb8	98	char* cpErr=NULL;
	99	int i=0;
	100	for (; i<argc && iResult>=0; i++) {
	101	cpErr=NULL;
3e6ec852	102	argument=argv[i];
	103	if (argument.IsNull()) continue;
	104
2a32beb8	105	// -skip-corrupted, just for backward compatibility, not announced
4ca5858d	106	if (argument.CompareTo("-skip-corrupted")==0) {
	107	fSkipCorrupted=true;
	108
2a32beb8	109	// -keep-corrupted
	110	} else if (argument.CompareTo("-keep-corrupted")==0) {
	111	fSkipCorrupted=false;
	112
4ca5858d	113	// -talkative
	114	} else if (argument.CompareTo("-talkative")==0) {
	115	fTalkative=true;
2a32beb8	116
	117	// -start-timebin
	118	} else if (argument.CompareTo("-start-timebin")==0) {
9d5f00ad	119	if ((bMissingParam=(++i>=argc))) break;
2a32beb8	120	fStartTimeBin = strtoul( argv[i], &cpErr ,0);
	121	if ( *cpErr ) break;
	122
	123	// -end-timebin
	124	} else if (argument.CompareTo("-end-timebin")==0) {
9d5f00ad	125	if ((bMissingParam=(++i>=argc))) break;
2a32beb8	126	fEndTimeBin = strtoul( argv[i], &cpErr ,0);
	127	if ( *cpErr ) break;
	128
	129	// -signal-threshold
	130	} else if (argument.CompareTo("-signal-threshold")==0) {
9d5f00ad	131	if ((bMissingParam=(++i>=argc))) break;
2a32beb8	132	fSignalThreshold = strtoul( argv[i], &cpErr ,0);
	133	if ( *cpErr ) break;
	134
	135	// -rms-threshold
	136	} else if (argument.CompareTo("-rms-threshold")==0) {
9d5f00ad	137	if ((bMissingParam=(++i>=argc))) break;
2a32beb8	138	fRMSThreshold = strtoul( argv[i], &cpErr ,0);
	139	if ( *cpErr ) break;
	140
3e6ec852	141	} else {
4ca5858d	142	HLTError("unknown argument %s", argument.Data());
3e6ec852	143	iResult=-EINVAL;
	144	}
	145	}
	146
2a32beb8	147	if (cpErr && *cpErr) {
	148	HLTError("Cannot convert specifier '%s' for argument '%s'", argv[i], argument.Data());
	149	iResult=-EINVAL;
	150	} else if (bMissingParam) {
e1440dab	151	HLTError("missing parameter for argument %s", argument.Data());
	152	iResult=-EINVAL;
	153	}
	154
3e6ec852	155	return iResult;
	156	}
	157
	158	int AliHLTAltroChannelSelectorComponent::DoDeinit()
	159	{
	160	// see header file for class documentation
	161	return 0;
	162	}
	163
	164	int AliHLTAltroChannelSelectorComponent::DoEvent(const AliHLTComponentEventData& evtData,
	165	const AliHLTComponentBlockData* blocks,
	166	AliHLTComponentTriggerData& /trigData/,
	167	AliHLTUInt8_t* outputPtr,
	168	AliHLTUInt32_t& size,
	169	AliHLTComponentBlockDataList& outputBlocks )
	170	{
	171	// see header file for class documentation
	172	int iResult=0;
1ea9f419	173	AliHLTUInt32_t iCapacity=size;
	174	size=0;
	175
	176	const int cdhSize=32;//8 32-bit words so 32 bytes
	177	const int rcuSize=36;//9 32-bit words
3e6ec852	178
8d682af4	179	if (!IsDataEvent()) {
	180	size=0;
	181	return 0;
	182	}
	183
e1440dab	184	// process the DLL input
e1440dab	185	int blockno=0;
1d337632	186	const AliHLTComponentBlockData* pDesc=NULL;
ff4a65b9	187	std::auto_ptr<AliRawReaderMemory> pRawReader(new AliRawReaderMemory);
fa7ed16b	188	if (!pRawReader.get()) return -ENOMEM;
1d337632	189
4ca5858d	190	for (pDesc=GetFirstInputBlock(kAliHLTDataTypeDDLRaw); pDesc!=NULL; pDesc=GetNextInputBlock(), blockno++) {
4ca5858d	191	iResult=0;
359ba057	192	if (pDesc->fSize<=32) {
	193	continue;
	194	}
e1440dab	195
e1440dab	196	// search for the active pad information
64defa03	197	AliHLTUInt16_t* pActiveHwAddressArray=NULL;
64defa03	198	int iArraySize=0;
2a32beb8	199	if (fSignalThreshold==0 && fRMSThreshold==0) {
e1440dab	200	for (int i=0; i<(int)evtData.fBlockCnt; i++ ) {
64defa03	201	// search for selection data of hw address type
4ca5858d	202	// which matches the data specification of the block
	203	if (blocks[i].fDataType==kAliHLTDataTypeHwAddr16 && blocks[i].fSpecification==pDesc->fSpecification) {
	204	pActiveHwAddressArray=reinterpret_cast<AliHLTUInt16_t*>(blocks[i].fPtr);
	205	iArraySize=blocks[i].fSize/sizeof(AliHLTUInt16_t);
	206	break;
e1440dab	207	}
e1440dab	208	}
4ca5858d	209	if (pActiveHwAddressArray==NULL) {
2a32beb8	210	HLTWarning("no block of type %s for specification 0x%08x available, data block skipped",
64defa03	211	DataType2Text(kAliHLTDataTypeHwAddr16).c_str(),
1d337632	212	pDesc->fSpecification);
2a32beb8	213	break;
2a32beb8	214	}
4ca5858d	215	}
4ca5858d	216
ff4a65b9	217	pRawReader->Reset();
	218	int ddlid=AliHLTDAQ::DdlIDFromHLTBlockData(pDesc->fDataType.fOrigin, pDesc->fSpecification);
	219	if (ddlid<0) {
	220	HLTError("unable to extract DDL Id for data block %s 0x%08x", DataType2Text(pDesc->fDataType).c_str(), pDesc->fSpecification);
	221	continue;
	222	}
	223
	224	if (!pRawReader->AddBuffer((UChar_t*)pDesc->fPtr,pDesc->fSize, ddlid)) {
	225	ALIHLTERRORGUARD(1, "can not set up AltroDecoder for data block %s 0x%08x,"
	226	" skipping data block and suppressing further messages",
	227	DataType2Text(pDesc->fDataType).c_str(), pDesc->fSpecification);
	228	continue;
	229	}
	230
	231	std::auto_ptr<AliAltroRawStreamV3> altroRawStream(new AliAltroRawStreamV3(pRawReader.get()));
	232
	233	if (!altroRawStream.get()) {
	234	iResult=-ENOMEM;
	235	break;
	236	}
	237
	238	altroRawStream->Reset();
	239	if (!altroRawStream->NextDDL()) {
	240	ALIHLTERRORGUARD(1, "internal error, can not read data from AliRawReaderMemory");
	241	continue;
4ca5858d	242	}
4ca5858d	243
1ea9f419	244	unsigned int rcuTrailerLength=rcuSize;
	245	if(rcuTrailerLength>pDesc->fSize-cdhSize) HLTWarning("corrupted data block: RCU trailer length exceeds buffer size");
	246
4ca5858d	247	if (iResult<0) {
8d682af4	248	// TODO: here the trigger has to come into play. It is up to
	249	// policy if a corrupted data block should be kept (original
	250	// DDL) or discarded. In any case, the block is not going to
	251	// be part of HLTOUT
	252	HLTWarning("skipping corrupted data block for event %lu, data block %s 0x%08x", evtData.fEventID,
	253	DataType2Text(pDesc->fDataType).c_str(), pDesc->fSpecification);
3dcae0cc	254	iResult=0;
e1440dab	255	continue;
	256	}
	257
1d337632	258	int iSelected=0;
1d337632	259	int iTotal=0;
4ca5858d	260	int iCorrupted=0;
e1440dab	261	AliHLTUInt32_t iOutputSize=0;
1ea9f419	262
	263	//First add the common data header (CDH)
	264	memcpy(outputPtr+size,pDesc->fPtr,cdhSize);
	265	iOutputSize+=cdhSize;
8d682af4	266
ff4a65b9	267	while (altroRawStream->NextChannel() && iResult>=0) {
1d337632	268	iTotal++;
ff4a65b9	269	int hwAddress=altroRawStream->GetHWAddress();
1ea9f419	270
2a32beb8	271	if (fSignalThreshold!=0) {
	272	// treshold by adc counts
	273	unsigned int sumSignals=0;
	274	unsigned int maxSignal=0;
	275	unsigned int nofSignals=0;
ff4a65b9	276	while(altroRawStream->NextBunch()){
	277	const UShort_t *bunchData=altroRawStream->GetSignals();
	278	unsigned int time=altroRawStream->GetStartTimeBin();
	279	for(Int_t i=0;i<altroRawStream->GetBunchLength();i++){
2a32beb8	280	if(bunchData[i]>0){// disregarding 0 data.
	281	if(time+i>=fStartTimeBin && time+i<=fEndTimeBin){
	282	sumSignals+=bunchData[i];
	283	if (maxSignal<bunchData[i]) maxSignal=bunchData[i];
	284	nofSignals++;
	285	}
	286	}
	287	}
	288	}
	289	if (nofSignals==0 \|\| maxSignal<=(sumSignals/nofSignals)+fSignalThreshold) {
	290	continue;
	291	}
	292
	293	} else if (fRMSThreshold!=0) {
	294	// treshold by adc counts
	295	unsigned int sumSignals=0;
	296	unsigned int maxSignal=0;
	297	unsigned int nofSignals=0;
ff4a65b9	298	while(altroRawStream->NextBunch()){
	299	const UShort_t *bunchData=altroRawStream->GetSignals();
	300	unsigned int time=altroRawStream->GetStartTimeBin();
	301	for(Int_t i=0;i<altroRawStream->GetBunchLength();i++){
2a32beb8	302	if(bunchData[i]>0){// disregarding 0 data.
	303	if(time+i>=fStartTimeBin && time+i<=fEndTimeBin){
	304	sumSignals+=bunchData[i]*bunchData[i];
	305	if (maxSignal<bunchData[i]) maxSignal=bunchData[i];
	306	nofSignals++;
	307	}
	308	}
	309	}
	310	}
	311	if (nofSignals==0 \|\| maxSignal<=TMath::Sqrt(sumSignals/nofSignals)*fRMSThreshold) {
	312	continue;
	313	}
	314
	315	} else {
e1440dab	316	int active=0;
4ca5858d	317	for (active=0; active<iArraySize; active++) {
	318	if (pActiveHwAddressArray[active]==(AliHLTUInt16_t)hwAddress) {
	319	break;
e1440dab	320	}
e1440dab	321	}
64defa03	322	if (active>=iArraySize) {
4ca5858d	323	HLTDebug("ALTRO block %#x (%d) discarded (inactive)", hwAddress, hwAddress);
e1440dab	324	continue;
e1440dab	325	}
2a32beb8	326	}
e1440dab	327
4ca5858d	328	// no of 10 bit words is without the fill words to fill complete 40 bit words
4ca5858d	329	// in addition, align to complete 40 bit words (the '+3')
1ea9f419	330	// also, the 4 bytes of the Altro trailer must be added to get the full size
ff4a65b9	331	int channelSize=((altroRawStream->GetChannelPayloadSize()+2)/3)*4;
8d682af4	332	if (channelSize==0) {
	333	if (fTalkative) HLTWarning("skipping zero length channel (hw address %d)", hwAddress);
	334	iCorrupted++;
	335	continue;
	336	}
ff4a65b9	337	channelSize+=4;
4ca5858d	338	HLTDebug("ALTRO block hwAddress 0x%08x (%d) selected (active), size %d", hwAddress, hwAddress, channelSize);
4ca5858d	339
1ea9f419	340	if (iCapacity>=(channelSize+iOutputSize)) {
	341	const UChar_t *ChannelData=altroRawStream->GetChannelPayload();
	342	iResult=channelSize;
	343	memcpy(outputPtr+size+iOutputSize,ChannelData,channelSize);
	344	if (channelSize == iResult){
ff4a65b9	345	if (channelSize%4 == 0) {
4ca5858d	346	iOutputSize+=channelSize;
e1440dab	347	} else {
1ea9f419	348	if (fTalkative) HLTWarning("corrupted ALTRO channel: incomplete 32 bit word (channel hw address %d)", hwAddress);
4ca5858d	349	iCorrupted++;
4ca5858d	350	continue;
e1440dab	351	}
4ca5858d	352	} else {
8d682af4	353	if (fTalkative) HLTWarning("internal error: failed to copy full channel: %d out of %d bytes (hw address %d)", iResult, channelSize, hwAddress);
4ca5858d	354	iCorrupted++;
4ca5858d	355	continue;
e1440dab	356	}
e1440dab	357	} else {
8d682af4	358	if (fTalkative) HLTError("failed to write ALTRO channel of length %d for block %d (hw address %d)", channelSize, blockno, hwAddress);
4ca5858d	359	// corrupted channel, but keep going
	360	iCorrupted++;
	361	iResult=0;
	362	continue;
1ea9f419	363	}
4ca5858d	364	iSelected++;
e1440dab	365	}
1ea9f419	366
e1440dab	367	if (iResult>=0) {
1ea9f419	368	//Write the RCU Trailer
	369	AliHLTUInt8_t* pSrc=reinterpret_cast<AliHLTUInt8_t*>(pDesc->fPtr);
	370	pSrc+=pDesc->fSize-rcuTrailerLength;
	371	memcpy(outputPtr+size+iOutputSize,pSrc,rcuTrailerLength);
	372	iOutputSize+=rcuTrailerLength;
	373
	374	//Set new payload length of the new data bloock
	375	AliHLTUInt32_t* pSize=reinterpret_cast<AliHLTUInt32_t*>(outputPtr+size+iOutputSize-rcuTrailerLength);
	376	(*pSize)&=~0x3FFFFFF;
	377	(*pSize)\|=(iOutputSize-rcuTrailerLength-cdhSize+3)/4;//# of 32 bit words in payload
	378
	379	//Write DDL length
	380	AliHLTUInt32_t* pCdhSize=reinterpret_cast<AliHLTUInt32_t*>(outputPtr+size);
	381	*pCdhSize=iOutputSize;
	382
	383	// insert block descriptor
	384	AliHLTComponentBlockData bd;
	385	FillBlockData(bd);
	386	bd.fOffset=size;
	387	bd.fSize=iOutputSize;
	388	bd.fDataType=pDesc->fDataType;
	389	bd.fSpecification=pDesc->fSpecification;
	390	outputBlocks.push_back(bd);
	391	size+=iOutputSize;
e1440dab	392	}
359ba057	393	if (fTalkative) HLTImportant("data block %d (0x%08x): selected %d out of %d ALTRO channel(s), %d corrupted channels skipped", blockno, pDesc->fSpecification, iSelected, iTotal, iCorrupted);
3e6ec852	394	}
3e6ec852	395
e1440dab	396	if (iResult<0) {
e1440dab	397	outputBlocks.clear();
3e6ec852	398	}
3e6ec852	399	return iResult;
3e6ec852	400	}