SBNDOpT0FinderAna_module.cc

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// Class:       SBNDOpT0FinderAna
// Plugin Type: analyzer (art v3_05_01)
// File:        SBNDOpT0FinderAna_module.cc
//
// Generated at Mon Oct 12 13:52:37 2020 by Marco Del Tutto using cetskelgen
// from cetlib version v3_10_00.
////////////////////////////////////////////////////////////////////////

#include "art/Framework/Core/EDAnalyzer.h"
#include "art/Framework/Core/ModuleMacros.h"
#include "art/Framework/Principal/Event.h"
#include "art/Framework/Principal/Handle.h"
#include "art/Framework/Principal/Run.h"
#include "art/Framework/Principal/SubRun.h"
#include "canvas/Utilities/InputTag.h"
#include "canvas/Persistency/Common/FindManyP.h"
#include "fhiclcpp/ParameterSet.h"
#include "messagefacility/MessageLogger/MessageLogger.h"

#include "lardataobj/RecoBase/Slice.h"
#include "lardataobj/RecoBase/OpFlash.h"
#include "lardataobj/AnalysisBase/T0.h"


class SBNDOpT0FinderAna;


class SBNDOpT0FinderAna : public art::EDAnalyzer {
public:
  explicit SBNDOpT0FinderAna(fhicl::ParameterSet const& p);
  // The compiler-generated destructor is fine for non-base
  // classes without bare pointers or other resource use.

  // Plugins should not be copied or assigned.
  SBNDOpT0FinderAna(SBNDOpT0FinderAna const&) = delete;
  SBNDOpT0FinderAna(SBNDOpT0FinderAna&&) = delete;
  SBNDOpT0FinderAna& operator=(SBNDOpT0FinderAna const&) = delete;
  SBNDOpT0FinderAna& operator=(SBNDOpT0FinderAna&&) = delete;

  // Required functions.
  void analyze(art::Event const& e) override;

private:

  std::string _t0_producer; ///< The T0 producer (to be set)

};


SBNDOpT0FinderAna::SBNDOpT0FinderAna(fhicl::ParameterSet const& p)
  : EDAnalyzer{p}
{
  _t0_producer = p.get<std::string>("T0Producer");
}

void SBNDOpT0FinderAna::analyze(art::Event const& e)
{

  // Get all the T0 objects
  ::art::Handle<std::vector<anab::T0>> t0_h;
  e.getByLabel(_t0_producer, t0_h);
  if(!t0_h.isValid() || t0_h->empty()) {
    mf::LogWarning("SBNDOpT0FinderAna") << "Don't have good T0s." << std::endl;
    return;
  }
  std::vector<art::Ptr<anab::T0>> t0_v;
  art::fill_ptr_vector(t0_v, t0_h);

  // Get the T0->Slice association
  art::FindManyP<recob::Slice> t0_to_slices(t0_h, e, _t0_producer);
  art::FindManyP<recob::OpFlash> t0_to_flashes(t0_h, e, _t0_producer);

  for (size_t n_t0 = 0; n_t0 < t0_v.size(); n_t0++) {

    // The T0 object
    auto const t0 =  t0_v[n_t0];

    // The associations from T0 to both slices and flashes
    std::vector<art::Ptr<recob::Slice>> slice_v = t0_to_slices.at(t0.key());
    std::vector<art::Ptr<recob::OpFlash>> flash_v = t0_to_flashes.at(t0.key());

    // One T0 object is always associated to one and only one Slice and Flash
    assert(slice_v.size() == 1);
    assert(flash_v.size() == 1);

    std::cout << "T0 obj: " << n_t0 << ", score: " << t0->TriggerConfidence() << std::endl;
    std::cout << "\t is associated with slice ID " << slice_v[0]->ID() << std::endl;
    std::cout << "\t is associated with flash with time " << flash_v[0]->Time()
                                         << ", total PE " << flash_v[0]->TotalPE() << std::endl;
  }

}

DEFINE_ART_MODULE(SBNDOpT0FinderAna)