122 auto t0Collection = std::make_unique<std::vector<anab::T0>>();
123 auto pfpCollection = std::make_unique<std::vector<recob::PFParticle>>();
124 auto clusterCollection = std::make_unique<std::vector<recob::Cluster>>();
125 auto spCollection = std::make_unique<std::vector<recob::SpacePoint>>();
126 auto vtxCollection = std::make_unique<std::vector<recob::Vertex>>();
127 auto sliceCollection = std::make_unique<std::vector<recob::Slice>>();
128 auto pfpMetaCollection = std::make_unique<std::vector<larpandoraobj::PFParticleMetadata>>();
131 auto t0PFPAssn = std::make_unique<art::Assns<anab::T0, recob::PFParticle>>();
132 auto sliceHitAssn = std::make_unique<art::Assns<recob::Slice, recob::Hit>>();
133 auto pfpSliceAssn = std::make_unique<art::Assns<recob::PFParticle, recob::Slice>>();
134 auto pfpVtxAssn = std::make_unique<art::Assns<recob::PFParticle, recob::Vertex>>();
135 auto pfpSPAssn = std::make_unique<art::Assns<recob::PFParticle, recob::SpacePoint>>();
136 auto pfpClusterAssn = std::make_unique<art::Assns<recob::PFParticle, recob::Cluster>>();
138 std::make_unique<art::Assns<recob::PFParticle, larpandoraobj::PFParticleMetadata>>();
139 auto spHitAssn = std::make_unique<art::Assns<recob::SpacePoint, recob::Hit>>();
140 auto clusterHitAssn = std::make_unique<art::Assns<recob::Cluster, recob::Hit>>();
152 std::vector<art::Ptr<recob::Slice>> allSlices;
157 std::vector<art::Ptr<recob::Cluster>> allClusters;
162 std::vector<art::Ptr<recob::SpacePoint>> allSpacePoints;
167 std::vector<art::Ptr<recob::PFParticle>> allPFParticles;
172 std::vector<art::Ptr<recob::Track>> allTracks;
187 if (!fmSlicePFP.isValid()) {
188 throw cet::exception(
"SCECorrection") <<
"FindMany Slice-PFP is not Valid" << std::endl;
190 if (!fmPFPSP.isValid()) {
191 throw cet::exception(
"SCECorrection") <<
"FindMany PFP-SpacePoint is not Valid" << std::endl;
193 if (!fmSPHit.isValid()) {
194 throw cet::exception(
"SCECorrection") <<
"FindMany SpacePoint-Hit is not Valid" << std::endl;
203 for (
auto const& slice : allSlices) {
207 sliceCollection->push_back(newSlice);
211 const std::vector<art::Ptr<recob::PFParticle>> slicePFPs = fmSlicePFP.at(slice.key());
213 const std::map<art::Ptr<anab::T0>,
bool> sliceT0CorrectMap =
214 getSliceT0s(evt, slicePFPs, pfpHandle, trackHandle, fmPFPTrack);
216 const std::pair<art::Ptr<anab::T0>,
bool> sliceT0CorrectPair =
219 if (sliceT0CorrectPair.first.isNull() && !
fCorrectNoT0Tag) {
continue; }
223 if (!sliceT0CorrectPair.first.isNull()) {
225 t0Offset = detProp.DriftVelocity() * sliceT0CorrectPair.first->Time() / 1e3;
227 t0Collection->push_back(*sliceT0CorrectPair.first);
228 newT0Ptr = t0PtrMaker(t0Collection->size() - 1);
233 if (fmSliceHit.isValid()) {
234 const std::vector<art::Ptr<recob::Hit>> sliceHits = fmSliceHit.at(slice.key());
236 sliceHitAssn->addSingle(newSlicePtr, hitPtr);
241 for (
auto const& pfp : slicePFPs) {
245 pfpCollection->push_back(newPFP);
247 pfpSliceAssn->addSingle(newPFPPtr, newSlicePtr);
249 if (!newT0Ptr.
isNull()) { t0PFPAssn->addSingle(newT0Ptr, newPFPPtr); }
251 std::vector<art::Ptr<recob::SpacePoint>> pfpSPs = fmPFPSP.at(pfp.key());
253 if (fmPFPVertex.isValid()) {
254 std::vector<art::Ptr<recob::Vertex>> pfpVertices = fmPFPVertex.at(pfp.key());
255 for (
auto const& pfpVertex : pfpVertices) {
260 std::vector<art::Ptr<recob::SpacePoint>> vtxSPs = pfpSPs.size() ? pfpSPs : allSpacePoints;
262 double minVtxSPDist = std::numeric_limits<double>::max();
264 for (
auto const& sp : vtxSPs) {
265 geo::Point_t spPos(sp->XYZ()[0], sp->XYZ()[1], sp->XYZ()[2]);
267 if (vtxSPDiff.Mag2() < minVtxSPDist) {
269 minVtxSPDist = vtxSPDiff.Mag2();
273 if (spPtr.
isNull())
continue;
279 if (!sliceT0CorrectPair.first.isNull() && sliceT0CorrectPair.second) {
292 vtxPos, pfpVertex->covariance(), pfpVertex->chi2(), pfpVertex->ndof(), pfpVertex->ID());
293 vtxCollection->push_back(newVtx);
295 pfpVtxAssn->addSingle(newPFPPtr, newVtxPtr);
299 for (
auto const& sp : pfpSPs) {
302 geo::Point_t spPos(sp->XYZ()[0], sp->XYZ()[1], sp->XYZ()[2]);
310 if (!sliceT0CorrectPair.first.isNull() && sliceT0CorrectPair.second) {
322 Double32_t spXYZ[3] = {spPos.X(), spPos.Y(), spPos.Z()};
325 spCollection->push_back(correctedSP);
327 pfpSPAssn->addSingle(newPFPPtr, spPtr);
328 spHitAssn->addSingle(spPtr, spHitPtr);
332 if (fmPFPCluster.isValid() && fmClusterHit.isValid()) {
333 std::vector<art::Ptr<recob::Cluster>> pfpClusters = fmPFPCluster.at(pfp.key());
334 for (
auto const& pfpCluster : pfpClusters) {
336 clusterCollection->push_back(newCluster);
339 std::vector<art::Ptr<recob::Hit>> clusterHits = fmClusterHit.at(pfpCluster.key());
340 pfpClusterAssn->addSingle(newPFPPtr, newClusterPtr);
341 for (
auto const& clusterHit : clusterHits) {
342 clusterHitAssn->addSingle(newClusterPtr, clusterHit);
348 if (fmPFPMeta.isValid()) {
349 const std::vector<art::Ptr<larpandoraobj::PFParticleMetadata>> pfpMetas =
350 fmPFPMeta.at(pfp.key());
353 pfpMetaCollection->push_back(newPFPMeta);
355 pfpMetaPtrMaker(pfpMetaCollection->size() - 1);
356 pfpMetaAssn->addSingle(newPFPPtr, newPFPMetaPtr);
365 evt.
put(std::move(t0Collection));
366 evt.
put(std::move(sliceCollection));
367 evt.
put(std::move(clusterCollection));
368 evt.
put(std::move(pfpCollection));
369 evt.
put(std::move(spCollection));
370 evt.
put(std::move(vtxCollection));
371 evt.
put(std::move(pfpMetaCollection));
373 evt.
put(std::move(t0PFPAssn));
375 evt.
put(std::move(sliceHitAssn));
376 evt.
put(std::move(pfpSPAssn));
377 evt.
put(std::move(spHitAssn));
378 evt.
put(std::move(pfpVtxAssn));
379 evt.
put(std::move(pfpSliceAssn));
380 evt.
put(std::move(pfpClusterAssn));
381 evt.
put(std::move(clusterHitAssn));
382 evt.
put(std::move(pfpMetaAssn));
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double >, ROOT::Math::GlobalCoordinateSystemTag > Vector_t
Type for representation of momenta in 3D space.
virtual geo::Vector_t GetCalPosOffsets(geo::Point_t const &point, int const &TPCid) const =0
geo::Vector_t applyT0Shift(const double &t0, const geo::TPCID &tpcId) const
Set of hits with a 2D structure.
geo::WireID const & WireID() const
Initial tdc tick for hit.
Definition of vertex object for LArSoft.
spacecharge::SpaceCharge const * fSCE
PutHandle< PROD > put(std::unique_ptr< PROD > &&edp, std::string const &instance={})
std::map< art::Ptr< anab::T0 >, bool > getSliceT0s(const art::Event &evt, const std::vector< art::Ptr< recob::PFParticle >> &slicePFPs, const art::Handle< std::vector< recob::PFParticle >> &pfpHandle, const art::Handle< std::vector< recob::Track >> &trackHandle, const art::FindManyP< recob::Track > &fmPFPTrack) const
Metadata associated to PFParticles.
const std::string fPFPLabel
key_type key() const noexcept
bool isNull() const noexcept
std::pair< art::Ptr< anab::T0 >, bool > getSliceBestT0(const std::map< art::Ptr< anab::T0 >, bool > &sliceT0CorrectMap) const
The data type to uniquely identify a TPC.
constexpr TPCID const & asTPCID() const
Conversion to TPCID (for convenience of notation).
const std::string fTrackLabel
ROOT::Math::PositionVector3D< ROOT::Math::Cartesian3D< double >, ROOT::Math::GlobalCoordinateSystemTag > Point_t
Type for representation of position in physical 3D space.
virtual bool EnableCalSpatialSCE() const =0
Hierarchical representation of particle flow.
bool getByLabel(std::string const &label, std::string const &instance, Handle< PROD > &result) const
void fill_ptr_vector(std::vector< Ptr< T >> &ptrs, H const &h)
TPCID_t TPC
Index of the TPC within its cryostat.
const bool fCorrectNoT0Tag
cet::coded_exception< error, detail::translate > exception
1.8.11 
