Inheritance diagram for cosmic::CosmicClusterTagger:

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

void	produce (art::Event &e) override

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 Attributes
std::string	fClusterModuleLabel

int	fDetectorWidthTicks

int	fTickLimit

int	fMinTickDrift

int	fMaxTickDrift

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 32 of file CosmicClusterTagger_module.cc.

Constructor & Destructor Documentation

cosmic::CosmicClusterTagger::CosmicClusterTagger ( fhicl::ParameterSet const & p )

explicit

Definition at line 46 of file CosmicClusterTagger_module.cc.

                                                                          : EDProducer{p}
 {
   auto const* geo = lar::providerFrom<geo::Geometry>();
 
   auto const clock_data = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataForJob();
   auto const detp =
     art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataForJob(clock_data);
   double const samplingRate = sampling_rate(clock_data);
   fClusterModuleLabel = p.get<std::string>("ClusterModuleLabel", "cluster");
   fTickLimit = p.get<int>("TickLimit", 0);
   const double driftVelocity = detp.DriftVelocity(detp.Efield(), detp.Temperature()); // cm/us
 
   fDetectorWidthTicks =
     2 * geo->DetHalfWidth() / (driftVelocity * samplingRate / 1000); // ~3200 for uB
   fMinTickDrift = p.get("MinTickDrift", 3200);
   fMaxTickDrift = fMinTickDrift + fDetectorWidthTicks;
 
   produces<std::vector<anab::CosmicTag>>();
   produces<art::Assns<recob::Cluster, anab::CosmicTag>>();
 }

Member Function Documentation

void cosmic::CosmicClusterTagger::produce ( art::Event & e )

overridevirtual

Implements art::EDProducer.

Definition at line 68 of file CosmicClusterTagger_module.cc.

 {
   std::unique_ptr<std::vector<anab::CosmicTag>> cosmicTagClusterVector(
     new std::vector<anab::CosmicTag>);
   std::unique_ptr<art::Assns<recob::Cluster, anab::CosmicTag>> assnOutCosmicTagCluster(
     new art::Assns<recob::Cluster, anab::CosmicTag>);
 
   art::Handle<std::vector<recob::Cluster>> Cluster_h;
   e.getByLabel(fClusterModuleLabel, Cluster_h);
   std::vector<art::Ptr<recob::Cluster>> ClusterVec;
   art::fill_ptr_vector(ClusterVec, Cluster_h);
 
   /////////////////////////////////
   // LOOP OVER CLUSTERS
   /////////////////////////////////
 
   for (unsigned int iCluster = 0; iCluster < Cluster_h->size(); iCluster++) {
 
     float cosmicScore = 0;
     anab::CosmicTagID_t tag_id = anab::CosmicTagID_t::kUnknown;
 
     art::Ptr<recob::Cluster> tCluster = ClusterVec.at(iCluster);
     // art::Ptr<recob::Track> tTrack; // unused
 
     std::vector<float> endPt1;
     std::vector<float> endPt2;
 
     // Doing some checks on the cluster to determine if it's a cosmic
     bool failClusterTickCheck = false;
 
     auto const [t0, t1] = std::minmax(tCluster->StartTick(), tCluster->EndTick());
     if (t0 + fTickLimit < fMinTickDrift) { // This is into the pre-spill window
       failClusterTickCheck = true;
     }
     if (t1 - fTickLimit > fMaxTickDrift) { // This is into the post-spill window
       failClusterTickCheck = true;
     }
 
     if (failClusterTickCheck) {
       cosmicScore = 1.;
       tag_id = anab::CosmicTagID_t::kOutsideDrift_Partial;
     }
 
     if (endPt1.size() < 1) {
       for (int s = 0; s < 3; s++) {
         endPt1.push_back(-999);
         endPt2.push_back(-999);
       }
     }
 
     // Making stuff to save!
     cosmicTagClusterVector->emplace_back(endPt1, endPt2, cosmicScore, tag_id);
 
     util::CreateAssn(e, *cosmicTagClusterVector, tCluster, *assnOutCosmicTagCluster);
   }
 
   /////////////////////////////////
   // END OF CLUSTER LOOP
   /////////////////////////////////
 
   e.put(std::move(cosmicTagClusterVector));
   e.put(std::move(assnOutCosmicTagCluster));
 
 } // end of produce