evd::TWireProjPad Class Reference

A drawing pad for time vs wire. More...

#include <TWireProjPad.h>

Inheritance diagram for evd::TWireProjPad:

Classes
struct	DrawOptions_t

Public Member Functions
	TWireProjPad (const char *nm, const char *ti, double x1, double y1, double x2, double y2, unsigned int plane)
	~TWireProjPad ()
void	Draw (const char *opt=0)
void	GetWireRange (int *i1, int *i2) const
void	SetWireRange (int i1, int i2)
void	SetZoomRange (int i1, int i2, int y1, int y2)
DrawOptions_t const &	GetDrawOptions () const
	Return the current draw options. More...
void	SetDrawOptions (DrawOptions_t const &opt)
	Receive the full set of draw options. More...
void	SetZoomToRoI (bool bZoomToRoI)
	Sets the draw option about zooming to the region of interest. More...
void	SetZoomFromView ()
	Sets the zoom parameters from the current histogram view. More...
void	SaveHitList (double i1, double i2, double y1, double y2, double distance, const char *zoom_opt, bool good_plane=true)
double	SaveSeedList (std::vector< util::PxLine > seedlines, double distance)
void	ClearHitList ()
void	SelectOneHit (double x, double y, const char *zoom_opt)
unsigned int	GetPlane () const
void	ClearandUpdatePad ()
void	UpdatePad ()
void	DrawLinesinView (std::vector< util::PxLine > lines, bool deleting=false, const char *zoom_opt=0)
void	ShowFull (int override=0)
evdb::View2D *	View () const
std::vector< double > const &	GetCurrentZoom () const
Public Member Functions inherited from evd::DrawingPad
	DrawingPad (const char *nm, const char *ti, double x1, double y1, double y2, double x2)
	~DrawingPad ()
TPad *	Pad ()
HeaderDrawer *	HeaderDraw ()
evd_tool::IExperimentDrawer *	GeometryDraw ()
SimulationDrawer *	SimulationDraw ()
RawDataDrawer *	RawDataDraw ()
RecoBaseDrawer *	RecoBaseDraw ()
AnalysisBaseDrawer *	AnalysisBaseDraw ()
HitSelector *	HitSelectorGet ()

Private Attributes
std::vector< double >	fCurrentZoom
DrawOptions_t	fDrawOpts
	set of current draw options More...
unsigned int	fPlane
	Which plane in the detector. More...
TH1F *	fHisto
	Histogram to draw object on. More...
evdb::View2D *	fView
	Collection of graphics objects to render. More...
double	fXLo
	Low value of x axis. More...
double	fXHi
	High value of x axis. More...
double	fYLo
	Low value of y axis. More...
double	fYHi
	High value of y axis. More...
int	fOri
	Orientation of the axes - see RawDrawingOptions for values. More...

Additional Inherited Members
Protected Types inherited from evd::DrawingPad
using	IExperimentDrawerPtr = std::unique_ptr< evd_tool::IExperimentDrawer >
Protected Attributes inherited from evd::DrawingPad
TPad *	fPad
	The ROOT graphics pad. More...
HeaderDrawer *	fHeaderDraw
	Drawer for event header info. More...
IExperimentDrawerPtr	fGeometryDraw
	Drawer for detector geometry. More...
SimulationDrawer *	fSimulationDraw
	Drawer for simulation objects. More...
RawDataDrawer *	fRawDataDraw
	Drawer for raw data. More...
RecoBaseDrawer *	fRecoBaseDraw
	Drawer for recobase objects. More...
AnalysisBaseDrawer *	fAnalysisBaseDraw
	Drawer for analysisbase objects. More...

Detailed Description

A drawing pad for time vs wire.

Definition at line 26 of file TWireProjPad.h.

Constructor & Destructor Documentation

evd::TWireProjPad::TWireProjPad	(	const char *	nm,
		const char *	ti,
		double	x1,
		double	x2,
		double	y1,
		double	y2,
		unsigned int	plane
	)

Create a pad showing a single X-Z or Y-Z projection of the detector

Parameters

nm	: Name of the pad
ti	: Title of the pad
x1	: Location of left edge of pad (0-1)
x2	: Location of right edge of pad (0-1)
y1	: Location of bottom edge of pad (0-1)
y2	: Location of top edge of pad (0-1)
plane	: plane number of view

Definition at line 106 of file TWireProjPad.cxx.

    : DrawingPad(nm, ti, x1, x2, y1, y2), fPlane(plane)
  {
    fCurrentZoom.resize(4);

    art::ServiceHandle<geo::Geometry const> geo;

    this->Pad()->cd();

    this->Pad()->SetLeftMargin(0.070);
    this->Pad()->SetRightMargin(0.010);

    // how many planes in the detector and
    // which plane is this one?

    unsigned int planes = geo->Nplanes();
    this->Pad()->SetTopMargin(0.005);
    this->Pad()->SetBottomMargin(0.110);

    // there has to be a better way of doing this that does
    // not have a case for each number of planes in a detector
    if (planes == 2 && fPlane > 0) {
      this->Pad()->SetTopMargin(0.110);
      this->Pad()->SetBottomMargin(0.005);
    }
    else if (planes > 2) {
      if (fPlane == 1) {
        this->Pad()->SetTopMargin(0.055);
        this->Pad()->SetBottomMargin(0.055);
      }
      else if (fPlane == 2) {
        this->Pad()->SetTopMargin(0.110);
        this->Pad()->SetBottomMargin(0.005);
      }
    }

    TString planeNo = "fTWirePlane";
    planeNo += fPlane;

    // picking the information from the current TPC
    art::ServiceHandle<evd::RawDrawingOptions const> rawopt;
    auto const signalType = geo->SignalType({rawopt->CurrentTPC(), fPlane});
    TString xtitle = ";Induction Wire;t (tdc)";
    if (signalType == geo::kCollection) xtitle = ";Collection Wire;t (tdc)";

    unsigned int const nWires = geo->Nwires(fPlane);
    unsigned int const nTicks = RawDataDraw()->TotalClockTicks();

    fXLo = -0.005 * (nWires - 1);
    fXHi = 1.005 * (nWires - 1);
    fYLo = 0.990 * (unsigned int)(this->RawDataDraw()->StartTick());
    fYHi = 1.005 * std::min((unsigned int)(this->RawDataDraw()->StartTick() + nTicks), nTicks);

    fOri = rawopt->fAxisOrientation;
    if (fOri > 0) {
      fYLo = -0.005 * (nWires - 1);
      fYHi = 1.005 * (nWires - 1);
      fYLo = 0.990 * (unsigned int)(this->RawDataDraw()->StartTick());
      fYHi = 1.005 * std::min((unsigned int)(this->RawDataDraw()->StartTick() + nTicks), nTicks);
      fXLo = -0.005 * nTicks;
      fXHi = 1.010 * nTicks;
      xtitle = ";t (tdc);InductionWire";
      if (signalType == geo::kCollection) xtitle = ";t (tdc);Collection Wire";
    }

    // make the range of the histogram be the biggest extent
    // in both directions and then use SetRangeUser() to shrink it down
    // that will allow us to change the axes on the fly
    double min = std::min(fXLo, fYLo);
    double max = std::max(fXHi, fYHi);

    fHisto = new TH1F(*(fPad->DrawFrame(min, min, max, max)));

    fHisto->SetTitleOffset(0.5, "Y");
    fHisto->SetTitleOffset(0.75, "X");
    SetZoomRange(fXLo, fXHi, fYLo, fYHi);
    fHisto->GetYaxis()->SetLabelSize(0.05);
    fHisto->GetYaxis()->CenterTitle();
    fHisto->GetXaxis()->SetLabelSize(0.05);
    fHisto->GetXaxis()->CenterTitle();
    fHisto->Draw("AB");

    fView = new evdb::View2D();
  }

evd::TWireProjPad::~TWireProjPad

(

)

Definition at line 198 of file TWireProjPad.cxx.

  {
    if (fHisto) {
      delete fHisto;
      fHisto = 0;
    }
    if (fView) {
      delete fView;
      fView = 0;
    }
  }

Member Function Documentation

void evd::TWireProjPad::ClearandUpdatePad

(

)

Definition at line 449 of file TWireProjPad.cxx.

  {
    fPad->Clear();
    this->UpdatePad();

    return;
  }

void evd::TWireProjPad::ClearHitList

(

)

Definition at line 298 of file TWireProjPad.cxx.

  {
    art::ServiceHandle<evd::RecoDrawingOptions const> recoOpt;
    if (recoOpt->fUseHitSelector) {
      this->HitSelectorGet()->ClearHitList(fPlane);
      this->Draw();
    }
  }

void evd::TWireProjPad::Draw

(

const char *

opt = 0

)

Todo:

: Why is kSelectedColor hard coded?

Definition at line 212 of file TWireProjPad.cxx.

  {
    // DumpPadsInCanvas(fPad, "TWireProjPad", "Draw()");
    MF_LOG_DEBUG("TWireProjPad") << "Started to draw plane " << fPlane;

    ///\todo: Why is kSelectedColor hard coded?
    int kSelectedColor = 4;
    fView->Clear();

    // grab the singleton holding the art::Event
    art::Event const* evtPtr = evdb::EventHolder::Instance()->GetEvent();
    if (evtPtr) {
      auto const& evt = *evtPtr;
      auto const clockData =
        art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
      auto const detProp =
        art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt, clockData);
      art::ServiceHandle<evd::RecoDrawingOptions const> recoOpt;

      this->SimulationDraw()->MCTruthVectors2D(evt, fView, fPlane);

      // the 2D pads have too much detail to be rendered on screen;
      // to act smarter, RawDataDrawer needs to know the range being plotted
      this->RawDataDraw()->ExtractRange(fPad, &GetCurrentZoom());
      this->RawDataDraw()->RawDigit2D(
        evt, detProp, fView, fPlane, GetDrawOptions().bZoom2DdrawToRoI);

      this->RecoBaseDraw()->Wire2D(evt, fView, fPlane);
      this->RecoBaseDraw()->Hit2D(evt, detProp, fView, fPlane);

      if (recoOpt->fUseHitSelector)
        this->RecoBaseDraw()->Hit2D(
          this->HitSelectorGet()->GetSelectedHits(fPlane), kSelectedColor, fView, true);

      this->RecoBaseDraw()->Slice2D(evt, detProp, fView, fPlane);
      this->RecoBaseDraw()->Cluster2D(evt, clockData, detProp, fView, fPlane);
      this->RecoBaseDraw()->EndPoint2D(evt, fView, fPlane);
      this->RecoBaseDraw()->Prong2D(evt, clockData, detProp, fView, fPlane);
      this->RecoBaseDraw()->Vertex2D(evt, detProp, fView, fPlane);
      this->RecoBaseDraw()->Seed2D(evt, detProp, fView, fPlane);
      this->RecoBaseDraw()->OpFlash2D(evt, clockData, detProp, fView, fPlane);
      this->RecoBaseDraw()->Event2D(evt, fView, fPlane);
      this->RecoBaseDraw()->DrawTrackVertexAssns2D(evt, clockData, detProp, fView, fPlane);

      UpdatePad();
    } // if (evt)

    ClearandUpdatePad();

    // check if we need to swap the axis ranges
    art::ServiceHandle<evd::RawDrawingOptions const> rawopt;
    if (fOri != rawopt->fAxisOrientation) {
      fOri = rawopt->fAxisOrientation;
      double max = fXHi;
      double min = fXLo;
      fXHi = fYHi;
      fXLo = fYLo;
      fYHi = max;
      fYLo = min;

      SetZoomRange(fXLo, fXHi, fYLo, fYHi);

      TString xtitle = fHisto->GetXaxis()->GetTitle();
      fHisto->GetXaxis()->SetTitle(fHisto->GetYaxis()->GetTitle());
      fHisto->GetYaxis()->SetTitle(xtitle);
    }

    if (fPlane > 0)
      fHisto->Draw("X+");
    else
      fHisto->Draw("");

    // Check if we should zoom the displays;
    // if there is no event, we have no clue about the region of interest
    // and therefore we don't touch anything
    if (opt == 0 && evtPtr) { this->ShowFull(); }

    MF_LOG_DEBUG("TWireProjPad") << "Started rendering plane " << fPlane;

    fView->Draw();

    MF_LOG_DEBUG("TWireProjPad") << "Drawing of plane " << fPlane << " completed";
  }

void evd::TWireProjPad::DrawLinesinView	(	std::vector< util::PxLine >	lines,
		bool	deleting = false,
		const char *	zoom_opt = 0
	)

Definition at line 472 of file TWireProjPad.cxx.

  {
    fPad->cd();
    if (deleting) {
      fPad->Clear();
      this->Draw(zoom_opt);
    }
    else {
      fView->Clear();
      fView->Draw();
    }

    mf::LogVerbatim("TWireProjPad") << "Drawing " << lines.size() << " lines";

    for (size_t is = 0; is < lines.size(); ++is) {
      if (fPlane != lines[is].plane) continue;

      TLine& l = fView->AddLine(lines[is].w0, lines[is].t0, lines[is].w1, lines[is].t1);

      fView->Draw();
      evd::Style::FromPDG(l, 11);
    }

    fView->Draw();
    UpdatePad();
    fView->Draw();

    return;
  }

std::vector<double> const& evd::TWireProjPad::GetCurrentZoom ( ) const

inline

Definition at line 72 of file TWireProjPad.h.

{return fCurrentZoom;}

DrawOptions_t const& evd::TWireProjPad::GetDrawOptions ( ) const

inline

Return the current draw options.

Definition at line 44 of file TWireProjPad.h.

{ return fDrawOpts; }

unsigned int evd::TWireProjPad::GetPlane ( ) const

inline

Definition at line 62 of file TWireProjPad.h.

{ return fPlane; }

void evd::TWireProjPad::GetWireRange	(	int *	i1,
		int *	i2
	)		const

Definition at line 337 of file TWireProjPad.cxx.

  {
    if (fOri < 1) {
      *i1 = fHisto->GetXaxis()->GetFirst();
      *i2 = fHisto->GetXaxis()->GetLast();
    }
    else {
      *i1 = fHisto->GetYaxis()->GetFirst();
      *i2 = fHisto->GetYaxis()->GetLast();
    }
  }

void evd::TWireProjPad::SaveHitList	(	double	i1,
		double	i2,
		double	y1,
		double	y2,
		double	distance,
		const char *	zoom_opt,
		bool	good_plane = true
	)

Definition at line 397 of file TWireProjPad.cxx.

  {
    const art::Event* evtPtr = evdb::EventHolder::Instance()->GetEvent();
    if (evtPtr) {
      auto const& evt = *evtPtr;
      art::ServiceHandle<evd::RecoDrawingOptions const> recoopt;
      if (recoopt->fUseHitSelector) {
        this->HitSelectorGet()->SaveHits(evt, fPlane, i1, i2, y1, y2, distance, good_plane);
        this->Draw(zoom_opt);
      }
    }
  }

double evd::TWireProjPad::SaveSeedList	(	std::vector< util::PxLine >	seedlines,
		double	distance
	)

Definition at line 420 of file TWireProjPad.cxx.

  {
    double KineticEnergy = util::kBogusD;
    const art::Event* evt = evdb::EventHolder::Instance()->GetEvent();
    if (evt) {
      art::ServiceHandle<evd::RecoDrawingOptions const> recoopt;
      if (recoopt->fUseHitSelector)
        KineticEnergy = this->HitSelectorGet()->SaveSeedLines(*evt, seedlines, distance);
    }
    return KineticEnergy;
  }

void evd::TWireProjPad::SelectOneHit	(	double	x,
		double	y,
		const char *	zoom_opt
	)

Definition at line 434 of file TWireProjPad.cxx.

  {

    const art::Event* evt = evdb::EventHolder::Instance()->GetEvent();
    if (evt) {
      art::ServiceHandle<evd::RecoDrawingOptions const> recoopt;
      if (recoopt->fUseHitSelector) {
        this->HitSelectorGet()->ChangeHit(*evt, fPlane, x, y);
        this->Draw(zoom_opt);
      }
    }
  }

void evd::TWireProjPad::SetDrawOptions ( DrawOptions_t const &  opt )

inline

Receive the full set of draw options.

Definition at line 46 of file TWireProjPad.h.

{ fDrawOpts = opt; }

void evd::TWireProjPad::SetWireRange	(	int	i1,
		int	i2
	)

Definition at line 353 of file TWireProjPad.cxx.

  {
    if (fOri < 1) { fHisto->GetXaxis()->SetRange(i1, i2); }
    else {
      fHisto->GetYaxis()->SetRange(i1, i2);
    }
    fCurrentZoom[0] = i1;
    fCurrentZoom[1] = i2;
  }

void evd::TWireProjPad::SetZoomFromView

(

)

Sets the zoom parameters from the current histogram view.

Definition at line 384 of file TWireProjPad.cxx.

  {
    TAxis const& xaxis = *(fHisto->GetXaxis());
    fCurrentZoom[0] = xaxis.GetBinLowEdge(xaxis.GetFirst());
    fCurrentZoom[1] = xaxis.GetBinUpEdge(xaxis.GetLast());
    fCurrentZoom[2] = fHisto->GetMinimum();
    fCurrentZoom[3] = fHisto->GetMaximum();
    MF_LOG_DEBUG("TWireProjPad") << "Zoom set to wires (" << fCurrentZoom[0] << "; "
                                 << fCurrentZoom[1] << " ), tick (" << fCurrentZoom[2] << "; "
                                 << fCurrentZoom[3] << ") for plane #" << fPlane;
  } // TWireProjPad::SetZoomFromView()

void evd::TWireProjPad::SetZoomRange	(	int	i1,
		int	i2,
		int	y1,
		int	y2
	)

Definition at line 367 of file TWireProjPad.cxx.

  {
    MF_LOG_DEBUG("TWireProjPad") << "SetZoomRange(" << i1 << ", " << i2 << ", " << y1 << ", " << y2
                                 << ") on plane #" << fPlane;

    fHisto->GetXaxis()->SetRangeUser(i1, i2);
    fHisto->GetYaxis()->SetRangeUser(y1, y2);
    fCurrentZoom[0] = i1;
    fCurrentZoom[1] = i2;
    fCurrentZoom[2] = y1;
    fCurrentZoom[3] = y2;
  }

void evd::TWireProjPad::SetZoomToRoI ( bool  bZoomToRoI )

inline

Sets the draw option about zooming to the region of interest.

Definition at line 49 of file TWireProjPad.h.

      { fDrawOpts.bZoom2DdrawToRoI = bZoomToRoI; }

void evd::TWireProjPad::ShowFull

(

int

override = 0

)

Definition at line 311 of file TWireProjPad.cxx.

  {
    // x values are wire numbers, y values are ticks of the clock
    int xmin = fXLo;
    int xmax = fXHi;
    int ymax = fYHi;
    int ymin = fYLo;

    art::ServiceHandle<evd::EvdLayoutOptions const> evdlayoutopt;
    art::ServiceHandle<evd::RawDrawingOptions const> rawopt;

    if (GetDrawOptions().bZoom2DdrawToRoI && !override) {
      int test = 0;
      if (rawopt->fDrawRawDataOrCalibWires == 0)
        test = RawDataDraw()->GetRegionOfInterest((int)fPlane, xmin, xmax, ymin, ymax);
      else
        test = RecoBaseDraw()->GetRegionOfInterest((int)fPlane, xmin, xmax, ymin, ymax);

      if (test != 0) return;
    }

    SetZoomRange(xmin, xmax, ymin, ymax);
  }

void evd::TWireProjPad::UpdatePad

(

)

Definition at line 459 of file TWireProjPad.cxx.

  {
    fPad->cd();
    fPad->Modified();
    fPad->Update();
    fPad->GetFrame()->SetBit(TPad::kCannotMove, true);
    fPad->SetBit(TPad::kCannotMove, true);

    return;
  }

evdb::View2D* evd::TWireProjPad::View ( ) const

inline

Definition at line 70 of file TWireProjPad.h.

{ return fView; }

Member Data Documentation

std::vector<double> evd::TWireProjPad::fCurrentZoom

private

Definition at line 80 of file TWireProjPad.h.

DrawOptions_t evd::TWireProjPad::fDrawOpts

private

set of current draw options

Definition at line 81 of file TWireProjPad.h.

TH1F* evd::TWireProjPad::fHisto

private

Histogram to draw object on.

Definition at line 85 of file TWireProjPad.h.

int evd::TWireProjPad::fOri

private

Orientation of the axes - see RawDrawingOptions for values.

Definition at line 92 of file TWireProjPad.h.

unsigned int evd::TWireProjPad::fPlane

private

Which plane in the detector.

Definition at line 84 of file TWireProjPad.h.

evdb::View2D* evd::TWireProjPad::fView

private

Collection of graphics objects to render.

Definition at line 86 of file TWireProjPad.h.

double evd::TWireProjPad::fXHi

private

High value of x axis.

Definition at line 89 of file TWireProjPad.h.

double evd::TWireProjPad::fXLo

private

Low value of x axis.

Definition at line 88 of file TWireProjPad.h.

double evd::TWireProjPad::fYHi

private

High value of y axis.

Definition at line 91 of file TWireProjPad.h.

double evd::TWireProjPad::fYLo

private

Low value of y axis.

Definition at line 90 of file TWireProjPad.h.

---

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

• TWireProjPad.h
• TWireProjPad.cxx