evdb_tool::DrawSkewHits Class Reference

Inheritance diagram for evdb_tool::DrawSkewHits:

Public Member Functions
	DrawSkewHits (const fhicl::ParameterSet &pset)
	~DrawSkewHits ()
void	configure (const fhicl::ParameterSet &pset) override
void	Draw (evdb::View2D &, raw::ChannelID_t &) const override

Private Member Functions
double	EvalExpoFit (double x, double tau1, double tau2, double amplitude, double peaktime) const
double	EvalMultiExpoFit (double x, int HitNumber, int NHits, std::vector< double > tau1, std::vector< double > tau2, std::vector< double > amplitude, std::vector< double > peaktime) const

Private Attributes
std::vector< TPolyLine * >	fPolyLineVec

Detailed Description

Definition at line 27 of file DrawSkewHits_tool.cc.

Constructor & Destructor Documentation

evdb_tool::DrawSkewHits::DrawSkewHits ( const fhicl::ParameterSet &  pset )

explicit

Definition at line 52 of file DrawSkewHits_tool.cc.

References configure().

  {
    configure(pset);
  }

evdb_tool::DrawSkewHits::~DrawSkewHits

(

)

Definition at line 57 of file DrawSkewHits_tool.cc.

{}

Member Function Documentation

void evdb_tool::DrawSkewHits::configure ( const fhicl::ParameterSet &  pset )

overridevirtual

Implements evdb_tool::IWFHitDrawer.

Definition at line 59 of file DrawSkewHits_tool.cc.

Referenced by DrawSkewHits().

  {
    return;
  }

void evdb_tool::DrawSkewHits::Draw ( evdb::View2D &  view2D, raw::ChannelID_t &  channel  ) const

overridevirtual

Implements evdb_tool::IWFHitDrawer.

Definition at line 64 of file DrawSkewHits_tool.cc.

References evdb::View2D::AddPolyLine(), recob::Hit::Channel(), anab::FVectorReader< T, N >::create(), art::PtrVector< T >::empty(), EvalExpoFit(), EvalMultiExpoFit(), evd::RecoDrawingOptions::fHitLabels, fPolyLineVec, evdb::EventHolder::GetEvent(), evdb::EventHolder::Instance(), art::PtrVector< T >::push_back(), art::PtrVector< T >::size(), x, and y.

  {
    art::ServiceHandle<evd::RawDrawingOptions const> rawOpt;
    art::ServiceHandle<evd::RecoDrawingOptions const> recoOpt;
    art::ServiceHandle<geo::Geometry const> geo;

    //grab the singleton with the event
    const art::Event* event = evdb::EventHolder::Instance()->GetEvent();
    if (!event) return;

    fPolyLineVec.clear();

    for (size_t imod = 0; imod < recoOpt->fHitLabels.size(); ++imod) {
      art::InputTag const which = recoOpt->fHitLabels[imod];

      art::Handle<std::vector<recob::Hit>> hitVecHandle;
      event->getByLabel(which, hitVecHandle);

      // Get a container for the subset of hits we are working with and setup to fill it
      art::PtrVector<recob::Hit> hitPtrVec;
      bool stillSearching(true);
      int fitParamsOffset(0);

      // Loop through all hits and find those on the channel in question, also count up to the start of these hits
      for (size_t hitIdx = 0; hitIdx < hitVecHandle->size(); hitIdx++) {
        art::Ptr<recob::Hit> hit(hitVecHandle, hitIdx);

        if (hit->Channel() == channel) {
          hitPtrVec.push_back(hit);
          stillSearching = false;
        }
        else if (!stillSearching)
          break;

        if (stillSearching) fitParamsOffset++;
      }

      // No hits no work
      if (hitPtrVec.empty()) continue;

      // The fit parameters will be returned in an auxiliary object
      auto hitResults = anab::FVectorReader<recob::Hit, 4>::create(*event, "dprawhit");
      const auto& fitParamVecs = hitResults->vectors();

      // Containers for the piecess...
      std::vector<double> hitPeakTimeVec;
      std::vector<double> hitTau1Vec;
      std::vector<double> hitTau2Vec;
      std::vector<double> hitPeakAmpVec;
      std::vector<int> hitStartTVec;
      std::vector<int> hitEndTVec;
      std::vector<int> hitNMultiHitVec;
      std::vector<int> hitLocalIdxVec;

      // Ok, loop through the hits for this channnel and recover the parameters
      for (size_t idx = 0; idx < hitPtrVec.size(); ++idx) {
        const auto& fitParams = fitParamVecs[fitParamsOffset + idx];
        const auto& hit = hitPtrVec[idx];

        hitPeakTimeVec.push_back(fitParams[0]);
        hitTau1Vec.push_back(fitParams[1]);
        hitTau2Vec.push_back(fitParams[2]);
        hitPeakAmpVec.push_back(fitParams[3]);
        hitStartTVec.push_back(hit->StartTick());
        hitEndTVec.push_back(hit->EndTick());
        hitNMultiHitVec.push_back(hit->Multiplicity());
        hitLocalIdxVec.push_back(hit->LocalIndex());
      }

      // Now we can go through these and start filling the polylines
      for (size_t idx = 0; idx < hitPeakTimeVec.size(); idx++) {
        if (hitNMultiHitVec[idx] > 1 && hitLocalIdxVec[idx] == 0) {
          TPolyLine& p2 = view2D.AddPolyLine(1001, kRed, 3, 1);

          // set coordinates of TPolyLine based fitted function
          for (int j = 0; j < 1001; ++j) {
            double x = hitStartTVec[idx] +
                       j * (hitEndTVec[idx + hitNMultiHitVec[idx] - 1] - hitStartTVec[idx]) / 1000;
            double y = EvalMultiExpoFit(
              x, idx, hitNMultiHitVec[idx], hitTau1Vec, hitTau2Vec, hitPeakAmpVec, hitPeakTimeVec);
            p2.SetPoint(j, x, y);
          }

          fPolyLineVec.push_back(&p2);
          p2.Draw("same");
        }

        // Always draw the single peaks in addition to the sum of all peaks
        // create TPolyLine that actually gets drawn
        TPolyLine& p1 = view2D.AddPolyLine(1001, kOrange + 7, 3, 1);

        // set coordinates of TPolyLine based fitted function
        for (int j = 0; j < 1001; ++j) {
          double x = hitStartTVec[idx - hitLocalIdxVec[idx]] +
                     j *
                       (hitEndTVec[idx + hitNMultiHitVec[idx] - hitLocalIdxVec[idx] - 1] -
                        hitStartTVec[idx - hitLocalIdxVec[idx]]) /
                       1000;
          double y = EvalExpoFit(
            x, hitTau1Vec[idx], hitTau2Vec[idx], hitPeakAmpVec[idx], hitPeakTimeVec[idx]);
          p1.SetPoint(j, x, y);
        }

        fPolyLineVec.push_back(&p1);

        p1.Draw("same");
      }
    }

    return;
  }

double evdb_tool::DrawSkewHits::EvalExpoFit ( double  x, double  tau1, double  tau2, double  amplitude, double  peaktime  ) const

private

Definition at line 176 of file DrawSkewHits_tool.cc.

Referenced by Draw().

  {
    return (amplitude * exp(0.4 * (x - peaktime) / tau1) / (1 + exp(0.4 * (x - peaktime) / tau2)));
  }

double evdb_tool::DrawSkewHits::EvalMultiExpoFit ( double  x, int  HitNumber, int  NHits, std::vector< double >  tau1, std::vector< double >  tau2, std::vector< double >  amplitude, std::vector< double >  peaktime  ) const

private

Definition at line 186 of file DrawSkewHits_tool.cc.

References DEFINE_ART_CLASS_TOOL.

Referenced by Draw().

  {
    double x_sum = 0.;

    for (int i = HitNumber; i < HitNumber + NHits; i++) {
      x_sum += (amplitude[i] * exp(0.4 * (x - peaktime[i]) / tau1[i]) /
                (1 + exp(0.4 * (x - peaktime[i]) / tau2[i])));
    }

    return x_sum;
  }

Member Data Documentation

std::vector<TPolyLine*> evdb_tool::DrawSkewHits::fPolyLineVec

mutableprivate

Definition at line 47 of file DrawSkewHits_tool.cc.

Referenced by Draw().

---

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

• lareventdisplay/v09_11_05/source/lareventdisplay/EventDisplay/wfHitDrawers/DrawSkewHits_tool.cc