**************************************************************************/
#include "AliEveEventManager.h"
+#include "AliEveEventSelector.h"
+#include "AliEveMacroExecutor.h"
+
+#include <THashList.h>
+#include <TEveElement.h>
#include <TEveManager.h>
+#include <TEveViewer.h>
+#include <AliLog.h>
#include <AliRunLoader.h>
#include <AliRun.h>
+#include <AliESDRun.h>
#include <AliESDEvent.h>
#include <AliESDfriend.h>
+#include <AliAODEvent.h>
+
+#include <AliRecoParam.h>
+#include <AliCentralTrigger.h>
+#include <AliCDBEntry.h>
+#include <AliTriggerClass.h>
+#include <AliTriggerConfiguration.h>
+#include <AliTriggerCluster.h>
+#include <AliDetectorRecoParam.h>
+
#include <AliDAQ.h>
#include <AliRawEventHeaderBase.h>
#include <AliRawReaderRoot.h>
#include <AliRawReaderFile.h>
#include <AliRawReaderDate.h>
-#include <AliMagFMaps.h>
+#include <AliMagF.h>
#include <AliCDBManager.h>
+#include <AliCDBStorage.h>
+#include <AliGRPObject.h>
#include <AliHeader.h>
#include <AliGeomManager.h>
+#include <AliGRPManager.h>
+#include <AliSysInfo.h>
#include <TFile.h>
#include <TTree.h>
#include <TGeoManager.h>
+#include <TGeoGlobalMagField.h>
#include <TSystem.h>
+#include <TTimeStamp.h>
+#include <TPRegexp.h>
+#include <TError.h>
+#include <TEnv.h>
+#include <TString.h>
+#include <TMap.h>
//==============================================================================
//==============================================================================
//______________________________________________________________________________
//
-// Provide interface for loading and navigating standard AliRoot data
-// (AliRunLoader) and ESDs.
+// Provides interface for loading and navigating standard AliRoot data
+// (AliRunLoader), ESD, AOD and RAW.
+//
+// ESDfriend is attached automatically, if the file is found.
//
-// Missing support for raw-data. For now this is handled individually
-// by each sub-detector.
+// AODfriends are not attached automatically as there are several
+// possible files involved. To have a specific AODfriend attached, call
+// static method
+// AliEveEventManager::AddAODfriend("AliAOD.VertexingHF.root");
+// before initializing the event-manager.
//
// Also provides interface to magnetic-field and geometry. Mostly
// intended as wrappers over standard AliRoot functionality for
// convenient use from visualizateion macros.
+//
+// There can be a single main event-manger, it is stored in private
+// data member fgMaster and can be accessed via static member function
+// GetMaster().
+//
+// For event overlaying and embedding one can instantiate additional
+// event-managers via static method AddDependentManager(const TString& path).
+// This interface is under development.
ClassImp(AliEveEventManager)
-AliEveEventManager* gAliEveEvent = 0;
-
Bool_t AliEveEventManager::fgAssertRunLoader = kFALSE;
Bool_t AliEveEventManager::fgAssertESD = kFALSE;
+Bool_t AliEveEventManager::fgAssertAOD = kFALSE;
Bool_t AliEveEventManager::fgAssertRaw = kFALSE;
TString AliEveEventManager::fgESDFileName("AliESDs.root");
+TString AliEveEventManager::fgAODFileName("AliAOD.root");
TString AliEveEventManager::fgRawFileName("raw.root");
-TString AliEveEventManager::fgCdbUri("local://$ALICE_ROOT");
+TString AliEveEventManager::fgCdbUri;
+
+TList* AliEveEventManager::fgAODfriends = 0;
+
+Bool_t AliEveEventManager::fgRawFromStandardLoc = kFALSE;
+
+Bool_t AliEveEventManager::fgGRPLoaded = kFALSE;
+AliMagF* AliEveEventManager::fgMagField = 0;
+AliRecoParam* AliEveEventManager::fgRecoParam = 0;
+Bool_t AliEveEventManager::fgUniformField = kFALSE;
+
+AliEveEventManager* AliEveEventManager::fgMaster = 0;
+AliEveEventManager* AliEveEventManager::fgCurrent = 0;
+
+void AliEveEventManager::InitInternals()
+{
+ // Initialize internal members.
+
+ static const TEveException kEH("AliEveEventManager::InitInternals ");
+
+ if (fgCurrent != 0)
+ {
+ throw(kEH + "Dependent event-managers should be created via static method AddDependentManager().");
+ }
+
+ if (fgMaster == 0)
+ {
+ fgMaster = this;
+ }
+
+ fgCurrent = this;
+
+ fAutoLoadTimer = new TTimer;
+ fAutoLoadTimer->Connect("Timeout()", "AliEveEventManager", this, "AutoLoadNextEvent()");
+ fAutoLoadTimer->Connect("Timeout()", "AliEveEventManager", this, "Timeout()");
+
+ fExecutor = new AliEveMacroExecutor;
+
+ fTransients = new TEveElementList("Transients", "Transient per-event elements.");
+ fTransients->IncDenyDestroy();
+ gEve->AddToListTree(fTransients, kFALSE);
-AliMagF* AliEveEventManager::fgMagField = 0;
+ fTransientLists = new TEveElementList("Transient Lists", "Containers of transient elements.");
+ fTransientLists->IncDenyDestroy();
+ gEve->AddToListTree(fTransientLists, kFALSE);
+ fPEventSelector = new AliEveEventSelector(this);
-AliEveEventManager::AliEveEventManager() :
- TEveEventManager(),
+ fGlobal = new TMap; fGlobal->SetOwnerKeyValue();
+}
+
+AliEveEventManager::AliEveEventManager(const TString& name) :
+ TEveEventManager(name),
fPath ( ), fEventId (-1),
fRunLoader (0),
fESDFile (0), fESDTree (0), fESD (0),
fESDfriend (0), fESDfriendExists(kFALSE),
- fRawReader (0),
- fAutoLoad(kFALSE),
- fAutoLoadTime(5.),
- fAutoLoadTimer(0),
- fIsOnline(kFALSE)
+ fAODFile (0), fAODTree (0), fAOD (0),
+ fRawReader (0), fEventInfo(),
+ fAutoLoad (kFALSE), fAutoLoadTime (5.), fAutoLoadTimer(0),
+ fIsOpen (kFALSE), fHasEvent (kFALSE), fExternalCtrl (kFALSE),
+ fGlobal (0), fGlobalReplace (kTRUE), fGlobalUpdate (kTRUE),
+ fExecutor (0), fTransients(0), fTransientLists(0),
+ fPEventSelector(0),
+ fSubManagers (0),
+ fAutoLoadTimerRunning(kFALSE)
{
// Default constructor.
+
+ InitInternals();
}
-AliEveEventManager::AliEveEventManager(TString path, Int_t ev) :
- TEveEventManager("AliEVE AliEveEventManager"),
+AliEveEventManager::AliEveEventManager(const TString& name, const TString& path, Int_t ev) :
+ TEveEventManager(name, path),
fPath (path), fEventId(-1),
fRunLoader (0),
fESDFile (0), fESDTree (0), fESD (0),
fESDfriend (0), fESDfriendExists(kFALSE),
- fRawReader (0),
- fAutoLoad(kFALSE),
- fAutoLoadTime(5.),
- fAutoLoadTimer(0),
- fIsOnline(kFALSE)
+ fAODFile (0), fAODTree (0), fAOD (0),
+ fRawReader (0), fEventInfo(),
+ fAutoLoad (kFALSE), fAutoLoadTime (5), fAutoLoadTimer(0),
+ fIsOpen (kFALSE), fHasEvent (kFALSE), fExternalCtrl (kFALSE),
+ fGlobal (0), fGlobalReplace (kTRUE), fGlobalUpdate (kTRUE),
+ fExecutor (0), fTransients(0), fTransientLists(0),
+ fPEventSelector(0),
+ fSubManagers (0),
+ fAutoLoadTimerRunning(kFALSE)
{
// Constructor with event-directory URL and event-id.
+ InitInternals();
+
Open();
- if (ev >= 0) GotoEvent(ev);
+ if (ev >= 0)
+ {
+ GotoEvent(ev);
+ }
}
AliEveEventManager::~AliEveEventManager()
{
// Destructor.
+ fAutoLoadTimer->Stop();
+ fAutoLoadTimer->Disconnect("Timeout");
+
+ delete fSubManagers;
+
+ if (fIsOpen)
+ {
+ Close();
+ }
- if (fAutoLoadTimer) delete fAutoLoadTimer;
- // Somewhat unclear what to do here.
- // In principle should close all data sources and deregister from
- // TEveManager.
+ fTransients->DecDenyDestroy();
+ fTransients->Destroy();
+
+ fTransientLists->DecDenyDestroy();
+ fTransientLists->Destroy();
+
+ //delete fExecutor;
}
/******************************************************************************/
-void AliEveEventManager::SetESDFileName(const Text_t* esd)
+void AliEveEventManager::SetESDFileName(const TString& esd)
{
// Set file-name for opening ESD, default "AliESDs.root".
- if (esd) fgESDFileName = esd;
+ if ( ! esd.IsNull()) fgESDFileName = esd;
+}
+
+void AliEveEventManager::SetAODFileName(const TString& aod)
+{
+ // Set file-name for opening AOD, default "AliAOD.root".
+
+ if ( ! aod.IsNull()) fgAODFileName = aod;
+}
+
+void AliEveEventManager::AddAODfriend(const TString& friendFileName)
+{
+ // Add new AOD friend file-name to be attached when opening AOD.
+ // This should include '.root', as in 'AliAOD.VertexingHF.root'.
+
+ if (fgAODfriends == 0)
+ {
+ fgAODfriends = new TList;
+ fgAODfriends->SetOwner(kTRUE);
+ }
+ if (fgAODfriends->FindObject(friendFileName) == 0)
+ {
+ fgAODfriends->Add(new TObjString(friendFileName));
+ }
}
-void AliEveEventManager::SetRawFileName(const Text_t* raw)
+void AliEveEventManager::SetRawFileName(const TString& raw)
{
// Set file-name for opening of raw-data, default "raw.root"
- if (raw) fgRawFileName = raw;
+ if ( ! raw.IsNull()) fgRawFileName = raw;
}
-void AliEveEventManager::SetCdbUri(const Text_t* cdb)
+void AliEveEventManager::SetCdbUri(const TString& cdb)
{
- // Set path to CDB, default "local://$ALICE_ROOT".
+ // Set path to CDB, there is no default.
- if (cdb) fgCdbUri = cdb;
+ if ( ! cdb.IsNull()) fgCdbUri = cdb;
}
-void AliEveEventManager::SetAssertElements(Bool_t assertRunloader,
- Bool_t assertEsd,
- Bool_t assertRaw)
+void AliEveEventManager::SetAssertElements(Bool_t assertRunloader, Bool_t assertEsd,
+ Bool_t assertAod, Bool_t assertRaw)
{
// Set global flags that detrmine which parts of the event-data must
// be present when the event is opened.
fgAssertRunLoader = assertRunloader;
fgAssertESD = assertEsd;
+ fgAssertAOD = assertAod;
fgAssertRaw = assertRaw;
}
+void AliEveEventManager::SearchRawForCentralReconstruction()
+{
+ // Enable searching of raw data in standard location. The path passed to
+ // Open() is expected to point to a centrally reconstructed run, e.g.:
+ // "alien:///alice/data/2009/LHC09c/000101134/ESDs/pass1/09000101134018.10".
+
+ fgRawFromStandardLoc = kTRUE;
+}
+
/******************************************************************************/
void AliEveEventManager::Open()
static const TEveException kEH("AliEveEventManager::Open ");
+ if (fExternalCtrl)
+ {
+ throw (kEH + "Event-loop is under external control.");
+ }
+ if (fIsOpen)
+ {
+ throw (kEH + "Event-files already opened.");
+ }
+
gSystem->ExpandPathName(fPath);
- // The following magick is required for ESDriends to be loaded properly
+ // The following magick is required for ESDfriends to be loaded properly
// from non-current directory.
if (fPath.IsNull() || fPath == ".")
{
Int_t runNo = -1;
+ // Open ESD and ESDfriends
+
+ TString esdPath(Form("%s/%s", fPath.Data(), fgESDFileName.Data()));
+ if (fPath.EndsWith(".zip")) esdPath.Form("%s#%s",fPath.Data(),fgESDFileName.Data());
+ if ((fESDFile = TFile::Open(esdPath)))
+ {
+ fESD = new AliESDEvent();
+ fESDTree = (TTree*) fESDFile->Get("esdTree");
+ if (fESDTree != 0)
+ {
+ // Check if ESDfriends exists and attach the branch.
+ // We use TFile::Open() instead of gSystem->AccessPathName
+ // as it seems to work better when attachine alieve to a
+ // running reconstruction process with auto-save on.
+ TString p(Form("%s/AliESDfriends.root", fPath.Data()));
+ if (fPath.EndsWith(".zip")) p.Form("%s#AliESDfriends.root",fPath.Data());
+ TFile *esdFriendFile = TFile::Open(p);
+ if (esdFriendFile)
+ {
+ if (!esdFriendFile->IsZombie())
+ {
+ esdFriendFile->Close();
+ fESDfriendExists = kTRUE;
+ fESDTree->SetBranchStatus ("ESDfriend*", 1);
+ }
+ delete esdFriendFile;
+ }
+
+ fESD->ReadFromTree(fESDTree);
+ if (fESDfriendExists)
+ {
+ fESDfriend = (AliESDfriend*) fESD->FindListObject("AliESDfriend");
+ Info(kEH, "found and attached ESD friend.");
+ }
+ else
+ {
+ Warning(kEH, "ESDfriend not found.");
+ }
+
+ if (fESDTree->GetEntry(0) <= 0)
+ {
+ delete fESDFile; fESDFile = 0;
+ delete fESD; fESD = 0;
+ Warning(kEH, "failed getting the first entry from esdTree.");
+ }
+ else
+ {
+ if (runNo < 0)
+ runNo = fESD->GetESDRun()->GetRunNumber();
+ }
+ }
+ else // esdtree == 0
+ {
+ delete fESDFile; fESDFile = 0;
+ delete fESD; fESD = 0;
+ Warning(kEH, "failed getting the esdTree.");
+ }
+ }
+ else // esd not readable
+ {
+ Warning(kEH, "can not read ESD file '%s'.", esdPath.Data());
+ }
+ if (fESDTree == 0)
+ {
+ if (fgAssertESD)
+ {
+ throw (kEH + "ESD not initialized. Its precence was requested.");
+ } else {
+ Warning(kEH, "ESD not initialized.");
+ }
+ }
+
+ // Open AOD and registered friends
+
+ TString aodPath(Form("%s/%s", fPath.Data(), fgAODFileName.Data()));
+ if (fPath.EndsWith(".zip")) aodPath.Form("%s#%s",fPath.Data(),fgAODFileName.Data());
+ if ((fAODFile = TFile::Open(aodPath)))
+ {
+ fAOD = new AliAODEvent();
+ fAODTree = (TTree*) fAODFile->Get("aodTree");
+ if (fAODTree != 0)
+ {
+ // Check if AODfriends exist and attach them.
+ TIter friends(fgAODfriends);
+ TObjString *name;
+ while ((name = (TObjString*) friends()) != 0)
+ {
+ TString p(Form("%s/%s", fPath.Data(), name->GetName()));
+ if (fPath.EndsWith(".zip")) p.Form("%s#%s",fPath.Data(),name->GetName());
+ if (gSystem->AccessPathName(p, kReadPermission) == kFALSE)
+ {
+ fAODTree->AddFriend("aodTree", name->GetName());
+ }
+ }
+
+ fAOD->ReadFromTree(fAODTree);
+
+ if (fAODTree->GetEntry(0) <= 0)
+ {
+ delete fAODFile; fAODFile = 0;
+ delete fAOD; fAOD = 0;
+ Warning(kEH, "failed getting the first entry from addTree.");
+ }
+ else
+ {
+ if (runNo < 0)
+ runNo = fAOD->GetRunNumber();
+ }
+ }
+ else // aodtree == 0
+ {
+ delete fAODFile; fAODFile = 0;
+ delete fAOD; fAOD = 0;
+ Warning(kEH, "failed getting the aodTree.");
+ }
+ }
+ else // aod not readable
+ {
+ Warning(kEH, "can not read AOD file '%s'.", aodPath.Data());
+ }
+ if (fAODTree == 0)
+ {
+ if (fgAssertAOD)
+ {
+ throw (kEH + "AOD not initialized. Its precence was requested.");
+ } else {
+ Warning(kEH, "AOD not initialized.");
+ }
+ }
+
+ // Open RunLoader from galice.root
+
TString gaPath(Form("%s/galice.root", fPath.Data()));
+ if (fPath.EndsWith(".zip")) gaPath.Form("%s#%s",fPath.Data(),"galice.root");
// If i use open directly, we get fatal.
- // Is this (AccessPathName check) ok for xrootd / alien?
- if (gSystem->AccessPathName(gaPath, kReadPermission) == kFALSE)
+ // Is AccessPathName check ok for xrootd / alien? Yes, not for http.
+ // Seems not to work for alien anymore.
+ // Fixed in ROOT on 27.10.2009, rev 30888.
+ // To revert after we move to root-5.26.
+ TFile *gafile = TFile::Open(gaPath);
+ if (gafile)
{
- fRunLoader = AliRunLoader::Open(gaPath);
+ gafile->Close();
+ delete gafile;
+ // if (gSystem->AccessPathName(gaPath, kReadPermission) == kFALSE)
+ // {
+ fRunLoader = AliRunLoader::Open(gaPath, GetName());
if (fRunLoader)
{
TString alicePath = fPath + "/";
if (fRunLoader->LoadHeader() == 0)
{
- runNo = fRunLoader->GetHeader()->GetRun();
+ if (runNo < 0)
+ runNo = fRunLoader->GetHeader()->GetRun();
}
else
{
Warning(kEH, "Bootstraping of run-loader failed.");
}
+ // Open raw-data file
- TString esdPath(Form("%s/%s", fPath.Data(), fgESDFileName.Data()));
- if ((fESDFile = TFile::Open(esdPath)))
- {
- fESD = new AliESDEvent();
- fESDTree = (TTree*) fESDFile->Get("esdTree");
- if (fESDTree != 0)
- {
- fESD->ReadFromTree(fESDTree);
- fESDTree->GetEntry(0);
- runNo = fESD->GetESDRun()->GetRunNumber();
-
- // Check if ESDfriends exists and attach the branch
- TString p(Form("%s/AliESDfriends.root", fPath.Data()));
- if (gSystem->AccessPathName(p, kReadPermission) == kFALSE)
- {
- fESDfriendExists = kTRUE;
- fESDTree->SetBranchStatus ("ESDfriend*", 1);
- fESDTree->SetBranchAddress("ESDfriend.", &fESDfriend);
- }
- }
- else // esdtree == 0
- {
- delete fESDFile; fESDFile = 0;
- Warning(kEH, "failed getting the esdTree.");
- }
- }
- else // esd not readable
+ TString rawPath;
+ if (fgRawFromStandardLoc)
{
- Warning(kEH, "can not read ESD file '%s'.", esdPath.Data());
+ if (!fPath.BeginsWith("alien:"))
+ throw kEH + "Standard raw search requested, but the directory is not in AliEn.";
+ if (!fPath.Contains("/ESDs/"))
+ throw kEH + "Standard raw search requested, but does not contain 'ESDs' directory.";
+
+ TPMERegexp chunk("/([\\d\\.])+/?$");
+ Int_t nm = chunk.Match(fPath);
+ if (nm != 2)
+ throw kEH + "Standard raw search requested, but the path does not end with chunk-id directory.";
+
+ TPMERegexp esdstrip("/ESDs/.*");
+ rawPath = fPath;
+ esdstrip.Substitute(rawPath, "/raw/");
+ rawPath += chunk[0];
+ rawPath += ".root";
+
+ Info(kEH, "Standard raw search requested, using the following path:\n %s\n", rawPath.Data());
}
- if (fESDTree == 0)
+ else
{
- if (fgAssertESD)
- {
- throw (kEH + "ESD not initialized. Its precence was requested.");
- } else {
- Warning(kEH, "ESD not initialized.");
- }
+ rawPath.Form("%s/%s", fPath.Data(), fgRawFileName.Data());
}
-
- TString rawPath(Form("%s/%s", fPath.Data(), fgRawFileName.Data()));
// If i use open directly, raw-reader reports an error but i have
// no way to detect it.
- // Is this (AccessPathName check) ok for xrootd / alien?
+ // Is this (AccessPathName check) ok for xrootd / alien? Yes, not for http.
+ AliLog::EType_t oldLogLevel = (AliLog::EType_t) AliLog::GetGlobalLogLevel();
+ if (fgAssertRaw == kFALSE)
+ {
+ AliLog::SetGlobalLogLevel(AliLog::kFatal);
+ }
if (gSystem->AccessPathName(rawPath, kReadPermission) == kFALSE)
{
- if (fgRawFileName.EndsWith("/"))
- {
- fRawReader = new AliRawReaderFile(rawPath);
- }
- else if (fgRawFileName.EndsWith(".root"))
- {
- fRawReader = new AliRawReaderRoot(rawPath);
- }
- else if (!fgRawFileName.IsNull())
- {
- fRawReader = new AliRawReaderDate(rawPath);
- }
+ fRawReader = AliRawReader::Create(rawPath);
+ }
+ else
+ {
+ fRawReader = AliRawReader::Create(fgRawFileName);
+ }
+ if (fgAssertRaw == kFALSE)
+ {
+ AliLog::SetGlobalLogLevel(oldLogLevel);
}
if (fRawReader == 0)
if (fgAssertRaw)
{
throw (kEH + "raw-data not initialized. Its precence was requested.");
- } else {
+ }
+ else
+ {
Warning(kEH, "raw-data not initialized.");
}
}
{
if (fRawReader)
{
- fRawReader->NextEvent();
+ if ( ! fRawReader->NextEvent())
+ {
+ throw (kEH + "can not go to first event in raw-reader to determine run-id.");
+ }
runNo = fRawReader->GetRunNumber();
- printf("Determining run-no from raw ... run=%d\n", runNo);
+ Info(kEH, "Determining run-no from raw ... run=%d.", runNo);
fRawReader->RewindEvents();
- } else {
+ }
+ else
+ {
throw (kEH + "unknown run number.");
}
}
+ // Initialize OCDB ... only in master event-manager
+
+ if (this == fgMaster)
{
AliCDBManager* cdb = AliCDBManager::Instance();
- cdb->SetDefaultStorage(fgCdbUri);
- if (cdb->IsDefaultStorageSet() == kFALSE)
- throw (kEH + "CDB initialization failed.");
- cdb->SetRun(runNo);
+ if (cdb->IsDefaultStorageSet() == kTRUE)
+ {
+ Warning(kEH, "CDB already set - using the old storage:\n '%s'",
+ cdb->GetDefaultStorage()->GetURI().Data());
+ }
+ else
+ {
+ if (fgCdbUri.IsNull())
+ {
+ gEnv->SetValue("Root.Stacktrace", "no");
+ Fatal("Open()", "OCDB path was not specified.");
+ }
+
+ // Handle some special cases for MC (should be in OCDBManager).
+ if (fgCdbUri == "mcideal://")
+ cdb->SetDefaultStorage("MC", "Ideal");
+ else if (fgCdbUri == "mcresidual://")
+ cdb->SetDefaultStorage("MC", "Residual");
+ else if (fgCdbUri == "mcfull://")
+ cdb->SetDefaultStorage("MC", "Full");
+ else if (fgCdbUri == "local://") {
+ fgCdbUri = "local://$ALICE_ROOT/OCDB";
+ cdb->SetDefaultStorage(fgCdbUri);
+ } else
+ cdb->SetDefaultStorage(fgCdbUri);
+
+ cdb->SetRun(runNo);
+
+ if (cdb->IsDefaultStorageSet() == kFALSE)
+ throw kEH + "CDB initialization failed for '" + fgCdbUri + "'.";
+ }
+
+ if (fgCdbUri.BeginsWith("local://"))
+ {
+ TString grp = "GRP/GRP/Data";
+ TString grppath = fPath + "/" + grp;
+ if (gSystem->AccessPathName(grppath, kReadPermission) == kFALSE)
+ {
+ if (cdb->GetSpecificStorage(grp))
+ {
+ Warning(kEH, "Local GRP exists, but the specific storage is already set.");
+ }
+ else
+ {
+ Info(kEH, "Setting CDB specific-storage for GRP from event directory.");
+ TString lpath("local://");
+ lpath += fPath;
+ cdb->SetSpecificStorage(grp, lpath);
+ }
+ }
+ }
}
- SetName(Form("Event %d", fEventId));
- SetTitle(fPath);
+ fIsOpen = kTRUE;
}
-void AliEveEventManager::SetEvent(AliRunLoader *runLoader, AliRawReader *rawReader, AliESDEvent *esd)
+void AliEveEventManager::SetEvent(AliRunLoader *runLoader, AliRawReader *rawReader, AliESDEvent *esd, AliESDfriend *esdf)
{
- // Set an event from an external source
- // The method is used in the online visualisation
+ // Set an event from an external source.
+ // The method is used in the online visualisation.
+ // AOD is not supported.
+
+ static const TEveException kEH("AliEveEventManager::SetEvent ");
+
+ if (fIsOpen)
+ {
+ Warning(kEH, "Event-files were open. Closing and switching to external control.");
+ Close();
+ }
+
fRunLoader = runLoader;
fRawReader = rawReader;
- fESD = esd;
- fIsOnline = kTRUE;
- SetTitle("Online event in memory");
- SetName("Online Event");
+ fESD = esd;
+ fESDfriend = esdf;
+ fAOD = 0;
+ fEventId++;
+ fHasEvent = kTRUE;
+ fExternalCtrl = kTRUE;
+
+ SetTitle("Online event in memory");
+ SetName ("Online Event");
ElementChanged();
+
AfterNewEventLoaded();
+
+ if (fAutoLoad) StartAutoLoadTimer();
+}
+
+Int_t AliEveEventManager::GetMaxEventId(Bool_t refreshESD) const
+{
+ // Returns maximum available event id.
+ // If under external control or event is not opened -1 is returned.
+ // If raw-data is the only data-source this can not be known
+ // and 10,000,000 is returned.
+ // If neither data-source is initialised an exception is thrown.
+ // If refresh_esd is true and ESD is the primary event-data source
+ // its header is re-read from disk.
+
+ static const TEveException kEH("AliEveEventManager::GetMaxEventId ");
+
+ if (fExternalCtrl || fIsOpen == kFALSE)
+ {
+ return -1;
+ }
+
+ if (fESDTree)
+ {
+ if (refreshESD)
+ {
+ fESDTree->Refresh();
+ fPEventSelector->Update();
+ }
+ return fESDTree->GetEntries() - 1;
+ }
+ else if (fAODTree)
+ {
+ return fAODTree->GetEntries() - 1;
+ }
+ else if (fRunLoader)
+ {
+ return fRunLoader->GetNumberOfEvents() - 1;
+ }
+ else if (fRawReader)
+ {
+ Int_t n = fRawReader->GetNumberOfEvents() - 1;
+ return n > -1 ? n : 10000000;
+ }
+ else
+ {
+ throw (kEH + "neither ESD, AOD, RunLoader nor Raw loaded.");
+ }
}
void AliEveEventManager::GotoEvent(Int_t event)
// After successful loading of event, the virtual function
// AfterNewEventLoaded() is called. This executes commands that
// were registered via TEveEventManager::AddNewEventCommand().
+ //
+ // If event is negative, it is subtracted from the number of
+ // available events, thus passing -1 will load the last event.
+ // This is not supported when raw-data is the only data-source
+ // as the number of events is not known.
static const TEveException kEH("AliEveEventManager::GotoEvent ");
- if (event < 0) {
- Error(kEH, "event must be non-negative.");
- return;
+ if (fAutoLoadTimerRunning)
+ {
+ throw (kEH + "Event auto-load timer is running.");
+ }
+ if (fExternalCtrl)
+ {
+ throw (kEH + "Event-loop is under external control.");
+ }
+ else if (!fIsOpen)
+ {
+ throw (kEH + "Event-files not opened.");
}
+ fEventInfo.Reset();
+
+ fHasEvent = kFALSE;
+
Int_t maxEvent = 0;
- if (fRunLoader) {
- maxEvent = fRunLoader->GetNumberOfEvents() - 1;
- } else if (fESDTree) {
+ if (fESDTree)
+ {
+ if (event >= fESDTree->GetEntries())
+ fESDTree->Refresh();
maxEvent = fESDTree->GetEntries() - 1;
- } else if (fRawReader) {
- maxEvent = 10000000;
- Info(kEH, "number of events unknown for raw-data, setting max-event id to 10M.");
- } else {
+ if (event < 0)
+ event = fESDTree->GetEntries() + event;
+ }
+ else if (fAODTree)
+ {
+ maxEvent = fAODTree->GetEntries() - 1;
+ if (event < 0)
+ event = fAODTree->GetEntries() + event;
+ }
+ else if (fRunLoader)
+ {
+ maxEvent = fRunLoader->GetNumberOfEvents() - 1;
+ if (event < 0)
+ event = fRunLoader->GetNumberOfEvents() + event;
+ }
+ else if (fRawReader)
+ {
+ maxEvent = fRawReader->GetNumberOfEvents() - 1;
+ if (maxEvent < 0)
+ {
+ maxEvent = 10000000;
+ if (event < 0) {
+ Error(kEH, "current raw-data source does not support direct event access.");
+ return;
+ }
+ Info(kEH, "number of events unknown for current raw-data source, setting max-event id to 10M.");
+ }
+ else
+ {
+ if (event < 0)
+ event = fRawReader->GetNumberOfEvents() + event;
+ }
+ }
+ else
+ {
throw (kEH + "neither RunLoader, ESD nor Raw loaded.");
}
if (event < 0 || event > maxEvent)
+ {
throw (kEH + Form("event %d not present, available range [%d, %d].",
event, 0, maxEvent));
+ }
+
+ TString sysInfoHeader;
+ sysInfoHeader.Form("AliEveEventManager::GotoEvent(%d) - ", event);
+ AliSysInfo::AddStamp(sysInfoHeader + "Start");
TEveManager::TRedrawDisabler rd(gEve);
gEve->Redraw3D(kFALSE, kTRUE); // Enforce drop of all logicals.
// !!! MT this is somewhat brutal; at least optionally, one could be
// a bit gentler, checking for objs owning their external refs and having
// additinal parents.
+ gEve->GetViewers()->DeleteAnnotations();
+ fTransients->DestroyElements();
+ for (TEveElement::List_i i = fTransientLists->BeginChildren();
+ i != fTransientLists->EndChildren(); ++i)
+ {
+ (*i)->DestroyElements();
+ }
DestroyElements();
- if (fRunLoader) {
- if (fRunLoader->GetEvent(event) != 0)
- throw (kEH + "failed getting required event.");
- }
+ AliSysInfo::AddStamp(sysInfoHeader + "PostDestroy");
if (fESDTree) {
if (fESDTree->GetEntry(event) <= 0)
fESD->SetESDfriend(fESDfriend);
}
- if (fRawReader)
- {
- Int_t rawEv = fEventId;
- if (event < rawEv)
- {
- fRawReader->RewindEvents();
- rawEv = -1;
- }
+ if (fAODTree) {
+ if (fAODTree->GetEntry(event) <= 0)
+ throw (kEH + "failed getting required event from AOD.");
+ }
- while (rawEv < event)
+ if (fRunLoader) {
+ if (fRunLoader->GetEvent(event) != 0)
+ throw (kEH + "failed getting required event.");
+ }
+
+ if (fRawReader)
+ {
+ // AliRawReader::GotoEvent(Int_t) works for AliRawReaderRoot/Chain.
+ if (fRawReader->GotoEvent(event) == kFALSE)
{
- if ( ! fRawReader->NextEvent())
+ // Use fallback method - iteration with NextEvent().
+ Int_t rawEv = fEventId;
+ if (event < rawEv)
{
fRawReader->RewindEvents();
- fEventId = -1;
- throw (kEH + Form("Error going to next raw-event from event %d.", rawEv));
+ rawEv = -1;
+ }
+
+ while (rawEv < event)
+ {
+ if ( ! fRawReader->NextEvent())
+ {
+ fRawReader->RewindEvents();
+ fEventId = -1;
+ throw (kEH + Form("Error going to next raw-event from event %d.", rawEv));
+ }
+ ++rawEv;
}
- ++rawEv;
+ Warning(kEH, "Loaded raw-event %d with fallback method.\n", rawEv);
}
+ }
- printf ("Loaded raw-event %d.\n", rawEv);
+ fHasEvent = kTRUE;
+ fEventId = event;
+ if (this == fgMaster)
+ {
+ SetName(Form("Event %d", fEventId));
+ ElementChanged();
}
- fEventId = event;
- SetName(Form("Event %d", fEventId));
- ElementChanged();
+ AliSysInfo::AddStamp(sysInfoHeader + "PostLoadEvent");
AfterNewEventLoaded();
+
+ AliSysInfo::AddStamp(sysInfoHeader + "PostUserActions");
+}
+
+void AliEveEventManager::Timeout()
+{
+ Emit("Timeout()");
+}
+
+void AliEveEventManager::NextEvent()
+{
+ // Loads next event.
+ // Does magick needed for online display when under external event control.
+
+ static const TEveException kEH("AliEveEventManager::NextEvent ");
+
+ if (fAutoLoadTimerRunning)
+ {
+ throw (kEH + "Event auto-load timer is running.");
+ }
+
+ if (fExternalCtrl)
+ {
+ // !!! This should really go somewhere else. It is done in GotoEvent(),
+ // so here we should do it in SetEvent().
+ DestroyElements();
+ gSystem->ExitLoop();
+
+ }
+ else if (fESDTree)
+ {
+ Int_t nextevent=0;
+ if (fPEventSelector->FindNext(nextevent))
+ {
+ GotoEvent(nextevent);
+ }
+ }
+ else if (fEventId < GetMaxEventId(kTRUE))
+ {
+ GotoEvent(fEventId + 1);
+ }
+}
+
+void AliEveEventManager::PrevEvent()
+{
+ // Loads previous event.
+
+ static const TEveException kEH("AliEveEventManager::PrevEvent ");
+
+ if (fAutoLoadTimerRunning)
+ {
+ throw (kEH + "Event auto-load timer is running.");
+ }
+ if (fExternalCtrl)
+ {
+ throw (kEH + "Event-loop is under external control.");
+ }
+
+ if (fESDTree)
+ {
+ Int_t nextevent=0;
+ if (fPEventSelector->FindPrev(nextevent))
+ {
+ GotoEvent(nextevent);
+ }
+ }
+ else if (fEventId > 0)
+ {
+ GotoEvent(fEventId - 1);
+ }
}
void AliEveEventManager::Close()
{
- // Close the event files.
- // For the moment only ESD is closed. Needs to be investigated for
- // AliRunLoader and Raw.
+ // Close the event data-files and delete ESD, ESDfriend, run-loader
+ // and raw-reader.
+
+ static const TEveException kEH("AliEveEventManager::Close ");
+
+ if (!fIsOpen)
+ {
+ throw (kEH + "Event-files not opened.");
+ }
+
+ if (fAutoLoadTimerRunning)
+ StopAutoLoadTimer();
if (fESDTree) {
delete fESD; fESD = 0;
- delete fESDfriend; fESDfriend = 0;
+ // delete fESDfriend; // friend tree is deleted with the tree
+ fESDfriend = 0;
+ fESDfriendExists = kFALSE;
+
+ delete fESDTree; fESDTree = 0;
+ delete fESDFile; fESDFile = 0;
+ }
+
+ if (fAODTree) {
+ delete fAOD; fAOD = 0;
- delete fESDTree; fESDTree = 0;
- delete fESDFile; fESDFile = 0;
+ delete fAODTree; fAODTree = 0;
+ delete fAODFile; fAODFile = 0;
}
+
+ if (fRunLoader) {
+ delete fRunLoader; fRunLoader = 0;
+ }
+
+ if (fRawReader) {
+ delete fRawReader; fRawReader = 0;
+ }
+
+ fEventId = -1;
+ fIsOpen = kFALSE;
+ fHasEvent = kFALSE;
}
-/******************************************************************************/
+//------------------------------------------------------------------------------
// Static convenience functions, mainly used from macros.
-/******************************************************************************/
+//------------------------------------------------------------------------------
+
+Int_t AliEveEventManager::CurrentEventId()
+{
+ // Return current event-id.
+
+ static const TEveException kEH("AliEveEventManager::CurrentEventId ");
+
+ if (fgCurrent == 0 || fgCurrent->fHasEvent == kFALSE)
+ throw (kEH + "ALICE event not ready.");
+ return fgCurrent->GetEventId();
+}
+
+Bool_t AliEveEventManager::HasRunLoader()
+{
+ // Check if AliRunLoader is initialized.
+
+ return fgCurrent && fgCurrent->fHasEvent && fgCurrent->fRunLoader;
+}
+
+Bool_t AliEveEventManager::HasESD()
+{
+ // Check if AliESDEvent is initialized.
+
+ return fgCurrent && fgCurrent->fHasEvent && fgCurrent->fESD;
+}
+
+Bool_t AliEveEventManager::HasESDfriend()
+{
+ // Check if AliESDfriend is initialized.
+
+ return fgCurrent && fgCurrent->fHasEvent && fgCurrent->fESDfriend;
+}
+
+Bool_t AliEveEventManager::HasAOD()
+{
+ // Check if AliESDEvent is initialized.
+
+ return fgCurrent && fgCurrent->fHasEvent && fgCurrent->fAOD;
+}
+
+Bool_t AliEveEventManager::HasRawReader()
+{
+ // Check if raw-reader is initialized.
+
+ return fgCurrent && fgCurrent->fHasEvent && fgCurrent->fRawReader;
+}
AliRunLoader* AliEveEventManager::AssertRunLoader()
{
static const TEveException kEH("AliEveEventManager::AssertRunLoader ");
- if (gAliEveEvent == 0)
+ if (fgCurrent == 0 || fgCurrent->fHasEvent == kFALSE)
throw (kEH + "ALICE event not ready.");
- if (gAliEveEvent->fRunLoader == 0)
+ if (fgCurrent->fRunLoader == 0)
throw (kEH + "AliRunLoader not initialised.");
- return gAliEveEvent->fRunLoader;
+ return fgCurrent->fRunLoader;
}
AliESDEvent* AliEveEventManager::AssertESD()
static const TEveException kEH("AliEveEventManager::AssertESD ");
- if (gAliEveEvent == 0)
+ if (fgCurrent == 0 || fgCurrent->fHasEvent == kFALSE)
throw (kEH + "ALICE event not ready.");
- if (gAliEveEvent->fESD == 0)
+ if (fgCurrent->fESD == 0)
throw (kEH + "AliESD not initialised.");
- return gAliEveEvent->fESD;
+ return fgCurrent->fESD;
}
AliESDfriend* AliEveEventManager::AssertESDfriend()
static const TEveException kEH("AliEveEventManager::AssertESDfriend ");
- if (gAliEveEvent == 0)
+ if (fgCurrent == 0 || fgCurrent->fHasEvent == kFALSE)
throw (kEH + "ALICE event not ready.");
- if (gAliEveEvent->fESDfriend == 0)
+ if (fgCurrent->fESDfriend == 0)
throw (kEH + "AliESDfriend not initialised.");
- return gAliEveEvent->fESDfriend;
+ return fgCurrent->fESDfriend;
+}
+
+AliAODEvent* AliEveEventManager::AssertAOD()
+{
+ // Make sure AliAODEvent is initialized and return it.
+ // Throws exception in case AOD is not available.
+ // Static utility for macros.
+
+ static const TEveException kEH("AliEveEventManager::AssertAOD ");
+
+ if (fgCurrent == 0 || fgCurrent->fHasEvent == kFALSE)
+ throw (kEH + "ALICE event not ready.");
+ if (fgCurrent->fAOD == 0)
+ throw (kEH + "AliAOD not initialised.");
+ return fgCurrent->fAOD;
}
AliRawReader* AliEveEventManager::AssertRawReader()
static const TEveException kEH("AliEveEventManager::AssertRawReader ");
- if (gAliEveEvent == 0)
+ if (fgCurrent == 0 || fgCurrent->fHasEvent == kFALSE)
throw (kEH + "ALICE event not ready.");
- if (gAliEveEvent->fRawReader == 0)
+ if (fgCurrent->fRawReader == 0)
throw (kEH + "RawReader not ready.");
- return gAliEveEvent->fRawReader;
+ return fgCurrent->fRawReader;
}
-AliMagF* AliEveEventManager::AssertMagField()
-{
- // Make sure AliMagF is initialized and return it.
- // Throws exception in case magnetic field is not available.
- // Static utility for macros.
+//==============================================================================
+
+AliMagF* AliEveEventManager::AssertMagField()
+{
+ // Make sure AliMagF is initialized and returns it.
+ // Throws exception in case magnetic field is not available.
+ // Static utility for macros.
+
+ static const TEveException kEH("AliEveEventManager::AssertMagField ");
+
+ if (fgMagField)
+ return fgMagField;
- if (fgMagField == 0)
+ if (TGeoGlobalMagField::Instance()->GetField())
{
- if (gAliEveEvent && gAliEveEvent->fRunLoader && gAliEveEvent->fRunLoader->GetAliRun())
- fgMagField = gAliEveEvent->fRunLoader->GetAliRun()->Field();
- else
- fgMagField = new AliMagFMaps("Maps","Maps", 1, 1., 10., AliMagFMaps::k5kG);
+ fgMagField = dynamic_cast<AliMagF*>(TGeoGlobalMagField::Instance()->GetField());
+ if (fgMagField == 0)
+ throw kEH + "Global field set, but it is not AliMagF.";
+ return fgMagField;
}
- return fgMagField;
+
+ if (!fgGRPLoaded)
+ {
+ InitGRP();
+ }
+
+ if (TGeoGlobalMagField::Instance()->GetField())
+ {
+ fgMagField = dynamic_cast<AliMagF*>(TGeoGlobalMagField::Instance()->GetField());
+ if (fgMagField == 0)
+ throw kEH + "Global field set, but it is not AliMagF.";
+ }
+ else
+ {
+ throw kEH + "Could not initialize magnetic field.";
+ }
+
+ return fgMagField;
}
TGeoManager* AliEveEventManager::AssertGeometry()
return gGeoManager;
}
-void AliEveEventManager::SetAutoLoad(Bool_t autoLoad)
+AliRecoParam* AliEveEventManager::AssertRecoParams()
{
- // Set the automatic event loading mode
- fAutoLoad = autoLoad;
- StartStopAutoLoadTimer();
+ if(!fgRecoParam)
+ InitRecoParam();
+
+ return fgRecoParam;
+}
+
+Bool_t AliEveEventManager::InitRecoParam()
+{
+// This is mostly a reap-off from reconstruction
+// The method accesses OCDB and retrieves all
+// the available reco-param objects from there.
+
+ fgRecoParam = new AliRecoParam;
+ const Int_t kNDetectors = 14;
+
+ static const TEveException kEH("AliEveEventManager::InitRecoParam ");
+
+ Bool_t isOK = kTRUE;
+
+ if (fgRecoParam->GetDetRecoParamArray(kNDetectors)) {
+ ::Info(kEH, "Using custom GRP reconstruction parameters");
+ }
+ else {
+ ::Info(kEH, "Loading GRP reconstruction parameter objects");
+
+ AliCDBPath path("GRP","Calib","RecoParam");
+ AliCDBEntry *entry=AliCDBManager::Instance()->Get(path.GetPath());
+ if(!entry){
+ ::Warning(kEH, "Couldn't find GRP RecoParam entry in OCDB");
+ isOK = kFALSE;
+ }
+ else {
+ TObject *recoParamObj = entry->GetObject();
+ if (dynamic_cast<TObjArray*>(recoParamObj)) {
+ // GRP has a normal TobjArray of AliDetectorRecoParam objects
+ // Registering them in AliRecoParam
+ fgRecoParam->AddDetRecoParamArray(kNDetectors,dynamic_cast<TObjArray*>(recoParamObj));
+ }
+ else if (dynamic_cast<AliDetectorRecoParam*>(recoParamObj)) {
+ // GRP has only onse set of reco parameters
+ // Registering it in AliRecoParam
+ ::Info(kEH, "Single set of GRP reconstruction parameters found");
+ dynamic_cast<AliDetectorRecoParam*>(recoParamObj)->SetAsDefault();
+ fgRecoParam->AddDetRecoParam(kNDetectors,dynamic_cast<AliDetectorRecoParam*>(recoParamObj));
+ }
+ else {
+ ::Error(kEH, "No valid GRP RecoParam object found in the OCDB");
+ isOK = kFALSE;
+ }
+ entry->SetOwner(0);
+ }
+ }
+
+ const char* fgkDetectorName[kNDetectors] = {"ITS", "TPC", "TRD", "TOF", "PHOS", "HMPID", "EMCAL", "MUON", "FMD", "ZDC", "PMD", "T0", "VZERO", "ACORDE" };
+
+
+ for (Int_t iDet = 0; iDet < kNDetectors; iDet++) {
+
+ if (fgRecoParam->GetDetRecoParamArray(iDet)) {
+ ::Info(kEH, Form("Using custom reconstruction parameters for detector %s",fgkDetectorName[iDet]));
+ continue;
+ }
+
+ ::Info(kEH, Form("Loading reconstruction parameter objects for detector %s",fgkDetectorName[iDet]));
+
+ AliCDBPath path(fgkDetectorName[iDet],"Calib","RecoParam");
+ AliCDBEntry *entry=AliCDBManager::Instance()->Get(path.GetPath());
+ if(!entry){
+ ::Warning(kEH, Form("Couldn't find RecoParam entry in OCDB for detector %s",fgkDetectorName[iDet]));
+ isOK = kFALSE;
+ }
+ else {
+ TObject *recoParamObj = entry->GetObject();
+ if (dynamic_cast<TObjArray*>(recoParamObj)) {
+ // The detector has a normal TobjArray of AliDetectorRecoParam objects
+ // Registering them in AliRecoParam
+ fgRecoParam->AddDetRecoParamArray(iDet,dynamic_cast<TObjArray*>(recoParamObj));
+ }
+ else if (dynamic_cast<AliDetectorRecoParam*>(recoParamObj)) {
+ // The detector has only onse set of reco parameters
+ // Registering it in AliRecoParam
+ ::Info(kEH, Form("Single set of reconstruction parameters found for detector %s",fgkDetectorName[iDet]));
+ dynamic_cast<AliDetectorRecoParam*>(recoParamObj)->SetAsDefault();
+ fgRecoParam->AddDetRecoParam(iDet,dynamic_cast<AliDetectorRecoParam*>(recoParamObj));
+ }
+ else {
+ ::Error(kEH, Form("No valid RecoParam object found in the OCDB for detector %s",fgkDetectorName[iDet]));
+ isOK = kFALSE;
+ }
+ entry->SetOwner(0);
+
+ }
+ }
+
+ if(!isOK) {
+ delete fgRecoParam;
+ fgRecoParam = 0;
+ }
+
+ return isOK;
}
-void AliEveEventManager::SetAutoLoadTime(Double_t time)
+
+//------------------------------------------------------------------------------
+
+AliEveEventManager* AliEveEventManager::AddDependentManager(const TString& name, const TString& path)
+{
+ // Create and attach a dependent event-manager.
+ // It is not added into eve list tree.
+
+ static const TEveException kEH("AliEveEventManager::AddDependentManager ");
+
+ if (fgMaster == 0)
+ throw(kEH + "Master event-manager must be instantiated first.");
+
+ if (fgMaster->fSubManagers == 0)
+ {
+ fgMaster->fSubManagers = new TList;
+ fgMaster->fSubManagers->SetOwner(kTRUE);
+ }
+
+ AliEveEventManager* new_mgr = 0;
+ fgCurrent = 0;
+ try
+ {
+ new_mgr = new AliEveEventManager(name, path, fgMaster->fEventId);
+ fgMaster->fSubManagers->Add(new_mgr);
+ }
+ catch (TEveException& exc)
+ {
+ ::Error(kEH, "Creation of new event-manager failed: '%s'.", exc.Data());
+ }
+ fgCurrent = fgMaster;
+
+ return new_mgr;
+}
+
+AliEveEventManager* AliEveEventManager::GetDependentManager(const TString& name)
+{
+ // Get a dependant manager by name.
+ // This will not change the current manager, use helper class
+ // AliEveEventManager::CurrentChanger for that.
+
+ static const TEveException kEH("AliEveEventManager::GetDependentManager ");
+
+ if (fgMaster == 0)
+ throw(kEH + "Master event-manager must be instantiated first.");
+
+ if (fgMaster->fSubManagers == 0)
+ return 0;
+
+ return dynamic_cast<AliEveEventManager*>(fgMaster->fSubManagers->FindObject(name));
+}
+
+AliEveEventManager* AliEveEventManager::GetMaster()
+{
+ // Get master event-manager.
+
+ return fgMaster;
+}
+
+AliEveEventManager* AliEveEventManager::GetCurrent()
+{
+ // Get current event-manager.
+
+ return fgCurrent;
+}
+
+void AliEveEventManager::RegisterTransient(TEveElement* element)
+{
+ GetCurrent()->fTransients->AddElement(element);
+}
+
+void AliEveEventManager::RegisterTransientList(TEveElement* element)
+{
+ GetCurrent()->fTransientLists->AddElement(element);
+}
+
+//------------------------------------------------------------------------------
+// Autoloading of events
+//------------------------------------------------------------------------------
+
+void AliEveEventManager::SetAutoLoadTime(Float_t time)
{
// Set the auto-load time in seconds
+
fAutoLoadTime = time;
- StartStopAutoLoadTimer();
}
-void AliEveEventManager::StartStopAutoLoadTimer()
+void AliEveEventManager::SetAutoLoad(Bool_t autoLoad)
{
- // Create if needed and start
- // the automatic event loading timer
- if (fAutoLoad) {
- if (!fAutoLoadTimer) {
- fAutoLoadTimer = new TTimer;
- fAutoLoadTimer->Connect("Timeout()","AliEveEventManager",this,"NextEvent()");
- }
- fAutoLoadTimer->Start((Long_t)fAutoLoadTime*1000,kTRUE);
+ // Set the automatic event loading mode
+
+ static const TEveException kEH("AliEveEventManager::SetAutoLoad ");
+
+ if (fAutoLoad == autoLoad)
+ {
+ Warning(kEH, "Setting autoload to the same value as before - %s. Ignoring.", fAutoLoad ? "true" : "false");
+ return;
}
- else {
- if (fAutoLoadTimer) fAutoLoadTimer->Stop();
+
+ fAutoLoad = autoLoad;
+ if (fAutoLoad)
+ {
+ StartAutoLoadTimer();
+ }
+ else
+ {
+ StopAutoLoadTimer();
}
}
-void AliEveEventManager::PrevEvent()
+void AliEveEventManager::SetTrigSel(Int_t trig)
{
- // Loads previous event
- // only in case of manual mode
- if (!fIsOnline) {
- GotoEvent(fEventId - 1);
- StartStopAutoLoadTimer();
+ static const TEveException kEH("AliEveEventManager::SetTrigSel ");
+
+ if (!fRawReader)
+ {
+ Warning(kEH, "No Raw-reader exists. Ignoring the call.");
+ return;
+ }
+ else
+ {
+ ULong64_t trigMask = 0;
+ if (trig >= 0) trigMask = (1ull << trig);
+ Info(kEH,"Trigger selection: 0x%llx",trigMask);
+ fRawReader->SelectEvents(-1,trigMask,NULL);
}
}
-void AliEveEventManager::NextEvent()
+void AliEveEventManager::StartAutoLoadTimer()
{
- // Loads next event
- // either in automatic (online) or
- // manual mode
-
- if (fIsOnline) {
- if (fAutoLoadTimer) fAutoLoadTimer->Stop();
+ // Start the auto-load timer.
- DestroyElements();
+ fAutoLoadTimer->SetTime((Long_t)(1000*fAutoLoadTime));
+ fAutoLoadTimer->Reset();
+ fAutoLoadTimer->TurnOn();
+ fAutoLoadTimerRunning = kTRUE;
+}
+
+void AliEveEventManager::StopAutoLoadTimer()
+{
+ // Stop the auto-load timer.
+
+ fAutoLoadTimerRunning = kFALSE;
+ fAutoLoadTimer->TurnOff();
+}
- gSystem->ExitLoop();
+void AliEveEventManager::AutoLoadNextEvent()
+{
+ // Called from auto-load timer, so it has to be public.
+ // Do NOT call it directly.
+
+ static const TEveException kEH("AliEveEventManager::AutoLoadNextEvent ");
+
+ if ( ! fAutoLoadTimerRunning || ! fAutoLoadTimer->HasTimedOut())
+ {
+ Warning(kEH, "Called unexpectedly - ignoring the call. Should ONLY be called from an internal timer.");
+ return;
}
- else {
- GotoEvent(fEventId + 1);
- StartStopAutoLoadTimer();
+
+ StopAutoLoadTimer();
+ NextEvent();
+ if (fAutoLoad && !fExternalCtrl)
+ StartAutoLoadTimer();
+}
+
+//------------------------------------------------------------------------------
+// Post event-loading functions
+//------------------------------------------------------------------------------
+
+void AliEveEventManager::AfterNewEventLoaded()
+{
+ // Execute registered macros and commands.
+ // At the end emit NewEventLoaded signal.
+ //
+ // Virtual from TEveEventManager.
+
+ static const TEveException kEH("AliEveEventManager::AfterNewEventLoaded ");
+
+ NewEventDataLoaded();
+
+ if (fExecutor)
+ fExecutor->ExecMacros();
+
+ TEveEventManager::AfterNewEventLoaded();
+
+ NewEventLoaded();
+
+ if (this == fgMaster && fSubManagers != 0)
+ {
+ TIter next(fSubManagers);
+ while ((fgCurrent = dynamic_cast<AliEveEventManager*>(next())) != 0)
+ {
+ gEve->SetCurrentEvent(fgCurrent);
+ try
+ {
+ fgCurrent->GotoEvent(fEventId);
+ }
+ catch (TEveException& exc)
+ {
+ // !!! Should somehow tag / disable / remove it?
+ Error(kEH, "Getting event %d for sub-event-manager '%s' failed: '%s'.",
+ fEventId, fgCurrent->GetName(), exc.Data());
+ }
+ }
+ fgCurrent = fgMaster;
+ gEve->SetCurrentEvent(fgMaster);
}
}
-const char* AliEveEventManager::GetEventInfo() const
+void AliEveEventManager::NewEventDataLoaded()
{
- // Dumps the event-header contents
+ // Emit NewEventDataLoaded signal.
- static TString eventInfo;
+ Emit("NewEventDataLoaded()");
+}
- if (!fRawReader) {
- eventInfo.Form("No raw-data event info is available!\n");
+void AliEveEventManager::NewEventLoaded()
+{
+ // Emit NewEventLoaded signal.
+
+ Emit("NewEventLoaded()");
+}
+
+
+//------------------------------------------------------------------------------
+// Event info dumpers
+//------------------------------------------------------------------------------
+
+const AliEventInfo* AliEveEventManager::GetEventInfo()
+{
+ // Fill the event info object
+
+ AliCentralTrigger *aCTP = NULL;
+ if (fRawReader) {
+ fEventInfo.SetEventType(fRawReader->GetType());
+
+ ULong64_t mask = fRawReader->GetClassMask();
+ fEventInfo.SetTriggerMask(mask);
+ UInt_t clmask = fRawReader->GetDetectorPattern()[0];
+ fEventInfo.SetTriggerCluster(AliDAQ::ListOfTriggeredDetectors(clmask));
+
+ aCTP = new AliCentralTrigger();
+ TString configstr("");
+ if (!aCTP->LoadConfiguration(configstr)) { // Load CTP config from OCDB
+ AliError("No trigger configuration found in OCDB! The trigger configuration information will not be used!");
+ delete aCTP;
+ return 0;
+ }
+ aCTP->SetClassMask(mask);
+ aCTP->SetClusterMask(clmask);
+
+ if (fRunLoader) {
+ AliCentralTrigger* rlCTP = fRunLoader->GetTrigger();
+ if (rlCTP) {
+ rlCTP->SetClassMask(mask);
+ rlCTP->SetClusterMask(clmask);
+ }
+ }
}
else {
- const UInt_t* attr = fRawReader->GetAttributes();
- eventInfo.Form("Raw-data event info:\nRun#: %d\nEvent type: %d (%s)\nPeriod: %x\nOrbit: %x BC: %x\nTrigger: %llx\nDetectors: %x (%s)\nAttributes:%x-%x-%x\nTimestamp: %x\n",
- fRawReader->GetRunNumber(),fRawReader->GetType(),AliRawEventHeaderBase::GetTypeName(fRawReader->GetType()),
- fRawReader->GetPeriod(),fRawReader->GetOrbitID(),fRawReader->GetBCID(),
- fRawReader->GetClassMask(),
- *fRawReader->GetDetectorPattern(),AliDAQ::ListOfTriggeredDetectors(*fRawReader->GetDetectorPattern()),
- attr[0],attr[1],attr[2],
- fRawReader->GetTimestamp());
+ fEventInfo.SetEventType(AliRawEventHeaderBase::kPhysicsEvent);
+
+ if (fRunLoader && (!fRunLoader->LoadTrigger())) {
+ aCTP = fRunLoader->GetTrigger();
+ fEventInfo.SetTriggerMask(aCTP->GetClassMask());
+ // get inputs from actp - just get
+ AliESDHeader* esdheader = fESD->GetHeader();
+ esdheader->SetL0TriggerInputs(aCTP->GetL0TriggerInputs());
+ esdheader->SetL1TriggerInputs(aCTP->GetL1TriggerInputs());
+ esdheader->SetL2TriggerInputs(aCTP->GetL2TriggerInputs());
+ fEventInfo.SetTriggerCluster(AliDAQ::ListOfTriggeredDetectors(aCTP->GetClusterMask()));
+ }
+ else {
+ AliWarning("No trigger can be loaded! The trigger information will not be used!");
+ return 0;
+ }
}
- if (!fESD) {
- eventInfo.Append(Form("\nNo ESD event info is available!\n"));
+
+ AliTriggerConfiguration *config = aCTP->GetConfiguration();
+ if (!config) {
+ AliError("No trigger configuration has been found! The trigger configuration information will not be used!");
+ if (fRawReader) delete aCTP;
+ return 0;
+ }
+
+ TString declTriggerClasses;
+
+ // Load trigger aliases and declare the trigger classes included in aliases
+ AliCDBEntry * entry = AliCDBManager::Instance()->Get("GRP/CTP/Aliases");
+ if (entry) {
+ THashList * lst = dynamic_cast<THashList*>(entry->GetObject());
+ if (lst) {
+ lst->Sort(kSortDescending); // to avoid problems with substrings
+ if (fRawReader) fRawReader->LoadTriggerAlias(lst);
+ // Now declare all the triggers present in the aliases
+ TIter iter(lst);
+ TNamed *nmd = 0;
+ while((nmd = dynamic_cast<TNamed*>(iter.Next()))){
+ declTriggerClasses += " ";
+ declTriggerClasses += nmd->GetName();
+ }
+ }
+ else {
+ AliError("Cannot cast the object with trigger aliases to THashList!");
+ }
}
else {
- TString acttrclasses = fESD->GetESDRun()->GetActiveTriggerClasses();
+ AliError("No OCDB entry for the trigger aliases!");
+ }
+
+ // Load trigger classes for this run
+ UChar_t clustmask = 0;
+ TString trclasses;
+ ULong64_t trmask = fEventInfo.GetTriggerMask();
+ const TObjArray& classesArray = config->GetClasses();
+ Int_t nclasses = classesArray.GetEntriesFast();
+ for( Int_t iclass=0; iclass < nclasses; iclass++ ) {
+ AliTriggerClass* trclass = (AliTriggerClass*)classesArray.At(iclass);
+ if (trclass && trclass->GetMask()>0) {
+ Int_t trindex = TMath::Nint(TMath::Log2(trclass->GetMask()));
+ if (fESD) fESD->SetTriggerClass(trclass->GetName(),trindex);
+ if (fRawReader) fRawReader->LoadTriggerClass(trclass->GetName(),trindex);
+ if (trmask & (1ull << trindex)) {
+ trclasses += " ";
+ trclasses += trclass->GetName();
+ trclasses += " ";
+ clustmask |= trclass->GetCluster()->GetClusterMask();
+ }
+ }
+ }
+ fEventInfo.SetTriggerClasses(trclasses);
+
+ if (!aCTP->CheckTriggeredDetectors()) {
+ if (fRawReader) delete aCTP;
+ return 0;
+ }
+
+ if (fRawReader) delete aCTP;
+
+// everything went ok, return pointer
+ return (&fEventInfo);
+}
+
+
+TString AliEveEventManager::GetEventInfoHorizontal() const
+{
+ // Dumps the event-header contents in vertical formatting.
+
+ TString rawInfo, esdInfo;
+
+ if (!fRawReader)
+ {
+ rawInfo = "No raw-data event info is available!\n";
+ }
+ else
+ {
+ const UInt_t* attr = fRawReader->GetAttributes();
+ TTimeStamp ts(fRawReader->GetTimestamp());
+ rawInfo.Form("RAW event info: Run#: %d Event type: %d (%s) Period: %x Orbit: %x BC: %x\n"
+ "Trigger: %llx\nDetectors: %x (%s)\nAttributes:%x-%x-%x Timestamp: %s\n",
+ fRawReader->GetRunNumber(),fRawReader->GetType(),AliRawEventHeaderBase::GetTypeName(fRawReader->GetType()),
+ fRawReader->GetPeriod(),fRawReader->GetOrbitID(),fRawReader->GetBCID(),
+ fRawReader->GetClassMask(),
+ *fRawReader->GetDetectorPattern(),AliDAQ::ListOfTriggeredDetectors(*fRawReader->GetDetectorPattern()),
+ attr[0],attr[1],attr[2], ts.AsString("s"));
+ }
+
+ if (!fESD)
+ {
+ esdInfo = "No ESD event info is available!";
+ }
+ else
+ {
+ TString acttrclasses = fESD->GetESDRun()->GetActiveTriggerClasses();
TString firedtrclasses = fESD->GetFiredTriggerClasses();
- eventInfo.Append(Form("\nESD event info:\nRun#: %d\nActive trigger classes: %s\nEvent type: %d (%s)\nPeriod: %x\nOrbit: %x BC: %x\nTrigger: %llx (%s)\nEvent# in file:%d\nTimestamp: %x\n",
- fESD->GetRunNumber(),
- acttrclasses.Data(),
- fESD->GetEventType(),AliRawEventHeaderBase::GetTypeName(fESD->GetEventType()),
- fESD->GetPeriodNumber(),fESD->GetOrbitNumber(),fESD->GetBunchCrossNumber(),
- fESD->GetTriggerMask(),firedtrclasses.Data(),
- fESD->GetEventNumberInFile(),
- fESD->GetTimeStamp()));
+ TTimeStamp ts(fESD->GetTimeStamp());
+ esdInfo.Form("ESD event info: Run#: %d Event type: %d (%s) Period: %x Orbit: %x BC: %x\n"
+ "Active trigger classes: %s\nTrigger: %llx (%s)\nEvent# in file: %d Timestamp: %s, MagField: %.2e",
+ fESD->GetRunNumber(),
+ fESD->GetEventType(),AliRawEventHeaderBase::GetTypeName(fESD->GetEventType()),
+ fESD->GetPeriodNumber(),fESD->GetOrbitNumber(),fESD->GetBunchCrossNumber(),
+ acttrclasses.Data(),
+ fESD->GetTriggerMask(),firedtrclasses.Data(),
+ fESD->GetEventNumberInFile(), ts.AsString("s"), fESD->GetMagneticField());
}
- return eventInfo.Data();
+ return rawInfo + esdInfo;
+}
+
+TString AliEveEventManager::GetEventInfoVertical() const
+{
+ // Dumps the event-header contents in vertical formatting.
+
+ TString rawInfo, esdInfo;
+
+ if (!fRawReader)
+ {
+ rawInfo = "No raw-data event info is available!\n";
+ }
+ else
+ {
+ const UInt_t* attr = fRawReader->GetAttributes();
+ rawInfo.Form("Raw-data event info:\nRun#: %d\nEvent type: %d (%s)\nPeriod: %x\nOrbit: %x BC: %x\nTrigger: %llx\nDetectors: %x (%s)\nAttributes:%x-%x-%x\nTimestamp: %x\n",
+ fRawReader->GetRunNumber(),fRawReader->GetType(),AliRawEventHeaderBase::GetTypeName(fRawReader->GetType()),
+ fRawReader->GetPeriod(),fRawReader->GetOrbitID(),fRawReader->GetBCID(),
+ fRawReader->GetClassMask(),
+ *fRawReader->GetDetectorPattern(),AliDAQ::ListOfTriggeredDetectors(*fRawReader->GetDetectorPattern()),
+ attr[0],attr[1],attr[2],
+ fRawReader->GetTimestamp());
+ }
+
+ if (!fESD)
+ {
+ esdInfo = "No ESD event info is available!\n";
+ }
+ else
+ {
+ TString acttrclasses = fESD->GetESDRun()->GetActiveTriggerClasses();
+ TString firedtrclasses = fESD->GetFiredTriggerClasses();
+ esdInfo.Form("ESD event info:\nRun#: %d\nActive trigger classes: %s\nEvent type: %d (%s)\nPeriod: %x\nOrbit: %x BC: %x\nTrigger: %llx (%s)\nEvent# in file:%d\nTimestamp: %x\n",
+ fESD->GetRunNumber(),
+ acttrclasses.Data(),
+ fESD->GetEventType(),AliRawEventHeaderBase::GetTypeName(fESD->GetEventType()),
+ fESD->GetPeriodNumber(),fESD->GetOrbitNumber(),fESD->GetBunchCrossNumber(),
+ fESD->GetTriggerMask(),firedtrclasses.Data(),
+ fESD->GetEventNumberInFile(),
+ fESD->GetTimeStamp());
+ }
+
+ return rawInfo + "\n" + esdInfo;
+}
+
+
+//==============================================================================
+// Reading of GRP and MagneticField.
+// This is a reap-off from reconstruction ... should really be a common
+// code to do this somewhere in STEER.
+//==============================================================================
+
+Bool_t AliEveEventManager::InitGRP()
+{
+ //------------------------------------
+ // Initialization of the GRP entry
+ //------------------------------------
+
+ static const TEveException kEH("AliEveEventManager::InitGRP ");
+
+ AliGRPManager grpMgr;
+ if (!grpMgr.ReadGRPEntry()) {
+ return kFALSE;
+ }
+ fgGRPLoaded = kTRUE;
+ if (!grpMgr.SetMagField()) {
+ throw kEH + "Setting of field failed!";
+ }
+
+ //*** Get the diamond profiles from OCDB
+ // Eventually useful.
+
+ /*
+ entry = AliCDBManager::Instance()->Get("GRP/Calib/MeanVertexSPD");
+ if (entry) {
+ fDiamondProfileSPD = dynamic_cast<AliESDVertex*> (entry->GetObject());
+ } else {
+ ::Error(kEH, "No SPD diamond profile found in OCDB!");
+ }
+
+ entry = AliCDBManager::Instance()->Get("GRP/Calib/MeanVertex");
+ if (entry) {
+ fDiamondProfile = dynamic_cast<AliESDVertex*> (entry->GetObject());
+ } else {
+ ::Error(kEH, "No diamond profile found in OCDB!");
+ }
+
+ entry = AliCDBManager::Instance()->Get("GRP/Calib/MeanVertexTPC");
+ if (entry) {
+ fDiamondProfileTPC = dynamic_cast<AliESDVertex*> (entry->GetObject());
+ } else {
+ ::Error(kEH, "No TPC diamond profile found in OCDB!");
+ }
+ */
+
+ return kTRUE;
+}
+
+//------------------------------------
+// Global variables management
+//------------------------------------
+
+Bool_t AliEveEventManager::InsertGlobal(const TString& tag, TEveElement* model)
+{
+ // Insert a new visualization-parameter database entry with the default
+ return InsertGlobal(tag, model, fGlobalReplace, fGlobalUpdate);
+}
+
+Bool_t AliEveEventManager::InsertGlobal(const TString& tag, TEveElement* model,
+ Bool_t replace, Bool_t update)
+{
+ TPair* pair = (TPair*) fGlobal->FindObject(tag);
+ if (pair)
+ {
+ if (replace)
+ {
+ model->IncDenyDestroy();
+ model->SetRnrChildren(kFALSE);
+
+ TEveElement* old_model = dynamic_cast<TEveElement*>(pair->Value());
+ if(!old_model) AliFatal("old_model == 0, dynamic cast failed\n");
+ while (old_model->HasChildren())
+ {
+ TEveElement *el = old_model->FirstChild();
+ el->SetVizModel(model);
+ if (update)
+ {
+ el->CopyVizParams(model);
+ el->PropagateVizParamsToProjecteds();
+ }
+ }
+ old_model->DecDenyDestroy();
+
+ pair->SetValue(dynamic_cast<TObject*>(model));
+ return kTRUE;
+ }
+ else
+ {
+ return kFALSE;
+ }
+ }
+ else
+ {
+ model->IncDenyDestroy();
+ model->SetRnrChildren(kFALSE);
+ fGlobal->Add(new TObjString(tag), dynamic_cast<TObject*>(model));
+ return kTRUE;
+ }
+}
+
+TEveElement* AliEveEventManager::FindGlobal(const TString& tag)
+{
+ return dynamic_cast<TEveElement*>(fGlobal->GetValue(tag));
}
-