trkf::TCTrack Class Reference

Inheritance diagram for trkf::TCTrack:

Public Member Functions
	TCTrack (fhicl::ParameterSet const &pset)

Private Member Functions
void	produce (art::Event &evt) override

Private Attributes
SpacePointAlg	fSptalg
art::InputTag	fPFPtag

Detailed Description

Definition at line 32 of file TCTrack_module.cc.

Constructor & Destructor Documentation

trkf::TCTrack::TCTrack ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 46 of file TCTrack_module.cc.

    : EDProducer{pset}
    , fSptalg{pset.get<fhicl::ParameterSet>("SpacePointAlg")}
    , fPFPtag{pset.get<std::string>("PFPModuleLabel")}
  {
    produces<std::vector<recob::SpacePoint>>();
    produces<art::Assns<recob::SpacePoint, recob::Hit>>();
  }

Member Function Documentation

void trkf::TCTrack::produce ( art::Event &  evt )

overrideprivate

Definition at line 57 of file TCTrack_module.cc.

  {
    // all data products are assumed to be produced by the same module that produced the PFParticles -> TrajCluster_module
    auto pfpHandle = evt.getValidHandle<std::vector<recob::PFParticle>>(fPFPtag);
    auto clsHandle = evt.getValidHandle<std::vector<recob::Cluster>>(fPFPtag);

    art::FindManyP<recob::Cluster> pfp_cls(pfpHandle, evt, fPFPtag);
    art::FindManyP<recob::Hit> cls_hit(clsHandle, evt, fPFPtag);

    auto sphitassn = std::make_unique<art::Assns<recob::SpacePoint, recob::Hit>>();
    auto spts = std::make_unique<std::vector<recob::SpacePoint>>();

    auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
    auto const detProp =
      art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt, clockData);

    art::PtrVector<recob::Hit> hits;
    for (unsigned short ipfp = 0; ipfp < pfpHandle->size(); ++ipfp) {
      // Get the clusters associated with this PFParticle - there should be one in each plane
      std::vector<art::Ptr<recob::Cluster>> clsList;
      pfp_cls.get(ipfp, clsList);
      hits.clear();
      std::cout << "PFP " << ipfp << "\n";
      for (unsigned short icl = 0; icl < clsList.size(); ++icl) {
        std::vector<art::Ptr<recob::Hit>> hitList;
        unsigned int clsIndex = clsList[icl]->ID() - 1;
        cls_hit.get(clsIndex, hitList);
        std::cout << " cls index " << clsIndex << " hits size " << hitList.size() << "  "
                  << (int)clsList[icl]->StartWire() << ":" << (int)clsList[icl]->StartTick()
                  << "  EndWire " << (int)clsList[icl]->EndWire() << ":"
                  << (int)clsList[icl]->EndTick() << "\n";
        hits.reserve(hits.size() + hitList.size());
        hits.insert(hits.end(), hitList.begin(), hitList.end());
      } // icl
      // make new space points using these hits
      std::vector<recob::SpacePoint> new_spts;
      fSptalg.makeSpacePoints(clockData, detProp, hits, new_spts);
      if (new_spts.empty()) continue;

      int nspt = spts->size();
      spts->insert(spts->end(), new_spts.begin(), new_spts.end());
      // associate the hits with the spacepoints
      art::PtrVector<recob::SpacePoint> sptvec;
      for (unsigned int ispt = nspt; ispt < spts->size(); ++ispt) {
        const recob::SpacePoint& spt = (*spts)[ispt];
        const art::PtrVector<recob::Hit>& hits = fSptalg.getAssociatedHits(spt);
        util::CreateAssn(evt, *spts, hits, *sphitassn, ispt);
      } // ispt
    }   // ipfp

    evt.put(std::move(spts));
    evt.put(std::move(sphitassn));

  } // TCTrack::produce()

Member Data Documentation

art::InputTag trkf::TCTrack::fPFPtag

private

Definition at line 42 of file TCTrack_module.cc.

SpacePointAlg trkf::TCTrack::fSptalg

private

Definition at line 40 of file TCTrack_module.cc.

---

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

• TCTrack_module.cc