**************************************************************************/
#include "AliEveEventManager.h"
+#include "AliEveEventSelector.h"
#include "AliEveMacroExecutor.h"
#include <TEveManager.h>
#include <TEveViewer.h>
#include <AliMagF.h>
#include <AliCDBManager.h>
#include <AliCDBStorage.h>
-#include <AliCDBEntry.h>
#include <AliGRPObject.h>
#include <AliHeader.h>
#include <AliGeomManager.h>
+#include <AliGRPManager.h>
+#include <AliSysInfo.h>
#include <TFile.h>
#include <TTree.h>
#include <TTimeStamp.h>
#include <TPRegexp.h>
#include <TError.h>
+#include <TEnv.h>
//==============================================================================
//==============================================================================
TString AliEveEventManager::fgESDFileName("AliESDs.root");
TString AliEveEventManager::fgAODFileName("AliAOD.root");
TString AliEveEventManager::fgRawFileName("raw.root");
-TString AliEveEventManager::fgCdbUri("local://$ALICE_ROOT/OCDB");
+TString AliEveEventManager::fgCdbUri;
TList* AliEveEventManager::fgAODfriends = 0;
-AliGRPObject* AliEveEventManager::fgGRPData = 0;
+Bool_t AliEveEventManager::fgGRPLoaded = kFALSE;
AliMagF* AliEveEventManager::fgMagField = 0;
Bool_t AliEveEventManager::fgUniformField = kFALSE;
fTransientLists = new TEveElementList("Transient Lists", "Containers of transient elements.");
fTransientLists->IncDenyDestroy();
gEve->AddToListTree(fTransientLists, kFALSE);
+
+ fPEventSelector = new AliEveEventSelector(this);
}
AliEveEventManager::AliEveEventManager(const TString& name) :
fRawReader (0),
fAutoLoad (kFALSE), fAutoLoadTime (5.), fAutoLoadTimer(0),
fIsOpen (kFALSE), fHasEvent (kFALSE), fExternalCtrl (kFALSE),
- fSelectOnTriggerType(kFALSE), fTriggerType(""),
fExecutor (0), fTransients(0), fTransientLists(0),
+ fPEventSelector(0),
fSubManagers (0),
fAutoLoadTimerRunning(kFALSE)
{
fRawReader (0),
fAutoLoad (kFALSE), fAutoLoadTime (5), fAutoLoadTimer(0),
fIsOpen (kFALSE), fHasEvent (kFALSE), fExternalCtrl (kFALSE),
- fSelectOnTriggerType(kFALSE), fTriggerType(""),
fExecutor (0), fTransients(0), fTransientLists(0),
+ fPEventSelector(0),
fSubManagers (0),
fAutoLoadTimerRunning(kFALSE)
{
fgAODfriends = new TList;
fgAODfriends->SetOwner(kTRUE);
}
- fgAODfriends->Add(new TObjString(friendFileName));
+ if (fgAODfriends->FindObject(friendFileName) == 0)
+ {
+ fgAODfriends->Add(new TObjString(friendFileName));
+ }
}
void AliEveEventManager::SetRawFileName(const TString& raw)
void AliEveEventManager::SetCdbUri(const TString& cdb)
{
- // Set path to CDB, default "local://$ALICE_ROOT/OCDB".
+ // Set path to CDB, there is no default.
if ( ! cdb.IsNull()) fgCdbUri = cdb;
}
// 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.
- // There was also a problem with TTree::Refresh() - it didn't
- // save the friend branch on a separate file, fixed in 5.22.2 -
- // so we might want to try the old way again soon.
TString p(Form("%s/AliESDfriends.root", fPath.Data()));
TFile *esdFriendFile = TFile::Open(p);
if (esdFriendFile)
TString gaPath(Form("%s/galice.root", fPath.Data()));
// If i use open directly, we get fatal.
// Is AccessPathName check ok for xrootd / alien? Yes, not for http.
- if (gSystem->AccessPathName(gaPath, kReadPermission) == kFALSE)
- {
+ // 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)
+ {
+ gafile->Close();
+ delete gafile;
+ // if (gSystem->AccessPathName(gaPath, kReadPermission) == kFALSE)
+ // {
fRunLoader = AliRunLoader::Open(gaPath, GetName());
if (fRunLoader)
{
}
else
{
- cdb->SetDefaultStorage(fgCdbUri);
+ 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.");
+ 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);
+ }
+ }
}
- cdb->SetRun(runNo);
}
fIsOpen = kTRUE;
if (fAutoLoad) StartAutoLoadTimer();
}
-Int_t AliEveEventManager::GetMaxEventId(Bool_t /*refreshESD*/) const
+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 (fESDTree)
{
- // Refresh crashes with root-5.21.1-alice.
- // Fixed by Philippe 5.8.2008 r25053, can be reactivated
- // when we move to a newer root.
- // if (refreshESD)
- // fESDTree->Refresh();
+ if (refreshESD)
+ {
+ fESDTree->Refresh();
+ fPEventSelector->Update();
+ }
return fESDTree->GetEntries() - 1;
}
else if (fAODTree)
Int_t maxEvent = 0;
if (fESDTree)
{
- // Refresh crashes with root-5.21.1-alice.
- // Fixed by Philippe 5.8.2008 r25053, can be reactivated
- // when we move to a newer root.
- // fESDTree->Refresh();
+ if (event >= fESDTree->GetEntries())
+ fESDTree->Refresh();
maxEvent = fESDTree->GetEntries() - 1;
if (event < 0)
event = fESDTree->GetEntries() + event;
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.
}
DestroyElements();
+ AliSysInfo::AddStamp(sysInfoHeader + "PostDestroy");
+
if (fESDTree) {
if (fESDTree->GetEntry(event) <= 0)
throw (kEH + "failed getting required event from ESD.");
ElementChanged();
}
+ AliSysInfo::AddStamp(sysInfoHeader + "PostLoadEvent");
+
AfterNewEventLoaded();
+
+ AliSysInfo::AddStamp(sysInfoHeader + "PostUserActions");
}
void AliEveEventManager::NextEvent()
gSystem->ExitLoop();
}
- else
+ else if (fESDTree)
{
- Int_t nexteventbytrigger=0;
- if (fSelectOnTriggerType)
+ Int_t nextevent=0;
+ if (fPEventSelector->FindNext(nextevent))
{
- if (FindNextByTrigger(nexteventbytrigger)) //if not found do nothing
- GotoEvent(nexteventbytrigger);
+ GotoEvent(nextevent);
}
- else if (fEventId < GetMaxEventId(kTRUE))
- GotoEvent(fEventId + 1);
- else
- GotoEvent(0);
+ }
+ else if (fEventId < GetMaxEventId(kTRUE))
+ {
+ GotoEvent(fEventId + 1);
}
}
{
throw (kEH + "Event-loop is under external control.");
}
- Int_t nexteventbytrigger=0;
- if (fSelectOnTriggerType)
+
+ if (fESDTree)
{
- if (FindPrevByTrigger(nexteventbytrigger))
- GotoEvent(nexteventbytrigger);
+ Int_t nextevent=0;
+ if (fPEventSelector->FindPrev(nextevent))
+ {
+ GotoEvent(nextevent);
+ }
}
- else
+ else if (fEventId > 0)
+ {
GotoEvent(fEventId - 1);
+ }
}
void AliEveEventManager::Close()
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;
if (fgMagField)
return fgMagField;
- if (fgGRPData == 0)
+ if (TGeoGlobalMagField::Instance()->GetField())
+ {
+ fgMagField = dynamic_cast<AliMagF*>(TGeoGlobalMagField::Instance()->GetField());
+ if (fgMagField == 0)
+ throw kEH + "Global field set, but it is not AliMagF.";
+ return fgMagField;
+ }
+
+ if (!fgGRPLoaded)
{
InitGRP();
}
- if (TGeoGlobalMagField::Instance())
+ if (TGeoGlobalMagField::Instance()->GetField())
{
fgMagField = dynamic_cast<AliMagF*>(TGeoGlobalMagField::Instance()->GetField());
if (fgMagField == 0)
}
}
+void AliEveEventManager::SetTrigSel(Int_t trig)
+{
+ 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::StartAutoLoadTimer()
{
// Start the auto-load timer.
StartAutoLoadTimer();
}
-
-//------------------------------------------------------------------------------
-// Event selection by trigger
-//------------------------------------------------------------------------------
-
-Bool_t AliEveEventManager::FindNextByTrigger(Int_t& event)
-{
- // Find next event that matches the trigger.
- // If a matching event is not found, we loop around and eventually
- // end up at the same event.
-
- static const TEveException kEH("AliEveEventManager::FindNextByTrigger ");
-
- if (!fESDTree) return kFALSE;
- TString firedtrclasses;
- for (Int_t i = fEventId+1; i<GetMaxEventId(kTRUE)+1; i++)
- {
- if (fESDTree->GetEntry(i) <= 0)
- throw (kEH + "failed getting required event from ESD.");
- firedtrclasses = fESD->GetFiredTriggerClasses();
- if (firedtrclasses.Contains(fTriggerType))
- {
- event=i;
- return kTRUE;
- }
- }
- for (Int_t i = 0; i<fEventId+1; i++)
- {
- if (fESDTree->GetEntry(i) <= 0)
- throw (kEH + "failed getting required event from ESD.");
- firedtrclasses = fESD->GetFiredTriggerClasses();
- if (firedtrclasses.Contains(fTriggerType))
- {
- event=i;
- return kTRUE;
- }
- }
- return kFALSE;
-}
-
-Bool_t AliEveEventManager::FindPrevByTrigger(Int_t& event)
-{
- // Find previous event that matches the trigger.
-
- static const TEveException kEH("AliEveEventManager::FindPrevByTrigger ");
-
- if (!fESDTree) return kFALSE;
- TString firedtrclasses;
- for (Int_t i = fEventId-1; i>=0; i--)
- {
- if (fESDTree->GetEntry(i) <= 0)
- throw (kEH + "failed getting required event from ESD.");
- firedtrclasses = fESD->GetFiredTriggerClasses();
- if (firedtrclasses.Contains(fTriggerType))
- {
- event=i;
- return kTRUE;
- }
- }
- for (Int_t i = GetMaxEventId(kTRUE); i>fEventId-1; i--)
- {
- if (fESDTree->GetEntry(i) <= 0)
- throw (kEH + "failed getting required event from ESD.");
- firedtrclasses = fESD->GetFiredTriggerClasses();
- if (firedtrclasses.Contains(fTriggerType))
- {
- event=i;
- return kTRUE;
- }
- }
- return kFALSE;
-}
-
-
//------------------------------------------------------------------------------
// Post event-loading functions
//------------------------------------------------------------------------------
static const TEveException kEH("AliEveEventManager::AfterNewEventLoaded ");
+ NewEventDataLoaded();
+
if (fExecutor)
fExecutor->ExecMacros();
}
}
+void AliEveEventManager::NewEventDataLoaded()
+{
+ // Emit NewEventDataLoaded signal.
+
+ Emit("NewEventDataLoaded()");
+}
+
void AliEveEventManager::NewEventLoaded()
{
// Emit NewEventLoaded signal.
static const TEveException kEH("AliEveEventManager::InitGRP ");
- AliCDBEntry* entry = AliCDBManager::Instance()->Get("GRP/GRP/Data");
-
- if (entry)
- {
- TMap* m = dynamic_cast<TMap*>(entry->GetObject()); // old GRP entry
-
- if (m)
- {
- ::Info(kEH, "Found a TMap in GRP/GRP/Data, converting it into an AliGRPObject");
- fgGRPData = new AliGRPObject();
- fgGRPData->ReadValuesFromMap(m);
- }
- else
- {
- ::Info(kEH, "Found an AliGRPObject in GRP/GRP/Data, reading it");
- fgGRPData = dynamic_cast<AliGRPObject*>(entry->GetObject()); // new GRP entry
- entry->SetOwner(0);
- }
-
- AliCDBManager::Instance()->UnloadFromCache("GRP/GRP/Data");
- }
-
- if (!fgGRPData)
- {
- ::Error(kEH, "No GRP entry found in OCDB!");
+ AliGRPManager grpMgr;
+ if (!grpMgr.ReadGRPEntry()) {
return kFALSE;
}
-
- TString lhcState = fgGRPData->GetLHCState();
- if (lhcState==AliGRPObject::GetInvalidString())
- {
- ::Error(kEH, "GRP/GRP/Data entry: missing value for the LHC state ! Using UNKNOWN");
- lhcState = "UNKNOWN";
- }
-
- TString beamType = fgGRPData->GetBeamType();
- if (beamType==AliGRPObject::GetInvalidString())
- {
- ::Error(kEH, "GRP/GRP/Data entry: missing value for the beam type ! Using UNKNOWN");
- beamType = "UNKNOWN";
- }
-
- Float_t beamEnergy = fgGRPData->GetBeamEnergy();
- if (beamEnergy==AliGRPObject::GetInvalidFloat())
- {
- ::Error(kEH, "GRP/GRP/Data entry: missing value for the beam energy ! Using 0");
- beamEnergy = 0;
+ fgGRPLoaded = kTRUE;
+ if (!grpMgr.SetMagField()) {
+ throw kEH + "Setting of field failed!";
}
- beamEnergy /= 120E3; // energy is provided in MeV*120
-
- TString runType = fgGRPData->GetRunType();
- if (runType==AliGRPObject::GetInvalidString())
- {
- ::Error(kEH, "GRP/GRP/Data entry: missing value for the run type ! Using UNKNOWN");
- runType = "UNKNOWN";
- }
-
- Int_t activeDetectors = fgGRPData->GetDetectorMask();
- if (activeDetectors==AliGRPObject::GetInvalidUInt())
- {
- ::Error(kEH, "GRP/GRP/Data entry: missing value for the detector mask ! Using 1074790399");
- activeDetectors = 1074790399;
- }
-
-
- // Might become useful.
- /*
- fRunInfo = new AliRunInfo(lhcState, beamType, beamEnergy, runType, activeDetectors);
- printf("AliRunInfo - %s %s %f %s %d\n", lhcState.Data(), beamType.Data(), beamEnergy, runType.Data(), activeDetectors);
- fRunInfo->Dump();
-
-
- // Process the list of active detectors
- if (activeDetectors) {
- UInt_t detMask = activeDetectors;
- fRunLocalReconstruction = MatchDetectorList(fRunLocalReconstruction,detMask);
- fRunTracking = MatchDetectorList(fRunTracking,detMask);
- fFillESD = MatchDetectorList(fFillESD,detMask);
- fQADetectors = MatchDetectorList(fQADetectors,detMask);
- fLoadCDB.Form("%s %s %s %s",
- fRunLocalReconstruction.Data(),
- fRunTracking.Data(),
- fFillESD.Data(),
- fQADetectors.Data());
- fLoadCDB = MatchDetectorList(fLoadCDB,detMask);
- if (!((detMask >> AliDAQ::DetectorID("ITSSPD")) & 0x1)) {
- // switch off the vertexer
- ::Info(kEH, "SPD is not in the list of active detectors. Vertexer switched off.");
- fRunVertexFinder = kFALSE;
- }
- if (!((detMask >> AliDAQ::DetectorID("TRG")) & 0x1)) {
- // switch off the reading of CTP raw-data payload
- if (fFillTriggerESD) {
- ::Info(kEH, "CTP is not in the list of active detectors. CTP data reading switched off.");
- fFillTriggerESD = kFALSE;
- }
- }
- }
-
- ::Info(kEH, "===================================================================================");
- ::Info(kEH, "Running local reconstruction for detectors: %s",fRunLocalReconstruction.Data());
- ::Info(kEH, "Running tracking for detectors: %s",fRunTracking.Data());
- ::Info(kEH, "Filling ESD for detectors: %s",fFillESD.Data());
- ::Info(kEH, "Quality assurance is active for detectors: %s",fQADetectors.Data());
- ::Info(kEH, "CDB and reconstruction parameters are loaded for detectors: %s",fLoadCDB.Data());
- ::Info(kEH, "===================================================================================");
- */
-
- //*** Dealing with the magnetic field map
- if (TGeoGlobalMagField::Instance()->IsLocked())
- {
- ::Info(kEH, "Running with externally locked B field!");
- }
- else
- {
- // Construct the field map out of the information retrieved from GRP.
-
- Float_t l3Current = fgGRPData->GetL3Current((AliGRPObject::Stats)0);
- if (l3Current == AliGRPObject::GetInvalidFloat())
- throw kEH + "GRPData entry: missing value for the L3 current!";
-
- Char_t l3Polarity = fgGRPData->GetL3Polarity();
- if (l3Polarity == AliGRPObject::GetInvalidChar())
- throw kEH + "GRPData entry: missing value for the L3 polarity!";
- Float_t diCurrent = fgGRPData->GetDipoleCurrent((AliGRPObject::Stats)0);
- if (diCurrent == AliGRPObject::GetInvalidFloat())
- throw kEH + "GRPData entry: missing value for the dipole current!";
-
- Char_t diPolarity = fgGRPData->GetDipolePolarity();
- if (diPolarity == AliGRPObject::GetInvalidChar())
- throw kEH + "GRPData entry: missing value for the dipole polarity!";
-
- if (!SetFieldMap(l3Current, diCurrent, l3Polarity ? -1:1, diPolarity ? -1:1))
- throw kEH + "Failed to create B field map!";
-
- ::Info(kEH, "Running with the B field constructed from GRP.");
- }
-
//*** Get the diamond profiles from OCDB
// Eventually useful.
return kTRUE;
}
-
-Bool_t AliEveEventManager::SetFieldMap(Float_t l3Cur, Float_t diCur,
- Float_t l3Pol, Float_t diPol,
- Float_t beamenergy,
- const Char_t *beamtype,
- const Char_t *path)
-{
- //------------------------------------------------
- // The magnetic field map, defined externally...
- // L3 current 30000 A -> 0.5 T
- // L3 current 12000 A -> 0.2 T
- // dipole current 6000 A
- // The polarities must be the same
- //------------------------------------------------
-
- static const TEveException kEH("AliEveEventManager::SetFieldMap ");
-
- const Float_t l3NominalCurrent1 = 30000.0; // (A)
- const Float_t l3NominalCurrent2 = 12000.0; // (A)
- const Float_t diNominalCurrent = 6000.0; // (A)
-
- const Float_t tolerance = 0.03; // relative current tolerance
- const Float_t zero = 77.0; // "zero" current (A)
-
- TString s = (l3Pol < 0) ? "L3: -" : "L3: +";
-
- AliMagF::BMap_t map = AliMagF::k5kG;
-
- double fcL3, fcDip;
-
- l3Cur = TMath::Abs(l3Cur);
- if (TMath::Abs(l3Cur-l3NominalCurrent1)/l3NominalCurrent1 < tolerance)
- {
- fcL3 = l3Cur/l3NominalCurrent1;
- map = AliMagF::k5kG;
- s += "0.5 T; ";
- }
- else if (TMath::Abs(l3Cur-l3NominalCurrent2)/l3NominalCurrent2 < tolerance)
- {
- fcL3 = l3Cur/l3NominalCurrent2;
- map = AliMagF::k2kG;
- s += "0.2 T; ";
- }
- else if (l3Cur <= zero)
- {
- fcL3 = 0;
- map = AliMagF::k5kGUniform;
- s += "0.0 T; ";
- fgUniformField = kTRUE; // track with the uniform (zero) B field
- }
- else
- {
- ::Error(kEH, "Wrong L3 current (%f A)!", l3Cur);
- return kFALSE;
- }
-
- diCur = TMath::Abs(diCur);
- if (TMath::Abs(diCur-diNominalCurrent)/diNominalCurrent < tolerance)
- {
- // 3% current tolerance...
- fcDip = diCur/diNominalCurrent;
- s += "Dipole ON";
- }
- else if (diCur <= zero)
- { // some small current..
- fcDip = 0.;
- s += "Dipole OFF";
- }
- else
- {
- ::Error(kEH, "Wrong dipole current (%f A)!", diCur);
- return kFALSE;
- }
-
- if (l3Pol != diPol && (map==AliMagF::k5kG || map==AliMagF::k2kG) && fcDip != 0)
- {
- ::Error(kEH, "L3 and Dipole polarities must be the same");
- return kFALSE;
- }
-
- if (l3Pol<0) fcL3 = -fcL3;
- if (diPol<0) fcDip = -fcDip;
-
- AliMagF::BeamType_t btype = AliMagF::kNoBeamField;
- TString btypestr = beamtype;
- btypestr.ToLower();
-
- TPRegexp protonBeam("(proton|p)\\s*-?\\s*\\1");
- TPRegexp ionBeam ("(lead|pb|ion|a)\\s*-?\\s*\\1");
-
- if (btypestr.Contains(ionBeam))
- btype = AliMagF::kBeamTypeAA;
- else if (btypestr.Contains(protonBeam))
- btype = AliMagF::kBeamTypepp;
- else
- ::Info(kEH, "Cannot determine the beam type from %s, assume no LHC magnet field",
- beamtype);
-
- AliMagF* fld = new AliMagF("MagneticFieldMap", s.Data(), 2, fcL3, fcDip, 10., map, path,
- btype,beamenergy);
- TGeoGlobalMagField::Instance()->SetField(fld);
- TGeoGlobalMagField::Instance()->Lock();
-
- return kTRUE;
-}