CosmicClusterTagger_module.cc

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// Class:       CosmicTagger
// Module Type: producer
// File:        CosmicTagger_module.cc
//
// Generated at Mon Sep 24 18:21:00 2012 by Sarah Lockwitz using artmod
// from art v1_02_02.
// artmod -e beginJob -e reconfigure -e endJob producer trkf::CosmicTagger
////////////////////////////////////////////////////////////////////////

#include "art/Framework/Core/EDProducer.h"
#include "art/Framework/Core/ModuleMacros.h"
#include "art/Framework/Principal/Event.h"

#include <algorithm> // std::minmax() ...
#include <memory>    // std::unique_ptr<>
#include <string>
#include <utility> // std::pair<>, std::move()

#include "larcore/Geometry/Geometry.h"
#include "lardata/DetectorInfoServices/DetectorClocksService.h"
#include "lardata/DetectorInfoServices/DetectorPropertiesService.h"
#include "lardata/Utilities/AssociationUtil.h"
#include "lardataobj/AnalysisBase/CosmicTag.h"
#include "lardataobj/RecoBase/Cluster.h"
#include "lardataobj/RecoBase/Track.h"

namespace cosmic {
  class CosmicClusterTagger;
}

class cosmic::CosmicClusterTagger : public art::EDProducer {
public:
  explicit CosmicClusterTagger(fhicl::ParameterSet const& p);

  void produce(art::Event& e) override;

private:
  std::string fClusterModuleLabel;
  int fDetectorWidthTicks;
  int fTickLimit;
  int fMinTickDrift;
  int fMaxTickDrift;
};

cosmic::CosmicClusterTagger::CosmicClusterTagger(fhicl::ParameterSet const& p) : 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>>();
}

void
cosmic::CosmicClusterTagger::produce(art::Event& e)
{
  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
//////////////////////////////////////////////////////////////////////////////////////////////////////

DEFINE_ART_MODULE(cosmic::CosmicClusterTagger)