PIDAAnalyzer_module.cc

Go to the documentation of this file.

// Class:       PIDAAnalyzer
// Module Type: analyzer
// File:        PIDAAnalyzer_module.cc
//
// Generated at Sat Nov  1 23:10:27 2014 by Wesley Ketchum using artmod
// from cetpkgsupport v1_07_01.

#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/Services/Registry/ServiceHandle.h"
#include "art_root_io/TFileService.h"
#include "fhiclcpp/ParameterSet.h"

#include "PIDAAlg.h"
#include "lardataobj/AnalysisBase/Calorimetry.h"

#include "TH1F.h"
#include "TTree.h"

namespace pid {
  class PIDAAnalyzer;
}

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

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

  void analyze(art::Event const& e) override;

  void beginJob() override;

private:
  std::string fCaloModuleLabel;
  PIDAAlg fPIDAAlg;
};

pid::PIDAAnalyzer::PIDAAnalyzer(fhicl::ParameterSet const& p)
  : EDAnalyzer(p)
  , fCaloModuleLabel(p.get<std::string>("CaloModuleLabel"))
  , fPIDAAlg(p.get<fhicl::ParameterSet>("PIDAAlg"))
{}

void pid::PIDAAnalyzer::beginJob()
{
  art::ServiceHandle<art::TFileService const> tfs;

  std::vector<TH1F*> kde_hists;
  for (size_t i_b = 0; i_b < fPIDAAlg.getNKDEBandwidths(); i_b++) {
    std::stringstream hname;
    hname << "hkde_" << i_b;
    kde_hists.push_back(tfs->make<TH1F>(hname.str().c_str(), "PIDA KDE Distribution", 100, 0, 30));
  }

  fPIDAAlg.SetPIDATree(tfs->make<TTree>("pida", "PIDAPropertiesTree"),
                       tfs->make<TH1F>("hvals", "PIDA Distribution", 100, 0, 30),
                       kde_hists);
}

void pid::PIDAAnalyzer::analyze(art::Event const& e)
{
  art::Handle<std::vector<anab::Calorimetry>> caloHandle;
  e.getByLabel(fCaloModuleLabel, caloHandle);
  std::vector<anab::Calorimetry> const& caloVector(*caloHandle);

  for (size_t i_calo = 0; i_calo < caloVector.size(); i_calo++) {
    fPIDAAlg.FillPIDATree(e.run(), e.event(), i_calo, caloVector[i_calo]);
  }
}

DEFINE_ART_MODULE(pid::PIDAAnalyzer)