Inheritance diagram for sbn::TrackStoppingChi2Fitter:

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

	TrackStoppingChi2Fitter (TrackStoppingChi2Fitter const &)=delete

	TrackStoppingChi2Fitter (TrackStoppingChi2Fitter &&)=delete

TrackStoppingChi2Fitter &	operator= (TrackStoppingChi2Fitter const &)=delete

TrackStoppingChi2Fitter &	operator= (TrackStoppingChi2Fitter &&)=delete

void	produce (art::Event &e) override

Private Attributes
const art::InputTag	fTrackLabel

const art::InputTag	fCaloLabel

const float	fMinTrackLength

sbn::TrackStoppingChi2Alg	fTrackStoppingChi2Alg

Detailed Description

Definition at line 30 of file TrackStoppingChi2Fitter_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 53 of file TrackStoppingChi2Fitter_module.cc.

     : EDProducer { p }
     , fTrackLabel(p.get<std::string>("TrackLabel"))
     , fCaloLabel(p.get<std::string>("CaloLabel"))
     , fMinTrackLength(p.get<float>("MinTrackLength"))
     , fTrackStoppingChi2Alg(p)
 {
   produces<std::vector<StoppingChi2Fit>>();
   produces<art::Assns<recob::Track, StoppingChi2Fit>>();
   produces<art::Assns<anab::Calorimetry, StoppingChi2Fit>>();
 }

sbn::TrackStoppingChi2Fitter::TrackStoppingChi2Fitter ( TrackStoppingChi2Fitter const & )

delete

sbn::TrackStoppingChi2Fitter::TrackStoppingChi2Fitter ( TrackStoppingChi2Fitter && )

delete

Member Function Documentation

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

delete

TrackStoppingChi2Fitter& sbn::TrackStoppingChi2Fitter::operator= ( TrackStoppingChi2Fitter && )

delete

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

override

Definition at line 65 of file TrackStoppingChi2Fitter_module.cc.

 {
   // Implementation of required member function here.
   auto const trackHandle(e.getValidHandle<std::vector<recob::Track>>(fTrackLabel));
 
   std::vector<art::Ptr<recob::Track>> tracks;
   art::fill_ptr_vector(tracks, trackHandle);
 
   art::FindManyP<anab::Calorimetry> fmTrackCalo(trackHandle, e, fCaloLabel);
 
   auto fitVec = std::make_unique<std::vector<StoppingChi2Fit>>();
   auto trackAssns = std::make_unique<art::Assns<recob::Track, StoppingChi2Fit>>();
   auto caloAssns = std::make_unique<art::Assns<anab::Calorimetry, StoppingChi2Fit>>();
 
   for (auto const& track : tracks) {
 
     if (track->Length() < fMinTrackLength)
       continue;
 
     const std::vector<art::Ptr<anab::Calorimetry>> caloVec(fmTrackCalo.at(track.key()));
 
     if (caloVec.size() != 3)
       continue;
 
     // Find the plane with the most hits: prefer collection > 1st induction > 2nd induction if multiple planes have the same number
     const unsigned int maxHits(std::max({ caloVec[0]->dEdx().size(), caloVec[1]->dEdx().size(), caloVec[2]->dEdx().size() }));
     const int bestPlane((caloVec[2]->dEdx().size() == maxHits) ? 2 : (caloVec[0]->dEdx().size() == maxHits) ? 0 : (caloVec[1]->dEdx().size() == maxHits) ? 1 : -1);
 
     if (bestPlane == -1)
       continue;
 
     StoppingChi2Fit thisFit(fTrackStoppingChi2Alg.RunFit(*caloVec.at(bestPlane)));
 
     if (thisFit.pol0Chi2 < 0.f || thisFit.expChi2 < 0.f)
       continue;
 
     fitVec->push_back(thisFit);
     util::CreateAssn(*this, e, *fitVec, track, *trackAssns);
     util::CreateAssn(*this, e, *fitVec, caloVec[bestPlane], *caloAssns);
   }
 
   e.put(std::move(fitVec));
   e.put(std::move(trackAssns));
   e.put(std::move(caloAssns));
 }