Producer for TrajectoryMCSFitterICARUS. More...

Inheritance diagram for trkf::MCSFitProducerICARUS:

Classes
struct	Config

struct	Inputs

Public Types
using	Parameters = art::EDProducer::Table< Config >

Public Member Functions
	MCSFitProducerICARUS (Parameters const &p)

	~MCSFitProducerICARUS ()

	MCSFitProducerICARUS (MCSFitProducerICARUS const &)=delete

	MCSFitProducerICARUS (MCSFitProducerICARUS &&)=delete

MCSFitProducerICARUS &	operator= (MCSFitProducerICARUS const &)=delete

MCSFitProducerICARUS &	operator= (MCSFitProducerICARUS &&)=delete

void	produce (art::Event &e) override

std::vector< recob::Hit >	projectHitsOnPlane (art::Event &e, const recob::Track &traj, unsigned int p) const

Private Attributes
Parameters	p_

art::InputTag	inputTag

TrajectoryMCSFitterICARUS	mcsfitter

Detailed Description

Producer for TrajectoryMCSFitterICARUS.

Producer for TrajectoryMCSFitterICARUS, which performs a Maximum Likelihood fit of Multiple Coulomb Scattering angles between segments within a Track or Trajectory. It reads a recob::Track collection and produces a collection of recob::MCSFitResult where the elements are in the same order as the input collection (no explicit association is written).

For configuration options see MCSFitProducer::Inputs and MCSFsitProducer::Config

Author: F. Varanini (Padova, ICARUS)

Date: 2018

Version: 1.0

Definition at line 41 of file MCSFitProducerICARUS_module.cc.

Member Typedef Documentation

using trkf::MCSFitProducerICARUS::Parameters = art::EDProducer::Table<Config>

Definition at line 62 of file MCSFitProducerICARUS_module.cc.

Constructor & Destructor Documentation

trkf::MCSFitProducerICARUS::MCSFitProducerICARUS ( Parameters const & p )

explicit

Definition at line 83 of file MCSFitProducerICARUS_module.cc.

   : EDProducer{p}, p_(p), mcsfitter(p_().fitter)
 {
   inputTag = art::InputTag(p_().inputs().inputLabel());
   produces<std::vector<recob::MCSFitResult> >();
 }

trkf::MCSFitProducerICARUS::~MCSFitProducerICARUS ( )

Definition at line 90 of file MCSFitProducerICARUS_module.cc.

90 {}

trkf::MCSFitProducerICARUS::MCSFitProducerICARUS ( MCSFitProducerICARUS const & )

delete

trkf::MCSFitProducerICARUS::MCSFitProducerICARUS ( MCSFitProducerICARUS && )

delete

Member Function Documentation

MCSFitProducerICARUS& trkf::MCSFitProducerICARUS::operator= ( MCSFitProducerICARUS const & )

delete

MCSFitProducerICARUS& trkf::MCSFitProducerICARUS::operator= ( MCSFitProducerICARUS && )

delete

void trkf::MCSFitProducerICARUS::produce ( art::Event & e )

override

Definition at line 92 of file MCSFitProducerICARUS_module.cc.

 {
   //std::cout << " MCSFitProducerICARUS produce " << std::endl;
   //
   auto output  = std::make_unique<std::vector<recob::MCSFitResult> >();
   //
   art::Handle<std::vector<recob::Track> > inputH;
   bool ok = e.getByLabel(inputTag,inputH);
   if (!ok) throw cet::exception("MCSFitProducerICARUS") << "Cannot find input art::Handle with inputTag " << inputTag;
   const auto& inputVec = *(inputH.product());
 
 //std::cout << " inputh size " << inputVec.size() << std::endl;
 
 for (const auto& element : inputVec) {
     //fit
     std::vector<recob::Hit> hits2d=projectHitsOnPlane(e,element,2);
     mcsfitter.set2DHits(hits2d);
     mcsfitter.ComputeD3P();
     try{
     recob::MCSFitResult result = mcsfitter.fitMcs(element);
     output->emplace_back(std::move(result));
     } catch(...)
     {
         continue;
     }
   }
 
   e.put(std::move(output));
 }

std::vector< recob::Hit > trkf::MCSFitProducerICARUS::projectHitsOnPlane	(	art::Event &	e,
		const recob::Track &	traj,
		unsigned int	p
	)		const

Definition at line 122 of file MCSFitProducerICARUS_module.cc.

 {
 std::vector<recob::Hit> v;
   // Get track collection proxy and parallel mcs fit data (associated hits loaded by default)
   // Note: if tracks were produced from a TrackTrajectory collection you could access the original trajectories adding ',proxy::withOriginalTrajectory()' to the list of arguments
 //std::cout << " before calling proxy " << std::endl;
 auto const& tracks   = proxy::getCollection<proxy::Tracks>(e,inputTag);
 //std::cout << " after calling proxy " << std::endl;
 const auto& track = tracks[0];
 //std::cout << " proxy nhits " << track.nHits() << std::endl;
  
         for (const art::Ptr<recob::Hit>& h : track.hits())
          if(h->WireID().Plane == p) {
           //std::cout << "collection hit wire=" << h->WireID() << " peak time=" << h->PeakTime() << std::endl;
           v.emplace_back(*h);
         }
 return v;
 }