Inheritance diagram for cluster::SmallClusterFinder:

Public Member Functions
	SmallClusterFinder (fhicl::ParameterSet const &pset)

Public Member Functions inherited from art::EDProducer
	EDProducer (fhicl::ParameterSet const &pset)

template<typename Config >
	EDProducer (Table< Config > const &config)

std::string	workerType () const

Public Member Functions inherited from art::detail::Producer
virtual	~Producer () noexcept

	Producer (fhicl::ParameterSet const &)

	Producer (Producer const &)=delete

	Producer (Producer &&)=delete

Producer &	operator= (Producer const &)=delete

Producer &	operator= (Producer &&)=delete

void	doBeginJob (SharedResources const &resources)

void	doEndJob ()

void	doRespondToOpenInputFile (FileBlock const &fb)

void	doRespondToCloseInputFile (FileBlock const &fb)

void	doRespondToOpenOutputFiles (FileBlock const &fb)

void	doRespondToCloseOutputFiles (FileBlock const &fb)

bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)

bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)

bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)

Public Member Functions inherited from art::Modifier
	~Modifier () noexcept

	Modifier ()

	Modifier (Modifier const &)=delete

	Modifier (Modifier &&)=delete

Modifier &	operator= (Modifier const &)=delete

Modifier &	operator= (Modifier &&)=delete

Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept

	ModuleBase ()

ModuleDescription const &	moduleDescription () const

void	setModuleDescription (ModuleDescription const &)

std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const

void	sortConsumables (std::string const &current_process_name)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Member Functions
void	beginJob () override

void	produce (art::Event &evt) override

int	GetPlaneAndTPC (art::Ptr< recob::Hit > a, unsigned int &p, unsigned int &cs, unsigned int &t, unsigned int &w)

Private Attributes
art::ServiceHandle< geo::Geometry const >	geom

std::string	fHitFinderModuleLabel

bool	verbose

SmallClusterFinderAlg	fSmallClusterFinderAlg

unsigned int	fNPlanes

Additional Inherited Members
Public Types inherited from art::EDProducer
using	ModuleType = EDProducer

using	WorkerType = WorkerT< EDProducer >

Public Types inherited from art::detail::Producer
template<typename UserConfig , typename KeysToIgnore = void>
using	Table = Modifier::Table< UserConfig, KeysToIgnore >

Public Types inherited from art::Modifier
template<typename UserConfig , typename UserKeysToIgnore = void>
using	Table = ProducerTable< UserConfig, detail::ModuleConfig, UserKeysToIgnore >

Static Public Member Functions inherited from art::EDProducer
static void	commitEvent (EventPrincipal &ep, Event &e)

Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 62 of file SmallClusterFinder_module.cc.

Constructor & Destructor Documentation

cluster::SmallClusterFinder::SmallClusterFinder ( fhicl::ParameterSet const & pset )

explicit

Definition at line 94 of file SmallClusterFinder_module.cc.

     : EDProducer{pset}, fSmallClusterFinderAlg(pset.get<fhicl::ParameterSet>("smallClustAlg"))
   {
     fHitFinderModuleLabel = pset.get<std::string>("HitFinderModuleLabel");
     verbose = pset.get<bool>("Verbose");
 
     produces<std::vector<recob::Cluster>>();            //This code makes clusters
     produces<art::Assns<recob::Cluster, recob::Hit>>(); //Matches clusters with hits
   }

Member Function Documentation

void cluster::SmallClusterFinder::beginJob ( )

overrideprivatevirtual

Reimplemented from art::EDProducer.

Definition at line 106 of file SmallClusterFinder_module.cc.

   {
     // this will not change on a run per run basis.
     fNPlanes = geom->Nplanes(); //get the number of planes in the TPC
   }

int cluster::SmallClusterFinder::GetPlaneAndTPC	(	art::Ptr< recob::Hit >	a,
		unsigned int &	p,
		unsigned int &	cs,
		unsigned int &	t,
		unsigned int &	w
	)

private

void cluster::SmallClusterFinder::produce ( art::Event & evt )

overrideprivatevirtual

Get Clusters

Implements art::EDProducer.

Definition at line 115 of file SmallClusterFinder_module.cc.

   {
     /**Get Clusters*/
 
     //Get the hits for this event:
     art::Handle<std::vector<recob::Hit>> HitListHandle;
     evt.getByLabel(fHitFinderModuleLabel, HitListHandle);
 
     //A vector to hold hits, not yet filled:
     std::vector<art::Ptr<recob::Hit>> hitlist;
 
     //How many hits in this event?  Tell user:
     if (verbose)
       std::cout << " ++++ Hitsreceived received " << HitListHandle->size() << " +++++ "
                 << std::endl;
     //Catch the case were there are no hits in the event:
     if (HitListHandle->size() == 0) {
       if (verbose) std::cout << "No hits received! Exiting." << std::endl;
       return;
     }
     hitlist.resize(HitListHandle->size());
 
     //wrap the hits in art::Ptrs to pass to the Alg
     for (unsigned int iHit = 0; iHit < hitlist.size(); iHit++) {
       hitlist[iHit] = art::Ptr<recob::Hit>(HitListHandle, iHit);
     }
 
     //std::cout << "Passing " << hitlist.size() << " hits to the alg." << std::endl;
 
     // Now run the alg to find the gammas:
     auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
     auto const detProp =
       art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt, clockData);
     util::GeometryUtilities const gser{*geom, clockData, detProp};
     fSmallClusterFinderAlg.FindSmallClusters(gser, clockData, detProp, hitlist);
 
     // make an art::PtrVector of the clusters
     auto SmallClusterFinder = std::make_unique<std::vector<recob::Cluster>>();
     auto assn = std::make_unique<art::Assns<recob::Cluster, recob::Hit>>();
 
     // prepare the algorithm to compute the cluster characteristics;
     // we use the "standard" one here; configuration would happen here,
     // but we are using the default configuration for that algorithm
     ClusterParamsImportWrapper<StandardClusterParamsAlg> ClusterParamAlgo;
 
     for (unsigned int iplane = 0; iplane < fNPlanes; iplane++) {
 
       auto const leftoverHits = fSmallClusterFinderAlg.GetLeftoversByPlane(iplane);
 
       //write the leftover hits as a cluster:
       if (leftoverHits.size() != 0) {
         // pick some information from the first hit
         geo::PlaneID planeID = leftoverHits.front()->WireID().planeID();
         if (verbose)
           std::cout << "Writing leftover hits to cluster ID: " << iplane * 100 << std::endl;
 
         ClusterParamAlgo.ImportHits(gser, leftoverHits);
 
         // create the recob::Cluster directly in the vector
         ClusterCreator leftover(gser,
                                 ClusterParamAlgo, // algo
                                 0.,               // start_wire
                                 0.,               // sigma_start_wire
                                 0.,               // start_tick
                                 0.,               // sigma_start_tick
                                 0.,               // end_wire
                                 0.,               // sigma_end_wire,
                                 0.,               // end_tick
                                 0.,               // sigma_end_tick
                                 iplane * 100,     // ID
                                 geom->Plane(iplane, planeID.TPC, planeID.Cryostat).View(),
                                 planeID,               // plane
                                 recob::Cluster::Sentry // sentry
         );
 
         SmallClusterFinder->emplace_back(leftover.move());
 
         util::CreateAssn(evt, *SmallClusterFinder, leftoverHits, *assn);
       } //leftovers are written for this plane, if they exist.
 
       auto const smallClusters = fSmallClusterFinderAlg.GetSmallClustersByPlane(iplane);
       for (unsigned int i = 0; i < smallClusters.size(); i++) {
         // pick some information from the first hit
         geo::PlaneID planeID; // invalid by default
         if (!smallClusters.empty()) planeID = smallClusters[i].front()->WireID().planeID();
 
         ClusterParamAlgo.ImportHits(gser, smallClusters[i]);
 
         // create the recob::Cluster directly in the vector
         ClusterCreator clust(gser,
                              ClusterParamAlgo,     // algo
                              0.,                   // start_wire
                              0.,                   // sigma_start_wire
                              0.,                   // start_tick
                              0.,                   // sigma_start_tick
                              0.,                   // end_wire
                              0.,                   // sigma_end_wire,
                              0.,                   // end_tick
                              0.,                   // sigma_end_tick
                              iplane * 100 + i + 1, // ID
                              geom->Plane(iplane, planeID.TPC, planeID.Cryostat).View(),
                              planeID,               // plane
                              recob::Cluster::Sentry // sentry
         );
 
         SmallClusterFinder->emplace_back(clust.move());
         // associate the hits to this cluster
         util::CreateAssn(evt, *SmallClusterFinder, smallClusters[i], *assn);
       }
     }
 
     evt.put(std::move(SmallClusterFinder));
     evt.put(std::move(assn));
   } //end produce