event_demo.C File Reference

Detailed Description

This example display geometry, tracks and hits in web browser.

 
#include <vector>
#include <string>
#include <iostream>
 
#include "TClass.h"
#include "TRandom.h"
#include "TGeoTube.h"
#include "TGeoSphere.h"
#include "TParticle.h"
#include "TApplication.h"
#include "TMatrixDSym.h"
#include "TVector.h"
#include "TMatrixDEigen.h"
 
#include <ROOT/REveGeoShape.hxx>
#include <ROOT/REveScene.hxx>
#include <ROOT/REveViewer.hxx>
#include <ROOT/REveElement.hxx>
#include <ROOT/REveCompound.hxx>
#include <ROOT/REveManager.hxx>
#include <ROOT/REveUtil.hxx>
#include <ROOT/REveGeoShape.hxx>
#include <ROOT/REveProjectionManager.hxx>
#include <ROOT/REveProjectionBases.hxx>
#include <ROOT/REvePointSet.hxx>
#include <ROOT/REveJetCone.hxx>
#include <ROOT/REveTrans.hxx>
#include <ROOT/REveTrack.hxx>
#include <ROOT/REveTrackPropagator.hxx>
 
namespace REX = ROOT::Experimental;
 
// globals
REX::REveManager *eveMng = nullptr;
REX::REveProjectionManager *mngRhoPhi = nullptr;
REX::REveProjectionManager *mngRhoZ = nullptr;
REX::REveScene *rPhiGeomScene = nullptr, *rPhiEventScene = nullptr;
REX::REveScene *rhoZGeomScene = nullptr, *rhoZEventScene = nullptr;
REX::REveViewer *rphiView = nullptr;
REX::REveViewer *rhoZView = nullptr;
 
const Double_t kR_min = 240;
const Double_t kR_max = 250;
const Double_t kZ_d = 300;
 
REX::REvePointSet *getPointSet(int npoints = 2, float s = 2, int color = 28)
{
   TRandom &r = *gRandom;
 
   auto ps = new REX::REvePointSet("fu", "", npoints);
 
   for (Int_t i = 0; i < npoints; ++i)
      ps->SetNextPoint(r.Uniform(-s, s), r.Uniform(-s, s), r.Uniform(-s, s));
 
   ps->SetMarkerColor(color);
   ps->SetMarkerSize(8);
   ps->SetMarkerStyle(4);
   return ps;
}
 
void addPoints()
{
   REX::REveElement *event = eveMng->GetEventScene();
 
   auto pntHolder = new REX::REveCompound("Hits");
 
   auto ps1 = getPointSet(20, 100);
   ps1->SetName("Points_1");
   ps1->SetTitle("Points_1 title"); // used as tooltip
 
   pntHolder->AddElement(ps1);
 
   auto ps2 = getPointSet(10, 200, 4);
   ps2->SetName("Points_2");
   ps2->SetTitle("Points_2 title"); // used as tooltip
   ps2->SetAlwaysSecSelect(true);
   pntHolder->AddElement(ps2);
 
   event->AddElement(pntHolder);
}
 
void addTracks()
{
   TRandom &r = *gRandom;
 
   REX::REveElement *event = eveMng->GetEventScene();
   auto prop = new REX::REveTrackPropagator();
   prop->SetMagFieldObj(new REX::REveMagFieldDuo(350, 3.5, -2.0));
   prop->SetMaxR(300);
   prop->SetMaxZ(600);
   prop->SetMaxOrbs(6);
 
   auto trackHolder = new REX::REveCompound("Tracks");
 
   double v = 0.2;
   double m = 5;
 
   int N_Tracks = 10 + r.Integer(20);
   for (int i = 0; i < N_Tracks; i++) {
      TParticle p;
 
      int pdg = 11 * (r.Integer(2) > 0 ? 1 : -1);
      p.SetPdgCode(pdg);
 
      p.SetProductionVertex(r.Uniform(-v, v), r.Uniform(-v, v), r.Uniform(-v, v), 1);
      p.SetMomentum(r.Uniform(-m, m), r.Uniform(-m, m), r.Uniform(-m, m) * r.Uniform(1, 3), 1);
      auto track = new REX::REveTrack(&p, 1, prop);
      track->MakeTrack();
      if (i % 4 == 3)
         track->SetLineStyle(kDashed); // enabled dashed style for some tracks
      track->SetMainColor(kBlue);
      track->SetName(Form("RandomTrack_%d", i));
      track->SetTitle(Form("RandomTrack_%d title", i)); // used as tooltip
      trackHolder->AddElement(track);
   }
 
   event->AddElement(trackHolder);
}
 
void addJets()
{
   TRandom &r = *gRandom;
 
   REX::REveElement *event = eveMng->GetEventScene();
   auto jetHolder = new REX::REveCompound("Jets");
 
   int N_Jets = 5 + r.Integer(5);
   for (int i = 0; i < N_Jets; i++) {
      auto jet = new REX::REveJetCone(Form("Jet_%d", i));
      jet->SetTitle(Form("Jet_%d\n  pT = %.2f", i, r.Uniform(1, 40))); // used as tooltip
      jet->SetCylinder(2 * kR_max, 2 * kZ_d);
      jet->AddEllipticCone(r.Uniform(-3.5, 3.5), r.Uniform(0, TMath::TwoPi()), r.Uniform(0.02, 0.2),
                           r.Uniform(0.02, 0.3));
      jet->SetFillColor(kPink - 8);
      jet->SetLineColor(kBlack);
 
      jetHolder->AddElement(jet);
   }
   event->AddElement(jetHolder);
}
 
void makeEventScene()
{
   addPoints();
   addTracks();
   addJets();
}
 
void makeGeometryScene()
{
   auto b1 = new REX::REveGeoShape("Barrel 1");
   b1->SetShape(new TGeoTube(kR_min, kR_max, kZ_d));
   b1->SetMainColor(kCyan);
   b1->SetNSegments(80);
   b1->SetMainTransparency(70);
   eveMng->GetGlobalScene()->AddElement(b1);
 
   // Debug of surface fill in RPhi (index buffer screwed).
   // b1->SetNSegments(3);
   b1->SetNSegments(40);
 
   // an example of axis guides
   eveMng->GetDefaultViewer()->SetAxesType(REX::REveViewer::EAxesType::kAxesEdge);
}
 
void createProjectionStuff()
{
   // project RhoPhi
   rPhiGeomScene = eveMng->SpawnNewScene("RPhi Geometry", "RPhi");
   rPhiEventScene = eveMng->SpawnNewScene("RPhi Event Data", "RPhi");
 
   mngRhoPhi = new REX::REveProjectionManager(REX::REveProjection::kPT_RPhi);
 
   rphiView = eveMng->SpawnNewViewer("RPhi View", "");
   rphiView->SetCameraType(REX::REveViewer::kCameraOrthoXOY);
   rphiView->AddScene(rPhiGeomScene);
   rphiView->AddScene(rPhiEventScene);
 
   // ----------------------------------------------------------------
 
   rhoZGeomScene = eveMng->SpawnNewScene("RhoZ Geometry", "RhoZ");
   rhoZEventScene = eveMng->SpawnNewScene("RhoZ Event Data", "RhoZ");
 
   mngRhoZ = new REX::REveProjectionManager(REX::REveProjection::kPT_RhoZ);
 
   rhoZView = eveMng->SpawnNewViewer("RhoZ View", "");
   rhoZView->SetCameraType(REX::REveViewer::kCameraOrthoXOY);
   rhoZView->AddScene(rhoZGeomScene);
   rhoZView->AddScene(rhoZEventScene);
}
 
void projectScenes(bool geomp, bool eventp)
{
   if (geomp) {
      for (auto &ie : eveMng->GetGlobalScene()->RefChildren()) {
         mngRhoPhi->SetCurrentDepth(0);
         mngRhoPhi->ImportElements(ie, rPhiGeomScene);
         mngRhoZ->SetCurrentDepth(0);
         mngRhoZ->ImportElements(ie, rhoZGeomScene);
      }
   }
   if (eventp) {
      int depth = 50;
      for (auto &ie : eveMng->GetEventScene()->RefChildren()) {
         mngRhoPhi->SetCurrentDepth(depth);
         mngRhoPhi->ImportElements(ie, rPhiEventScene);
         mngRhoZ->SetCurrentDepth(depth);
         mngRhoZ->ImportElements(ie, rhoZEventScene);
         depth -= 10;
      }
   }
 
   // auto t0 = eveMng->GetEventScene()->FindChild("Tracks")->FirstChild();
   // printf("t0=%p, %s %s\n", t0, t0->GetElementName(), t0->IsA()->GetName());
   // dynamic_cast<REX::REveTrack*>(t0)->Print("all");
 
   // auto t1 = rPhiEventScene->FindChild("Tracks [P]")->FirstChild();
   // printf("t1=%p, %s %s\n", t1, t1->GetElementName(), t1->IsA()->GetName());
   // dynamic_cast<REX::REveTrack*>(t1)->Print("all");
}
 
//==============================================================================
 
class EventManager : public REX::REveElement {
private:
   bool fAutoplay{false};
   int fPlayDelay{10};
   int fCount{0};
   std::chrono::duration<double> fDeltaTime{1};
 
   std::thread *fTimerThread{nullptr};
   std::mutex fMutex;
   std::condition_variable fCV;
 
public:
   EventManager()
   {
      std::chrono::milliseconds ms(100);
      fDeltaTime = ms;
   }
 
   ~EventManager() override {}
 
   void NextEvent()
   {
      auto scene = eveMng->GetEventScene();
      scene->DestroyElements();
      makeEventScene();
      projectScenes(false, true);
      // if (++fCount % 10 == 0) printf("At event %d\n", fCount);
   }
 
   void autoplay_scheduler()
   {
      while (true) {
         bool autoplay;
         {
            std::unique_lock<std::mutex> lock{fMutex};
            if (!fAutoplay) {
               // printf("exit thread pre wait\n");
               return;
            }
            if (fCV.wait_for(lock, fDeltaTime) != std::cv_status::timeout) {
               printf("autoplay not timed out \n");
               if (!fAutoplay) {
                  printf("exit thread post wait\n");
                  return;
               } else {
                  continue;
               }
            }
            autoplay = fAutoplay;
         }
         if (autoplay) {
            REX::REveManager::ChangeGuard ch;
            NextEvent();
         } else {
            return;
         }
      }
   }
 
   void Autoplay()
   {
      static std::mutex autoplay_mutex;
      std::unique_lock<std::mutex> aplock{autoplay_mutex};
      {
         std::unique_lock<std::mutex> lock{fMutex};
         fAutoplay = !fAutoplay;
         if (fAutoplay) {
            if (fTimerThread) {
               fTimerThread->join();
               delete fTimerThread;
               fTimerThread = nullptr;
            }
            NextEvent();
            fTimerThread = new std::thread{[this] { autoplay_scheduler(); }};
         } else {
            fCV.notify_all();
         }
      }
   }
 
   virtual void QuitRoot()
   {
      printf("Quit ROOT\n");
      REX::REveManager::QuitRoot();
   }
};
 
void event_demo()
{
   // disable browser cache - all scripts and html files will be loaded every time, useful for development
   // gEnv->SetValue("WebGui.HttpMaxAge", 0);
 
   gRandom->SetSeed(0); // make random seed
 
   eveMng = REX::REveManager::Create();
   eveMng->AllowMultipleRemoteConnections(false, false);
 
   auto eventMng = new EventManager();
   eventMng->SetName("EventManager");
   eveMng->GetWorld()->AddElement(eventMng);
 
   eveMng->GetWorld()->AddCommand("QuitRoot", "sap-icon://log", eventMng, "QuitRoot()");
 
   eveMng->GetWorld()->AddCommand("NextEvent", "sap-icon://step", eventMng, "NextEvent()");
 
   eveMng->GetWorld()->AddCommand("Autoplay", "sap-icon://refresh", eventMng, "Autoplay()");
 
   makeGeometryScene();
   makeEventScene();
 
   if (true) {
      createProjectionStuff();
      projectScenes(true, true);
   }
 
   eveMng->Show();
}

Author

Alja Mrak-Tadel

Definition in file event_demo.C.