Inheritance diagram for sbn::ShowerCosmicDistance:

Public Member Functions
	ShowerCosmicDistance (fhicl::ParameterSet const &p)

	ShowerCosmicDistance (ShowerCosmicDistance const &)=delete

	ShowerCosmicDistance (ShowerCosmicDistance &&)=delete

ShowerCosmicDistance &	operator= (ShowerCosmicDistance const &)=delete

ShowerCosmicDistance &	operator= (ShowerCosmicDistance &&)=delete

void	produce (art::Event &e) override

Private Member Functions
const std::vector< art::Ptr < recob::PFParticle > >	GetCosmicPFPs (const std::vector< art::Ptr< recob::PFParticle >> &pfps, const art::FindManyP< larpandoraobj::PFParticleMetadata > fmPFPMeta) const

const float	FindShowerResidual (const recob::Shower &shower, const std::vector< art::Ptr< recob::PFParticle >> &cosmicPFPs, const art::FindManyP< recob::SpacePoint > &fmPFPSP) const

Private Attributes
const art::InputTag	fPandoraLabel

const art::InputTag	fShowerLabel

const float	fMinShowerEnergy

Detailed Description

Definition at line 37 of file ShowerCosmicDistance_module.cc.

Constructor & Destructor Documentation

sbn::ShowerCosmicDistance::ShowerCosmicDistance ( fhicl::ParameterSet const & p )

explicit

Definition at line 67 of file ShowerCosmicDistance_module.cc.

     : EDProducer { p }
     , fPandoraLabel(p.get<art::InputTag>("PandoraLabel"))
     , fShowerLabel(p.get<art::InputTag>("ShowerLabel"))
     , fMinShowerEnergy(p.get<float>("MinShowerEnergy"))
 {
   produces<std::vector<float>>();
   produces<art::Assns<recob::Shower, float>>();
 }

sbn::ShowerCosmicDistance::ShowerCosmicDistance ( ShowerCosmicDistance const & )

delete

sbn::ShowerCosmicDistance::ShowerCosmicDistance ( ShowerCosmicDistance && )

delete

Member Function Documentation

const float sbn::ShowerCosmicDistance::FindShowerResidual	(	const recob::Shower &	shower,
		const std::vector< art::Ptr< recob::PFParticle >> &	cosmicPFPs,
		const art::FindManyP< recob::SpacePoint > &	fmPFPSP
	)		const

private

Definition at line 146 of file ShowerCosmicDistance_module.cc.

 {
   float res(std::numeric_limits<float>::max());
   const TVector3 showerStart3(shower.ShowerStart());
   const geo::Point_t showerStart(showerStart3.X(), showerStart3.Y(), showerStart3.Z());
 
   // Loop over the cosmic PFPs and find the closest
   for (auto const& cosmicPFP : cosmicPFPs) {
     auto const& pfpSPs(fmPFPSP.at(cosmicPFP.key()));
     // Loop over all of the spacepoints in the PFP and find the closest
     for (auto const& sp : pfpSPs) {
       const float dist((showerStart - sp->position()).r());
       res = std::min(res, dist);
     }
   }
   return res;
 }

const std::vector< art::Ptr< recob::PFParticle > > sbn::ShowerCosmicDistance::GetCosmicPFPs	(	const std::vector< art::Ptr< recob::PFParticle >> &	pfps,
		const art::FindManyP< larpandoraobj::PFParticleMetadata >	fmPFPMeta
	)		const

private

Definition at line 121 of file ShowerCosmicDistance_module.cc.

 {
   std::vector<art::Ptr<recob::PFParticle>> cosmicPFPs;
   for (auto const& pfp : pfps) {
 
     if (!pfp->IsPrimary())
       continue;
 
     auto const& pfpMetaVec(fmPFPMeta.at(pfp.key()));
     if (pfpMetaVec.size() != 1)
       throw cet::exception("ShowerCosmicDistance") << "Wrong metadata entries for PFP: " << pfpMetaVec.size() << ". Stopping";
 
     // Ignore anything that is not a clear cosmic
     if (!pfpMetaVec.front()->GetPropertiesMap().count("IsClearCosmic"))
       continue;
 
     // TODO: Should expand to check for T0 tags from CRT/Flash matching
 
     cosmicPFPs.push_back(pfp);
   }
   return cosmicPFPs;
 }

ShowerCosmicDistance& sbn::ShowerCosmicDistance::operator= ( ShowerCosmicDistance const & )

delete

ShowerCosmicDistance& sbn::ShowerCosmicDistance::operator= ( ShowerCosmicDistance && )

delete

void sbn::ShowerCosmicDistance::produce ( art::Event & e )

override

Definition at line 77 of file ShowerCosmicDistance_module.cc.

 {
   //Get the showers
   auto const showerHandle = e.getValidHandle<std::vector<recob::Shower>>(fShowerLabel);
   auto const pfpHandle = e.getValidHandle<std::vector<recob::PFParticle>>(fPandoraLabel);
 
   std::vector<art::Ptr<recob::Shower>> showers;
   art::fill_ptr_vector(showers, showerHandle);
 
   std::vector<art::Ptr<recob::PFParticle>> pfps;
   art::fill_ptr_vector(pfps, pfpHandle);
 
   const art::FindManyP<larpandoraobj::PFParticleMetadata> fmPFPMeta(pfpHandle, e, fPandoraLabel);
   if (!fmPFPMeta.isValid()) {
     throw cet::exception("ShowerCosmicDistance") << "PFP-Meta association is somehow not valid. Stopping";
     return;
   }
   const art::FindManyP<recob::SpacePoint> fmPFPSP(pfpHandle, e, fPandoraLabel);
   if (!fmPFPSP.isValid()) {
     throw cet::exception("ShowerCosmicDistance") << "PFP-SP association is somehow not valid. Stopping";
     return;
   }
 
   const std::vector<art::Ptr<recob::PFParticle>> cosmicPFPs(GetCosmicPFPs(pfps, fmPFPMeta));
 
   std::unique_ptr<std::vector<float>> residualCol(std::make_unique<std::vector<float>>());
   std::unique_ptr<art::Assns<recob::Shower, float>> residualAssns(std::make_unique<art::Assns<recob::Shower, float>>());
 
   for (auto const& shower : showers) {
 
     // To reduce the combinatorics, put a minumum energy cut on the showers we consider
     if (shower->best_plane() < 0 || shower->Energy().at(shower->best_plane()) < fMinShowerEnergy)
       continue;
 
     const float res(FindShowerResidual(*shower, cosmicPFPs, fmPFPSP));
 
     residualCol->push_back(res);
     util::CreateAssn(*this, e, *residualCol, shower, *residualAssns);
   }
 
   e.put(std::move(residualCol));
   e.put(std::move(residualAssns));
 }