some quick hacks to improve speed for ActiveChannelSelection: terminate ZeroSuppressi...
[u/mrichter/AliRoot.git] / HLT / RCU / AliHLTAltroChannelSelectorComponent.cxx
CommitLineData
3e6ec852 1// $Id$
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
19/** @file AliHLTAltroChannelSelectorComponent.cxx
20 @author Matthias Richter
21 @date
13398559 22 @brief A filter/selective readout component for Altro data.
23*/
3e6ec852 24
25// see header file for class documentation
26// or
27// refer to README to build package
28// or
29// visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
30
e1440dab 31#include <cassert>
3e6ec852 32#include "AliHLTAltroChannelSelectorComponent.h"
4ca5858d 33#include "AliAltroDecoder.h"
34#include "AliAltroData.h"
3e6ec852 35
36/** ROOT macro for the implementation of ROOT specific class methods */
37ClassImp(AliHLTAltroChannelSelectorComponent)
38
39AliHLTAltroChannelSelectorComponent::AliHLTAltroChannelSelectorComponent()
40 :
e1440dab 41 AliHLTProcessor(),
4ca5858d 42 fSkipCorrupted(false),
43 fTalkative(false)
3e6ec852 44{
45 // see header file for class documentation
46 // or
47 // refer to README to build package
48 // or
49 // visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
50}
51
52AliHLTAltroChannelSelectorComponent::~AliHLTAltroChannelSelectorComponent()
53{
54 // see header file for class documentation
55}
56
57const char* AliHLTAltroChannelSelectorComponent::GetComponentID()
58{
59 // see header file for class documentation
60 return "AltroChannelSelector";
61}
62
63void AliHLTAltroChannelSelectorComponent::GetInputDataTypes(AliHLTComponentDataTypeList& list)
64{
65 // see header file for class documentation
66 list.clear();
67 list.push_back(kAliHLTDataTypeDDLRaw|kAliHLTDataOriginTPC);
64defa03 68 list.push_back(kAliHLTDataTypeHwAddr16);
3e6ec852 69}
70
71AliHLTComponentDataType AliHLTAltroChannelSelectorComponent::GetOutputDataType()
72{
73 // see header file for class documentation
74 return kAliHLTDataTypeDDLRaw|kAliHLTDataOriginTPC;
75}
76
77void AliHLTAltroChannelSelectorComponent::GetOutputDataSize(unsigned long& constBase, double& inputMultiplier)
78{
79 // see header file for class documentation
80 constBase=0;
81 inputMultiplier=1.0;
82}
83
84AliHLTComponent* AliHLTAltroChannelSelectorComponent::Spawn()
85{
86 // see header file for class documentation
87 return new AliHLTAltroChannelSelectorComponent;
88}
89
90int AliHLTAltroChannelSelectorComponent::DoInit(int argc, const char** argv)
91{
92 // see header file for class documentation
93 int iResult=0;
94 TString argument="";
e1440dab 95 bool bMissingParam=0;
3e6ec852 96 for (int i=0; i<argc && iResult>=0; i++) {
97 argument=argv[i];
98 if (argument.IsNull()) continue;
99
4ca5858d 100 // -skip-corrupted
101 if (argument.CompareTo("-skip-corrupted")==0) {
102 fSkipCorrupted=true;
103
104 // -talkative
105 } else if (argument.CompareTo("-talkative")==0) {
106 fTalkative=true;
3e6ec852 107 } else {
4ca5858d 108 HLTError("unknown argument %s", argument.Data());
3e6ec852 109 iResult=-EINVAL;
110 }
111 }
112
e1440dab 113 if (bMissingParam) {
114 HLTError("missing parameter for argument %s", argument.Data());
115 iResult=-EINVAL;
116 }
117
3e6ec852 118 return iResult;
119}
120
121int AliHLTAltroChannelSelectorComponent::DoDeinit()
122{
123 // see header file for class documentation
124 return 0;
125}
126
127int AliHLTAltroChannelSelectorComponent::DoEvent(const AliHLTComponentEventData& evtData,
128 const AliHLTComponentBlockData* blocks,
129 AliHLTComponentTriggerData& /*trigData*/,
130 AliHLTUInt8_t* outputPtr,
131 AliHLTUInt32_t& size,
132 AliHLTComponentBlockDataList& outputBlocks )
133{
134 // see header file for class documentation
135 int iResult=0;
4ca5858d 136 const int cdhSize=32;
3e6ec852 137
8d682af4 138 if (!IsDataEvent()) {
139 size=0;
140 return 0;
141 }
142
e1440dab 143 // process the DLL input
144 int blockno=0;
1d337632 145 const AliHLTComponentBlockData* pDesc=NULL;
146
4ca5858d 147 AliAltroDecoder* decoder=NULL;
148 for (pDesc=GetFirstInputBlock(kAliHLTDataTypeDDLRaw); pDesc!=NULL; pDesc=GetNextInputBlock(), blockno++) {
149 iResult=0;
740e48ad 150 if (pDesc->fSize<=32) continue;
e1440dab 151
152 // search for the active pad information
64defa03 153 AliHLTUInt16_t* pActiveHwAddressArray=NULL;
154 int iArraySize=0;
e1440dab 155 for (int i=0; i<(int)evtData.fBlockCnt; i++ ) {
64defa03 156 // search for selection data of hw address type
4ca5858d 157 // which matches the data specification of the block
158 if (blocks[i].fDataType==kAliHLTDataTypeHwAddr16 && blocks[i].fSpecification==pDesc->fSpecification) {
159 pActiveHwAddressArray=reinterpret_cast<AliHLTUInt16_t*>(blocks[i].fPtr);
160 iArraySize=blocks[i].fSize/sizeof(AliHLTUInt16_t);
161 break;
e1440dab 162 }
163 }
4ca5858d 164 if (pActiveHwAddressArray==NULL) {
165 HLTWarning("no block of type %s for specification 0x%08x available, data block unchanged",
64defa03 166 DataType2Text(kAliHLTDataTypeHwAddr16).c_str(),
1d337632 167 pDesc->fSpecification);
4ca5858d 168 iResult=-EFAULT;
169 }
170
171 if (decoder) delete decoder;
172 decoder=new AliAltroDecoder;
173 if (decoder->SetMemory(reinterpret_cast<UChar_t*>(pDesc->fPtr), pDesc->fSize)<0) {
174 HLTWarning("corrupted data block: initialization of decoder failed for block: %s specification %#x size %d",
175 DataType2Text(pDesc->fDataType).c_str(), pDesc->fSpecification, pDesc->fSize);
176 iResult=-EFAULT;
177 } else {
178 if (decoder->Decode()) {
179 HLTDebug("init decoder %p size %d", pDesc->fPtr,pDesc->fSize);
180 } else {
181 HLTWarning("corrupted data block: decoding failed for raw data block: %s specification %#x size %d",
182 DataType2Text(pDesc->fDataType).c_str(), pDesc->fSpecification, pDesc->fSize);
183 iResult=-EFAULT;
184 }
185 }
186
8d682af4 187 unsigned int rcuTrailerLength=0;
188 if (iResult>=0 &&
189 ((rcuTrailerLength=decoder->GetRCUTrailerSize())==0 ||
190 rcuTrailerLength>pDesc->fSize-cdhSize)) {
191 if (rcuTrailerLength>0) {
192 HLTWarning("corrupted data block: RCU trailer length exceeds buffer size");
193 } else {
194 HLTWarning("corrupted data block: RCU trailer of zero length");
195 }
4ca5858d 196 iResult=-EFAULT;
197 }
198
199 if (iResult<0) {
8d682af4 200 // TODO: here the trigger has to come into play. It is up to
201 // policy if a corrupted data block should be kept (original
202 // DDL) or discarded. In any case, the block is not going to
203 // be part of HLTOUT
204 HLTWarning("skipping corrupted data block for event %lu, data block %s 0x%08x", evtData.fEventID,
205 DataType2Text(pDesc->fDataType).c_str(), pDesc->fSpecification);
3dcae0cc 206 iResult=0;
e1440dab 207 continue;
208 }
209
1d337632 210 int iSelected=0;
211 int iTotal=0;
4ca5858d 212 int iCorrupted=0;
e1440dab 213 AliHLTUInt32_t iOutputSize=0;
4ca5858d 214 AliHLTUInt32_t iNofAltro40=0;
e1440dab 215 AliHLTUInt32_t iCapacity=size;
4ca5858d 216 AliAltroData channel;
8d682af4 217
218 // first add the RCU trailer
219 AliHLTUInt8_t* pSrc=reinterpret_cast<AliHLTUInt8_t*>(pDesc->fPtr);
220 pSrc+=pDesc->fSize-rcuTrailerLength;
221 if ((iResult=CopyBlockToEnd(outputPtr, iCapacity, iOutputSize, pSrc, rcuTrailerLength))>=0) {
222 assert(iResult==(int)rcuTrailerLength);
223 iOutputSize+=rcuTrailerLength;
224 } else {
225 HLTError("failed to write RCU trailer of length %d for block %d, too little space in output buffer?", rcuTrailerLength, blockno);
226 iResult=-ENOSPC;
227 break;
228 }
229
4ca5858d 230 while (decoder->NextChannel(&channel) && iResult>=0) {
1d337632 231 iTotal++;
64defa03 232
4ca5858d 233 int hwAddress=channel.GetHadd();
e1440dab 234 int active=0;
4ca5858d 235 for (active=0; active<iArraySize; active++) {
236 if (pActiveHwAddressArray[active]==(AliHLTUInt16_t)hwAddress) {
237 break;
e1440dab 238 }
239 }
64defa03 240 if (active>=iArraySize) {
4ca5858d 241 HLTDebug("ALTRO block %#x (%d) discarded (inactive)", hwAddress, hwAddress);
e1440dab 242 continue;
243 }
244
4ca5858d 245 // no of 10 bit words is without the fill words to fill complete 40 bit words
246 // in addition, align to complete 40 bit words (the '+3')
247 // also, the 5 bytes of the Altro trailer must be added to get the full size
248 int channelSize=((channel.GetDataSize()+3)/4)*5;
8d682af4 249 if (channelSize==0) {
250 if (fTalkative) HLTWarning("skipping zero length channel (hw address %d)", hwAddress);
251 iCorrupted++;
252 continue;
253 }
4ca5858d 254 channelSize+=5;
255 HLTDebug("ALTRO block hwAddress 0x%08x (%d) selected (active), size %d", hwAddress, hwAddress, channelSize);
4ca5858d 256
257 if ((iResult=decoder->CopyBackward(outputPtr, iCapacity-iOutputSize))>=0) {
258 if (channelSize == iResult) {
259 if (channelSize%5 == 0) {
260 iNofAltro40+=channelSize/5;
261 iOutputSize+=channelSize;
e1440dab 262 } else {
8d682af4 263 if (fTalkative) HLTWarning("corrupted ALTRO channel: incomplete 40 bit word (channel hw address %d)", hwAddress);
4ca5858d 264 iCorrupted++;
265 continue;
e1440dab 266 }
4ca5858d 267 } else {
8d682af4 268 if (fTalkative) HLTWarning("internal error: failed to copy full channel: %d out of %d bytes (hw address %d)", iResult, channelSize, hwAddress);
4ca5858d 269 iCorrupted++;
270 continue;
e1440dab 271 }
e1440dab 272 } else {
8d682af4 273 if (fTalkative) HLTError("failed to write ALTRO channel of length %d for block %d (hw address %d)", channelSize, blockno, hwAddress);
4ca5858d 274 // corrupted channel, but keep going
275 iCorrupted++;
276 iResult=0;
277 continue;
e1440dab 278 }
4ca5858d 279 iSelected++;
e1440dab 280 }
281 if (iResult>=0) {
282 // write the common data header
1d337632 283 if ((iResult=CopyBlockToEnd(outputPtr, iCapacity, iOutputSize, pDesc->fPtr, cdhSize))>=0) {
e1440dab 284 assert(iResult==cdhSize);
285 iOutputSize+=cdhSize;
286
287 // set new length of the data block
4ca5858d 288 AliHLTUInt32_t* pCdhSize=reinterpret_cast<AliHLTUInt32_t*>(outputPtr+iCapacity-iOutputSize);
e1440dab 289 *pCdhSize=iOutputSize;
290
4ca5858d 291 // set number of Altro words
292 AliHLTUInt32_t* pNofAltro40=reinterpret_cast<AliHLTUInt32_t*>(outputPtr+iCapacity-rcuTrailerLength);
293 *pNofAltro40=iNofAltro40;
294
e1440dab 295 // insert block descriptor
296 AliHLTComponentBlockData bd;
297 FillBlockData(bd);
298 bd.fOffset=iCapacity-iOutputSize;
299 bd.fSize=iOutputSize;
1d337632 300 bd.fDataType=pDesc->fDataType;
301 bd.fSpecification=pDesc->fSpecification;
e1440dab 302 outputBlocks.push_back(bd);
303 iCapacity-=iOutputSize;
304 } else {
305 HLTError("failed to write CDH of length %d for block %d", cdhSize, blockno);
306 break;
307 }
308 }
4ca5858d 309 HLTInfo("data block %d (0x%08x): selected %d out of %d ALTRO channel(s), %d corrupted channels skipped", blockno, pDesc->fSpecification, iSelected, iTotal, iCorrupted);
3e6ec852 310 }
4ca5858d 311 if (decoder) delete decoder;
3e6ec852 312
e1440dab 313 if (iResult<0) {
314 outputBlocks.clear();
3e6ec852 315 }
316
e1440dab 317 // !!! do not change the size since the output buffer is filled from the end !!!
318
319 return iResult;
320}
321
322int AliHLTAltroChannelSelectorComponent::CopyBlockToEnd(AliHLTUInt8_t* pTgt, unsigned capacity, unsigned position, void* pSrc, unsigned size)
323{
324 int iResult=0;
325 if (pTgt==NULL || pSrc==NULL) return -EINVAL;
326 if (capacity-position<size) return -ENOSPC;
327
328 memcpy(pTgt+(capacity-position-size), pSrc, size);
329 iResult=size;
330
3e6ec852 331 return iResult;
332}