2 /**************************************************************************
3 * This file is property of and copyright by the ALICE HLT Project *
4 * ALICE Experiment at CERN, All rights reserved. *
6 * Primary Authors: Artur Szostak <artursz@iafrica.com> *
7 * for The ALICE HLT Project. *
9 * Permission to use, copy, modify and distribute this software and its *
10 * documentation strictly for non-commercial purposes is hereby granted *
11 * without fee, provided that the above copyright notice appears in all *
12 * copies and that both the copyright notice and this permission notice *
13 * appear in the supporting documentation. The authors make no claims *
14 * about the suitability of this software for any purpose. It is *
15 * provided "as is" without express or implied warranty. *
16 **************************************************************************/
18 /// @file AliHLTGlobalTriggerComponent.cxx
19 /// @author Artur Szostak <artursz@iafrica.com>
21 /// @brief Implementation of the AliHLTGlobalTriggerComponent component class.
23 /// The AliHLTGlobalTriggerComponent class applies the global HLT trigger to all
24 /// trigger information produced by components deriving from AliHLTTrigger.
26 #include "AliHLTGlobalTriggerComponent.h"
27 #include "AliHLTGlobalTriggerDecision.h"
28 #include "AliHLTGlobalTrigger.h"
29 #include "AliHLTGlobalTriggerConfig.h"
30 #include "AliHLTTriggerMenu.h"
31 #include "AliHLTCTPData.h"
32 #include "AliCDBManager.h"
33 #include "AliCDBStorage.h"
34 #include "AliCDBEntry.h"
35 #include "AliRawDataHeader.h"
40 #include "TClonesArray.h"
41 #include "TObjString.h"
43 #include "TInterpreter.h"
54 ClassImp(AliHLTGlobalTriggerComponent)
56 const char* AliHLTGlobalTriggerComponent::fgkTriggerMenuCDBPath = "HLT/ConfigHLT/HLTGlobalTrigger";
62 * This method is used as a comparison functor with the STL sort routines.
64 bool AliHLTDescendingNumbers(UInt_t a, UInt_t b)
71 AliHLTGlobalTriggerComponent::AliHLTGlobalTriggerComponent() :
75 fRuntimeCompile(true),
76 fDeleteCodeFile(false),
77 fMakeSoftwareTriggers(true),
81 fBufferSizeConst(2*(sizeof(AliHLTGlobalTriggerDecision) + sizeof(AliHLTReadoutList))),
82 fBufferSizeMultiplier(1.),
83 fIncludePaths(TObjString::Class()),
84 fIncludeFiles(TObjString::Class()),
87 fDataEventsOnly(true),
90 fSoftwareTrigger(true, "SOFTWARE"),
91 fTotalEventCounter(0),
94 // Default constructor.
96 ClearInfoForNewEvent(false);
100 AliHLTGlobalTriggerComponent::~AliHLTGlobalTriggerComponent()
102 // Default destructor.
104 if (fTrigger != NULL) delete fTrigger;
107 fCTPDecisions->Delete();
108 delete fCTPDecisions;
113 void AliHLTGlobalTriggerComponent::GetOutputDataTypes(AliHLTComponentDataTypeList& list) const
115 // Returns the kAliHLTDataTypeGlobalTrigger type as output.
116 list.push_back(kAliHLTDataTypeGlobalTrigger);
120 void AliHLTGlobalTriggerComponent::GetOutputDataSize(unsigned long& constBase, double& inputMultiplier)
122 // Returns the output data size estimate.
124 constBase = fBufferSizeConst;
125 inputMultiplier = fBufferSizeMultiplier;
129 Int_t AliHLTGlobalTriggerComponent::DoInit(int argc, const char** argv)
131 // Initialises the global trigger component.
136 fDeleteCodeFile = false;
137 fMakeSoftwareTriggers = true;
138 const char* configFileName = NULL;
139 const char* codeFileName = NULL;
140 fIncludePaths.Clear();
141 fIncludeFiles.Clear();
142 SetBit(kIncludeInput);
143 fDataEventsOnly = true;
144 UInt_t randomSeed = 0;
145 bool randomSeedSet = false;
148 for (int i = 0; i < argc; i++)
150 if (strcmp(argv[i], "-config") == 0)
152 if (configFileName != NULL)
154 HLTWarning("Trigger configuration macro was already specified."
155 " Will replace previous value given by -config."
160 HLTError("The trigger configuration macro filename was not specified." );
163 configFileName = argv[i+1];
168 if (strcmp(argv[i], "-includepath") == 0)
172 HLTError("The include path was not specified." );
177 new (fIncludePaths[fIncludePaths.GetEntriesFast()]) TObjString(argv[i+1]);
179 catch (const std::bad_alloc&)
181 HLTError("Could not allocate more memory for the fIncludePaths array.");
188 if (strcmp(argv[i], "-include") == 0)
192 HLTError("The include file name was not specified." );
197 new (fIncludeFiles[fIncludeFiles.GetEntriesFast()]) TObjString(argv[i+1]);
199 catch (const std::bad_alloc&)
201 HLTError("Could not allocate more memory for the fIncludeFiles array.");
208 if (strcmp(argv[i], "-debug") == 0)
210 if (fDebugMode == true)
212 HLTWarning("The debug flag was already specified. Ignoring this instance.");
218 if (strcmp(argv[i], "-cint") == 0)
220 fRuntimeCompile = false;
224 if (strcmp(argv[i], "-usecode") == 0)
226 if (codeFileName != NULL)
228 HLTWarning("Custom trigger code file was already specified."
229 " Will replace previous value given by -usecode."
234 HLTError("The custom trigger code filename was not specified." );
237 codeFileName = argv[i+1];
240 HLTError("The custom trigger class name was not specified." );
243 fClassName = argv[i+2];
248 if (strcmp(argv[i], "-skipctp") == 0)
250 HLTInfo("Skipping CTP counters in trigger decision");
255 if (strcmp(argv[i], "-forward-input") == 0)
257 HLTInfo("Forwarding input objects and trigger decisions");
258 SetBit(kForwardInput);
259 SetBit(kIncludeShort);
260 SetBit(kIncludeInput, false);
264 if (strstr(argv[i], "-include-input") == argv[i])
266 SetBit(kForwardInput,false);
267 TString param=argv[i];
268 param.ReplaceAll("-include-input", "");
269 if (param.CompareTo("=none")==0)
271 HLTInfo("skipping objects and trigger decisions");
272 SetBit(kIncludeShort, false);
273 SetBit(kIncludeInput, false);
275 else if (param.CompareTo("=short")==0)
277 HLTInfo("including short info on objects and trigger decisions");
278 SetBit(kIncludeShort);
279 SetBit(kIncludeInput, false);
281 else if (param.CompareTo("=both")==0)
283 HLTInfo("including input objects, trigger decisions and short info");
284 SetBit(kIncludeShort);
285 SetBit(kIncludeInput);
287 else if (param.CompareTo("=objects")==0 || param.IsNull())
289 HLTInfo("including input objects and trigger decisions");
290 SetBit(kIncludeShort, false);
291 SetBit(kIncludeInput);
295 HLTError("unknown parameter '%s' for argument '-include-input'", param.Data());
300 if (strcmp(argv[i], "-process-all-events") == 0)
302 fDataEventsOnly = false;
306 if (strstr(argv[i], "-monitoring") == argv[i])
308 TString param=argv[i];
309 param.ReplaceAll("-monitoring", "");
313 // enable monitoring trigger for all events
317 // enable monitoring trigger once every n seconds
318 param.ReplaceAll("=", "");
319 if (param.IsDigit()) {
320 fMonitorPeriod=param.Atoi();
322 HLTError("expecting number as parameter for argument '-monitoring=', got %s, skipping monitoring trigger", param.Data());
328 if (strcmp(argv[i], "-dont-make-software-triggers") == 0)
330 fMakeSoftwareTriggers = false;
333 if (strcmp(argv[i], "-randomseed") == 0)
337 HLTWarning("The random seed was already specified previously with the option -randomseed."
338 "Will replace the previous value of %d with the new one.",
344 HLTError("The number to use as the seed was not specified for the -randomseed option." );
347 TString numstr = argv[i+1];
348 if (not numstr.IsDigit())
350 HLTError("The number specified in the -randomseed option is not a valid decimal integer." );
353 randomSeed = numstr.Atoi();
354 randomSeedSet = true;
359 HLTError("Unknown option '%s'.", argv[i]);
363 const AliHLTTriggerMenu* menu = NULL;
364 if (configFileName != NULL)
367 cmd += configFileName;
368 gROOT->ProcessLine(cmd);
369 menu = AliHLTGlobalTriggerConfig::Menu();
372 // Try load the trigger menu from the CDB if it is not loaded yet with the
374 int result = -ENOENT;
377 result = LoadTriggerMenu(fgkTriggerMenuCDBPath, menu);
381 HLTError("No trigger menu configuration found or specified.");
385 if (codeFileName == NULL)
387 result = GenerateTrigger(menu, fClassName, fCodeFileName, fIncludePaths, fIncludeFiles);
388 if (result == 0) fDeleteCodeFile = true;
392 result = LoadTriggerClass(codeFileName, fIncludePaths);
393 if (result == 0) fCodeFileName = codeFileName;
395 if (result != 0) return result;
399 fTrigger = AliHLTGlobalTrigger::CreateNew(fClassName.Data());
401 catch (const std::bad_alloc&)
403 HLTError("Could not allocate memory for the AliHLTGlobalTrigger instance.");
406 if (fTrigger == NULL)
408 HLTError("Could not create a new instance of '%s'.", fClassName.Data());
412 fTrigger->FillFromMenu(*menu);
413 if (fTrigger->CallFailed()) return -EPROTO;
415 // setup the CTP accounting in AliHLTComponent
418 // Set the default values from the trigger menu.
419 SetDescription(menu->DefaultDescription());
420 SetTriggerDomain(menu->DefaultTriggerDomain());
422 // Initialise the random number generator seed value.
423 // NOTE: The GenerateTrigger method called above will set the fUniqueID value
424 // with a random value based on a GUID that should be unique across the system.
425 // This is then used as the seed to the random number generator if the -randomseed
426 // option is not used.
427 if (not randomSeedSet) randomSeed = fUniqueID;
428 gRandom->SetSeed(randomSeed);
430 fTotalEventCounter = 0;
435 Int_t AliHLTGlobalTriggerComponent::DoDeinit()
437 // Cleans up the global trigger component.
439 if (fTrigger != NULL)
446 fCTPDecisions->Delete();
447 delete fCTPDecisions;
452 Int_t result = UnloadTriggerClass(fCodeFileName);
453 if (result != 0) return result;
455 if (fDeleteCodeFile and !fCodeFileName.IsNull() && gSystem->AccessPathName(fCodeFileName)==0 && !fDebugMode) {
456 fCodeFileName.ReplaceAll(".cxx", "*");
457 TString command="rm "; command+=fCodeFileName;
458 gSystem->Exec(command);
461 fDeleteCodeFile=false;
467 AliHLTComponent* AliHLTGlobalTriggerComponent::Spawn()
469 // Creates a new object instance.
470 AliHLTComponent* comp = NULL;
473 comp = new AliHLTGlobalTriggerComponent;
475 catch (const std::bad_alloc&)
477 HLTError("Could not allocate memory for a new instance of AliHLTGlobalTriggerComponent.");
484 int AliHLTGlobalTriggerComponent::DoTrigger()
486 // This method will apply the global trigger decision.
488 if (fTrigger == NULL)
490 HLTFatal("Global trigger implementation object is NULL!");
494 AliHLTUInt32_t eventType=0;
495 if (!IsDataEvent(&eventType)) {
498 if (eventType == gkAliEventTypeEndOfRun) PrintStatistics(fTrigger, kHLTLogImportant);
499 IgnoreEvent(); // dont generate any trigger decision.
504 fCDH = NULL; // have to reset this in case ExtractTriggerData fails.
505 ExtractTriggerData(*GetTriggerData(), NULL, NULL, &fCDH, NULL);
507 // Copy the trigger counters in case we need to set them back to their original
508 // value because the PushBack method fails with ENOSPC.
509 TArrayL64 originalCounters = fTrigger->GetCounters();
510 if (fTrigger->CallFailed()) return -EPROTO;
512 fTrigger->NewEvent();
513 if (fTrigger->CallFailed()) return -EPROTO;
515 // Fill in the input data.
516 const TObject* obj = GetFirstInputObject();
519 fTrigger->Add(obj, GetDataType(), GetSpecification());
520 if (fTrigger->CallFailed()) return -EPROTO;
521 obj = GetNextInputObject();
524 // add trigger decisions for every CTP class
525 const AliHLTCTPData* pCTPData=CTPData();
527 AddCTPDecisions(fTrigger, pCTPData, GetTriggerData());
530 bool softwareTriggerIsValid = FillSoftwareTrigger();
531 if (softwareTriggerIsValid)
533 fTrigger->Add(&fSoftwareTrigger, kAliHLTDataTypeTriggerDecision, kAliHLTVoidDataSpec);
536 // Calculate the global trigger result and trigger domain, then create and push
537 // back the new global trigger decision object.
539 AliHLTTriggerDomain triggerDomain;
540 bool triggerResult = false;
541 bool matchedItems = fTrigger->CalculateTriggerDecision(triggerResult, triggerDomain, description);
542 if (fTrigger->CallFailed()) return -EPROTO;
543 AliHLTGlobalTriggerDecision decision(
545 // The following will cause the decision to be generated with default values
546 // (set in fTriggerDomain and fDescription) if the trigger result is false.
547 (matchedItems == true) ? triggerDomain : GetTriggerDomain(),
548 (matchedItems == true) ? description.Data() : GetDescription()
551 decision.SetUniqueID(fUniqueID);
552 decision.SetCounters(fTrigger->GetCounters(), fTotalEventCounter+1);
553 if (fTrigger->CallFailed()) return -EPROTO;
555 TClonesArray shortInfo(TNamed::Class(), GetNumberOfInputBlocks());
557 // Add the input objects used to make the global decision.
558 obj = GetFirstInputObject();
561 const AliHLTTriggerDecision* intrig = dynamic_cast<const AliHLTTriggerDecision*>(obj);
563 if (TestBit(kForwardInput)) Forward(obj);
565 if (TestBit(kIncludeInput))
569 decision.AddTriggerInput(*intrig);
573 // The const_cast should be safe in this case because the list of inputObjects
574 // will be added to the global decision with AddInputObjectRef, which only
575 // modifies the kCanDelete bit and nothing else.
576 // This is necessary since GetFirstInputObject only returns const objects.
577 decision.AddInputObjectRef( const_cast<TObject*>(obj) );
581 if (TestBit(kIncludeShort))
583 int entries = shortInfo.GetEntriesFast();
586 new (shortInfo[entries]) TNamed(obj->GetName(), obj->GetTitle());
588 catch (const std::bad_alloc&)
590 HLTError("Could not allocate more memory for the short list of input objects.");
595 shortInfo[entries]->SetBit(BIT(16)); // indicate that this is a trigger decision
596 shortInfo[entries]->SetBit(BIT(15), intrig->Result());
600 obj = GetNextInputObject();
602 if (TestBit(kIncludeShort)) decision.AddInputObjectRef(&shortInfo);
604 // The const_cast should be safe in this case because AddInputObjectRef just
605 // modifies the kCanDelete bit and nothing else.
606 if (!TestBit(kSkipCTP) && CTPData()) decision.AddInputObjectRef(const_cast<AliHLTCTPData*>(CTPData()));
608 if (softwareTriggerIsValid and TestBit(kIncludeInput)) decision.AddTriggerInput(fSoftwareTrigger);
610 static UInt_t lastTime=0;
612 if (time.Get()-lastTime>60)
615 PrintStatistics(fTrigger, kHLTLogImportant);
617 else if (eventType==gkAliEventTypeEndOfRun)
619 PrintStatistics(fTrigger, kHLTLogImportant);
622 // add readout filter to event done data
623 CreateEventDoneReadoutFilter(decision.TriggerDomain(), 3);
625 // add monitoring filter to event done data if enabled by setting
626 // a monitoring period >=0: -1 means off, 0 means for every event
627 // configured by argument '-monitoring[=n]'
628 if (fMonitorPeriod>=0) {
629 static UInt_t lastMonitorEvent=0;
631 AliHLTTriggerDomain monitoringFilter(decision.TriggerDomain());
632 if (decision.Result() &&
633 int(time.Get()-lastMonitorEvent)>fMonitorPeriod) {
634 lastMonitorEvent=time.Get();
635 // add monitoring event command for triggered events
636 CreateEventDoneReadoutFilter(decision.TriggerDomain(), 5);
638 // empty filter list if events are not triggered
639 // or within the monitoring interval
640 monitoringFilter.Clear();
642 // add monitoring filter list
643 CreateEventDoneReadoutFilter(monitoringFilter, 4);
646 // Mask the readout list according to the CTP trigger
647 // if the classes have been initialized (mask non-zero).
648 // If we are dealing with a software trigger on the other hand then
649 // mask with the participating detector list.
650 // In both cases we must make sure that HLT is part of the readout mask.
651 if (CTPData() != NULL and CTPData()->Mask() != 0x0)
653 AliHLTReadoutList readoutlist = decision.ReadoutList();
654 AliHLTReadoutList ctpreadout = CTPData()->ReadoutList(*GetTriggerData());
655 ctpreadout.Enable(AliHLTReadoutList::kHLT);
656 readoutlist.AndEq(ctpreadout);
657 decision.ReadoutList(readoutlist); // override the readout list with the masked one.
659 else if (softwareTriggerIsValid)
661 assert(fCDH != NULL);
662 AliHLTReadoutList readoutlist = decision.ReadoutList();
663 UInt_t detectors = fCDH->GetSubDetectors();
664 AliHLTReadoutList softwareReadout(Int_t(detectors | AliHLTReadoutList::kHLT));
665 readoutlist.AndEq(softwareReadout);
666 decision.ReadoutList(readoutlist); // override the readout list with the masked one.
669 if (TriggerEvent(&decision, kAliHLTDataTypeGlobalTrigger) == -ENOSPC)
671 // Increase the estimated buffer space required if the PushBack methods in TriggerEvent
672 // returned the "no buffer space" error code. Also remember to set the trigger counters
673 // back to what they were, otherwise triggers will be double counted when we try to reprocess
674 // this event with larger buffers.
675 fBufferSizeConst += 1024*1024;
676 fBufferSizeMultiplier *= 2.;
677 fTrigger->SetCounters(originalCounters);
678 if (fTrigger->CallFailed()) return -EPROTO;
682 ++fTotalEventCounter;
687 int AliHLTGlobalTriggerComponent::Reconfigure(const char* cdbEntry, const char* chainId)
689 // Reconfigure the component by loading the trigger menu and recreating the
690 // trigger logic class.
692 const char* path = fgkTriggerMenuCDBPath;
693 const char* id = "(unknown)";
694 if (cdbEntry != NULL) path = cdbEntry;
695 if (chainId != NULL and chainId[0] != '\0') id = chainId;
696 HLTInfo("Reconfiguring from '%s' for chain component '%s'.", path, id);
698 const AliHLTTriggerMenu* menu = NULL;
699 int result = LoadTriggerMenu(path, menu);
700 if (result != 0) return result;
703 TString codeFileName;
704 result = GenerateTrigger(menu, className, codeFileName, fIncludePaths, fIncludeFiles);
705 if (result != 0) return result;
707 AliHLTGlobalTrigger* trigger = NULL;
710 trigger = AliHLTGlobalTrigger::CreateNew(className.Data());
712 catch (const std::bad_alloc&)
714 HLTError("Could not allocate memory for the AliHLTGlobalTrigger instance.");
719 HLTError("Could not create a new instance of '%s'.", className.Data());
720 // Make sure to cleanup after the new code file.
721 UnloadTriggerClass(codeFileName);
722 if (not codeFileName.IsNull() and gSystem->AccessPathName(codeFileName)==0 and not fDebugMode)
724 codeFileName.ReplaceAll(".cxx", "*");
725 TString command="rm "; command+=codeFileName;
726 gSystem->Exec(command);
731 if (fTrigger != NULL)
738 fTrigger->FillFromMenu(*menu);
739 if (fTrigger->CallFailed()) return -EPROTO;
741 // Set the default values from the trigger menu.
742 SetDescription(menu->DefaultDescription());
743 SetTriggerDomain(menu->DefaultTriggerDomain());
745 // Cleanup the old class code.
746 UnloadTriggerClass(fCodeFileName);
747 if (fDeleteCodeFile and not fCodeFileName.IsNull() and gSystem->AccessPathName(fCodeFileName)==0 and not fDebugMode)
749 fCodeFileName.ReplaceAll(".cxx", "*");
750 TString command="rm "; command+=fCodeFileName;
751 gSystem->Exec(command);
753 fCodeFileName = codeFileName;
754 fDeleteCodeFile = true; // Since we generated a new class.
760 int AliHLTGlobalTriggerComponent::LoadTriggerMenu(const char* cdbPath, const AliHLTTriggerMenu*& menu)
762 // Loads the trigger menu object from the CDB path.
764 HLTDebug("Trying to load trigger menu from '%s'.", cdbPath);
765 if (AliCDBManager::Instance() == NULL)
767 HLTError("CDB manager object not found.");
770 AliCDBStorage* store = AliCDBManager::Instance()->GetDefaultStorage();
773 HLTError("Could not get the the default storage for the CDB.");
776 Int_t version = store->GetLatestVersion(cdbPath, GetRunNo());
779 HLTError("Could not find an entry in the CDB for \"%s\".", cdbPath);
782 Int_t subVersion = store->GetLatestSubVersion(cdbPath, GetRunNo(), version);
783 AliCDBEntry* entry = AliCDBManager::Instance()->Get(cdbPath, GetRunNo(), version, subVersion);
786 HLTError("Could not get the CDB entry for \"%s\".", cdbPath);
789 TObject* obj = entry->GetObject();
792 HLTError("Configuration object for \"%s\" is missing.", cdbPath);
795 if (obj->IsA() != AliHLTTriggerMenu::Class())
797 HLTError("Wrong type for configuration object in \"%s\". Found a %s but we expect a AliHLTTriggerMenu.",
798 cdbPath, obj->ClassName()
802 menu = static_cast<AliHLTTriggerMenu*>(obj);
807 void AliHLTGlobalTriggerComponent::GenerateFileName(TString& name, TString& filename)
809 // Creates a unique file name for the generated code.
811 TUUID guid = GenerateGUID();
818 fUniqueID = bufAsInt[0];
819 TString guidstr = guid.AsString();
820 // Replace '-' with '_' in string.
821 for (int i = 0; i < guidstr.Length(); ++i)
823 if (guidstr[i] == '-') guidstr[i] = '_';
825 name = "AliHLTGlobalTriggerImpl_";
827 filename = name + ".cxx";
829 // For good measure, check that the file names are not used. If they are then regenerate.
830 fstream file(filename.Data(), ios_base::in);
834 GenerateFileName(name, filename);
839 int AliHLTGlobalTriggerComponent::GenerateTrigger(
840 const AliHLTTriggerMenu* menu, TString& name, TString& filename,
841 const TClonesArray& includePaths, const TClonesArray& includeFiles
844 // Generates the global trigger class that will implement the specified trigger menu.
845 // See header for more details.
847 GenerateFileName(name, filename);
848 HLTDebug("Generating custom HLT trigger class named %s, in file %s, using trigger menu %p.",
849 name.Data(), filename.Data(), ((void*)menu)
852 // Open a text file to write the code and generate the new class.
853 fstream code(filename.Data(), ios_base::out | ios_base::trunc);
856 HLTError("Could not open file '%s' for writing.", filename.Data());
860 TClonesArray symbols(AliHLTTriggerMenuSymbol::Class());
861 int result = BuildSymbolList(menu, symbols);
862 if (result != 0) return result;
864 code << "#if !defined(__CINT__) || defined(__MAKECINT__)" << endl;
865 code << "#include <cstring>" << endl;
866 code << "#include \"TClass.h\"" << endl;
867 code << "#include \"TString.h\"" << endl;
868 code << "#include \"TClonesArray.h\"" << endl;
869 code << "#include \"TRandom3.h\"" << endl;
870 code << "#include \"AliHLTLogging.h\"" << endl;
871 code << "#include \"AliHLTGlobalTrigger.h\"" << endl;
872 code << "#include \"AliHLTGlobalTriggerDecision.h\"" << endl;
873 code << "#include \"AliHLTDomainEntry.h\"" << endl;
874 code << "#include \"AliHLTTriggerDomain.h\"" << endl;
875 code << "#include \"AliHLTReadoutList.h\"" << endl;
876 code << "#include \"AliHLTTriggerMenu.h\"" << endl;
877 code << "#include \"AliHLTTriggerMenuItem.h\"" << endl;
878 code << "#include \"AliHLTTriggerMenuSymbol.h\"" << endl;
880 // Add any include files that were specified on the command line.
881 for (Int_t i = 0; i < includeFiles.GetEntriesFast(); i++)
883 TString file = static_cast<TObjString*>(includeFiles.UncheckedAt(i))->String();
884 code << "#include \"" << file.Data() << "\"" << endl;
889 code << "#else" << endl;
890 code << "const char* gFunctionName = \"???\";" << endl;
891 code << "#define HLTDebug(msg) if (CheckFilter(kHLTLogDebug) && CheckGroup(Class_Name())) SendMessage(kHLTLogDebug, Class_Name(), ::gFunctionName, __FILE__, __LINE__, msg)" << endl;
893 code << "#endif" << endl;
895 code << "class " << name << " :" << endl;
896 // Add appropriate #ifdef sections since we need to prevent inheritance from
897 // AliHLTGlobalTrigger. CINT does not seem to support multiple inheritance nor
898 // multiple levels of inheritance. Neither of the following schemes worked:
900 // 1) class AliHLTGlobalTrigger : public AliHLTLogging {};
901 // class AliHLTGlobalTriggerImpl_xyz : public AliHLTGlobalTrigger {};
903 // 2) class AliHLTGlobalTrigger {};
904 // class AliHLTGlobalTriggerImpl_xyz : public AliHLTGlobalTrigger, public AliHLTLogging {};
906 // Thus, we are forced to just inherit from AliHLTLogging when running in the CINT
907 // interpreter. But we anyway have to call the global trigger implementation class
908 // through the AliHLTGlobalTriggerWrapper so this is not such a problem.
909 code << "#if !defined(__CINT__) || defined(__MAKECINT__)" << endl;
910 code << " public AliHLTGlobalTrigger," << endl;
911 code << "#endif" << endl;
912 code << " public AliHLTLogging" << endl;
914 code << "public:" << endl;
916 // Generate constructor method.
917 code << " " << name << "() :" << endl;
918 code << "#if !defined(__CINT__) || defined(__MAKECINT__)" << endl;
919 code << " AliHLTGlobalTrigger()," << endl;
920 code << "#endif" << endl;
921 code << " AliHLTLogging()";
922 // Write the symbols in the trigger menu in the initialisation list.
923 for (Int_t i = 0; i < symbols.GetEntriesFast(); i++)
926 AliHLTTriggerMenuSymbol* symbol = static_cast<AliHLTTriggerMenuSymbol*>( symbols.UncheckedAt(i) );
927 code << " " << symbol->Name() << "()," << endl;
928 if (strcmp(symbol->ObjectClass(), "AliHLTTriggerDecision") == 0)
930 code << " " << symbol->Name() << "TriggerDomain()," << endl;
932 code << " " << symbol->Name() << "DomainEntry(kAliHLTAnyDataType)";
934 for (UInt_t i = 0; i < menu->NumberOfItems(); i++)
936 code << "," << endl << " fMenuItemDescription" << i << "()";
938 code << endl << " {" << endl;
941 code << "#ifdef __CINT__" << endl;
942 code << " gFunctionName = \"" << name.Data() <<"\";" << endl;
943 code << "#endif" << endl;
944 code << " SetLocalLoggingLevel(kHLTLogAll);" << endl;
945 code << " HLTDebug(Form(\"Creating new instance at %p.\", this));" << endl;
947 code << " }" << endl;
949 code << " virtual ~" << name << "() {" << endl;
952 code << "#ifdef __CINT__" << endl;
953 code << " gFunctionName = \"~" << name.Data() << "\";" << endl;
954 code << "#endif" << endl;
955 code << " HLTDebug(Form(\"Deleting instance at %p.\", this));" << endl;
957 code << " }" << endl;
959 // Generate the FillFromMenu method.
960 code << " virtual void FillFromMenu(const AliHLTTriggerMenu& menu) {" << endl;
963 code << "#ifdef __CINT__" << endl;
964 code << " gFunctionName = \"FillFromMenu\";" << endl;
965 code << "#endif" << endl;
966 code << " HLTDebug(Form(\"Filling description entries from trigger menu for global trigger %p.\", this));" << endl;
968 code << " fCounter.Set(menu.NumberOfItems());" << endl;
969 for (UInt_t i = 0; i < menu->NumberOfItems(); i++)
971 code << " fMenuItemDescription" << i << " = (menu.Item(" << i
972 << ") != NULL) ? menu.Item(" << i << ")->Description() : \"\";" << endl;
976 code << " HLTDebug(Form(\"Finished filling description entries from trigger menu.\"));" << endl;
977 code << " HLTDebug(Form(\"Filling domain entries from trigger menu symbols for global trigger %p.\", this));" << endl;
979 code << " for (Int_t i = 0; i < menu.SymbolArray().GetEntriesFast(); i++) {" << endl;
980 // 30 Oct 2009 - CINT sometimes evaluates the dynamic_cast incorrectly.
981 // Have to use the TClass system for extra protection.
982 code << " if (menu.SymbolArray().UncheckedAt(i) == NULL) continue;" << endl;
983 code << " if (menu.SymbolArray().UncheckedAt(i)->IsA() != AliHLTTriggerMenuSymbol::Class()) continue;" << endl;
984 code << " const AliHLTTriggerMenuSymbol* symbol = dynamic_cast<const"
985 " AliHLTTriggerMenuSymbol*>(menu.SymbolArray().UncheckedAt(i));" << endl;
986 code << " if (symbol == NULL) continue;" << endl;
987 for (Int_t i = 0; i < symbols.GetEntriesFast(); i++)
989 AliHLTTriggerMenuSymbol* symbol = static_cast<AliHLTTriggerMenuSymbol*>( symbols.UncheckedAt(i) );
990 code << " if (strcmp(symbol->RealName(), \"" << symbol->RealName() << "\") == 0) {" << endl;
993 code << " HLTDebug(Form(\"Assinging domain entry value corresponding with symbol '%s' to '%s'.\","
994 " symbol->RealName(), symbol->BlockType().AsString().Data()));" << endl;
996 code << " " << symbol->Name() << "DomainEntry = symbol->BlockType();" << endl;
997 code << " continue;" << endl;
998 code << " }" << endl;
1000 code << " }" << endl;
1001 // The following is an optimisation where symbols without any assignment operators
1002 // are treated as constant and only initialised in FillFromMenu rather than reseting
1003 // them in the NewEvent method.
1004 // Note: we putting this initialisation into the constructor can lead to seg faults
1005 // under CINT interpretation. Thus we must put it into the FillFromMenu method instead.
1006 for (Int_t i = 0; i < symbols.GetEntriesFast(); i++)
1008 AliHLTTriggerMenuSymbol* symbol = static_cast<AliHLTTriggerMenuSymbol*>( symbols.UncheckedAt(i) );
1009 if (TString(symbol->AssignExpression()) != "") continue;
1010 if (strcmp(symbol->ObjectClass(), "AliHLTTriggerDecision") == 0) continue;
1011 // CINT has problems with the implicit equals operator for complex types, so if
1012 // the type has an equals operater we need to write the operator call explicitly.
1013 TClass* clas = TClass::GetClass(symbol->Type());
1014 if (clas != NULL and clas->GetMethodAny("operator=") != NULL)
1016 code << " " << symbol->Name() << ".operator = (" << symbol->DefaultValue() << ");" << endl;
1020 code << " " << symbol->Name() << " = " << symbol->DefaultValue() << ";" << endl;
1025 code << " HLTDebug(Form(\"Finished filling domain entries from trigger menu symbols.\"));" << endl;
1027 code << " }" << endl;
1029 // Generate the NewEvent method.
1030 code << " virtual void NewEvent() {" << endl;
1033 code << "#ifdef __CINT__" << endl;
1034 code << " gFunctionName = \"NewEvent\";" << endl;
1035 code << "#endif" << endl;
1036 code << " HLTDebug(Form(\"New event for global trigger object %p, initialising variables to default values.\", this));" << endl;
1038 // Write code to initialise the symbols in the trigger menu to their default values.
1039 for (Int_t i = 0; i < symbols.GetEntriesFast(); i++)
1041 AliHLTTriggerMenuSymbol* symbol = static_cast<AliHLTTriggerMenuSymbol*>( symbols.UncheckedAt(i) );
1042 // The following is an optimisation. If the symbol does not have an assignment expression
1043 // then it is effectively a constant symbol and can be initialised earlier and only once.
1044 // In this case we initialise it in the FillFromMenu method instead.
1045 if (TString(symbol->AssignExpression()) == "") continue;
1046 // CINT has problems with the implicit equals operator for complex types, so if
1047 // the type has an equals operater we need to write the operator call explicitly.
1048 TClass* clas = TClass::GetClass(symbol->Type());
1049 if (clas != NULL and clas->GetMethodAny("operator=") != NULL)
1051 code << " " << symbol->Name() << ".operator = (" << symbol->DefaultValue() << ");" << endl;
1055 code << " " << symbol->Name() << " = " << symbol->DefaultValue() << ";" << endl;
1057 if (strcmp(symbol->ObjectClass(), "AliHLTTriggerDecision") == 0)
1059 code << " " << symbol->Name() << "TriggerDomain.Clear();" << endl;
1064 code << " HLTDebug(Form(\"Finished initialising variables.\"));" << endl;
1066 code << " }" << endl;
1068 // Generate the Add method.
1069 bool haveAssignments = false;
1070 for (Int_t i = 0; i < symbols.GetEntriesFast(); i++)
1072 // First check if we have any symbols with assignment expressions.
1073 // This is needed to get rid of the on the fly compilation warning about '_object_' not being used.
1074 AliHLTTriggerMenuSymbol* symbol = static_cast<AliHLTTriggerMenuSymbol*>( symbols.UncheckedAt(i) );
1075 TString expr = symbol->AssignExpression();
1076 if (expr == "") continue; // Skip entries that have no assignment expression.
1077 haveAssignments = true;
1080 if (haveAssignments or fDebugMode)
1082 code << " virtual void Add(const TObject* _object_, const AliHLTComponentDataType& _type_, AliHLTUInt32_t _spec_) {" << endl;
1086 code << " virtual void Add(const TObject* /*_object_*/, const AliHLTComponentDataType& _type_, AliHLTUInt32_t _spec_) {" << endl;
1090 code << "#ifdef __CINT__" << endl;
1091 code << " gFunctionName = \"Add\";" << endl;
1092 code << "#endif" << endl;
1094 code << " AliHLTDomainEntry _type_spec_(_type_, _spec_);" << endl;
1097 code << " HLTDebug(Form(\"Adding TObject %p, with class name '%s' from data block"
1098 " '%s', to global trigger object %p\", _object_, _object_->ClassName(),"
1099 " _type_spec_.AsString().Data(), this));" << endl;
1100 code << " _object_->Print();" << endl;
1101 code << " bool _object_assigned_ = false;" << endl;
1103 for (Int_t i = 0; i < symbols.GetEntriesFast(); i++)
1105 // Write code to check if the block type, specification and class name is correct.
1106 // Then write code to assign the variable from the input object.
1107 AliHLTTriggerMenuSymbol* symbol = static_cast<AliHLTTriggerMenuSymbol*>( symbols.UncheckedAt(i) );
1108 TString expr = symbol->AssignExpression();
1109 if (expr == "") continue; // Skip entries that have no assignment expression.
1110 bool isTrigDecision = strcmp(symbol->ObjectClass(), "AliHLTTriggerDecision") == 0;
1115 code << " HLTDebug(Form(\"Trying to match input object to class '"
1116 << symbol->ObjectClass() << "', trigger name '" << symbol->RealName()
1117 << "' and block type '%s'\", " << symbol->Name()
1118 << "DomainEntry.AsString().Data()));" << endl;
1122 code << " HLTDebug(Form(\"Trying to match input object to class '"
1123 << symbol->ObjectClass() << "' and block type '%s'\", "
1124 << symbol->Name() << "DomainEntry.AsString().Data()));" << endl;
1127 // 30 Oct 2009 - CINT sometimes evaluates the dynamic_cast incorrectly.
1128 // Have to use the TClass system for extra protection.
1129 code << " const " << symbol->ObjectClass() << "* " << symbol->Name() << "_object_ = NULL;" << endl;
1130 code << " if (_object_->InheritsFrom(" << symbol->ObjectClass() << "::Class())) " << symbol->Name()
1131 << "_object_ = dynamic_cast<const " << symbol->ObjectClass()
1132 << "*>(_object_);" << endl;
1133 code << " if (" << symbol->Name() << "_object_ != NULL && ";
1136 code << "strcmp(" << symbol->Name() << "_object_->Name(), \""
1137 << symbol->RealName() << "\") == 0 && ";
1139 code << symbol->Name() << "DomainEntry == _type_spec_) {" << endl;
1140 TString fullname = symbol->Name();
1141 fullname += "_object_";
1142 expr.ReplaceAll("this", fullname);
1143 code << " this->" << symbol->Name() << " = " << expr.Data() << ";" << endl;
1146 code << " this->" << symbol->Name() << "TriggerDomain = "
1147 << fullname.Data() << "->TriggerDomain();" << endl;
1151 code << " HLTDebug(Form(\"Added TObject %p with class name '%s' to variable "
1152 << symbol->Name() << "\", _object_, _object_->ClassName()));" << endl;
1153 code << " _object_assigned_ = true;" << endl;
1155 code << " }" << endl;
1159 code << " if (! _object_assigned_) {" << endl;
1160 code << " HLTDebug(Form(\"Did not assign TObject %p"
1161 " with class name '%s' to any variable.\", _object_, _object_->ClassName()));"
1163 code << " }" << endl;
1165 code << " }" << endl;
1167 // Generate the CalculateTriggerDecision method.
1168 // This requires code to be generated that checks which items in the trigger menu
1169 // have their conditions asserted and then the trigger domain is generated from
1170 // those fired items.
1171 // The processing will start from the highest priority trigger group and stop
1172 // after at least one trigger from the current priority group being processed
1173 // is positive. For each priority group all the trigger menu items are checked.
1174 // Their combined trigger condition expression must be true for the trigger priority
1175 // group to be triggered positive. The full condition expression is formed by
1176 // concatenating the individual condition expressions. If no trailing operators are
1177 // used in the individual expressions then the default condition operator is placed
1178 // between two concatenated condition expressions.
1179 // If a trigger priority group has at least one trigger fired then the trigger domain
1180 // is calculated such that it will give the same result as the concatenated trigger
1181 // domain merging expressions for all the individual trigger menu items with
1182 // positive results. Again, if no trailing operators are used in the individual
1183 // merging expressions then the default domain operator is placed between two
1184 // expression fragments.
1185 code << " virtual bool CalculateTriggerDecision(bool& _trigger_result_, AliHLTTriggerDomain& _domain_, TString& _description_) {" << endl;
1188 code << "#ifdef __CINT__" << endl;
1189 code << " gFunctionName = \"CalculateTriggerDecision\";" << endl;
1190 code << "#endif" << endl;
1191 code << " HLTDebug(Form(\"Calculating global HLT trigger result with trigger object at %p.\", this));" << endl;
1194 // Build a list of priorities used in the trigger menu.
1195 std::vector<UInt_t> priorities;
1196 for (UInt_t i = 0; i < menu->NumberOfItems(); i++)
1198 const AliHLTTriggerMenuItem* item = menu->Item(i);
1199 bool priorityNotInList = std::find(priorities.begin(), priorities.end(), item->Priority()) == priorities.end();
1200 if (priorityNotInList) priorities.push_back(item->Priority());
1202 std::sort(priorities.begin(), priorities.end(), AliHLTDescendingNumbers);
1203 // From the priority list, build the priority groups in the correct order,
1204 // i.e. highest priority first.
1205 // The priority group is a list of vectors of integers. The integers are the
1206 // index numbers into the trigger menu item list for the items which form part
1207 // of the priority group.
1208 std::vector<std::vector<Int_t> > priorityGroup;
1209 priorityGroup.insert(priorityGroup.begin(), priorities.size(), std::vector<Int_t>());
1210 for (size_t n = 0; n < priorities.size(); n++)
1212 UInt_t priority = priorities[n];
1213 for (UInt_t i = 0; i < menu->NumberOfItems(); i++)
1215 const AliHLTTriggerMenuItem* item = menu->Item(i);
1216 if (item->Priority() == priority) priorityGroup[n].push_back(i);
1220 for (size_t n = 0; n < priorityGroup.size(); n++)
1224 code << " HLTDebug(Form(\"Processing trigger priority group " << priorities[n] << "\"));" << endl;
1227 if (n == 0) code << "UInt_t ";
1228 code << "_previous_match_ = 0xFFFFFFFF;" << endl;
1230 if (n == 0) code << "bool ";
1231 code << "_trigger_matched_ = false;" << endl;
1233 if (n == 0) code << "bool ";
1234 code << "_group_result_ = false;" << endl;
1235 std::vector<TString> conditionOperator;
1236 conditionOperator.insert(conditionOperator.begin(), priorityGroup[n].size(), TString(""));
1237 std::vector<TString> domainOperator;
1238 domainOperator.insert(domainOperator.begin(), priorityGroup[n].size(), TString(""));
1239 for (size_t m = 0; m < priorityGroup[n].size(); m++)
1241 UInt_t i = priorityGroup[n][m];
1242 const AliHLTTriggerMenuItem* item = menu->Item(i);
1243 TString triggerCondition = item->TriggerCondition();
1244 TString mergeExpr = item->MergeExpression();
1245 // Replace the symbols found in the trigger condition and merging expressions
1246 // with appropriate compilable versions.
1247 for (Int_t j = 0; j < symbols.GetEntriesFast(); j++)
1249 AliHLTTriggerMenuSymbol* symbol = static_cast<AliHLTTriggerMenuSymbol*>( symbols.UncheckedAt(j) );
1250 bool symbolNamesDifferent = strcmp(symbol->RealName(), symbol->Name()) != 0;
1251 if (strcmp(symbol->ObjectClass(), "AliHLTTriggerDecision") == 0)
1253 TString newname = symbol->Name();
1254 newname += "TriggerDomain";
1255 mergeExpr.ReplaceAll(symbol->RealName(), newname);
1259 if (symbolNamesDifferent) mergeExpr.ReplaceAll(symbol->RealName(), symbol->Name());
1261 if (symbolNamesDifferent) triggerCondition.ReplaceAll(symbol->RealName(), symbol->Name());
1263 // We allow the trigger conditions and merging expressions to have trailing operators.
1264 // Thus, we need to extract the operators and cleanup the expressions so that they will
1265 // compile. This means that we silently ignore the trailing operator if not needed.
1266 if (ExtractedOperator(triggerCondition, conditionOperator[m]))
1268 // If the trailing operator is the same as the default operator then reset
1269 // the value in the operator list so that the default is used in the generated
1270 // code. This creates more compact code.
1271 if (conditionOperator[m] == menu->DefaultConditionOperator()) conditionOperator[m] = "";
1273 if (ExtractedOperator(mergeExpr, domainOperator[m]))
1275 if (domainOperator[m] == menu->DefaultDomainOperator()) domainOperator[m] = "";
1279 code << " HLTDebug(Form(\"Trying trigger condition " << i
1280 << " (Description = '%s').\", fMenuItemDescription" << i << ".Data()));"
1284 if (n == 0 and m == 0) code << "bool ";
1285 code << "_item_result_ = false;" << endl;
1286 code << " if (" << triggerCondition << ") {" << endl;
1287 code << " ++fCounter[" << i << "];" << endl;
1288 const char* indentation = "";
1289 // Generate the code to handle the prescalar and scale-down
1290 bool havePrescalar = item->PreScalar() != 0;
1291 bool haveScaledown = item->ScaleDown() < 1;
1292 if (havePrescalar or haveScaledown)
1296 if (havePrescalar) code << "(fCounter[" << i << "] % " << item->PreScalar() << ") == 1";
1297 if (havePrescalar and haveScaledown) code << " && ";
1300 std::streamsize oldprecision = code.precision(17);
1301 code << "gRandom->Rndm() < " << item->ScaleDown();
1302 code.precision(oldprecision);
1304 code << ") {" << endl;
1306 code << indentation << " _item_result_ = true;" << endl;
1309 code << indentation << " HLTDebug(Form(\"Matched trigger condition " << i
1310 << " (Description = '%s').\", fMenuItemDescription" << i << ".Data()));" << endl;
1312 if (havePrescalar or haveScaledown)
1314 code << " }" << endl;
1316 code << " }" << endl;
1319 // Since this is the first item of the trigger group,
1320 // the generated trigger logic can be simplified a little.
1321 code << " _group_result_ = _item_result_;" << endl;
1322 code << " if (_item_result_) {" << endl;
1323 code << " _domain_ = " << mergeExpr.Data() << ";" << endl;
1324 code << " _description_ = fMenuItemDescription" << i << ";" << endl;
1325 code << " _previous_match_ = " << i << ";" << endl;
1326 code << " _trigger_matched_ = true;" << endl;
1327 code << " }" << endl;
1331 if (conditionOperator[m-1] == "")
1333 code << " _group_result_ = _group_result_ "
1334 << menu->DefaultConditionOperator() << " _item_result_;" << endl;
1338 code << " _group_result_ = _group_result_ "
1339 << conditionOperator[m-1] << " _item_result_;" << endl;
1341 code << " if (_item_result_) {" << endl;
1342 code << " if (_trigger_matched_) {" << endl;
1343 bool switchWillBeEmpty = true;
1344 for (size_t k = 0; k < m; k++)
1346 if (domainOperator[k] == "") continue;
1347 switchWillBeEmpty = false;
1349 if (switchWillBeEmpty)
1351 code << " _domain_ = _domain_ " << menu->DefaultDomainOperator() << " "
1352 << mergeExpr.Data() << ";" << endl;
1356 code << " switch(_previous_match_) {" << endl;
1357 for (size_t k = 0; k < m; k++)
1359 if (domainOperator[k] == "") continue;
1360 code << " case " << k << ": _domain_ = _domain_ "
1361 << domainOperator[k] << " " << mergeExpr.Data() << "; break;" << endl;
1363 code << " default: _domain_ = _domain_ "
1364 << menu->DefaultDomainOperator() << " " << mergeExpr.Data() << ";" << endl;
1365 code << " }" << endl;
1367 code << " _description_ += \",\";" << endl;
1368 code << " _description_ += fMenuItemDescription" << i << ";" << endl;
1369 code << " } else {" << endl;
1370 code << " _domain_ = " << mergeExpr.Data() << ";" << endl;
1371 code << " _description_ = fMenuItemDescription" << i << ";" << endl;
1372 code << " }" << endl;
1373 code << " _previous_match_ = " << i << ";" << endl;
1374 code << " _trigger_matched_ = true;" << endl;
1375 code << " }" << endl;
1378 code << " if (_group_result_) {" << endl;
1381 if (n < priorities.size() - 1)
1383 code << " HLTDebug(Form(\"Matched triggers in trigger priority group " << priorities[n]
1384 << ". Stopping processing here because all other trigger groups have lower priority.\"));" << endl;
1388 code << " HLTDebug(Form(\"Matched triggers in trigger priority group " << priorities[n] << ".\"));" << endl;
1391 bool methodReturnResult = true;
1392 if (priorityGroup[n].size() > 0)
1394 const AliHLTTriggerMenuItem* item = menu->Item(priorityGroup[n][0]);
1395 methodReturnResult = item->DefaultResult();
1397 // Check to see if the items of the group all have the same default result.
1398 // If not then warn the user since only the first item's value will be used.
1399 for (size_t m = 1; m < priorityGroup[n].size(); m++)
1401 const AliHLTTriggerMenuItem* item = menu->Item(priorityGroup[n][m]);
1402 if (item->DefaultResult() != methodReturnResult)
1404 HLTWarning("Found items with different default results set for priority group %d."
1405 "Will only use the value from the first item.",
1411 code << " _trigger_result_ = " << (methodReturnResult ? "true" : "false") << ";" << endl;
1412 code << " return true;" << endl;
1413 code << " }" << endl;
1415 code << " _domain_.Clear();" << endl;
1416 code << " _description_ = \"\";" << endl;
1417 code << " _trigger_result_ = " << (menu->DefaultResult() ? "true" : "false") << ";" << endl;
1418 code << " return false;" << endl;
1419 code << " }" << endl;
1421 // Generate getter and setter methods for the counters.
1422 code << " const TArrayL64& GetCounters() const { return fCounter; }" << endl;
1423 code << " void SetCounters(const TArrayL64& counters) { fCounter = counters; }" << endl;
1425 code << "private:" << endl;
1426 // Add the symbols in the trigger menu to the list of private variables.
1427 for (Int_t i = 0; i < symbols.GetEntriesFast(); i++)
1429 AliHLTTriggerMenuSymbol* symbol = static_cast<AliHLTTriggerMenuSymbol*>( symbols.UncheckedAt(i) );
1430 code << " " << symbol->Type() << " " << symbol->Name() << ";" << endl;
1431 if (strcmp(symbol->ObjectClass(), "AliHLTTriggerDecision") == 0)
1433 code << " AliHLTTriggerDomain " << symbol->Name() << "TriggerDomain;" << endl;
1435 code << " AliHLTDomainEntry " << symbol->Name() << "DomainEntry;" << endl;
1437 for (UInt_t i = 0; i < menu->NumberOfItems(); i++)
1439 code << " TString fMenuItemDescription" << i << ";" << endl;
1441 code << " TArrayL64 fCounter;" << endl;
1442 code << "#if !defined(__CINT__) || defined(__MAKECINT__)" << endl;
1443 code << " ClassDef(" << name.Data() << ", 0)" << endl;
1444 code << "#else" << endl;
1445 code << " virtual const char* Class_Name() const { return \"" << name.Data() << "\"; }" << endl;
1446 code << "#endif" << endl;
1447 code << "};" << endl;
1448 code << "#if !defined(__CINT__) || defined(__MAKECINT__)" << endl;
1449 code << "ClassImp(" << name.Data() << ")" << endl;
1450 code << "#endif" << endl;
1454 // Now we need to compile and load the new class.
1455 result = LoadTriggerClass(filename, includePaths);
1460 int AliHLTGlobalTriggerComponent::LoadTriggerClass(
1461 const char* filename, const TClonesArray& includePaths
1464 // Loads the code for a custom global trigger class implementation on the fly.
1466 HLTDebug("Loading HLT trigger class from file '%s'.", filename);
1468 TString compiler = gSystem->GetBuildCompilerVersion();
1469 if (fRuntimeCompile && (compiler.Contains("gcc") or compiler.Contains("icc")))
1471 TString includePath;
1472 #if defined(PKGINCLUDEDIR)
1473 // this is especially for the HLT build system where the package is installed
1474 // in a specific directory including proper treatment of include files
1475 includePath.Form("-I%s", PKGINCLUDEDIR);
1477 // the default AliRoot behavior, all include files can be found in the
1478 // $ALICE_ROOT subfolders
1479 includePath = "-I${ALICE_ROOT}/include -I${ALICE_ROOT}/HLT/BASE -I${ALICE_ROOT}/HLT/trigger";
1481 // Add any include paths that were specified on the command line.
1482 for (Int_t i = 0; i < includePaths.GetEntriesFast(); i++)
1484 TString path = static_cast<TObjString*>(includePaths.UncheckedAt(i))->String();
1485 includePath += " -I";
1486 includePath += path;
1488 HLTDebug("using include settings: %s", includePath.Data());
1489 gSystem->SetIncludePath(includePath);
1490 gSystem->SetFlagsOpt("-O3 -DNDEBUG");
1491 gSystem->SetFlagsDebug("-g3 -DDEBUG -D__DEBUG");
1496 result = gSystem->CompileMacro(filename, "g");
1500 result = gSystem->CompileMacro(filename, "O");
1502 if (result != kTRUE)
1504 HLTFatal("Could not compile and load global trigger menu implementation.");
1510 // Store the library state to be checked later in UnloadTriggerClass.
1511 fLibStateAtLoad = gSystem->GetLibraries();
1513 // If we do not support the compiler then try interpret the class instead.
1514 TString cmd = ".L ";
1516 Int_t errorcode = TInterpreter::kNoError;
1517 gROOT->ProcessLine(cmd, &errorcode);
1518 if (errorcode != TInterpreter::kNoError)
1520 HLTFatal("Could not load interpreted global trigger menu implementation"
1521 " (Interpreter error code = %d).",
1532 int AliHLTGlobalTriggerComponent::UnloadTriggerClass(const char* filename)
1534 // Unloads the code previously loaded by LoadTriggerClass.
1536 HLTDebug("Unloading HLT trigger class in file '%s'.", filename);
1538 TString compiler = gSystem->GetBuildCompilerVersion();
1539 if (fRuntimeCompile && (compiler.Contains("gcc") or compiler.Contains("icc")))
1541 // Generate the library name.
1542 TString libname = filename;
1543 Ssiz_t dotpos = libname.Last('.');
1544 if (0 <= dotpos and dotpos < libname.Length()) libname[dotpos] = '_';
1546 libname += gSystem->GetSoExt();
1548 // This is a workaround for a problem with unloading shared libraries in ROOT.
1549 // If the trigger logic library is loaded before the libAliHLTHOMER.so library
1550 // or any other library is loaded afterwards, then during the gInterpreter->UnloadFile
1551 // call all the subsequent libraries get unloded. This means that any objects created
1552 // from classes implemented in the libAliHLTHOMER.so library will generate segfaults
1553 // since the executable code has been unloaded.
1554 // We need to check if there are any more libraries loaded after the class we
1555 // are unloading and in that case don't unload the class.
1556 TString libstring = gSystem->GetLibraries();
1557 TString token, lastlib;
1559 Int_t numOfLibs = 0, posOfLib = -1;
1560 while (libstring.Tokenize(token, from, " "))
1564 if (token.Contains(libname)) posOfLib = numOfLibs;
1566 if (numOfLibs != posOfLib)
1568 HLTWarning(Form("ROOT limitation! Cannot properly cleanup and unload the shared"
1569 " library '%s' since another library '%s' was loaded afterwards. Trying to"
1570 " unload this library will remove the others and lead to serious memory faults.",
1571 libname.Data(), lastlib.Data()
1578 if ((path = gSystem->DynamicPathName(libname)) != NULL)
1580 result = gInterpreter->UnloadFile(path);
1583 if (result != TInterpreter::kNoError) return -ENOENT;
1587 // This is again a workaround for the problem with unloading files in ROOT.
1588 // If the trigger logic class is loaded before the libAliHLTHOMER.so library
1589 // or any other library is loaded afterwards, then during the gInterpreter->UnloadFile
1590 // call all the subsequent libraries get unloded.
1591 // We need to check if the list of loaded libraries has changed since the last
1592 // call to LoadTriggerClass. If it has then don't unload the class.
1593 if (fLibStateAtLoad != gSystem->GetLibraries())
1595 TString libstring = gSystem->GetLibraries();
1598 while (libstring.Tokenize(token, from, " "))
1600 if (not fLibStateAtLoad.Contains(token)) break;
1602 HLTWarning(Form("ROOT limitation! Cannot properly cleanup and unload the file"
1603 " '%s' since another library '%s' was loaded afterwards. Trying to unload"
1604 " this file will remove the other library and lead to serious memory faults.",
1605 filename, token.Data()
1610 // If we did not compile the trigger logic then remove the interpreted class.
1611 TString cmd = ".U ";
1613 Int_t errorcode = TInterpreter::kNoError;
1614 gROOT->ProcessLine(cmd, &errorcode);
1615 if (errorcode != TInterpreter::kNoError)
1617 HLTFatal("Could not unload interpreted global trigger menu implementation"
1618 " (Interpreter error code = %d).",
1629 int AliHLTGlobalTriggerComponent::FindSymbol(const char* name, const TClonesArray& list)
1631 // Searches for the named symbol in the given list.
1632 // See header for more details.
1634 for (int i = 0; i < list.GetEntriesFast(); i++)
1636 const AliHLTTriggerMenuSymbol* symbol = dynamic_cast<const AliHLTTriggerMenuSymbol*>( list.UncheckedAt(i) );
1637 if (symbol == NULL) continue;
1638 if (strcmp(symbol->Name(), name) == 0) return i;
1647 * Helper routine to compare two trigger menu symbols to see if b is a subset of a.
1648 * \returns true if b is a subset of a or they are unrelated.
1650 bool AliHLTCheckForContainment(const AliHLTTriggerMenuSymbol* a, const AliHLTTriggerMenuSymbol* b)
1652 TString bstr = b->Name();
1653 return bstr.Contains(a->Name());
1656 } // end of namespace
1659 int AliHLTGlobalTriggerComponent::BuildSymbolList(const AliHLTTriggerMenu* menu, TClonesArray& list)
1661 // Builds the list of symbols to use in the custom global trigger menu
1662 // implementation class.
1663 // See header for more details.
1665 // Note: when we build the symbol list we must use the symbol name as returned
1666 // by the Name() method and not the RealName() method when using FindSymbol.
1667 // This is so that we avoid problems with the generated code not compiling
1668 // because names like "abc-xyz" and "abc_xyz" are synonymous.
1669 // Name() returns the converted C++ symbol name as used in the generated code.
1671 for (UInt_t i = 0; i < menu->NumberOfSymbols(); i++)
1673 const AliHLTTriggerMenuSymbol* symbol = menu->Symbol(i);
1674 if (FindSymbol(symbol->Name(), list) != -1)
1676 HLTError("Multiple symbols with the name '%s' defined in the trigger menu.", symbol->Name());
1681 new (list[list.GetEntriesFast()]) AliHLTTriggerMenuSymbol(*symbol);
1683 catch (const std::bad_alloc&)
1685 HLTError("Could not allocate more memory for the symbols list when adding a trigger menu symbol.");
1689 Int_t initialEntryCount = list.GetEntriesFast();
1691 // Note: the \\. must not be the first element in the character class, otherwise
1692 // it is interpreted as an "any character" dot symbol.
1693 TRegexp exp("[_a-zA-Z][-\\._a-zA-Z0-9]*");
1694 for (UInt_t i = 0; i < menu->NumberOfItems(); i++)
1696 const AliHLTTriggerMenuItem* item = menu->Item(i);
1697 TString str = item->TriggerCondition();
1702 Ssiz_t pos = exp.Index(str, &length, start);
1703 if (pos == kNPOS) break;
1706 // Check if there is a numerical character before the found
1707 // regular expression. If so, then the symbol is not a valid one
1708 // and should be skipped.
1711 bool notValid = false;
1714 case '0': case '1': case '2': case '3': case '4':
1715 case '5': case '6': case '7': case '8': case '9':
1722 if (notValid) continue;
1724 TString s = str(pos, length);
1726 if (s == "and" or s == "and_eq" or s == "bitand" or s == "bitor" or
1727 s == "compl" or s == "not" or s == "not_eq" or s == "or" or
1728 s == "or_eq" or s == "xor" or s == "xor_eq" or s == "true" or
1732 // Ignore iso646.h and other keywords.
1736 // We need to handle the special case where the symbol contains a dot.
1737 // In C++ this is a dereferencing operator. So we need to check if the
1738 // current symbol we are handling starts with the same string as any of
1739 // the existing symbols defined manually in the symbols table.
1740 // If we do find such a case then revert to treating the dot as an operator
1741 // rather than part of the symbol name. i.e. skip adding the automatic symbol.
1742 bool dereferencedSymbol = false;
1743 for (int j = 0; j < initialEntryCount; j++)
1745 const AliHLTTriggerMenuSymbol* symbol = dynamic_cast<const AliHLTTriggerMenuSymbol*>( list.UncheckedAt(j) );
1746 if (symbol == NULL) continue;
1747 TString symstr = symbol->Name();
1749 if (s.BeginsWith(symstr))
1751 dereferencedSymbol = true;
1755 if (dereferencedSymbol) continue;
1757 // Need to create the symbols first and check if its name is in the list
1758 // before actually adding it to the symbols list.
1759 AliHLTTriggerMenuSymbol newSymbol;
1760 newSymbol.Name(s.Data());
1761 newSymbol.Type("bool");
1762 newSymbol.ObjectClass("AliHLTTriggerDecision");
1763 newSymbol.AssignExpression("this->Result()");
1764 newSymbol.DefaultValue("false");
1765 if (FindSymbol(newSymbol.Name(), list) == -1)
1769 new (list[list.GetEntriesFast()]) AliHLTTriggerMenuSymbol(newSymbol);
1771 catch (const std::bad_alloc&)
1773 HLTError("Could not allocate more memory for the symbols list when adding a trigger name symbol.");
1778 while (start < str.Length());
1781 // This last part is necessary to make sure that symbols are replaced in the
1782 // trigger condition and domain merging expressions in a greedy manner.
1783 // I.e. we need to make sure that if one symbol's string representation is
1784 // contained inside another (string subset) that the longer symbol name is
1785 // always first in the symbols list.
1786 // This will work because the symbol table is traversed from first to last
1787 // element and TString::ReplaceAll is used to replace the substrings inside
1788 // the AliHLTGlobalTriggerComponent::GenerateTrigger method.
1789 std::vector<AliHLTTriggerMenuSymbol*> orderedList;
1790 for (Int_t i = 0; i < list.GetEntriesFast(); i++)
1792 orderedList.push_back( static_cast<AliHLTTriggerMenuSymbol*>(list.UncheckedAt(i)) );
1794 std::sort(orderedList.begin(), orderedList.end(), AliHLTCheckForContainment);
1795 //std::sort(orderedList.begin(), orderedList.end());
1796 // Now swap values around according to the orderedList.
1797 for (Int_t i = 0; i < list.GetEntriesFast(); i++)
1799 AliHLTTriggerMenuSymbol* target = static_cast<AliHLTTriggerMenuSymbol*>(list.UncheckedAt(i));
1800 AliHLTTriggerMenuSymbol tmp = *target;
1801 *target = *orderedList[i];
1802 *orderedList[i] = tmp;
1809 bool AliHLTGlobalTriggerComponent::ExtractedOperator(TString& expr, TString& op)
1811 // Extracts the trailing operator from the expression.
1814 // First skip the trailing whitespace.
1815 bool whitespace = true;
1816 for (i = expr.Length()-1; i >= 0 and whitespace; i--)
1820 case ' ': case '\t': case '\r': case '\n':
1827 if (i < 0 or whitespace) return false;
1829 // Now find the first whitespace character before the trailing symbol.
1830 bool nonwhitespace = true;
1831 for (; i >= 0 and nonwhitespace; i--)
1835 case ' ': case '\t': case '\r': case '\n':
1836 nonwhitespace = false;
1839 nonwhitespace = true;
1842 if (i < 0 or nonwhitespace) return false;
1844 // Extract the last symbols and check if it is a valid operator.
1847 if (s == "and" or s == "and_eq" or s == "bitand" or s == "bitor" or
1848 s == "compl" or s == "not" or s == "not_eq" or s == "or" or
1849 s == "or_eq" or s == "xor" or s == "xor_eq" or s == "&&" or
1850 s == "&=" or s == "&" or s == "|" or s == "~" or s == "!" or
1851 s == "!=" or s == "||" or s == "|=" or s == "^" or s == "^=" or
1852 s == "==" or s == "+" or s == "-" or s == "*" or s == "/" or
1853 s == "%" or s == ">" or s == "<" or s == ">=" or s == "<="
1865 bool AliHLTGlobalTriggerComponent::FillSoftwareTrigger()
1867 // Fills the fSoftwareTrigger structure.
1869 if (fCDH == NULL) return false;
1870 UChar_t l1msg = fCDH->GetL1TriggerMessage();
1871 if ((l1msg & 0x1) == 0x0) return false; // skip physics events.
1872 // From here on everything must be a software trigger.
1873 if (((l1msg >> 2) & 0xF) == 0xE)
1875 fSoftwareTrigger.Name("START_OF_DATA");
1876 fSoftwareTrigger.Description("Generated internal start of data trigger.");
1878 else if (((l1msg >> 2) & 0xF) == 0xF)
1880 fSoftwareTrigger.Name("END_OF_DATA");
1881 fSoftwareTrigger.Description("Generated internal end of data trigger.");
1883 else if (((l1msg >> 6) & 0x1) == 0x1)
1885 fSoftwareTrigger.Name("CALIBRATION");
1886 fSoftwareTrigger.Description("Generated internal calibration trigger.");
1890 fSoftwareTrigger.Name("SOFTWARE");
1891 fSoftwareTrigger.Description("Generated internal software trigger.");
1893 UInt_t detectors = fCDH->GetSubDetectors();
1894 fSoftwareTrigger.ReadoutList( AliHLTReadoutList(Int_t(detectors)) );
1899 int AliHLTGlobalTriggerComponent::PrintStatistics(const AliHLTGlobalTrigger* pTrigger, AliHLTComponentLogSeverity level, int offset) const
1901 // print some statistics
1902 int totalEvents=fTotalEventCounter+offset;
1903 const TArrayL64& counters = pTrigger->GetCounters();
1904 if (pTrigger->CallFailed()) return -EPROTO;
1905 for (int i = 0; i < counters.GetSize(); i++) {
1906 ULong64_t count = counters[i];
1908 if (totalEvents>0) ratio=100*(float)count/totalEvents;
1909 HLTLog(level, "Item %d: total events: %d - counted events: %llu (%.1f%%)", i, totalEvents, count, ratio);
1914 int AliHLTGlobalTriggerComponent::AddCTPDecisions(AliHLTGlobalTrigger* pTrigger, const AliHLTCTPData* pCTPData, const AliHLTComponentTriggerData* trigData)
1916 // add trigger decisions for the valid CTP classes
1917 if (!pCTPData || !pTrigger) return 0;
1919 AliHLTUInt64_t triggerMask=pCTPData->Mask();
1920 AliHLTUInt64_t bit0=0x1;
1921 if (!fCTPDecisions) {
1924 fCTPDecisions=new TClonesArray(AliHLTTriggerDecision::Class(), gkNCTPTriggerClasses);
1926 catch (const std::bad_alloc&)
1928 HLTError("Could not allocate memory for the CTP decisions array.");
1931 if (!fCTPDecisions) return -ENOMEM;
1935 fCTPDecisions->ExpandCreate(gkNCTPTriggerClasses);
1937 catch (const std::bad_alloc&)
1939 HLTError("Could not allocate more memory for the CTP decisions array.");
1942 for (int i=0; i<gkNCTPTriggerClasses; i++) {
1943 const char* name=pCTPData->Name(i);
1944 if (triggerMask&(bit0<<i) && name) {
1945 AliHLTTriggerDecision* pDecision=dynamic_cast<AliHLTTriggerDecision*>(fCTPDecisions->At(i));
1948 delete fCTPDecisions;
1952 pDecision->Name(name);
1957 for (int i=0; i<gkNCTPTriggerClasses; i++) {
1958 const char* name=pCTPData->Name(i);
1959 if ((triggerMask&(bit0<<i))==0 || name==NULL) continue;
1960 AliHLTTriggerDecision* pDecision=dynamic_cast<AliHLTTriggerDecision*>(fCTPDecisions->At(i));
1961 HLTDebug("updating CTP trigger decision %d %s (%p casted %p)", i, name, fCTPDecisions->At(i), pDecision);
1962 if (!pDecision) return -ENOENT;
1965 // 13 March 2010 - Optimisation:
1966 // Dont use the EvaluateCTPTriggerClass method, which uses slow TFormula objects.
1967 AliHLTUInt64_t triggers = 0;
1968 if (trigData) triggers = pCTPData->ActiveTriggers(*trigData);
1969 else triggers = pCTPData->Triggers();
1970 result = (triggers&((AliHLTUInt64_t)0x1<<i)) ? true : false;
1971 //if (trigData) result=pCTPData->EvaluateCTPTriggerClass(name, *trigData);
1972 //else result=pCTPData->EvaluateCTPTriggerClass(name);
1973 pDecision->Result(result);
1974 pDecision->TriggerDomain().Clear();
1975 if (trigData) pDecision->TriggerDomain().Add(pCTPData->ReadoutList(*trigData));
1976 else pDecision->TriggerDomain().Add(pCTPData->ReadoutList());
1978 pTrigger->Add(fCTPDecisions->At(i), kAliHLTDataTypeTriggerDecision, kAliHLTVoidDataSpec);
1979 if (pTrigger->CallFailed()) return -EPROTO;