//-*- Mode: C++ -*- // ** USED macros : // *************************************************** // - hlt_alieve_init.C // - VizDB_scan.C // - geom_gentle_hlt.C // - geom_gentle_muon.C // *************************************************** #if !defined(__CINT__) || defined(__MAKECINT__) //****************** ROOT ****************************************** #include "TRandom.h" #include "TVirtualPad.h" #include "TGLViewer.h" #include "TThread.h" #include "TGFileBrowser.h" #include "TStyle.h" #include "TList.h" #include "TDirectory.h" #include "TFile.h" #include "TH2.h" #include "TCanvas.h" #include "TColor.h" #include "TVector3.h" //****************** ROOT/EVE ************************************** #include "TEveManager.h" #include "AliEveHOMERManager.h" #include "AliEveHLTEventManagerOffline.h" #include "geom_gentle_hlt.C" //*************************************************************** #include "HLT/rec/AliHLTReconstructor.h" #endif // ----------------------------------------------------------------- // -- Geometry / Scenes -- // ----------------------------------------------------------------- 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; // ----------------------------------------------------------------- // -- Geometry / Scenes Parameters -- // ----------------------------------------------------------------- // -- Parameters to show different geometries Bool_t gShowMUON = kTRUE; Bool_t gShowMUONRPhi = kFALSE; Bool_t gShowMUONRhoZ = kTRUE; Bool_t gShowTRD = kFALSE; // ----------------------------------------------------------------- // -- Members -- // ----------------------------------------------------------------- // -- Timer for automatic event loop TTimer eventTimer; TTimer eventTimerFast; // -- HOMERManager AliEveHLTEventManagerOffline* geventManager = 0; // -- Geometry Manager TGeoManager* gGeoManager = 0; AliPHOSGeometry* gPHOSGeom = 0; // -- Cluster members TEvePointSet* gSPDClusters = 0; TEvePointSet* gSSDClusters = 0; TEvePointSet* gSDDClusters = 0; TEvePointSet* gTRDClusters = 0; TEvePointSetArray* gTRDColClusters = 0; TEvePointSet* gTPCClusters = 0; TEvePointSet* gTPCTestClusters = 0; TEvePointSetArray* gTPCColClusters = 0; TEveBoxSet* gPHOSBoxSet[5] = {0, 0, 0, 0, 0}; TEveBoxSet* gEMCALBoxSet[13] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}; TEvePointSet* gMUONClusters = 0; TEveStraightLineSet* gMUONTracks = 0; // -- Text output members TEveText* gHLTText = 0; // -- Tracks members TEveTrackList* gTPCTrack = 0; // -- Canvas for histograms TCanvas* gTRDCanvas = 0; TCanvas* gTPCCanvas = 0; TCanvas* gTPCClustCanvas = 0; TCanvas* gTRDCalibCanvas = 0; TCanvas* gTRDEORCanvas = 0; TCanvas* gPrimVertexCanvas = 0; TCanvas* gSPDVertexCanvas = 0; TCanvas* gITSCanvas = 0; TCanvas* gSSDCanvas0 = 0; TCanvas* gSSDCanvas1 = 0; TCanvas* gV0Canvas = 0; TCanvas* gPHOSCanvas = NULL; TCanvas* gEMCALCanvas = 0; // -- vertex -- Int_t gSPDVertexHistoCount = 0; // -- TRD -- Int_t gTRDHistoCount = 0; Int_t gTRDEvents = 0; Int_t gTRDBins = 12; // -- TPC -- Int_t gTPCBins = 15; TH1F* gTPCCharge = 0; TH1F* gTPCQMax = 0; TH1F* gTPCQMaxOverCharge = 0; TH1F* gTPCPt = 0; // KK TH1F* gTPCEta = 0; TH1F* gTPCPsi = 0; TH1F* gTPCnClusters = 0; TH1F* gTPCMult = 0; // -- PHOS -- TEveElementList* gPHOSElementList = 0; Int_t gPHOSHistoCount =0; // -- EMCAL TEveElementList* gEMCALElementList = 0; TGeoNode* gEMCALNode = 0; Int_t gEMCALHistoCount =0; // --- Flag if eventloop is running Bool_t gEventLoopStarted = kFALSE; //Container for gGeoManager till it is broken TGeoManager *fGeoManager = 0; // ----------------------------------------------------------------- // -- Methods -- // ----------------------------------------------------------------- Int_t initializeEveViewer( Bool_t showBarrel, Bool_t showMuon ); void writeToFile(); // ################################################################# // ################################################################# // ################################################################# // ----------------------------------------------------------------- void off ( Bool_t showBarrel = kTRUE, Bool_t showMuon = kFALSE ) { // -- Loading Geometry // --------------------- Int_t run = 67179; AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB"); AliCDBManager::Instance()->SetRun(run); AliGeomManager::LoadGeometry(); // The default in the simulation is the following line // TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1, AliMagF::k5kG)); // However for the current setting of +ve L3 and +ve Dipole magnetic field // the following setting creates the field close to real field with currect polarity if(showMuon) TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1, AliMagF::k5kG)); // Get the pointer to gGeoManager before it's broken (bug in alieve) fGeoManager = gGeoManager; // -- Initialize pointsets and add macros // ---------------------------------------- //TEveUtil::LoadMacro("hlt_alieve_init.C"); //hlt_alieve_init(".", -1); // -- Initialize Eve // ------------------- initializeEveViewer( showBarrel, showMuon ); // -- Reset gGeoManager to the original pointer // ---------------------------------------------- // -- Finalize Eve // ----------------- gSystem->ProcessEvents(); gEve->Redraw3D(kTRUE); // -- Create new hM object // ------------------------- gEventManager = new AliEveHLTEventManagerOffline("/home/slindal/alice/data/PbPb/AliESDs.root"); gEventManager->SetEveManager(gEve); gEventManager->SetGeoManager(gGeoManager); gEventManager->SetRPhiManager(gRPhiMgr); gEventManager->SetRPhiEventScene(gRPhiEventScene); gEventManager->SetRPhiViewer(gRPhiView); gEventManager->SetRhoZManager(gRhoZMgr); gEventManager->SetRhoZEventScene(gRhoZEventScene); gEventManager->SetRhoZViewer(gRhoZView); //gEventManager->SetBarrelFlag(showBarrel); //gEventManager->SetMuonFlag(showMuon); // -- Add hM to EveTree // ---------------------- gEve->AddToListTree(gEventManager, kTRUE); // -- Create SourceList // ---------------------- gGeoManager = fGeoManager; } // ------------------------------------------------------------------------- Int_t initializeEveViewer( Bool_t showBarrel, Bool_t showMuon ) { //============================================================================= // Visualization database //============================================================================ TEveUtil::AssertMacro("VizDB_scan.C"); // alieve_vizdb(); //============================================================================== // -- Geometry, scenes, projections and viewers //============================================================================== TEveBrowser *browser = gEve->GetBrowser(); browser->ShowCloseTab(kFALSE); // -- Disable extra geometry // --------------------------- if (!showMuon) gShowMUON = gShowMUONRPhi = gShowMUONRhoZ = kFALSE; // -- Load Geometry // ------------------ TEveUtil::LoadMacro("geom_gentle_hlt.C"); gGeomGentle = geom_gentle_hlt(); gGeomGentleRPhi = geom_gentle_rphi(); gGeomGentleRPhi->IncDenyDestroy(); gGeomGentleRhoZ = geom_gentle_rhoz(); gGeomGentleRhoZ->IncDenyDestroy(); gGeomGentleTRD = geom_gentle_trd(); gGeoManager = fGeoManager; gEMCALNode = gGeoManager->GetTopVolume()->FindNode("XEN1_1"); TEveGeoTopNode* emcal_re = new TEveGeoTopNode(gGeoManager, gEMCALNode); gEve->AddGlobalElement(emcal_re); gEve->Redraw3D(); if (gShowMUON) gGeomGentleMUON = geom_gentle_muon(kFALSE); // -- 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->SetCurrentDepth(-10); gRPhiMgr->ImportElements(gGeomGentleRPhi, gRPhiGeomScene); gRPhiMgr->SetCurrentDepth(0); 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->SetCurrentDepth(-10); gRhoZMgr->ImportElements(gGeomGentleRhoZ, gRhoZGeomScene); gRhoZMgr->SetCurrentDepth(0); gRhoZMgr->ImportElements(gGeomGentleTRD, gRhoZGeomScene); if (gShowMUONRhoZ) gRhoZMgr->ImportElements(gGeomGentleMUON, gRhoZGeomScene); // -- Viewers // ------------ TEveWindowSlot *slot = TEveWindow::CreateWindowInTab(browser->GetTabRight()); TEveWindowPack *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); //Add HLT Text to windows TGLOverlayButton *ob1 = new TGLOverlayButton(g3DView->GetGLViewer(), "HLT", 0, 20, 110, 60); ob1->SetAlphaValues(0.8, 0.8); // cout << "color" << ob1->GetBackColor() << endl; //ob1->SetBackColor(8421631); //ob1->SetBackColor(10492431); TGLOverlayButton *ob2 = new TGLOverlayButton(g3DView->GetGLViewer(), "ALICE", 0, 0, 110, 20); ob2->SetAlphaValues(0.8, 0.8); //ob2->SetBackColor(0.2); TGLOverlayButton *ob3 = new TGLOverlayButton(gEve->GetDefaultGLViewer(), "HLT", 0, 20, 110, 60); ob3->SetAlphaValues(0.8, 0.8); TGLOverlayButton *ob4 = new TGLOverlayButton(gEve->GetDefaultGLViewer(), "ALICE", 0, 0, 110, 20); ob4->SetAlphaValues(0.8, 0.8); TGLOverlayButton *ne = new TGLOverlayButton(gEve->GetDefaultGLViewer(), "Next Event", 110, 0, 210, 20); ne->SetAlphaValues(0.0, 0.8); // -- List of Viewers // -------------------- TEveViewerList *viewerlist = new TEveViewerList(); viewerlist->AddElement(gEve->GetDefaultViewer()); viewerlist->AddElement(g3DView); viewerlist->AddElement(gRhoZView); viewerlist->AddElement(gRPhiView); viewerlist->SwitchColorSet(); //============================================================================== // -- Macros / QA histograms //============================================================================== // -- Registration of per-event macros // ------------------------------------- AliEveMacroExecutor *exec = new AliEveMacroExecutor(); gStyle->SetPalette(1, 0); gStyle->SetOptFit(1); return 0; }