RootParticleReader.hpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file RootParticleReader.hpp

// This file is part of the Acts project.
//
// Copyright (C) 2017-2021 CERN for the benefit of the Acts project
//
// This Source Code Form is subject to the terms of the Mozilla Public
// License, v. 2.0. If a copy of the MPL was not distributed with this
// file, You can obtain one at http://mozilla.org/MPL/2.0/.

#pragma once

#include "ActsExamples/EventData/SimParticle.hpp"
#include "ActsExamples/Framework/DataHandle.hpp"
#include "ActsExamples/Framework/IReader.hpp"
#include "ActsExamples/Framework/ProcessCode.hpp"
#include <Acts/Definitions/Algebra.hpp>
#include <Acts/Propagator/MaterialInteractor.hpp>
#include <Acts/Utilities/Logger.hpp>

#include <cstddef>
#include <cstdint>
#include <memory>
#include <mutex>
#include <string>
#include <utility>
#include <vector>

class TChain;

namespace ActsExamples {
struct AlgorithmContext;

class RootParticleReader : public IReader {
 public:
  struct Config {
    std::string particleCollection =
        "particleCollection";             
    std::string vertexPrimaryCollection;  
    std::string
        vertexSecondaryCollection;  
    std::string treeName = "particles";  
    std::string filePath;                

    bool orderedEvents = true;
  };

  RootParticleReader(const Config& config, Acts::Logging::Level level);

  ~RootParticleReader() override;

  std::string name() const override { return "RootParticleReader"; }

  std::pair<size_t, size_t> availableEvents() const override;

  ProcessCode read(const ActsExamples::AlgorithmContext& context) override;

  const Config& config() const { return m_cfg; }

 private:
  const Acts::Logger& logger() const { return *m_logger; }

  Config m_cfg;

  WriteDataHandle<SimParticleContainer> m_outputParticles{this,
                                                          "OutputParticles"};

  WriteDataHandle<std::vector<uint32_t>> m_outputPrimaryVertices{
      this, "OutputPrimaryVertices"};
  WriteDataHandle<std::vector<uint32_t>> m_outputSecondaryVertices{
      this, "OutputSecondaryVertices"};

  std::unique_ptr<const Acts::Logger> m_logger;

  std::mutex m_read_mutex;

  size_t m_events = 0;

  TChain* m_inputChain = nullptr;

  uint32_t m_eventId = 0;

  std::vector<long long> m_entryNumbers = {};

  std::vector<uint64_t>* m_particleId = new std::vector<uint64_t>;
  std::vector<int32_t>* m_particleType = new std::vector<int32_t>;
  std::vector<uint32_t>* m_process = new std::vector<uint32_t>;
  std::vector<float>* m_vx = new std::vector<float>;
  std::vector<float>* m_vy = new std::vector<float>;
  std::vector<float>* m_vz = new std::vector<float>;
  std::vector<float>* m_vt = new std::vector<float>;
  std::vector<float>* m_px = new std::vector<float>;
  std::vector<float>* m_py = new std::vector<float>;
  std::vector<float>* m_pz = new std::vector<float>;
  std::vector<float>* m_m = new std::vector<float>;
  std::vector<float>* m_q = new std::vector<float>;
  std::vector<float>* m_eta = new std::vector<float>;
  std::vector<float>* m_phi = new std::vector<float>;
  std::vector<float>* m_pt = new std::vector<float>;
  std::vector<float>* m_p = new std::vector<float>;
  std::vector<uint32_t>* m_vertexPrimary = new std::vector<uint32_t>;
  std::vector<uint32_t>* m_vertexSecondary = new std::vector<uint32_t>;
  std::vector<uint32_t>* m_particle = new std::vector<uint32_t>;
  std::vector<uint32_t>* m_generation = new std::vector<uint32_t>;
  std::vector<uint32_t>* m_subParticle = new std::vector<uint32_t>;
};

}  // namespace ActsExamples