108 std::unique_ptr<std::vector<recob::Track>> outputTracks(
new std::vector<recob::Track>);
109 std::unique_ptr<art::Assns<recob::PFParticle, recob::Track>> outputParticlesToTracks(
111 std::unique_ptr<art::Assns<recob::Track, recob::Hit>> outputTracksToHits(
113 std::unique_ptr<art::Assns<recob::Track, recob::Hit, recob::TrackHitMeta>>
118 const float wirePitchW(detType->WirePitchW());
142 pfParticlesToSpacePoints.find(pPFParticle));
144 if (pfParticlesToSpacePoints.end() == particleToSpacePointIter) {
145 mf::LogDebug(
"LArPandoraTrackCreation") <<
"No spacepoints associated to particle ";
151 pfParticlesToClusters.find(pPFParticle));
153 if (pfParticlesToClusters.end() == particleToClustersIter) {
154 mf::LogDebug(
"LArPandoraShowerCreation") <<
"No clusters associated to particle ";
160 pfParticlesToVertices.find(pPFParticle));
162 if ((pfParticlesToVertices.end() == particleToVertexIter) ||
163 (1 != particleToVertexIter->second.size())) {
164 mf::LogDebug(
"LArPandoraTrackCreation") <<
"Unexpected number of vertices for particle ";
169 pandora::CartesianPointVector cartesianPointVector;
171 cartesianPointVector.emplace_back(pandora::CartesianVector(
172 spacePoint->
XYZ()[0], spacePoint->
XYZ()[1], spacePoint->
XYZ()[2]));
174 double vertexXYZ[3] = {0., 0., 0.};
175 particleToVertexIter->second.front()->XYZ(vertexXYZ);
176 const pandora::CartesianVector vertexPosition(vertexXYZ[0], vertexXYZ[1], vertexXYZ[2]);
189 catch (
const pandora::StatusCodeException&) {
190 mf::LogDebug(
"LArPandoraTrackCreation") <<
"Unable to extract sliding fit trajectory";
196 <<
"Insufficient input trajectory points to build track: " << trackStateVector.size();
200 HitVector hitsFromSpacePoints, hitsFromClusters, hitsInParticle;
205 particleToSpacePointIter->second,
211 for (
unsigned int hitIndex = 0; hitIndex < hitsFromSpacePoints.size(); hitIndex++) {
212 hitsInParticle.push_back(hitsFromSpacePoints.at(hitIndex));
213 (void)hitsInParticleSet.insert(hitsFromSpacePoints.at(hitIndex));
216 for (
unsigned int hitIndex = 0; hitIndex < hitsFromClusters.size(); hitIndex++) {
217 if (hitsInParticleSet.count(hitsFromClusters.at(hitIndex)) == 0)
218 hitsInParticle.push_back(hitsFromClusters.at(hitIndex));
222 if (trackStateVector.size() > hitsFromSpacePoints.size()) {
224 <<
"trackStateVector.size() is greater than hitsFromSpacePoints.size()";
226 const unsigned int nInvalidPoints = hitsInParticle.size() - trackStateVector.size();
227 for (
unsigned int i = 0; i < nInvalidPoints; ++i) {
235 outputTracks->emplace_back(
240 util::CreateAssn(evt, pTrack, pPFParticle, *(outputParticlesToTracks.get()));
241 util::CreateAssn(evt, *(outputTracks.get()), hitsInParticle, *(outputTracksToHits.get()));
244 for (
unsigned int hitIndex = 0; hitIndex < hitsInParticle.size(); hitIndex++) {
246 const int index((hitIndex < hitsFromSpacePoints.size()) ? hitIndex :
247 std::numeric_limits<int>::max());
249 outputTracksToHitsWithMeta->addSingle(pTrack, pHit, metadata);
254 <<
"Number of new tracks: " << outputTracks->size() << std::endl;
256 evt.
put(std::move(outputTracks));
257 evt.
put(std::move(outputTracksToHits));
258 evt.
put(std::move(outputTracksToHitsWithMeta));
259 evt.
put(std::move(outputParticlesToTracks));
std::unordered_set< art::Ptr< recob::Hit > > HitSet
std::map< art::Ptr< recob::PFParticle >, ClusterVector > PFParticlesToClusters
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
unsigned int m_slidingFitHalfWindow
The sliding fit half window.
std::map< art::Ptr< recob::PFParticle >, VertexVector > PFParticlesToVertices
PutHandle< PROD > put(std::unique_ptr< PROD > &&edp, std::string const &instance={})
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.
std::vector< art::Ptr< recob::PFParticle > > PFParticleVector
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.
const Double32_t * XYZ() const
std::map< art::Ptr< recob::PFParticle >, SpacePointVector > PFParticlesToSpacePoints
static void CollectPFParticles(const art::Event &evt, const std::string &label, PFParticleVector &particleVector)
Collect the reconstructed PFParticles from the ART event record.
bool CreateAssn(art::Event &evt, std::vector< T > const &a, art::Ptr< U > const &b, art::Assns< U, T > &assn, std::string a_instance, size_t index=UINT_MAX)
Creates a single one-to-one association.
std::vector< art::Ptr< recob::Hit > > HitVector
LArPandoraDetectorType * GetDetectorType()
Factory class that returns the correct detector type interface.
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::vector< art::Ptr< recob::Vertex > > VertexVector
std::vector< LArTrackState > LArTrackStateVector
unsigned int m_minTrajectoryPoints
The minimum number of trajectory points.
cet::coded_exception< error, detail::translate > exception