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
Public Member Functions inherited from evdb_tool::IWFHitDrawer
virtual	~IWFHitDrawer () noexcept=default

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 57 of file DrawSkewHits_tool.cc.

{
    configure(pset);
}

evdb_tool::DrawSkewHits::~DrawSkewHits

(

)

Definition at line 62 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 66 of file DrawSkewHits_tool.cc.

{
    return;
}

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

overridevirtual

Implements evdb_tool::IWFHitDrawer.

Definition at line 72 of file DrawSkewHits_tool.cc.

{
    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 184 of file DrawSkewHits_tool.cc.

{
    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 194 of file DrawSkewHits_tool.cc.

{
    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 52 of file DrawSkewHits_tool.cc.

---

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

• DrawSkewHits_tool.cc