trkf::VertexFitter Class Reference

Module for fitting a vertex using the Geometric3DVertexFitter. More...

Inheritance diagram for trkf::VertexFitter:

Classes
struct	Config
struct	Inputs

Public Types
using	Parameters = art::EDProducer::Table< Config >

Public Member Functions
	VertexFitter (Parameters const &p)
	VertexFitter (VertexFitter const &)=delete
	VertexFitter (VertexFitter &&)=delete
VertexFitter &	operator= (VertexFitter const &)=delete
VertexFitter &	operator= (VertexFitter &&)=delete

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

Private Attributes
art::InputTag	pfParticleInputTag
art::InputTag	trackInputTag
Geometric3DVertexFitter	fitter

Detailed Description

Module for fitting a vertex using the Geometric3DVertexFitter.

It selects primary PFParticles, and then collects all tracks associated to its daughters; if at least 2 tracks are found, they are passed to the vertex fitter.

Inputs are: a PFParticle collection, and the associated tracks.

Outputs are: vector of recob::Vertex, Assns of (neutrino) recob::PFParticle to recob::Vertex, Assns of recob::Vertex and recob::Track with recob::VertexAssnMeta.

For configuration options see Geometric3DVertexFitter::Parameters

Author

G. Cerati (FNAL, MicroBooNE)

Date

2017

Version

1.0

Definition at line 42 of file VertexFitter_module.cc.

Member Typedef Documentation

using trkf::VertexFitter::Parameters = art::EDProducer::Table<Config>

Definition at line 63 of file VertexFitter_module.cc.

Constructor & Destructor Documentation

trkf::VertexFitter::VertexFitter ( Parameters const &  p )

explicit

Definition at line 82 of file VertexFitter_module.cc.

  : EDProducer{p}
  , pfParticleInputTag(p().inputs().inputPFParticleLabel())
  , trackInputTag(p().inputs().inputTracksLabel())
  , fitter(p().geom3dvtxfit, p().propagator)
{
  produces<std::vector<recob::Vertex>>();
  produces<art::Assns<recob::PFParticle, recob::Vertex>>();
  produces<art::Assns<recob::Vertex, recob::Track, recob::VertexAssnMeta>>();
}

trkf::VertexFitter::VertexFitter ( VertexFitter const &  )

delete

trkf::VertexFitter::VertexFitter ( VertexFitter &&  )

delete

Member Function Documentation

VertexFitter& trkf::VertexFitter::operator= ( VertexFitter const &  )

delete

VertexFitter& trkf::VertexFitter::operator= ( VertexFitter &&  )

delete

void trkf::VertexFitter::produce ( art::Event &  e )

overrideprivate

Definition at line 94 of file VertexFitter_module.cc.

{
  using namespace std;

  auto outputVertices = make_unique<vector<recob::Vertex>>();
  auto outputPFVxAssn = make_unique<art::Assns<recob::PFParticle, recob::Vertex>>();
  auto outputVxTkMtAssn =
    make_unique<art::Assns<recob::Vertex, recob::Track, recob::VertexAssnMeta>>();

  const auto& inputPFParticle = e.getValidHandle<vector<recob::PFParticle>>(pfParticleInputTag);
  art::FindManyP<recob::Track> const assocTracks{inputPFParticle, e, trackInputTag};

  auto const detProp = art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(e);

  // PtrMakers for Assns
  art::PtrMaker<recob::Vertex> vtxPtrMaker(e);

  for (size_t iPF = 0; iPF < inputPFParticle->size(); ++iPF) {

    art::Ptr<recob::PFParticle> pfp(inputPFParticle, iPF);
    if (pfp->IsPrimary() == false || pfp->NumDaughters() < 2) continue;
    vector<art::Ptr<recob::Track>> tracks;
    auto& pfd = pfp->Daughters();
    for (auto ipfd : pfd) {
      // Daugthers returns the id as in pfp->Self() not the key
      // so need to find the key for the pfp with Self==ipfd
      for (size_t jPF = 0; jPF < inputPFParticle->size(); ++jPF) {
        art::Ptr<recob::PFParticle> pfpd(inputPFParticle, jPF);
        if (pfpd->Self() != ipfd) continue;
        vector<art::Ptr<recob::Track>> pftracks = assocTracks.at(jPF);
        for (auto t : pftracks) {
          tracks.push_back(t);
        }
        break;
      }
    }
    if (tracks.size() < 2) continue;

    VertexWrapper vtx = fitter.fitTracks(detProp, tracks);
    if (vtx.isValid() == false) continue;
    vtx.setVertexId(outputVertices->size());

    auto meta = fitter.computeMeta(detProp, vtx, tracks);

    // Fill the output collections

    outputVertices->emplace_back(vtx.vertex());
    const art::Ptr<recob::Vertex> aptr = vtxPtrMaker(outputVertices->size() - 1);
    outputPFVxAssn->addSingle(art::Ptr<recob::PFParticle>(inputPFParticle, iPF), aptr);

    size_t itt = 0;
    for (auto t : tracks) {
      outputVxTkMtAssn->addSingle(aptr, t, meta[itt]);
      itt++;
    }
  }

  e.put(std::move(outputVertices));
  e.put(std::move(outputPFVxAssn));
  e.put(std::move(outputVxTkMtAssn));
}

Member Data Documentation

Geometric3DVertexFitter trkf::VertexFitter::fitter

private

Definition at line 78 of file VertexFitter_module.cc.

art::InputTag trkf::VertexFitter::pfParticleInputTag

private

Definition at line 76 of file VertexFitter_module.cc.

art::InputTag trkf::VertexFitter::trackInputTag

private

Definition at line 77 of file VertexFitter_module.cc.

---

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

• VertexFitter_module.cc