106 std::unique_ptr< std::vector<recob::Track> > outputTracks(
new std::vector<recob::Track> );
137 if (pfParticlesToSpacePoints.end() == particleToSpacePointIter)
139 mf::LogDebug(
"LArPandoraTrackCreation") <<
"No spacepoints associated to particle ";
146 if (pfParticlesToClusters.end() == particleToClustersIter)
148 mf::LogDebug(
"LArPandoraShowerCreation") <<
"No clusters associated to particle ";
155 if ((pfParticlesToVertices.end() == particleToVertexIter) || (1 != particleToVertexIter->second.size()))
157 mf::LogDebug(
"LArPandoraTrackCreation") <<
"Unexpected number of vertices for particle ";
162 pandora::CartesianPointVector cartesianPointVector;
164 cartesianPointVector.emplace_back(pandora::CartesianVector(spacePoint->
XYZ()[0], spacePoint->
XYZ()[1], spacePoint->
XYZ()[2]));
166 double vertexXYZ[3] = {0., 0., 0.};
167 particleToVertexIter->second.front()->XYZ(vertexXYZ);
168 const pandora::CartesianVector vertexPosition(vertexXYZ[0], vertexXYZ[1], vertexXYZ[2]);
177 catch (
const pandora::StatusCodeException &)
179 mf::LogDebug(
"LArPandoraTrackCreation") <<
"Unable to extract sliding fit trajectory";
185 mf::LogDebug(
"LArPandoraTrackCreation") <<
"Insufficient input trajectory points to build track: " << trackStateVector.size();
193 if (trackStateVector.size()>hitsInParticle.size())
195 throw cet::exception(
"LArPandoraTrackCreation") <<
"trackStateVector.size() is greater than hitsInParticle.size()";
197 const unsigned int nInvalidPoints = hitsInParticle.size()-trackStateVector.size();
198 for (
unsigned int i=0;i<nInvalidPoints;++i) {
208 util::CreateAssn(*
this, evt, pTrack, pPFParticle, *(outputParticlesToTracks.get()));
209 util::CreateAssn(*
this, evt, *(outputTracks.get()), hitsInParticle, *(outputTracksToHits.get()));
211 for (
unsigned int hitIndex = 0; hitIndex < hitsInParticle.size(); hitIndex++)
215 outputTracksToHitsWithMeta->addSingle(pTrack, pHit, metadata);
219 mf::LogDebug(
"LArPandoraTrackCreation") <<
"Number of new tracks: " << outputTracks->size() << std::endl;
221 evt.
put(std::move(outputTracks));
222 evt.
put(std::move(outputTracksToHits));
223 evt.
put(std::move(outputTracksToHitsWithMeta));
224 evt.
put(std::move(outputParticlesToTracks));
std::map< art::Ptr< recob::PFParticle >, ClusterVector > PFParticlesToClusters
std::vector< art::Ptr< recob::PFParticle > > PFParticleVector
static void GetSlidingFitTrajectory(const pandora::CartesianPointVector &pointVector, const pandora::CartesianVector &vertexPosition, const unsigned int layerWindow, const float layerPitch, LArTrackStateVector &trackStateVector, pandora::IntVector *const pIndexVector=nullptr)
Apply 3D sliding fit to a set of 3D points and return track trajectory.
std::vector< int > IntVector
geo::Length_t WirePitch(geo::PlaneID const &planeid) const
Returns the distance between two consecutive wires.
unsigned int m_slidingFitHalfWindow
The sliding fit half window.
ProductID put(std::unique_ptr< PROD > &&product)
unsigned int MaxPlanes() const
Returns the largest number of planes among all TPCs in this detector.
static void GetAssociatedHits(const art::Event &evt, const std::string &label, const std::vector< art::Ptr< T > > &inputVector, HitVector &associatedHits, const pandora::IntVector *const indexVector=nullptr)
Get all hits associated with input clusters.
bool m_useAllParticles
Build a recob::Track for every recob::PFParticle.
std::string m_pfParticleLabel
The pf particle label.
static void CollectVertices(const art::Event &evt, const std::string &label, VertexVector &vertexVector, PFParticlesToVertices &particlesToVertices)
Collect the reconstructed PFParticles and associated Vertices from the ART event record.
std::map< art::Ptr< recob::PFParticle >, VertexVector > PFParticlesToVertices
bool CreateAssn(PRODUCER const &prod, art::Event &evt, std::vector< T > const &a, art::Ptr< U > const &b, art::Assns< U, T > &assn, std::string a_instance, size_t indx=UINT_MAX)
Creates a single one-to-one association.
static void CollectPFParticles(const art::Event &evt, const std::string &label, PFParticleVector &particleVector)
Collect the reconstructed PFParticles from the ART event record.
std::vector< art::Ptr< recob::Hit > > HitVector
const double * XYZ() const
constexpr float kBogusF
obviously bogus float value
MaybeLogger_< ELseverityLevel::ELsev_success, false > LogDebug
static bool IsTrack(const art::Ptr< recob::PFParticle > particle)
Determine whether a particle has been reconstructed as track-like.
recob::Track BuildTrack(const int id, const lar_content::LArTrackStateVector &trackStateVector) const
Build a recob::Track object.
std::map< art::Ptr< recob::PFParticle >, SpacePointVector > PFParticlesToSpacePoints
std::vector< LArTrackState > LArTrackStateVector
unsigned int m_minTrajectoryPoints
The minimum number of trajectory points.
Planes which measure W (third view for Bo, MicroBooNE, etc).
cet::coded_exception< error, detail::translate > exception
std::vector< art::Ptr< recob::Vertex > > VertexVector