FWTrackHitsDetailView.cc

Go to the documentation of this file.

// ROOT includes
#include "TGLFontManager.h"
#include "TEveScene.h"
#include "TEveManager.h"
#include "TEveStraightLineSet.h"
#include "TEveTrack.h"
#include "TEveTrackPropagator.h"
#include "TEveTrans.h"
#include "TEveText.h"
#include "TEveGeoShape.h"
#include "TGSlider.h"
#include "TGButton.h"
#include "TGLabel.h"
#include "TGLViewer.h"
#include "TCanvas.h"
#include "TLatex.h"
#include "TLegend.h"

// CMSSW includes
#include "DataFormats/SiStripDetId/interface/SiStripDetId.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/MuonDetId/interface/MuonSubdetId.h"
#include "DataFormats/MuonDetId/interface/DTChamberId.h"

// Fireworks includes
#include "Fireworks/Core/interface/FWModelId.h"
#include "Fireworks/Core/interface/FWColorManager.h"
#include "Fireworks/Core/interface/FWGeometry.h"
#include "Fireworks/Core/interface/FWEventItem.h"
#include "Fireworks/Core/interface/CSGAction.h"
#include "Fireworks/Core/interface/FWIntValueListener.h"
#include "Fireworks/Core/interface/FWMagField.h"
#include "Fireworks/Core/interface/fwLog.h"

#include "Fireworks/Tracks/plugins/FWTrackHitsDetailView.h"
#include "Fireworks/Tracks/interface/TrackUtils.h"

#include <functional>

FWTrackHitsDetailView::FWTrackHitsDetailView()
    : m_modules(nullptr),
      m_moduleLabels(nullptr),
      m_hits(nullptr),
      m_slider(nullptr),
      m_sliderListener(),
      m_legend(nullptr) {}

FWTrackHitsDetailView::~FWTrackHitsDetailView() {}

void FWTrackHitsDetailView::build(const FWModelId& id, const reco::Track* track) {
  {
    TGCompositeFrame* f = new TGVerticalFrame(m_guiFrame);
    m_guiFrame->AddFrame(f);
    f->AddFrame(new TGLabel(f, "Module Transparency:"), new TGLayoutHints(kLHintsLeft, 2, 2, 0, 0));
    m_slider = new TGHSlider(f, 120, kSlider1 | kScaleNo);
    f->AddFrame(m_slider, new TGLayoutHints(kLHintsTop | kLHintsLeft, 2, 2, 1, 4));
    m_slider->SetRange(0, 100);
    m_slider->SetPosition(75);

    m_sliderListener = new FWIntValueListener();
    TQObject::Connect(
        m_slider, "PositionChanged(Int_t)", "FWIntValueListenerBase", m_sliderListener, "setValue(Int_t)");
    m_sliderListener->valueChanged_.connect(
        std::bind(&FWTrackHitsDetailView::transparencyChanged, this, std::placeholders::_1));
  }

  {
    CSGAction* action = new CSGAction(this, "Show Module Labels");
    TGCheckButton* b = new TGCheckButton(m_guiFrame, "Show Module Labels");
    b->SetState(kButtonUp, false);
    m_guiFrame->AddFrame(b, new TGLayoutHints(kLHintsNormal, 2, 3, 1, 4));
    TQObject::Connect(b, "Clicked()", "CSGAction", action, "activate()");
    action->activated.connect(sigc::mem_fun(*this, &FWTrackHitsDetailView::rnrLabels));
  }
  {
    CSGAction* action = new CSGAction(this, " Pick Camera Center ");
    action->createTextButton(m_guiFrame, new TGLayoutHints(kLHintsNormal, 2, 0, 1, 4));
    action->setToolTip("Click on object in viewer to set camera center.");
    action->activated.connect(sigc::mem_fun(*this, &FWTrackHitsDetailView::pickCameraCenter));
  }
  makeLegend();

  TGCompositeFrame* p = (TGCompositeFrame*)m_guiFrame->GetParent();
  p->MapSubwindows();
  p->Layout();

  m_modules = new TEveElementList("Modules");
  m_eveScene->AddElement(m_modules);
  m_moduleLabels = new TEveElementList("Modules");
  m_eveScene->AddElement(m_moduleLabels);
  m_hits = new TEveElementList("Hits");
  m_eveScene->AddElement(m_hits);
  if (track->extra().isAvailable()) {
    addModules(*track, id.item(), m_modules, true);
    addHits(*track, id.item(), m_hits, true);
  }
  for (TEveElement::List_i i = m_modules->BeginChildren(), end = m_modules->EndChildren(); i != end; ++i) {
    TEveGeoShape* gs = dynamic_cast<TEveGeoShape*>(*i);
    const auto& rhs = *(*(i));
    if (gs == nullptr && (*i != nullptr)) {
      std::cerr << "Got a " << typeid(rhs).name() << ", expecting TEveGeoShape. ignoring (it must be the clusters)."
                << std::endl;
      continue;
    }
    gs->SetMainTransparency(75);
    gs->SetPickable(kFALSE);

    TString name = gs->GetElementTitle();
    if (!name.Contains("BAD") && !name.Contains("INACTIVE") && !name.Contains("LOST")) {
      gs->SetMainColor(kBlue);
    }
    TEveText* text = new TEveText(name.Data());
    text->PtrMainTrans()->SetFrom(gs->RefMainTrans().Array());
    text->SetFontMode(TGLFont::kPixmap);
    text->SetFontSize(12);
    m_moduleLabels->AddElement(text);
  }
  m_moduleLabels->SetRnrChildren(false);

  TEveTrackPropagator* prop = new TEveTrackPropagator();
  prop->SetMagFieldObj(item()->context().getField(), false);
  prop->SetStepper(TEveTrackPropagator::kRungeKutta);
  prop->SetMaxR(123);
  prop->SetMaxZ(300);
  prop->SetMaxStep(1);
  TEveTrack* trk = fireworks::prepareTrack(*track, prop);
  trk->MakeTrack();
  trk->SetLineWidth(2);
  trk->SetTitle("Track and its ref states");
  prop->SetRnrDaughters(kTRUE);
  prop->SetRnrReferences(kTRUE);
  prop->SetRnrDecay(kTRUE);
  prop->SetRnrFV(kTRUE);
  trk->SetMainColor(id.item()->defaultDisplayProperties().color());
  m_eveScene->AddElement(trk);

  viewerGL()->SetStyle(TGLRnrCtx::kOutline);
  viewerGL()->SetDrawCameraCenter(kTRUE);
  viewerGL()->ResetCamerasAfterNextUpdate();
  viewerGL()->UpdateScene(kFALSE);
  gEve->Redraw3D();

  setTextInfo(id, track);
}

void FWTrackHitsDetailView::setBackgroundColor(Color_t col) {
  // Callback for cmsShow change of background

  FWColorManager::setColorSetViewer(viewerGL(), col);

  // adopt label colors to background, this should be implemneted in TEveText
  if (m_moduleLabels) {
    Color_t x = viewerGL()->GetRnrCtx()->ColorSet().Foreground().GetColorIndex();
    for (TEveElement::List_i i = m_moduleLabels->BeginChildren(); i != m_moduleLabels->EndChildren(); ++i)
      (*i)->SetMainColor(x);
  }
}

void FWTrackHitsDetailView::pickCameraCenter() { viewerGL()->PickCameraCenter(); }

void FWTrackHitsDetailView::transparencyChanged(int x) {
  for (TEveElement::List_i i = m_modules->BeginChildren(); i != m_modules->EndChildren(); ++i) {
    (*i)->SetMainTransparency(x);
  }
  gEve->Redraw3D();
}

void FWTrackHitsDetailView::setTextInfo(const FWModelId& id, const reco::Track* track) {
  m_infoCanvas->cd();

  float_t x = 0.02;
  float y = 0.95;

  TLatex* latex = new TLatex(x, y, "");
  const double textsize(0.07);
  latex->SetTextSize(2 * textsize);

  latex->DrawLatex(x, y, id.item()->modelName(id.index()).c_str());
  y -= latex->GetTextSize() * 0.6;

  latex->SetTextSize(textsize);
  float lineH = latex->GetTextSize() * 0.6;

  latex->DrawLatex(
      x, y, Form(" P_{T} = %.1f GeV, #eta = %0.2f, #varphi = %0.2f", track->pt(), track->eta(), track->phi()));
  y -= lineH;

  if (track->charge() > 0)
    latex->DrawLatex(x, y, " charge = +1");
  else
    latex->DrawLatex(x, y, " charge = -1");
  y -= lineH;
  y -= lineH;

  latex->DrawLatex(x, y, "Track modules:");
  y -= lineH;

  Double_t pos[4];
  pos[0] = x + 0.05;
  pos[2] = x + 0.20;
  Double_t boxH = 0.25 * textsize;

  pos[1] = y;
  pos[3] = pos[1] + boxH;
  FWDetailViewBase::drawCanvasBox(pos, kBlue);
  latex->DrawLatex(x + 0.25, y, "Module");
  y -= lineH;

  pos[1] = y;
  pos[3] = pos[1] + boxH;
  FWDetailViewBase::drawCanvasBox(pos, kRed);
  latex->DrawLatex(x + 0.25, y, "LOST Module");
  y -= lineH;

  pos[1] = y;
  pos[3] = pos[1] + boxH;
  FWDetailViewBase::drawCanvasBox(pos, 28);
  latex->DrawLatex(x + 0.25, y, "INACTIVE Module");
  y -= lineH;

  pos[1] = y;
  pos[3] = pos[1] + boxH;
  FWDetailViewBase::drawCanvasBox(pos, 218);
  latex->DrawLatex(x + 0.25, y, "BAD Module");
  y -= lineH;

  Float_t r = 0.01;
  Float_t r2 = 0.02;
  y -= lineH;
  drawCanvasDot(x + r2, y, r2, kGreen);
  y -= r;
  latex->DrawLatex(x + 3 * r2, y, "Pixel Hits");
  y -= lineH;

  drawCanvasDot(x + r2, y, r2, kRed);
  y -= r;
  latex->DrawLatex(x + 3 * r2, y, "Extra Pixel Hits");
  y -= lineH;

  m_legend->SetY2(y);
  m_legend->Draw();
  m_legend = nullptr;  // Deleted together with TPad.
}

void FWTrackHitsDetailView::makeLegend(void) {
  m_legend = new TLegend(0.01, 0.01, 0.99, 0.99, nullptr, "NDC");
  m_legend->SetFillColor(kWhite);
  m_legend->SetTextSize(0.07);
  m_legend->SetBorderSize(0);
  m_legend->SetMargin(0.15);
  m_legend->SetEntrySeparation(0.01);

  TEveStraightLineSet* legend = new TEveStraightLineSet("siStripCluster");
  legend->SetLineWidth(3);
  legend->SetLineColor(kGreen);
  m_legend->AddEntry(legend, "Exact SiStripCluster", "l");

  TEveStraightLineSet* legend2 = new TEveStraightLineSet("siStripCluster2");
  legend2->SetLineWidth(3);
  legend2->SetLineColor(kRed);
  m_legend->AddEntry(legend2, "Extra SiStripCluster", "l");

  TEveStraightLineSet* legend3 = new TEveStraightLineSet("refStates");
  legend3->SetDepthTest(kFALSE);
  legend3->SetMarkerColor(kYellow);
  legend3->SetMarkerStyle(kPlus);
  legend3->SetMarkerSize(2);
  m_legend->AddEntry(legend3, "Inner/Outermost States", "p");

  TEveStraightLineSet* legend4 = new TEveStraightLineSet("vertex");
  legend4->SetDepthTest(kFALSE);
  legend4->SetMarkerColor(kRed);
  legend4->SetMarkerStyle(kFullDotLarge);
  legend4->SetMarkerSize(2);
  m_legend->AddEntry(legend4, "Vertex", "p");

  TEveStraightLineSet* legend5 = new TEveStraightLineSet("cameraCenter");
  legend5->SetDepthTest(kFALSE);
  legend5->SetMarkerColor(kCyan);
  legend5->SetMarkerStyle(kFullDotLarge);
  legend5->SetMarkerSize(2);
  m_legend->AddEntry(legend5, "Camera center", "p");
}

void FWTrackHitsDetailView::addTrackerHits3D(std::vector<TVector3>& points,
                                             class TEveElementList* tList,
                                             Color_t color,
                                             int size) {
  // !AT this is  detail view specific, should move to track hits
  // detail view

  TEvePointSet* pointSet = new TEvePointSet();
  pointSet->SetMarkerSize(size);
  pointSet->SetMarkerStyle(4);
  pointSet->SetPickable(kTRUE);
  pointSet->SetTitle("Pixel Hits");
  pointSet->SetMarkerColor(color);

  for (std::vector<TVector3>::const_iterator it = points.begin(), itEnd = points.end(); it != itEnd; ++it) {
    pointSet->SetNextPoint(it->x(), it->y(), it->z());
  }
  tList->AddElement(pointSet);
}

void FWTrackHitsDetailView::addHits(const reco::Track& track,
                                    const FWEventItem* iItem,
                                    TEveElement* trkList,
                                    bool addNearbyHits) {
  std::vector<TVector3> pixelPoints;
  fireworks::pushPixelHits(pixelPoints, *iItem, track);
  TEveElementList* pixels = new TEveElementList("Pixels");
  trkList->AddElement(pixels);
  if (addNearbyHits) {
    // get the extra hits
    std::vector<TVector3> pixelExtraPoints;
    fireworks::pushNearbyPixelHits(pixelExtraPoints, *iItem, track);
    // draw first the others
    addTrackerHits3D(pixelExtraPoints, pixels, kRed, 1);
    // then the good ones, so they're on top
    addTrackerHits3D(pixelPoints, pixels, kGreen, 1);
  } else {
    // just add those points with the default color
    addTrackerHits3D(pixelPoints, pixels, iItem->defaultDisplayProperties().color(), 1);
  }

  // strips
  TEveElementList* strips = new TEveElementList("Strips");
  trkList->AddElement(strips);
  fireworks::addSiStripClusters(iItem, track, strips, addNearbyHits, false);
}

//______________________________________________________________________________

void FWTrackHitsDetailView::addModules(const reco::Track& track,
                                       const FWEventItem* iItem,
                                       TEveElement* trkList,
                                       bool addLostHits) {
  std::set<unsigned int> ids;
  for (trackingRecHit_iterator recIt = track.recHitsBegin(), recItEnd = track.recHitsEnd(); recIt != recItEnd;
       ++recIt) {
    DetId detid = (*recIt)->geographicalId();
    if (!addLostHits && !(*recIt)->isValid())
      continue;
    if (detid.rawId() != 0) {
      TString name("");
      switch (detid.det()) {
        case DetId::Tracker:
          name = iItem->getGeom()->getTrackerTopology()->print(detid);
          break;

        case DetId::Muon:
          switch (detid.subdetId()) {
            case MuonSubdetId::DT:
              name = "DT";
              detid = DetId(DTChamberId(detid));  // get rid of layer bits
              break;
            case MuonSubdetId::CSC:
              name = "CSC";
              break;
            case MuonSubdetId::RPC:
              name = "RPC";
              break;
            case MuonSubdetId::GEM:
              name = "GEM";
              break;
            case MuonSubdetId::ME0:
              name = "ME0";
              break;
            default:
              break;
          }
          break;
        default:
          break;
      }
      if (!ids.insert(detid.rawId()).second)
        continue;
      if (iItem->getGeom()) {
        TEveGeoShape* shape = iItem->getGeom()->getEveShape(detid);
        if (nullptr != shape) {
          shape->SetMainTransparency(65);
          shape->SetPickable(kTRUE);
          switch ((*recIt)->type()) {
            case TrackingRecHit::valid:
              shape->SetMainColor(iItem->defaultDisplayProperties().color());
              break;
            case TrackingRecHit::missing:
            case TrackingRecHit::missing_inner:
            case TrackingRecHit::missing_outer:
              name += "LOST ";
              shape->SetMainColor(kRed);
              break;
            case TrackingRecHit::inactive_inner:
            case TrackingRecHit::inactive_outer:
            case TrackingRecHit::inactive:
              name += "INACTIVE ";
              shape->SetMainColor(28);
              break;
            case TrackingRecHit::bad:
              name += "BAD ";
              shape->SetMainColor(218);
              break;
          }
          shape->SetTitle(name + ULong_t(detid.rawId()));
          trkList->AddElement(shape);
        } else {
          fwLog(fwlog::kInfo) << "Failed to get shape extract for a tracking rec hit: "
                              << "\n"
                              << fireworks::info(detid) << std::endl;
        }
      }
    }
  }
}

void FWTrackHitsDetailView::rnrLabels() {
  m_moduleLabels->SetRnrChildren(!m_moduleLabels->GetRnrChildren());
  gEve->Redraw3D();
}

REGISTER_FWDETAILVIEW(FWTrackHitsDetailView, Hits);