RootTrackParameterWriter.hpp

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

// This file is part of the Acts project.
//
// Copyright (C) 2017-2018 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/Utilities/Logger.hpp"
#include "ActsExamples/EventData/Index.hpp"
#include "ActsExamples/EventData/ProtoTrack.hpp"
#include "ActsExamples/EventData/SimHit.hpp"
#include "ActsExamples/EventData/SimParticle.hpp"
#include "ActsExamples/EventData/Track.hpp"
#include "ActsExamples/Framework/DataHandle.hpp"
#include "ActsExamples/Framework/ProcessCode.hpp"
#include "ActsExamples/Framework/WriterT.hpp"

#include <mutex>
#include <string>

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

namespace ActsExamples {
struct AlgorithmContext;

using TrackParameterWriter = WriterT<TrackParametersContainer>;

class RootTrackParameterWriter final : public TrackParameterWriter {
 public:
  using HitParticlesMap = IndexMultimap<ActsFatras::Barcode>;
  using HitSimHitsMap = IndexMultimap<Index>;

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

  RootTrackParameterWriter(const Config& config,
                           Acts::Logging::Level level = Acts::Logging::INFO);

  ~RootTrackParameterWriter() override;

  ProcessCode finalize() override;

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

 protected:
  ProcessCode writeT(const AlgorithmContext& ctx,
                     const TrackParametersContainer& trackParams) override;

 private:
  Config m_cfg;  

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

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

  float m_loc0{NaNfloat};   
  float m_loc1{NaNfloat};   
  float m_phi{NaNfloat};    
  float m_theta{NaNfloat};  
  float m_qop{NaNfloat};    
  float m_time{NaNfloat};   
  float m_p{NaNfloat};      
  float m_pt{NaNfloat};     
  float m_eta{NaNfloat};    

  int m_t_charge{0};            
  float m_t_loc0{NaNfloat};     
  float m_t_loc1{NaNfloat};     
  float m_t_phi{NaNfloat};      
  float m_t_theta{NaNfloat};    
  float m_t_qop{NaNfloat};      
  float m_t_time{NaNfloat};     
  bool m_truthMatched = false;  
};

}  // namespace ActsExamples