VtxScore_module.cc

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
// Class:       VtxScore
// Plugin Type: analyzer (art v3_02_06)
// File:        VtxScore_module.cc
//
// Generated at Tue Oct  1 11:26:39 2019 by Tingjun Yang using cetskelgen
// from cetlib version v3_07_02.
////////////////////////////////////////////////////////////////////////

#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 "art_root_io/TFileService.h"
#include "canvas/Utilities/InputTag.h"
#include "fhiclcpp/ParameterSet.h"
#include "messagefacility/MessageLogger/MessageLogger.h"
#include "lardata/ArtDataHelper/MVAReader.h"
#include "lardataobj/RecoBase/Hit.h"

#include "TTree.h"

namespace pdsp {
  class VtxScore;
}


class pdsp::VtxScore : public art::EDAnalyzer {
public:
  explicit VtxScore(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.
  VtxScore(VtxScore const&) = delete;
  VtxScore(VtxScore&&) = delete;
  VtxScore& operator=(VtxScore const&) = delete;
  VtxScore& operator=(VtxScore&&) = delete;

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

  // Selected optional functions.
  void beginJob() override;

private:

  // Declare member data here.
  // Input parameters
  art::InputTag fNNetModuleLabel;   // label of the module used for CNN tagging
  art::InputTag fHitsModuleLabel;   // label of hit finder module

  TTree *ftree;
  int run;
  int subrun;
  int event;
  std::vector<short> channel;
  std::vector<short> tpc;
  std::vector<short> plane;
  std::vector<short> wire;
  std::vector<double> charge;
  std::vector<double> peakt;
  std::vector<double> score_inel;
  std::vector<double> score_el;
  std::vector<double> score_none;
};


pdsp::VtxScore::VtxScore(fhicl::ParameterSet const& p)
  : EDAnalyzer{p},
  fNNetModuleLabel(p.get<art::InputTag>("NNetModuleLabel")),
  fHitsModuleLabel(p.get<art::InputTag>("HitsModuleLabel"))
{
}

void pdsp::VtxScore::analyze(art::Event const& e)
{

  run = e.run();
  subrun = e.subRun();
  event = e.id().event();
  channel.clear();
  tpc.clear();
  plane.clear();
  wire.clear();
  charge.clear();
  peakt.clear();
  score_inel.clear();
  score_el.clear();
  score_none.clear();

  anab::MVAReader<recob::Hit,3> hitResults(e, fNNetModuleLabel);
  
  std::vector<art::Ptr<recob::Hit> > hitlist;
  auto hitListHandle = e.getHandle< std::vector<recob::Hit> >(fHitsModuleLabel);
  if (hitListHandle) art::fill_ptr_vector(hitlist, hitListHandle);

  // loop over hits 
  //  for (size_t h = 0; h < hitResults.size(); ++h) {
  for (auto & hit : hitlist){
    
    // Get cnn output for hit h 
    std::array<float,3> cnn_out = hitResults.getOutput(hit);

    if (hit->WireID().Plane == 2){
      channel.push_back(hit->Channel());
      tpc.push_back(hit->WireID().TPC);
      plane.push_back(hit->WireID().Plane);
      wire.push_back(hit->WireID().Wire);
      charge.push_back(hit->Integral());
      peakt.push_back(hit->PeakTime());     
      score_inel.push_back(cnn_out[hitResults.getIndex("inel")]);
      score_el.push_back(cnn_out[hitResults.getIndex("el")]);
      score_none.push_back(cnn_out[hitResults.getIndex("none")]);
//      std::cout<<hit->WireID().TPC<<" "
//               <<hit->WireID().Wire<<" "
//               <<hit->PeakTime()<<" "
//               <<cnn_out[hitResults.getIndex("el")]<<" "
//               <<cnn_out[hitResults.getIndex("inel")]<<" "
//               <<cnn_out[hitResults.getIndex("none")]<<std::endl;
    }
  }
  if (!channel.empty()) ftree->Fill();
}

void pdsp::VtxScore::beginJob()
{
  art::ServiceHandle<art::TFileService> fileServiceHandle;
  ftree = fileServiceHandle->make<TTree>("ftree", "hit info");
  ftree->Branch("run", &run, "run/I");
  ftree->Branch("event", &event, "event/I");
  ftree->Branch("channel", &channel);
  ftree->Branch("tpc", &tpc);
  ftree->Branch("plane", &plane);
  ftree->Branch("wire", &wire);
  ftree->Branch("charge", &charge);
  ftree->Branch("peakt", &peakt);
  ftree->Branch("score_inel", &score_inel);
  ftree->Branch("score_el", &score_el);
  ftree->Branch("score_none", &score_none);
}

DEFINE_ART_MODULE(pdsp::VtxScore)