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