SimDriftedElectronClusterAna_module.cc

Go to the documentation of this file.

// Class:        SimDriftedElectronClusterAna
// Plugin Type: analyzer (art v2_05_00)
// File:         DriftedElectronsClustersSAna_module.cc
//

#include "lardataobj/Simulation/SimDriftedElectronCluster.h"

#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 "canvas/Utilities/InputTag.h"
#include "fhiclcpp/ParameterSet.h"

#include "TNtuple.h"

#include <iostream>
#include <vector>

namespace detsim {

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

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

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

  private:
    art::InputTag fEDepTag;

    TNtuple* fNtuple;
  };

  SimDriftedElectronClusterAna::SimDriftedElectronClusterAna(fhicl::ParameterSet const& p)
    : EDAnalyzer(p), fEDepTag{p.get<art::InputTag>("EDepModuleLabel")}
  {}

  void SimDriftedElectronClusterAna::analyze(art::Event const& e)
  {
    auto const& edep_handle =
      e.getValidHandle<std::vector<sim::SimDriftedElectronCluster>>(fEDepTag);
    auto const& edep_vec(*edep_handle);
    std::cout << "=====================edep" << edep_vec.size() << std::endl;
    for (auto const& edep : edep_vec) {
      fNtuple->Fill(e.run(), e.event(), edep.NumberOfElectrons(), edep.Time());
    }
  }

  void SimDriftedElectronClusterAna::beginJob()
  {
    art::ServiceHandle<art::TFileService const> tfs;
    fNtuple = tfs->make<TNtuple>("nt_is", "EDep IS Calc Ntuple", "run:event:ne:t");
  }
} // namespace detsim

DEFINE_ART_MODULE(detsim::SimDriftedElectronClusterAna)