Inheritance diagram for cluster::SmallClusterFinder:

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

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

Detailed Description

Definition at line 63 of file SmallClusterFinder_module.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 95 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 ( )

overrideprivate

Definition at line 107 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 )

overrideprivate

Get Clusters

Definition at line 116 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