RootTrajectorySummaryReader.hpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file RootTrajectorySummaryReader.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 "Acts/EventData/TrackParameters.hpp"
#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 RootTrajectorySummaryReader : public IReader {
 public:
  struct Config {
    std::string outputTracks = "outputTracks";  
    std::string outputParticles =
        "outputParticles";                  
    std::string treeName = "tracksummary";  
    std::string filePath;                   

    bool orderedEvents = true;
  };

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

  ~RootTrajectorySummaryReader() override;

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

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

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

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

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

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

  Config m_cfg;

  WriteDataHandle<std::vector<Acts::BoundTrackParameters>>
      m_outputTrackParameters{this, "OutputTrackParameters"};

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

  std::mutex m_read_mutex;

  size_t m_events = 0;

  TChain* m_inputChain = nullptr;

  uint32_t m_eventNr{0};  
  std::vector<uint32_t>* m_multiTrajNr =
      new std::vector<uint32_t>;  
  std::vector<unsigned int>* m_subTrajNr =
      new std::vector<unsigned int>;  


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

  std::vector<unsigned int>* m_nStates =
      new std::vector<unsigned int>;  
  std::vector<unsigned int>* m_nMeasurements =
      new std::vector<unsigned int>;  
  std::vector<unsigned int>* m_nOutliers =
      new std::vector<unsigned int>;  
  std::vector<unsigned int>* m_nHoles =
      new std::vector<unsigned int>;  
  std::vector<float>* m_chi2Sum = new std::vector<float>;  
  std::vector<unsigned int>* m_NDF =
      new std::vector<unsigned int>;  

  std::vector<std::vector<double>>* m_measurementChi2 =
      new std::vector<std::vector<double>>;  

  std::vector<std::vector<double>>* m_outlierChi2 =
      new std::vector<std::vector<double>>;  
  std::vector<std::vector<double>>* m_measurementVolume =
      new std::vector<std::vector<double>>;  

  std::vector<std::vector<double>>* m_measurementLayer =
      new std::vector<std::vector<double>>;  

  std::vector<std::vector<double>>* m_outlierVolume =
      new std::vector<std::vector<double>>;  
  std::vector<std::vector<double>>* m_outlierLayer =
      new std::vector<std::vector<double>>;  

  // The majority truth particle info
  std::vector<unsigned int>* m_nMajorityHits =
      new std::vector<unsigned int>;  

  std::vector<uint64_t>* m_majorityParticleId =
      new std::vector<uint64_t>;  
  std::vector<int>* m_t_charge =
      new std::vector<int>;  
  std::vector<float>* m_t_time =
      new std::vector<float>;  
  std::vector<float>* m_t_vx =
      new std::vector<float>;  
  std::vector<float>* m_t_vy =
      new std::vector<float>;  
  std::vector<float>* m_t_vz =
      new std::vector<float>;  
  std::vector<float>* m_t_px =
      new std::vector<float>;  
  std::vector<float>* m_t_py =
      new std::vector<float>;  
  std::vector<float>* m_t_pz =
      new std::vector<float>;  
  std::vector<float>* m_t_theta =
      new std::vector<float>;  
  std::vector<float>* m_t_phi =
      new std::vector<float>;  
  std::vector<float>* m_t_pT =
      new std::vector<float>;  
  std::vector<float>* m_t_eta =
      new std::vector<float>;  

  std::vector<bool>* m_hasFittedParams =
      new std::vector<bool>;  
  std::vector<float>* m_eLOC0_fit =
      new std::vector<float>;  
  std::vector<float>* m_eLOC1_fit =
      new std::vector<float>;  
  std::vector<float>* m_ePHI_fit =
      new std::vector<float>;  
  std::vector<float>* m_eTHETA_fit =
      new std::vector<float>;  
  std::vector<float>* m_eQOP_fit =
      new std::vector<float>;  
  std::vector<float>* m_eT_fit =
      new std::vector<float>;  
  std::vector<float>* m_err_eLOC0_fit =
      new std::vector<float>;  
  std::vector<float>* m_err_eLOC1_fit =
      new std::vector<float>;  
  std::vector<float>* m_err_ePHI_fit =
      new std::vector<float>;  
  std::vector<float>* m_err_eTHETA_fit =
      new std::vector<float>;  
  std::vector<float>* m_err_eQOP_fit =
      new std::vector<float>;  
  std::vector<float>* m_err_eT_fit =
      new std::vector<float>;  
};

}  // namespace ActsExamples