-// $Id: NLT_trackcount_init.C 24927 2008-04-04 13:46:04Z mtadel $
+// $Id$
// Main authors: Matevz Tadel & Alja Mrak-Tadel: 2006, 2007
/**************************************************************************
* full copyright notice. *
**************************************************************************/
+class AliEveMacroExecutor;
+
class TEveProjectionManager;
class TEveGeoShape;
class TEveUtil;
-R__EXTERN TEveProjectionManager *gRPhiMgr;
-R__EXTERN TEveProjectionManager *gRhoZMgr;
+TEveGeoShape *gGeomGentle = 0;
+TEveGeoShape *gGeomGentleRPhi = 0;
+TEveGeoShape *gGeomGentleRhoZ = 0;
+TEveGeoShape *gGeomGentleTRD = 0;
+TEveGeoShape *gGeomGentleMUON = 0;
+
+TEveScene *gRPhiGeomScene = 0;
+TEveScene *gRhoZGeomScene = 0;
+TEveScene *gRPhiEventScene = 0;
+TEveScene *gRhoZEventScene = 0;
+
+TEveProjectionManager *gRPhiMgr = 0;
+TEveProjectionManager *gRhoZMgr = 0;
+
+TEveViewer *g3DView = 0;
+TEveViewer *gRPhiView = 0;
+TEveViewer *gRhoZView = 0;
-TEveGeoShape *gGeomGentle = 0;
-TEveGeoShape *gGeomGentleTRD = 0;
+Bool_t gShowTRD = kFALSE;
+Bool_t gShowMUON = kTRUE;
+Bool_t gShowMUONRPhi = kFALSE;
+Bool_t gShowMUONRhoZ = kTRUE;
-Bool_t gShowTRD = kFALSE;
+Bool_t gCenterProjectionsAtPrimaryVertex = kFALSE;
-void visscan_init()
+void visscan_init(Bool_t show_extra_geo=kFALSE)
{
+ if (!show_extra_geo)
+ {
+ gShowTRD = gShowMUON = gShowMUONRPhi = gShowMUONRhoZ = kFALSE;
+ }
+
TEveUtil::LoadMacro("alieve_init.C");
alieve_init(".", -1);
- TEveBrowser* b = gEve->GetBrowser();
- b->StartEmbedding(TRootBrowser::kBottom);
- new AliEveEventManagerWindow;
- b->StopEmbedding("EventCtrl");
+ // TEveLine::SetDefaultSmooth(1);
+ TEveUtil::AssertMacro("VizDB_scan.C");
+
+ AliEveMacroExecutor *exec = AliEveEventManager::GetMaster()->GetExecutor();
+ TEveBrowser *browser = gEve->GetBrowser();
+
+ //==============================================================================
+ // Geometry, scenes, projections and viewers
+ //==============================================================================
+
+ browser->ShowCloseTab(kFALSE);
+
+ // Geometry
TEveUtil::LoadMacro("geom_gentle.C");
- if (gShowTRD) TEveUtil::LoadMacro("geom_gentle_trd.C");
-
- TEveUtil::LoadMacro("primary_vertex.C");
- TEveUtil::LoadMacro("esd_V0_points.C");
- TEveUtil::LoadMacro("esd_V0.C");
- TEveUtil::LoadMacro("esd_tracks.C");
- TEveUtil::LoadMacro("its_clusters.C+");
- TEveUtil::LoadMacro("tpc_clusters.C+");
- TEveUtil::LoadMacro("trd_clusters.C+");
- TEveUtil::LoadMacro("tof_clusters.C+");
-
- TEveLine::SetDefaultSmooth(1);
-
- /*
- TEveViewer* nv = gEve->SpawnNewViewer("NLT Projected");
- TEveScene* ns = gEve->SpawnNewScene("NLT");
- nv->AddScene(ns);
- TGLViewer* v = nv->GetGLViewer();
- v->SetCurrentCamera(TGLViewer::kCameraOrthoXOY);
- TGLCameraMarkupStyle* mup = v->GetCameraMarkup();
- if(mup) mup->SetShow(kFALSE);
- */
+ gGeomGentle = geom_gentle();
+ gGeomGentleRPhi = geom_gentle_rphi(); gGeomGentleRPhi->IncDenyDestroy();
+ gGeomGentleRhoZ = geom_gentle_rhoz(); gGeomGentleRhoZ->IncDenyDestroy();
+ if (gShowTRD) {
+ TEveUtil::LoadMacro("geom_gentle_trd.C");
+ gGeomGentleTRD = geom_gentle_trd();
+ }
+ if (gShowMUON) {
+ TEveUtil::LoadMacro("geom_gentle_muon.C");
+ gGeomGentleMUON = geom_gentle_muon();
+ }
+
+ // Scenes
+
+ gRPhiGeomScene = gEve->SpawnNewScene("RPhi Geometry",
+ "Scene holding projected geometry for the RPhi view.");
+ gRhoZGeomScene = gEve->SpawnNewScene("RhoZ Geometry",
+ "Scene holding projected geometry for the RhoZ view.");
+ gRPhiEventScene = gEve->SpawnNewScene("RPhi Event Data",
+ "Scene holding projected geometry for the RPhi view.");
+ gRhoZEventScene = gEve->SpawnNewScene("RhoZ Event Data",
+ "Scene holding projected geometry for the RhoZ view.");
+
+ // Projection managers
+
+ gRPhiMgr = new TEveProjectionManager();
+ gRPhiMgr->SetProjection(TEveProjection::kPT_RPhi);
+ gEve->AddToListTree(gRPhiMgr, kFALSE);
+ {
+ TEveProjectionAxes* a = new TEveProjectionAxes(gRPhiMgr);
+ a->SetMainColor(kWhite);
+ a->SetTitle("R-Phi");
+ a->SetTitleSize(0.05);
+ a->SetTitleFont(102);
+ a->SetLabelSize(0.025);
+ a->SetLabelFont(102);
+ gRPhiGeomScene->AddElement(a);
+ }
+ gRPhiMgr->ImportElements(gGeomGentleRPhi, gRPhiGeomScene);
+ if (gShowTRD) gRPhiMgr->ImportElements(gGeomGentleTRD, gRPhiGeomScene);
+ if (gShowMUONRPhi) gRPhiMgr->ImportElements(gGeomGentleMUON, gRPhiGeomScene);
+
+ gRhoZMgr = new TEveProjectionManager();
+ gRhoZMgr->SetProjection(TEveProjection::kPT_RhoZ);
+ gEve->AddToListTree(gRhoZMgr, kFALSE);
+ {
+ TEveProjectionAxes* a = new TEveProjectionAxes(gRhoZMgr);
+ a->SetMainColor(kWhite);
+ a->SetTitle("Rho-Z");
+ a->SetTitleSize(0.05);
+ a->SetTitleFont(102);
+ a->SetLabelSize(0.025);
+ a->SetLabelFont(102);
+ gRhoZGeomScene->AddElement(a);
+ }
+ gRhoZMgr->ImportElements(gGeomGentleRhoZ, gRhoZGeomScene);
+ if (gShowTRD) gRhoZMgr->ImportElements(gGeomGentleTRD, gRhoZGeomScene);
+ if (gShowMUONRhoZ) gRhoZMgr->ImportElements(gGeomGentleMUON, gRhoZGeomScene);
+
+ // Viewers
+
+ TEveWindowSlot *slot = 0;
+ TEveWindowPack *pack = 0;
+
+ slot = TEveWindow::CreateWindowInTab(browser->GetTabRight());
+ pack = slot->MakePack();
+ pack->SetElementName("Multi View");
+ pack->SetHorizontal();
+ pack->SetShowTitleBar(kFALSE);
+ pack->NewSlot()->MakeCurrent();
+ g3DView = gEve->SpawnNewViewer("3D View", "");
+ g3DView->AddScene(gEve->GetGlobalScene());
+ g3DView->AddScene(gEve->GetEventScene());
+
+ pack = pack->NewSlot()->MakePack();
+ pack->SetShowTitleBar(kFALSE);
+ pack->NewSlot()->MakeCurrent();
+ gRPhiView = gEve->SpawnNewViewer("RPhi View", "");
+ gRPhiView->GetGLViewer()->SetCurrentCamera(TGLViewer::kCameraOrthoXOY);
+ gRPhiView->AddScene(gRPhiGeomScene);
+ gRPhiView->AddScene(gRPhiEventScene);
+
+ pack->NewSlot()->MakeCurrent();
+ gRhoZView = gEve->SpawnNewViewer("RhoZ View", "");
+ gRhoZView->GetGLViewer()->SetCurrentCamera(TGLViewer::kCameraOrthoXOY);
+ gRhoZView->AddScene(gRhoZGeomScene);
+ gRhoZView->AddScene(gRhoZEventScene);
+
+
+ //==============================================================================
+ // Registration of per-event macros
+ //==============================================================================
+
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRunLoader, "SIM Track", "kine_tracks.C", "kine_tracks", "", kFALSE));
+
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRunLoader, "SIM Hit ITS", "its_hits.C", "its_hits", "", kFALSE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRunLoader, "SIM Hit TPC", "tpc_hits.C", "tpc_hits", "", kFALSE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRunLoader, "SIM Hit T0", "t0_hits.C", "t0_hits", "", kFALSE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRunLoader, "SIM Hit FMD", "fmd_hits.C", "fmd_hits", "", kFALSE));
+
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRunLoader, "DIG FMD", "fmd_digits.C", "fmd_digits", "", kFALSE));
+
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRawReader, "RAW TPC", "tpc_raw.C", "tpc_raw", "", kFALSE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRawReader, "RAW T0", "t0_raw.C", "t0_raw", "", kFALSE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRawReader, "RAW FMD", "fmd_raw.C", "fmd_raw", "", kFALSE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRawReader, "RAW VZERO", "vzero_raw.C", "vzero_raw", "", kFALSE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRawReader, "RAW ACORDE", "acorde_raw.C", "acorde_raw", "", kFALSE));
+
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC PVTX", "primary_vertex.C", "primary_vertex", "", kTRUE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC PVTX Ellipse", "primary_vertex.C", "primary_vertex_ellipse", "", kTRUE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC PVTX Box", "primary_vertex.C", "primary_vertex_box", "kFALSE, 3, 3, 3", kFALSE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC PVTX", "primary_vertex.C", "primary_vertex_spd", "", kTRUE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC PVTX Ellipse", "primary_vertex.C", "primary_vertex_ellipse_spd", "", kTRUE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC PVTX Box", "primary_vertex.C", "primary_vertex_box_spd", "kFALSE, 3, 3, 3", kFALSE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC PVTX", "primary_vertex.C", "primary_vertex_tpc", "", kFALSE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC PVTX Ellipse", "primary_vertex.C", "primary_vertex_ellipse_tpc", "", kFALSE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC PVTX Box", "primary_vertex.C", "primary_vertex_box_tpc", "kFALSE, 3, 3, 3", kFALSE));
+
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC V0", "esd_V0_points.C", "esd_V0_points_onfly"));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC V0", "esd_V0_points.C", "esd_V0_points_offline"));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC V0", "esd_V0.C", "esd_V0"));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC CSCD", "esd_cascade_points.C", "esd_cascade_points"));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC CSCD", "esd_cascade.C", "esd_cascade"));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC KINK", "esd_kink_points.C", "esd_kink_points"));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC KINK", "esd_kink.C", "esd_kink"));
+
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC Track", "esd_tracks.C", "esd_tracks", "", kFALSE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC Track", "esd_tracks.C", "esd_tracks_MI", "", kFALSE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC Track", "esd_tracks.C", "esd_tracks_by_category", "", kTRUE));
+
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC Tracklet", "esd_spd_tracklets.C", "esd_spd_tracklets", "", kFALSE));
+
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kESD, "REC ZDC", "esd_zdc.C", "esd_zdc", "", kFALSE));
+
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRunLoader, "REC Clus", "clusters.C+", "clusters", "", kFALSE));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRunLoader, "REC Clus ITS", "its_clusters.C+", "its_clusters"));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRunLoader, "REC Clus TPC", "tpc_clusters.C+", "tpc_clusters"));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRunLoader, "REC Clus TRD", "trd_clusters.C+", "trd_clusters"));
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRunLoader, "REC Clus TOF", "tof_clusters.C+", "tof_clusters"));
+
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kRunLoader, "REC Clus TPC", "vplot_tpc.C+", "vplot_tpc", "", kFALSE));
+
+ exec->AddMacro(new AliEveMacro(AliEveMacro::kAOD, "ANA Jets", "jetplane.C", "jetplane", "", kFALSE));
+
+ //==============================================================================
+ // Additional GUI components
+ //==============================================================================
+
+ slot = TEveWindow::CreateWindowInTab(browser->GetTabRight());
+ slot->StartEmbedding();
+ AliEveMacroExecutorWindow* exewin = new AliEveMacroExecutorWindow(exec);
+ slot->StopEmbedding("DataSelection");
+ exewin->PopulateMacros();
+
+ slot = TEveWindow::CreateWindowInTab(browser->GetTabRight());
+ slot->StartEmbedding();
+ new AliQAHistViewer(gClient->GetRoot(), 600, 400, kTRUE);
+ slot->StopEmbedding("QA histograms");
+
+ browser->GetTabRight()->SetTab(1);
+
+ browser->StartEmbedding(TRootBrowser::kBottom);
+ new AliEveEventManagerWindow(AliEveEventManager::GetMaster());
+ browser->StopEmbedding("EventCtrl");
+
+ slot = TEveWindow::CreateWindowInTab(browser->GetTabRight());
+ TEveWindowTab *store_tab = slot->MakeTab();
+ store_tab->SetElementNameTitle("WindowStore",
+ "Undocked windows whose previous container is not known\n"
+ "are placed here when the main-frame is closed.");
+ gEve->GetWindowManager()->SetDefaultContainer(store_tab);
+
+ //==============================================================================
+ // AliEve objects - global tools
+ //==============================================================================
AliEveTrackFitter* fitter = new AliEveTrackFitter();
gEve->AddToListTree(fitter, 1);
gEve->AddToListTree(g_trkcnt, kFALSE);
- // geometry
- gGeomGentle = geom_gentle();
- if (gShowTRD) gGeomGentleTRD = geom_gentle_trd();
-
-
- gROOT->ProcessLine(".L SplitGLView.C+");
- TEveBrowser* browser = gEve->GetBrowser();
- browser->ExecPlugin("SplitGLView", 0, "new SplitGLView(gClient->GetRoot(), 600, 450, kTRUE)");
-
- if (gRPhiMgr) {
- TEveProjectionAxes* a = new TEveProjectionAxes(gRPhiMgr);
- a->SetNumTickMarks(3);
- a->SetText("R-Phi");
- a->SetFontFile("comicbd");
- a->SetFontSize(10);
- gEve->GetScenes()->FindChild("R-Phi Projection")->AddElement(a);
- }
- if (gRhoZMgr) {
- TEveProjectionAxes* a = new TEveProjectionAxes(gRhoZMgr);
- a->SetNumTickMarks(3);
- a->SetText("Rho-Z");
- a->SetFontFile("comicbd");
- a->SetFontSize(10);
- gEve->GetScenes()->FindChild("Rho-Z Projection")->AddElement(a);
- }
-
- // event
- gAliEveEvent->AddNewEventCommand("on_new_event();");
- gAliEveEvent->GotoEvent(0);
+ //==============================================================================
+ // Final stuff
+ //==============================================================================
+
+ // A refresh to show proper window.
+ gEve->Redraw3D(kTRUE);
+ gSystem->ProcessEvents();
+
+ // Register command to call on each event.
+ AliEveEventManager::GetMaster()->AddNewEventCommand("on_new_event();");
+ AliEveEventManager::GetMaster()->GotoEvent(0);
gEve->EditElement(g_trkcnt);
void on_new_event()
{
- try {
- TEvePointSet* itsc = its_clusters();
- if (itsc) {
- itsc->SetMarkerColor(5);
- }
-
- TEvePointSet* tpcc = tpc_clusters();
- if (tpcc) {
- tpcc->SetMarkerColor(4);
- }
+ AliEveTrackCounter* g_trkcnt = AliEveTrackCounter::fgInstance;
+ g_trkcnt->Reset();
+ g_trkcnt->SetEventId(AliEveEventManager::GetMaster()->GetEventId());
- TEvePointSet* trdc = trd_clusters();
- if (trdc) {
- trdc->SetMarkerColor(7);
- trdc->SetMarkerStyle(4);
- trdc->SetMarkerSize(0.5);
+ if (g_esd_tracks_by_category_container != 0)
+ {
+ TEveElementList* cont = g_esd_tracks_by_category_container;
+
+ // Here we expect several TEveTrackList containers.
+ // First two have reasonable primaries (sigma-to-prim-vertex < 5).
+ // Others are almost certainly secondaries.
+ Int_t count = 1;
+ TEveElement::List_i i = cont->BeginChildren();
+ while (i != cont->EndChildren())
+ {
+ TEveTrackList* l = dynamic_cast<TEveTrackList*>(*i);
+ if (l != 0)
+ {
+ g_trkcnt->RegisterTracks(l, (count <= 2));
+ ++count;
+ }
+ ++i;
}
- TEvePointSet* tofc = tof_clusters();
- if (tofc) {
- tofc->SetMarkerColor(kOrange);
- tofc->SetMarkerStyle(4);
- tofc->SetMarkerSize(0.5);
- }
+ // Set it to zero, so that we do not reuse an old one.
+ g_esd_tracks_by_category_container = 0;
}
- catch(TEveException& exc) {
- printf("Exception loading ITS/TPC clusters: %s\n", exc.Data());
+ else
+ {
+ Warning("on_new_event", "g_esd_tracks_by_category_container not initialized.");
}
- primary_vertex(1, 1);
- esd_V0_points();
- esd_V0();
-
- TEveElementList* cont = esd_tracks_vertex_cut();
-
- // Here we expect five TEveTrackList containers.
- // First two have reasonable primaries (sigma-to-prim-vertex < 5).
- // Other three are almost certainly secondaries.
- Int_t count = 1;
- AliEveTrackCounter* g_trkcnt = AliEveTrackCounter::fgInstance;
- g_trkcnt->Reset();
- g_trkcnt->SetEventId(gAliEveEvent->GetEventId());
- TEveElement::List_i i = cont->BeginChildren();
- while (i != cont->EndChildren()) {
- TEveTrackList* l = dynamic_cast<TEveTrackList*>(*i);
- if (l != 0) {
- // l->SetLineWidth(2);
- g_trkcnt->RegisterTracks(l, (count <= 2));
- ++count;
- }
- ++i;
- }
+ Double_t x[3] = { 0, 0, 0 };
- AliESDEvent* esd = AliEveEventManager::AssertESD();
+ if (AliEveEventManager::HasESD())
{
+ AliESDEvent* esd = AliEveEventManager::AssertESD();
+ esd->GetPrimaryVertex()->GetXYZ(x);
+
TTimeStamp ts(esd->GetTimeStamp());
TString win_title("Eve Main Window -- Timestamp: ");
win_title += ts.AsString("s");
win_title += esd->GetEventNumberInFile();
gEve->GetBrowser()->SetWindowName(win_title);
}
- Double_t x[3];
- esd->GetPrimaryVertex()->GetXYZ(x);
TEveElement* top = gEve->GetCurrentEvent();
- if (gRPhiMgr && top) {
- gRPhiMgr->DestroyElements();
- gRPhiMgr->SetCenter(x[0], x[1], x[2]);
- gRPhiMgr->ImportElements(gGeomGentle);
- if (gShowTRD) gRPhiMgr->ImportElements(gGeomGentleTRD);
- gRPhiMgr->ImportElements(top);
- }
- if (gRhoZMgr && top) {
- gRhoZMgr->DestroyElements();
- gRhoZMgr->SetCenter(x[0], x[1], x[2]);
- gRhoZMgr->ImportElements(gGeomGentle);
- if (gShowTRD) gRhoZMgr->ImportElements(gGeomGentleTRD);
- gRhoZMgr->ImportElements(top);
- }
-
- gROOT->ProcessLine("SplitGLView::UpdateSummary()");
-}
-
-/******************************************************************************/
-
-TParticle* id(Int_t label=0, Bool_t showParents=kTRUE)
-{
- AliRunLoader* rl = AliEveEventManager::AssertRunLoader();
- rl->LoadKinematics();
- AliStack* stack = rl->Stack();
-
- printf("Number primaries %d, all particles %d, label %d\n",
- stack->GetNprimary(), stack->GetNtrack(), label);
- if (label < 0 || label >= stack->GetNtrack()) {
- printf(" Label exceeds available range.\n");
- return 0;
+ if (gRPhiMgr && top)
+ {
+ gRPhiEventScene->DestroyElements();
+ if (gCenterProjectionsAtPrimaryVertex)
+ gRPhiMgr->SetCenter(x[0], x[1], x[2]);
+ gRPhiMgr->ImportElements(top, gRPhiEventScene);
}
-
- TParticle* part = stack->Particle(label);
- if (part != 0) {
- part->Print();
- if (showParents) {
- while (part->GetMother(0) >= 0) {
- part = stack->Particle(part->GetMother(0));
- part->Print();
- }
- }
+ if (gRhoZMgr && top)
+ {
+ gRhoZEventScene->DestroyElements();
+ if (gCenterProjectionsAtPrimaryVertex)
+ gRhoZMgr->SetCenter(x[0], x[1], x[2]);
+ gRhoZMgr->ImportElements(top, gRhoZEventScene);
}
- return stack->Particle(label);
}