evdb_tool::OpHit3DDrawer Class Reference

Inheritance diagram for evdb_tool::OpHit3DDrawer:

Public Member Functions
	OpHit3DDrawer (const fhicl::ParameterSet &)
void	Draw (const art::Event &, evdb::View3D *) const override

Private Member Functions
void	DrawRectangularBox (evdb::View3D *, const Eigen::Vector3f &, const Eigen::Vector3f &, int, int, int) const

Detailed Description

Definition at line 26 of file OpHit3DDrawer_tool.cc.

Constructor & Destructor Documentation

evdb_tool::OpHit3DDrawer::OpHit3DDrawer ( const fhicl::ParameterSet &  )

inlineexplicit

Definition at line 28 of file OpHit3DDrawer_tool.cc.

References Draw(), and DrawRectangularBox().

{}

Member Function Documentation

void evdb_tool::OpHit3DDrawer::Draw ( const art::Event &  event, evdb::View3D *  view  ) const

overridevirtual

Implements evdb_tool::I3DDrawer.

Definition at line 41 of file OpHit3DDrawer_tool.cc.

References evd::ColorDrawingOptions::CalQ(), DrawRectangularBox(), evd::RecoDrawingOptions::fDrawOpHits, evd::RecoDrawingOptions::fFlashMinPE, evd::RecoDrawingOptions::fFlashTMax, evd::RecoDrawingOptions::fFlashTMin, evd::RecoDrawingOptions::fOpHitLabels, evd::ColorDrawingOptions::fRecoQHigh, evd::ColorDrawingOptions::fRecoQLow, Get, geo::OpDetGeo::GetCenter(), evdb::ColorScale::GetColor(), geo::OpDetGeo::HalfH(), geo::OpDetGeo::HalfW(), art::Handle< T >::isValid(), and geo::kCollection.

Referenced by OpHit3DDrawer().

  {
    art::ServiceHandle<evd::RecoDrawingOptions> recoOpt;

    if (recoOpt->fDrawOpHits == 0) return;

    // Service recovery
    art::ServiceHandle<evd::ColorDrawingOptions> cst;

    art::Handle<std::vector<recob::OpHit>> opHitHandle;

    // This seems like a time waster but we want to get the full color scale for all OpHits... so loops away...
    std::vector<float> opHitPEVec;

    // This is almost identically the same for loop we will re-excute below... sigh...
    // But the idea is to get a min/max range for the drawing colors...
    for (size_t idx = 0; idx < recoOpt->fOpHitLabels.size(); idx++) {
      art::InputTag opHitProducer = recoOpt->fOpHitLabels[idx];

      event.getByLabel(opHitProducer, opHitHandle);

      if (!opHitHandle.isValid()) continue;
      if (opHitHandle->size() == 0) continue;

      // Start the loop over flashes
      for (const auto& opHit : *opHitHandle) {
        // Make some selections...
        if (opHit.PE() < recoOpt->fFlashMinPE) continue;
        if (opHit.PeakTime() < recoOpt->fFlashTMin) continue;
        if (opHit.PeakTime() > recoOpt->fFlashTMax) continue;

        opHitPEVec.push_back(opHit.PE());
      }
    }

    // Do we have any flashes and hits?
    if (!opHitPEVec.empty()) {
      // Sorting is good for mind and body...
      std::sort(opHitPEVec.begin(), opHitPEVec.end());

      float minTotalPE = opHitPEVec.front();
      float maxTotalPE = opHitPEVec[0.9 * opHitPEVec.size()];

      // Now we can set the scaling factor for PE
      float opHitPEScale((cst->fRecoQHigh[geo::kCollection] - cst->fRecoQLow[geo::kCollection]) /
                         (maxTotalPE - minTotalPE));

      // We are meant to draw the flashes/hits, so loop over the list of input flashes
      auto const& wireReadoutGeom = art::ServiceHandle<geo::WireReadout const>()->Get();
      for (size_t idx = 0; idx < recoOpt->fOpHitLabels.size(); idx++) {
        art::InputTag opHitProducer = recoOpt->fOpHitLabels[idx];

        event.getByLabel(opHitProducer, opHitHandle);

        if (!opHitHandle.isValid()) continue;
        if (opHitHandle->size() == 0) continue;

        // Start the loop over flashes
        for (const auto& opHit : *opHitHandle) {
          // Make some selections...
          if (opHit.PE() < recoOpt->fFlashMinPE) continue;
          if (opHit.PeakTime() < recoOpt->fFlashTMin) continue;
          if (opHit.PeakTime() > recoOpt->fFlashTMax) continue;

          unsigned int opChannel = opHit.OpChannel();
          const geo::OpDetGeo& opHitGeo = wireReadoutGeom.OpDetGeoFromOpChannel(opChannel);
          const geo::Point_t& opHitPos = opHitGeo.GetCenter();
          float xWidth = opHit.Width();
          float zWidth = opHitGeo.HalfW();
          float yWidth = opHitGeo.HalfH();

          Eigen::Vector3f opHitLo(
            opHitPos.X() - xWidth, opHitPos.Y() - yWidth, opHitPos.Z() - zWidth);
          Eigen::Vector3f opHitHi(
            opHitPos.X() + xWidth, opHitPos.Y() + yWidth, opHitPos.Z() + zWidth);

          float peFactor = cst->fRecoQLow[geo::kCollection] +
                           opHitPEScale * std::min(maxTotalPE, float(opHit.PE()));

          int chargeColorIdx = cst->CalQ(geo::kCollection).GetColor(peFactor);

          DrawRectangularBox(view, opHitLo, opHitHi, chargeColorIdx, 2, 1);
        }
      }
    }
  }

void evdb_tool::OpHit3DDrawer::DrawRectangularBox ( evdb::View3D *  view, const Eigen::Vector3f &  coordsLo, const Eigen::Vector3f &  coordsHi, int  color, int  width, int  style  ) const

private

Definition at line 128 of file OpHit3DDrawer_tool.cc.

References evdb::View3D::AddPolyLine3D(), and DEFINE_ART_CLASS_TOOL.

Referenced by Draw(), and OpHit3DDrawer().

  {
    TPolyLine3D& top = view->AddPolyLine3D(5, color, width, style);
    top.SetPoint(0, coordsLo[0], coordsHi[1], coordsLo[2]);
    top.SetPoint(1, coordsHi[0], coordsHi[1], coordsLo[2]);
    top.SetPoint(2, coordsHi[0], coordsHi[1], coordsHi[2]);
    top.SetPoint(3, coordsLo[0], coordsHi[1], coordsHi[2]);
    top.SetPoint(4, coordsLo[0], coordsHi[1], coordsLo[2]);

    TPolyLine3D& side = view->AddPolyLine3D(5, color, width, style);
    side.SetPoint(0, coordsHi[0], coordsHi[1], coordsLo[2]);
    side.SetPoint(1, coordsHi[0], coordsLo[1], coordsLo[2]);
    side.SetPoint(2, coordsHi[0], coordsLo[1], coordsHi[2]);
    side.SetPoint(3, coordsHi[0], coordsHi[1], coordsHi[2]);
    side.SetPoint(4, coordsHi[0], coordsHi[1], coordsLo[2]);

    TPolyLine3D& side2 = view->AddPolyLine3D(5, color, width, style);
    side2.SetPoint(0, coordsLo[0], coordsHi[1], coordsLo[2]);
    side2.SetPoint(1, coordsLo[0], coordsLo[1], coordsLo[2]);
    side2.SetPoint(2, coordsLo[0], coordsLo[1], coordsHi[2]);
    side2.SetPoint(3, coordsLo[0], coordsHi[1], coordsHi[2]);
    side2.SetPoint(4, coordsLo[0], coordsHi[1], coordsLo[2]);

    TPolyLine3D& bottom = view->AddPolyLine3D(5, color, width, style);
    bottom.SetPoint(0, coordsLo[0], coordsLo[1], coordsLo[2]);
    bottom.SetPoint(1, coordsHi[0], coordsLo[1], coordsLo[2]);
    bottom.SetPoint(2, coordsHi[0], coordsLo[1], coordsHi[2]);
    bottom.SetPoint(3, coordsLo[0], coordsLo[1], coordsHi[2]);
    bottom.SetPoint(4, coordsLo[0], coordsLo[1], coordsLo[2]);
  }

---

The documentation for this class was generated from the following file:

• lareventdisplay/v10_00_10/source/lareventdisplay/EventDisplay/3DDrawers/OpHit3DDrawer_tool.cc