RootTrajectoryStatesWriter.hpp

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

// This file is part of the Acts project.
//
// Copyright (C) 2019-2023 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/Definitions/TrackParametrization.hpp"
#include "Acts/Utilities/Logger.hpp"
#include "ActsExamples/EventData/Index.hpp"
#include "ActsExamples/EventData/Measurement.hpp"
#include "ActsExamples/EventData/SimHit.hpp"
#include "ActsExamples/EventData/SimParticle.hpp"
#include "ActsExamples/EventData/Trajectories.hpp"
#include "ActsExamples/Framework/DataHandle.hpp"
#include "ActsExamples/Framework/ProcessCode.hpp"
#include "ActsExamples/Framework/WriterT.hpp"

#include <array>
#include <cstdint>
#include <mutex>
#include <string>
#include <vector>

class TFile;
class TTree;
namespace ActsFatras {
class Barcode;
}  // namespace ActsFatras

namespace ActsExamples {
struct AlgorithmContext;

class RootTrajectoryStatesWriter final : public WriterT<TrajectoriesContainer> {
 public:
  using HitParticlesMap = IndexMultimap<ActsFatras::Barcode>;
  using HitSimHitsMap = IndexMultimap<Index>;

  struct Config {
    std::string inputTrajectories;
    std::string inputParticles;
    std::string inputSimHits;
    std::string inputMeasurementParticlesMap;
    std::string inputMeasurementSimHitsMap;
    std::string filePath = "trackstates.root";
    std::string treeName = "trackstates";
    std::string fileMode = "RECREATE";
  };

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

  ~RootTrajectoryStatesWriter() override;

  ProcessCode finalize() override;

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

 protected:
  ProcessCode writeT(const AlgorithmContext& ctx,
                     const TrajectoriesContainer& trajectories) override;

 private:
  Config m_cfg;  

  ReadDataHandle<SimParticleContainer> m_inputParticles{this, "InputParticles"};
  ReadDataHandle<SimHitContainer> m_inputSimHits{this, "InputSimHits"};
  ReadDataHandle<HitParticlesMap> m_inputMeasurementParticlesMap{
      this, "InputMeasurementParticlesMaps"};
  ReadDataHandle<HitSimHitsMap> m_inputMeasurementSimHitsMap{
      this, "InputMeasurementSimHitsMap"};

  std::mutex m_writeMutex;  
  TFile* m_outputFile{nullptr};  
  TTree* m_outputTree{nullptr};  
  uint32_t m_eventNr{0};         
  uint32_t m_multiTrajNr{0};     
  unsigned int m_subTrajNr{0};   

  std::vector<float> m_t_x;  
  std::vector<float> m_t_y;  
  std::vector<float> m_t_z;  
  std::vector<float> m_t_r;  
  std::vector<float>
      m_t_dx;  
  std::vector<float>
      m_t_dy;  
  std::vector<float>
      m_t_dz;  

  std::vector<float> m_t_eLOC0;   
  std::vector<float> m_t_eLOC1;   
  std::vector<float> m_t_ePHI;    
  std::vector<float> m_t_eTHETA;  
  std::vector<float> m_t_eQOP;    
  std::vector<float> m_t_eT;      

  unsigned int m_nStates{0};        
  unsigned int m_nMeasurements{0};  
  std::vector<int> m_volumeID;      
  std::vector<int> m_layerID;       
  std::vector<int> m_moduleID;      
  std::vector<float> m_pathLength;  
  std::vector<float> m_lx_hit;      
  std::vector<float> m_ly_hit;      
  std::vector<float> m_x_hit;       
  std::vector<float> m_y_hit;       
  std::vector<float> m_z_hit;       
  std::vector<float> m_res_x_hit;   
  std::vector<float> m_res_y_hit;   
  std::vector<float> m_err_x_hit;   
  std::vector<float> m_err_y_hit;   
  std::vector<float> m_pull_x_hit;  
  std::vector<float> m_pull_y_hit;  
  std::vector<int> m_dim_hit;       

  std::array<int, 4>
      m_nParams{};  

  std::array<std::vector<bool>, 4>
      m_hasParams;  

  std::array<std::vector<float>, 4>
      m_eLOC0;  
  std::array<std::vector<float>, 4>
      m_eLOC1;  
  std::array<std::vector<float>, 4>
      m_ePHI;  
  std::array<std::vector<float>, 4>
      m_eTHETA;  
  std::array<std::vector<float>, 4>
      m_eQOP;  
  std::array<std::vector<float>, 4>
      m_eT;  
  std::array<std::vector<float>, 4>
      m_res_eLOC0;  

  std::array<std::vector<float>, 4>
      m_res_eLOC1;  

  std::array<std::vector<float>, 4>
      m_res_ePHI;  

  std::array<std::vector<float>, 4>
      m_res_eTHETA;  

  std::array<std::vector<float>, 4>
      m_res_eQOP;  

  std::array<std::vector<float>, 4>
      m_res_eT;  
  std::array<std::vector<float>, 4>
      m_err_eLOC0;  

  std::array<std::vector<float>, 4>
      m_err_eLOC1;  

  std::array<std::vector<float>, 4>
      m_err_ePHI;  

  std::array<std::vector<float>, 4>
      m_err_eTHETA;  

  std::array<std::vector<float>, 4>
      m_err_eQOP;  

  std::array<std::vector<float>, 4>
      m_err_eT;  
  std::array<std::vector<float>, 4>
      m_pull_eLOC0;  

  std::array<std::vector<float>, 4>
      m_pull_eLOC1;  

  std::array<std::vector<float>, 4>
      m_pull_ePHI;  

  std::array<std::vector<float>, 4>
      m_pull_eTHETA;  

  std::array<std::vector<float>, 4>
      m_pull_eQOP;  

  std::array<std::vector<float>, 4>
      m_pull_eT;  
  std::array<std::vector<float>, 4>
      m_x;  
  std::array<std::vector<float>, 4>
      m_y;  
  std::array<std::vector<float>, 4>
      m_z;  
  std::array<std::vector<float>, 4>
      m_px;  
  std::array<std::vector<float>, 4>
      m_py;  
  std::array<std::vector<float>, 4>
      m_pz;  
  std::array<std::vector<float>, 4>
      m_eta;  
  std::array<std::vector<float>, 4>
      m_pT;  

  std::vector<float> m_chi2;  
};

}  // namespace ActsExamples