Memory leak fixed + some comments updated (Alexei)
[u/mrichter/AliRoot.git] / HLT / BASE / AliHLTCTPData.cxx
CommitLineData
1dbbd625 1// $Id$
2
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> *
1dbbd625 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//**************************************************************************
18
19/** @file AliHLTCTPData.cxx
20 @author Matthias Richter
21 @date 2009-08-20
22 @brief Container for CTP trigger classes and counters
23*/
24
25#include "AliHLTCTPData.h"
1462df14 26#include "AliHLTReadoutList.h"
1dbbd625 27#include "TClass.h"
1dbbd625 28#include "TObjString.h"
29#include "TFormula.h"
30
31/** ROOT macro for the implementation of ROOT specific class methods */
32ClassImp(AliHLTCTPData)
33
34AliHLTCTPData::AliHLTCTPData()
1462df14 35 : TNamed("AliHLTCTPData", "HLT counters for the CTP")
1dbbd625 36 , AliHLTLogging()
37 , fMask(0)
1462df14 38 , fClassIds(AliHLTReadoutList::Class(), gkNCTPTriggerClasses)
59be6254 39 , fCounters(gkNCTPTriggerClasses)
1dbbd625 40{
41 // see header file for class documentation
42 // or
43 // refer to README to build package
44 // or
45 // visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
46}
47
48AliHLTCTPData::AliHLTCTPData(const char* parameter)
49 : TNamed("AliHLTCTPData", "HLT counters for the CTP")
50 , AliHLTLogging()
51 , fMask(0)
1462df14 52 , fClassIds(AliHLTReadoutList::Class(), gkNCTPTriggerClasses)
59be6254 53 , fCounters(gkNCTPTriggerClasses)
1dbbd625 54{
55 // see header file for class documentation
56 InitCTPTriggerClasses(parameter);
57}
58
59AliHLTCTPData::~AliHLTCTPData()
60{
61 // see header file for class documentation
62 fClassIds.Delete();
63}
64
59be6254 65AliHLTCTPData::AliHLTCTPData(const AliHLTCTPData& src)
66 : TNamed(src.GetName(), src.GetTitle())
67 , AliHLTLogging()
68 , fMask(src.Mask())
69 , fClassIds(src.fClassIds)
70 , fCounters(src.Counters())
71{
72 // see header file for class documentation
73}
74
75AliHLTCTPData& AliHLTCTPData::operator=(const AliHLTCTPData& src)
76{
77 // see header file for class documentation
78 if (this!=&src) {
79 SetName(src.GetName());
80 SetTitle(src.GetTitle());
81 fMask=src.Mask();
82 fClassIds.Delete();
83 fClassIds.ExpandCreate(gkNCTPTriggerClasses);
84 for (int i=0; i<gkNCTPTriggerClasses; i++) {
85 ((TNamed*)fClassIds.At(i))->SetName(src.fClassIds.At(i)->GetName());
86 ((TNamed*)fClassIds.At(i))->SetTitle(src.fClassIds.At(i)->GetTitle());
87 }
88 fCounters=src.Counters();
89 }
90
91 return *this;
92}
93
94int AliHLTCTPData::Add(const AliHLTCTPData& src, int factor, int &skipped)
95{
96 // see header file for class documentation
97
98 skipped=0;
99 for (int i=0; i<gkNCTPTriggerClasses; i++) {
100 TString c;
101 c=fClassIds.At(i)->GetName();
102 if (c.IsNull()) continue;
103 if (c.CompareTo(src.fClassIds.At(i)->GetName())==0) {
104 fCounters[i]+=factor*src.Counter(i);
105 } else {
106 skipped++;
107 }
108 }
109 return 0;
110}
111
112AliHLTCTPData& AliHLTCTPData::operator += (const AliHLTCTPData& src)
113{
114 // see header file for class documentation
115
116 int nofInconsistencies=0;
117 Add(src, 1, nofInconsistencies);
118 if (nofInconsistencies>0) {
119 HLTError("Inconsistent operants: skipping %d of %d CTP classes for operation", nofInconsistencies, gkNCTPTriggerClasses);
120 }
121 return *this;
122}
123
124AliHLTCTPData& AliHLTCTPData::operator -= (const AliHLTCTPData& src)
125{
126 // see header file for class documentation
127
128 int nofInconsistencies=0;
129 Add(src, -1, nofInconsistencies);
130 if (nofInconsistencies>0) {
131 HLTError("Inconsistent operants: skipping %d of %d CTP classes for operation", nofInconsistencies, gkNCTPTriggerClasses);
132 }
133 return *this;
134}
135
136AliHLTCTPData AliHLTCTPData::operator + (const AliHLTCTPData& src) const
137{
138 // see header file for class documentation
139
140 AliHLTCTPData result(*this);
141 result+=src;
142 return result;
143}
144
145AliHLTCTPData AliHLTCTPData::operator - (const AliHLTCTPData& src) const
146{
147 // see header file for class documentation
148
149 AliHLTCTPData result(*this);
150 result-=src;
151 return result;
152}
153
1dbbd625 154int AliHLTCTPData::InitCTPTriggerClasses(const char* ctpString)
155{
156 // see header file for function documentation
157 if (!ctpString) return -EINVAL;
158
5be104b8 159 HLTImportant("Parameter: %s", ctpString);
160
1dbbd625 161 fMask=0;
162 fClassIds.Delete();
163 fClassIds.ExpandCreate(gkNCTPTriggerClasses);
164
165 // general format of the CTP_TRIGGER_CLASS parameter
166 // <bit position>:<Trigger class identifier string>:<detector-id-nr>-<detector-id-nr>-...,<bit position>:<Trigger class identifier string>:<detector-id-nr>-<detector-id-nr>-...,...
1dbbd625 167 HLTDebug(": %s", ctpString);
168 TString string=ctpString;
1462df14 169 if (string.BeginsWith("CTP_TRIGGER_CLASS=")) string.ReplaceAll("CTP_TRIGGER_CLASS=", "");
1dbbd625 170 TObjArray* classEntries=string.Tokenize(",");
171 if (classEntries) {
1462df14 172 enum {kBit=0, kName, kDetectors};
1dbbd625 173 for (int i=0; i<classEntries->GetEntries(); i++) {
174 TString entry=((TObjString*)classEntries->At(i))->GetString();
175 TObjArray* entryParams=entry.Tokenize(":");
176 if (entryParams) {
177 if (entryParams->GetEntries()==3 &&
1462df14 178 (((TObjString*)entryParams->At(kBit))->GetString()).IsDigit()) {
179 int index=(((TObjString*)entryParams->At(kBit))->GetString()).Atoi();
1dbbd625 180 if (index<gkNCTPTriggerClasses) {
1462df14 181 AliHLTReadoutList* pCTPClass=dynamic_cast<AliHLTReadoutList*>(fClassIds.At(index));
182 if (pCTPClass) {
183 fMask|=(AliHLTUInt64_t)0x1 << index;
184 pCTPClass->SetTitle("CTP Class");
185 pCTPClass->SetName((((TObjString*)entryParams->At(kName))->GetString()).Data());
186 TObjArray* detectors=(((TObjString*)entryParams->At(kDetectors))->GetString()).Tokenize("-");
187 if (detectors) {
188 for (int dix=0; dix<detectors->GetEntriesFast(); dix++) {
189 if (!(((TObjString*)detectors->At(dix))->GetString()).IsDigit()) {
190 HLTError("invalid detector list format: trigger class entry %s", entry.Data());
191 break;
192 }
193 // see AliHLTReadoutList::EDetectorId for defines of detectors
194 pCTPClass->Enable(0x1<<(((TObjString*)detectors->At(dix))->GetString()).Atoi());
195 }
196 delete detectors;
197 }
198 } else {
199 }
1dbbd625 200 } else {
201 // the trigger bitfield is fixed to 50 bits (gkNCTPTriggerClasses)
1462df14 202 HLTError("invalid trigger class entry %s, index width of trigger bitfield exceeded (%d)", entry.Data(), gkNCTPTriggerClasses);
1dbbd625 203 }
204 } else {
205 HLTError("invalid trigger class entry %s", entry.Data());
206 }
207 delete entryParams;
208 }
209 }
210 delete classEntries;
211 }
212
213 ResetCounters();
214
215 return 0;
216}
217
218bool AliHLTCTPData::EvaluateCTPTriggerClass(const char* expression, AliHLTComponentTriggerData& trigData) const
219{
220 // see header file for function documentation
221 if (trigData.fDataSize != sizeof(AliHLTEventTriggerData)) {
222 HLTError("invalid trigger data size: %d expected %d", trigData.fDataSize, sizeof(AliHLTEventTriggerData));
223 return false;
224 }
225
226 // trigger mask is 50 bit wide and is stored in word 5 and 6 of the CDH
227 AliHLTEventTriggerData* evtData=reinterpret_cast<AliHLTEventTriggerData*>(trigData.fData);
5be104b8 228 AliHLTUInt64_t triggerMask=evtData->fCommonHeader[6]&0x3ffff;
1dbbd625 229 triggerMask<<=32;
230 triggerMask|=evtData->fCommonHeader[5];
231
232 if (fMask!=0 && (triggerMask & fMask)==0) {
5be104b8 233 HLTWarning("invalid trigger mask 0x%llx, unknown CTP trigger, initialized 0x%llx", triggerMask, fMask);
234 for (int i=0; i<8; i++) HLTWarning("\t CDH[%d]=0x%lx", i, evtData->fCommonHeader[i]);
1dbbd625 235 return false;
236 }
237
238 // use a TFormula to interprete the expression
239 // all classname are replaced by '[n]' which means the n'th parameter in the formula
240 // the parameters are set to 0 or 1 depending on the bit in the trigger mask
241 //
242 // TODO: this will most likely fail for class names like 'base', 'baseA', 'baseB'
243 // the class names must be fully unique, none must be contained as substring in
244 // another class name. Probably not needed for the moment but needs to be extended.
245 vector<Double_t> par;
246 TString condition=expression;
247 for (int i=0; i<gkNCTPTriggerClasses; i++) {
248 const char* className=Name(i);
249 if (className && strlen(className)>0) {
250 //HLTDebug("checking trigger class %s", className.Data());
251 if (condition.Contains(className)) {
252 TString replace; replace.Form("[%d]", par.size());
253 //HLTDebug("replacing %s with %s in \"%s\"", className.Data(), replace.Data(), condition.Data());
254 condition.ReplaceAll(className, replace);
255 if (triggerMask&((AliHLTUInt64_t)0x1<<i)) par.push_back(1.0);
256 else par.push_back(0.0);
257 }
258 }
259 }
260
261 TFormula form("trigger expression", condition);
262 if (form.Compile()!=0) {
263 HLTError("invalid expression %s", expression);
264 return false;
265 }
266 if (form.EvalPar(&par[0], &par[0])>0.5) return true;
267 return false;
268}
269
270void AliHLTCTPData::ResetCounters()
271{
272 // see header file for function documentation
273 fCounters.Set(gkNCTPTriggerClasses);
274 fCounters.Reset();
275}
276
277int AliHLTCTPData::Index(const char* name) const
278{
279 // see header file for function documentation
280 TObject* obj=fClassIds.FindObject(name);
281 return obj!=NULL?fClassIds.IndexOf(obj):-1;
282}
283
284void AliHLTCTPData::Increment(const char* classIds)
285{
286 // see header file for function documentation
287 TString string=classIds;
288 TObjArray* classEntries=string.Tokenize(",");
289 if (classEntries) {
290 for (int i=0; i<classEntries->GetEntries(); i++) {
291 int index=Index(((TObjString*)classEntries->At(i))->GetString().Data());
292 if (index>=0 && index<fCounters.GetSize()) fCounters[index]++;
293 }
294 delete classEntries;
295 }
296}
297
298void AliHLTCTPData::Increment(AliHLTUInt64_t triggerPattern)
299{
300 // see header file for function documentation
301 AliHLTUInt64_t pattern=triggerPattern&fMask;
302 for (int i=0; i<fCounters.GetSize(); i++) {
303 if ((pattern&((AliHLTUInt64_t)0x1<<i))==0) continue;
304 fCounters[i]++;
305 }
306}
307
308void AliHLTCTPData::Increment(int classIdx)
309{
310 // see header file for function documentation
311 if (classIdx<fCounters.GetSize() &&
312 (fMask&((AliHLTUInt64_t)0x1<<classIdx))) {
313 fCounters[classIdx]++;
314 }
315
316}
317
318int AliHLTCTPData::Increment(AliHLTComponentTriggerData& trigData)
319{
320 // see header file for function documentation
321 if (trigData.fDataSize != sizeof(AliHLTEventTriggerData)) {
322 HLTError("invalid trigger data size: %d expected %d", trigData.fDataSize, sizeof(AliHLTEventTriggerData));
323 return -EBADF;
324 }
325
326 // trigger mask is 50 bit wide and is stored in word 5 and 6 of the CDH
327 AliHLTEventTriggerData* evtData=reinterpret_cast<AliHLTEventTriggerData*>(trigData.fData);
5be104b8 328 AliHLTUInt64_t triggerMask=evtData->fCommonHeader[6]&0x3ffff;
1dbbd625 329 triggerMask<<=32;
330 triggerMask|=evtData->fCommonHeader[5];
331
332 if (fMask!=0 && (triggerMask & fMask)==0) {
333 HLTWarning("invalid trigger mask 0x%llx, unknown CTP trigger, initialized 0x%llx", triggerMask, fMask);
5be104b8 334 for (int i=0; i<8; i++) HLTWarning("\t CDH[%d]=0x%lx", i, evtData->fCommonHeader[i]);
1dbbd625 335 }
336 Increment(triggerMask);
337 return 0;
338}
339
340AliHLTUInt64_t AliHLTCTPData::Counter(int index) const
341{
342 // see header file for function documentation
343 if (index>=0 && index<Counters().GetSize()) return Counters()[index];
344 return 0;
345}
346
347AliHLTUInt64_t AliHLTCTPData::Counter(const char* classId) const
348{
349 // see header file for function documentation
350 return Counter(Index(classId));
351}
352
353const char* AliHLTCTPData::Name(int index) const
354{
355 // see header file for function documentation
356 if (index>fClassIds.GetLast()) return NULL;
357 return fClassIds.At(index)->GetName();
358}
359
1462df14 360AliHLTEventDDL AliHLTCTPData::ReadoutList(const AliHLTComponentTriggerData& trigData) const
361{
362 // see header file for function documentation
363 if (trigData.fDataSize != sizeof(AliHLTEventTriggerData)) {
364 HLTError("invalid trigger data size: %d expected %d", trigData.fDataSize, sizeof(AliHLTEventTriggerData));
365 AliHLTEventDDL dummy;
366 memset(&dummy, 0, sizeof(AliHLTEventDDL));
367 return dummy;
368 }
369
370 // trigger mask is 50 bit wide and is stored in word 5 and 6 of the CDH
371 AliHLTEventTriggerData* evtData=reinterpret_cast<AliHLTEventTriggerData*>(trigData.fData);
5be104b8 372 AliHLTUInt64_t triggerMask=evtData->fCommonHeader[6]&0x3ffff;
1462df14 373 triggerMask<<=32;
374 triggerMask|=evtData->fCommonHeader[5];
375
5be104b8 376 if (fMask!=0 && (triggerMask & fMask)==0) {
377 HLTWarning("invalid trigger mask 0x%llx, unknown CTP trigger, initialized 0x%llx", triggerMask, fMask);
378 for (int i=0; i<8; i++) HLTWarning("\t CDH[%d]=0x%lx", i, evtData->fCommonHeader[i]);
379 }
380
1462df14 381 // take an 'OR' of all active trigger classes
382 AliHLTReadoutList list;
383 for (int i=0; i<gkNCTPTriggerClasses; i++) {
384 if (i>fClassIds.GetLast()) break;
385 if ((triggerMask&((AliHLTUInt64_t)0x1<<i))==0) continue;
5be104b8 386 AliHLTReadoutList* tcrl=(AliHLTReadoutList*)fClassIds.At(i);
387 // 2009-08-27: this is a temorary bugfix:
388 // the operator functions of the AliHLTReadoutList class did not work
389 // when running online on the HLT cluster. The fix for the moment is
390 // to send out the readout list only for the first matching trigger
391 // class instead of merging the list. This is sufficient for the
392 // current trigger setups but needs to be corrected
393 return *tcrl;
394 list|=*tcrl;
1462df14 395 }
396
397 return list;
398}
399
1dbbd625 400void AliHLTCTPData::Print(Option_t* /*option*/) const
401{
402 // see header file for function documentation
59be6254 403 cout << GetTitle() << endl;
1dbbd625 404 int count=0;
405 for (int i=0; i<gkNCTPTriggerClasses; i++) {
406 if (i>=Counters().GetSize()) break;
407 if (i>fClassIds.GetLast()) break;
408 if ((fMask&((AliHLTUInt64_t)0x1<<i))==0) continue;
409 count++;
410 cout << "\t" << i << "\t" << Name(i) << "\t" << Counter(i) << endl;
411 }
412 if (count==0) cout << "\t(none)" << endl;
413}