3 //**************************************************************************
4 //* This file is property of and copyright by the ALICE HLT Project *
5 //* ALICE Experiment at CERN, All rights reserved. *
7 //* Primary Authors: Matthias Richter <Matthias.Richter@ift.uib.no> *
8 //* for The ALICE HLT Project. *
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 //**************************************************************************
19 // @file AliHLTCTPData.cxx
20 // @author Matthias Richter
22 // @brief Container for CTP trigger classes and counters
25 #include "AliHLTCTPData.h"
27 #include "TObjString.h"
29 #include "AliHLTComponent.h"
30 #include "AliRawDataHeader.h"
32 /** ROOT macro for the implementation of ROOT specific class methods */
33 ClassImp(AliHLTCTPData)
35 AliHLTCTPData::AliHLTCTPData()
36 : TNamed("AliHLTCTPData", "HLT counters for the CTP")
40 , fClassIds(AliHLTReadoutList::Class(), gkNCTPTriggerClasses)
41 , fCounters(gkNCTPTriggerClasses)
45 // see header file for class documentation
47 // refer to README to build package
49 // visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
52 AliHLTCTPData::AliHLTCTPData(const char* parameter)
53 : TNamed("AliHLTCTPData", "HLT counters for the CTP")
57 , fClassIds(AliHLTReadoutList::Class(), gkNCTPTriggerClasses)
58 , fCounters(gkNCTPTriggerClasses)
61 // constructor, init the CTP trigger classes
62 InitCTPTriggerClasses(parameter);
65 AliHLTCTPData::~AliHLTCTPData()
71 AliHLTCTPData::AliHLTCTPData(const AliHLTCTPData& src)
72 : TNamed(src.GetName(), src.GetTitle())
75 , fTriggers(src.fTriggers)
76 , fClassIds(src.fClassIds)
77 , fCounters(src.Counters())
84 AliHLTCTPData& AliHLTCTPData::operator=(const AliHLTCTPData& src)
86 // assignment operator, clone content
88 SetName(src.GetName());
89 SetTitle(src.GetTitle());
92 fClassIds.ExpandCreate(gkNCTPTriggerClasses);
93 for (int i=0; i<gkNCTPTriggerClasses; i++) {
94 if (i>src.fClassIds.GetLast()) break;
95 ((TNamed*)fClassIds.At(i))->SetName(src.fClassIds.At(i)->GetName());
96 ((TNamed*)fClassIds.At(i))->SetTitle(src.fClassIds.At(i)->GetTitle());
98 fCounters=src.Counters();
105 int AliHLTCTPData::Add(const AliHLTCTPData& src, int factor, int &skipped)
107 // see header file for class documentation
110 for (int i=0; i<gkNCTPTriggerClasses; i++) {
112 c=fClassIds.At(i)->GetName();
113 if (c.IsNull()) continue;
114 if (c.CompareTo(src.fClassIds.At(i)->GetName())==0) {
115 fCounters[i]+=factor*src.Counter(i);
123 AliHLTCTPData& AliHLTCTPData::operator += (const AliHLTCTPData& src)
125 // see header file for class documentation
127 int nofInconsistencies=0;
128 Add(src, 1, nofInconsistencies);
129 if (nofInconsistencies>0) {
130 HLTError("Inconsistent operants: skipping %d of %d CTP classes for operation", nofInconsistencies, gkNCTPTriggerClasses);
135 AliHLTCTPData& AliHLTCTPData::operator -= (const AliHLTCTPData& src)
137 // see header file for class documentation
139 int nofInconsistencies=0;
140 Add(src, -1, nofInconsistencies);
141 if (nofInconsistencies>0) {
142 HLTError("Inconsistent operants: skipping %d of %d CTP classes for operation", nofInconsistencies, gkNCTPTriggerClasses);
147 AliHLTCTPData AliHLTCTPData::operator + (const AliHLTCTPData& src) const
149 // see header file for class documentation
151 AliHLTCTPData result(*this);
156 AliHLTCTPData AliHLTCTPData::operator - (const AliHLTCTPData& src) const
158 // see header file for class documentation
160 AliHLTCTPData result(*this);
165 int AliHLTCTPData::InitCTPTriggerClasses(const char* ctpString)
167 // see header file for function documentation
168 if (!ctpString) return -EINVAL;
170 HLTImportant("Parameter: %s", ctpString);
174 fClassIds.ExpandCreate(gkNCTPTriggerClasses);
176 // general format of the CTP_TRIGGER_CLASS parameter
177 // <bit position>:<Trigger class identifier string>:<detector-id-nr>-<detector-id-nr>-...,<bit position>:<Trigger class identifier string>:<detector-id-nr>-<detector-id-nr>-...,...
178 HLTDebug(": %s", ctpString);
179 TString string=ctpString;
180 if (string.BeginsWith("CTP_TRIGGER_CLASS=")) string.ReplaceAll("CTP_TRIGGER_CLASS=", "");
181 TObjArray* classEntries=string.Tokenize(",");
183 enum {kBit=0, kName, kDetectors};
184 for (int i=0; i<classEntries->GetEntriesFast(); i++) {
185 TString entry=((TObjString*)classEntries->At(i))->GetString();
186 TObjArray* entryParams=entry.Tokenize(":");
188 if (entryParams->GetEntriesFast()==3 &&
189 (((TObjString*)entryParams->At(kBit))->GetString()).IsDigit()) {
190 int index=(((TObjString*)entryParams->At(kBit))->GetString()).Atoi();
191 if (index<gkNCTPTriggerClasses) {
192 AliHLTReadoutList* pCTPClass=dynamic_cast<AliHLTReadoutList*>(fClassIds.At(index));
194 fMask|=(AliHLTUInt64_t)0x1 << index;
195 pCTPClass->SetTitle("CTP Class");
196 pCTPClass->SetName((((TObjString*)entryParams->At(kName))->GetString()).Data());
197 TObjArray* detectors=(((TObjString*)entryParams->At(kDetectors))->GetString()).Tokenize("-");
199 for (int dix=0; dix<detectors->GetEntriesFast(); dix++) {
200 if (!(((TObjString*)detectors->At(dix))->GetString()).IsDigit()) {
201 HLTError("invalid detector list format: trigger class entry %s", entry.Data());
204 // see AliHLTReadoutList::EDetectorId for defines of detectors
205 pCTPClass->Enable(0x1<<(((TObjString*)detectors->At(dix))->GetString()).Atoi());
212 // the trigger bitfield is fixed to 50 bits (gkNCTPTriggerClasses)
213 HLTError("invalid trigger class entry %s, index width of trigger bitfield exceeded (%d)", entry.Data(), gkNCTPTriggerClasses);
216 HLTError("invalid trigger class entry %s", entry.Data());
230 AliHLTUInt64_t AliHLTCTPData::ActiveTriggers(const AliHLTComponentTriggerData& trigData)
232 // extract active triggers from the trigger data
234 const AliRawDataHeader* cdh = NULL;
235 if (AliHLTComponent::ExtractTriggerData(trigData, NULL, NULL, &cdh, NULL) != 0) return (AliHLTUInt64_t)0;
236 // trigger mask is 50 bit wide and is stored in word 5 and 6 of the CDH
237 AliHLTUInt64_t triggerMask = cdh->GetTriggerClasses();
241 bool AliHLTCTPData::EvaluateCTPTriggerClass(const char* expression, const AliHLTComponentTriggerData& trigData) const
243 // see header file for function documentation
245 const AliRawDataHeader* cdh = NULL;
246 if (AliHLTComponent::ExtractTriggerData(trigData, NULL, NULL, &cdh, NULL, true) != 0) return false;
247 // trigger mask is 50 bit wide and is stored in word 5 and 6 of the CDH
248 AliHLTUInt64_t triggerMask = cdh->GetTriggerClasses();
250 if (fMask!=0 && (triggerMask & fMask)==0) {
251 AliHLTEventTriggerData* evtData=reinterpret_cast<AliHLTEventTriggerData*>(trigData.fData);
252 HLTWarning("invalid trigger mask 0x%llx, unknown CTP trigger, initialized 0x%llx", triggerMask, fMask);
253 for (int i=0; i<8; i++) HLTWarning("\t CDH[%d]=0x%lx", i, evtData->fCommonHeader[i]);
257 return EvaluateCTPTriggerClass(expression, triggerMask);
260 bool AliHLTCTPData::EvaluateCTPTriggerClass(const char* expression, AliHLTUInt64_t triggerMask) const
262 // see header file for function documentation
264 // use a TFormula to interprete the expression
265 // all classname are replaced by '[n]' which means the n'th parameter in the formula
266 // the parameters are set to 0 or 1 depending on the bit in the trigger mask
267 const vector<unsigned> *pMap=&fMap;
268 vector<unsigned> tmp;
269 if (fMap.size()==0 && fClassIds.GetLast()>=0) {
270 // read map into temporary array and use it
273 static bool suppressWarning=false;
274 if (!suppressWarning) HLTWarning("map not yet initialized, creating local map (slow), suppressing further warnings");
275 suppressWarning=true;
277 vector<Double_t> par;
278 TString condition=expression;
279 for (unsigned index=0; index<pMap->size(); index++) {
280 const char* className=Name((*pMap)[index]);
281 if (className && strlen(className)>0) {
282 //HLTDebug("checking trigger class %s", className.Data());
283 if (condition.Contains(className)) {
284 TString replace; replace.Form("[%d]", par.size());
285 //HLTDebug("replacing %s with %s in \"%s\"", className.Data(), replace.Data(), condition.Data());
286 condition.ReplaceAll(className, replace);
287 if (triggerMask&((AliHLTUInt64_t)0x1<<(*pMap)[index])) par.push_back(1.0);
288 else par.push_back(0.0);
293 TFormula form("trigger expression", condition);
294 if (form.Compile()!=0) {
295 HLTError("invalid expression %s", expression);
298 if (form.EvalPar(&par[0], &par[0])>0.5) return true;
302 void AliHLTCTPData::ResetCounters()
304 // see header file for function documentation
305 fCounters.Set(gkNCTPTriggerClasses);
309 int AliHLTCTPData::Index(const char* name) const
311 // see header file for function documentation
312 TObject* obj=fClassIds.FindObject(name);
313 return obj!=NULL?fClassIds.IndexOf(obj):-1;
316 int AliHLTCTPData::CheckTrigger(const char* name) const
318 // check status of a trigger class
319 int index=Index(name);
320 if (index<0) return index;
321 return (fTriggers&(0x1<<index))>0?1:0;
324 void AliHLTCTPData::Increment(const char* classIds)
326 // see header file for function documentation
327 TString string=classIds;
328 TObjArray* classEntries=string.Tokenize(",");
330 for (int i=0; i<classEntries->GetEntriesFast(); i++) {
331 int index=Index(((TObjString*)classEntries->At(i))->GetString().Data());
332 if (index>=0 && index<fCounters.GetSize()) fCounters[index]++;
338 void AliHLTCTPData::Increment(AliHLTUInt64_t triggerPattern)
340 // see header file for function documentation
341 AliHLTUInt64_t pattern=triggerPattern&fMask;
342 for (int i=0; i<fCounters.GetSize(); i++) {
343 if ((pattern&((AliHLTUInt64_t)0x1<<i))==0) continue;
348 void AliHLTCTPData::Increment(int classIdx)
350 // see header file for function documentation
351 if (classIdx<fCounters.GetSize() &&
352 (fMask&((AliHLTUInt64_t)0x1<<classIdx))) {
353 fCounters[classIdx]++;
358 int AliHLTCTPData::Increment(AliHLTComponentTriggerData& trigData)
360 // see header file for function documentation
361 const AliRawDataHeader* cdh = NULL;
362 int result = AliHLTComponent::ExtractTriggerData(trigData, NULL, NULL, &cdh, NULL, true);
363 if (result != 0) return result;
364 // trigger mask is 50 bit wide and is stored in word 5 and 6 of the CDH
365 AliHLTUInt64_t triggerMask = cdh->GetTriggerClasses();
367 if (fMask!=0 && (triggerMask & fMask)==0) {
368 AliHLTEventTriggerData* evtData=reinterpret_cast<AliHLTEventTriggerData*>(trigData.fData);
369 HLTWarning("invalid trigger mask 0x%llx, unknown CTP trigger, initialized 0x%llx", triggerMask, fMask);
370 for (int i=0; i<8; i++) HLTWarning("\t CDH[%d]=0x%lx", i, evtData->fCommonHeader[i]);
372 Increment(triggerMask);
376 AliHLTUInt64_t AliHLTCTPData::Counter(int index) const
378 // see header file for function documentation
379 if (index>=0 && index<Counters().GetSize()) return Counters()[index];
383 AliHLTUInt64_t AliHLTCTPData::Counter(const char* classId) const
385 // see header file for function documentation
386 return Counter(Index(classId));
389 const char* AliHLTCTPData::Name(int index) const
391 // see header file for function documentation
392 if (index>fClassIds.GetLast()) return NULL;
393 return fClassIds.At(index)->GetName();
396 int AliHLTCTPData::ReadMap(vector<unsigned> &map) const
398 // read the index map for class names
399 // for nested class names (e.g. 'myclass' is contained in
400 // 'myclassA') the longer names is added first to the map.
401 for (int index=0; index<=fClassIds.GetLast(); index++) {
402 vector<unsigned>::iterator element=map.begin();
403 for (; element!=map.end(); element++) {
404 TString name=Name(index);
405 if (name.Contains(Name(*element))) {
406 // current name contains another one already in the map
407 // -> add before and go to next entry
408 element=map.insert(element, index);
413 if (element==map.end()) {
414 // unique class name, append to map
415 map.push_back(index);
422 AliHLTReadoutList AliHLTCTPData::ReadoutList(const AliHLTComponentTriggerData& trigData) const
424 // see header file for function documentation
426 const AliRawDataHeader* cdh = NULL;
427 if (AliHLTComponent::ExtractTriggerData(trigData, NULL, NULL, &cdh, NULL, true) != 0) return AliHLTReadoutList();
428 // trigger mask is 50 bit wide and is stored in word 5 and 6 of the CDH
429 AliHLTUInt64_t triggerMask = cdh->GetTriggerClasses();
431 if (fMask!=0 && (triggerMask & fMask)==0) {
432 AliHLTEventTriggerData* evtData=reinterpret_cast<AliHLTEventTriggerData*>(trigData.fData);
433 HLTWarning("invalid trigger mask 0x%llx, unknown CTP trigger, initialized 0x%llx", triggerMask, fMask);
434 for (int i=0; i<8; i++) HLTWarning("\t CDH[%d]=0x%lx", i, evtData->fCommonHeader[i]);
437 return ReadoutList(triggerMask);
440 AliHLTReadoutList AliHLTCTPData::ReadoutList(AliHLTUInt64_t triggerMask) const
442 // take an 'OR' of all active trigger classes
443 AliHLTReadoutList list;
444 for (int i=0; i<gkNCTPTriggerClasses; i++) {
445 if (i>fClassIds.GetLast()) break;
446 if ((triggerMask&((AliHLTUInt64_t)0x1<<i))==0) continue;
447 AliHLTReadoutList* tcrl=(AliHLTReadoutList*)fClassIds.At(i);
455 void AliHLTCTPData::Print(Option_t* /*option*/) const
457 // see header file for function documentation
458 cout << GetTitle() << endl;
459 cout << "\tactive trigger mask: 0x" << hex << fTriggers << dec << endl;
461 for (int i=0; i<gkNCTPTriggerClasses; i++) {
462 if (i>=Counters().GetSize()) break;
463 if (i>fClassIds.GetLast()) break;
464 if ((fMask&((AliHLTUInt64_t)0x1<<i))==0) continue;
466 cout << "\t" << i << "\t" << Name(i) << "\t" << Counter(i) << endl;
468 if (count==0) cout << "\t(none)" << endl;