RootPlanarClusterWriter.hpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file RootPlanarClusterWriter.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/Digitization/PlanarModuleCluster.hpp"
#include "Acts/Geometry/TrackingGeometry.hpp"
#include "Acts/Utilities/Logger.hpp"
#include "ActsExamples/EventData/GeometryContainers.hpp"
#include "ActsExamples/EventData/SimHit.hpp"
#include "ActsExamples/Framework/DataHandle.hpp"
#include "ActsExamples/Framework/ProcessCode.hpp"
#include "ActsExamples/Framework/WriterT.hpp"

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

class TFile;
class TTree;
namespace Acts {
class PlanarModuleCluster;
class TrackingGeometry;
}  // namespace Acts

namespace ActsExamples {
struct AlgorithmContext;

class RootPlanarClusterWriter
    : public WriterT<GeometryIdMultimap<Acts::PlanarModuleCluster>> {
 public:
  struct Config {
    std::string inputClusters;
    std::string inputSimHits;
    std::string filePath = "";          
    std::string fileMode = "RECREATE";  
    std::string treeName = "clusters";  

    std::shared_ptr<const Acts::TrackingGeometry> trackingGeometry;
  };

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

  ~RootPlanarClusterWriter() override;

  ProcessCode finalize() override;

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

 protected:
  ProcessCode writeT(
      const AlgorithmContext& ctx,
      const GeometryIdMultimap<Acts::PlanarModuleCluster>& clusters) override;

 private:
  Config m_cfg;                    
  std::mutex m_writeMutex;         
  TFile* m_outputFile{nullptr};    
  TTree* m_outputTree{nullptr};    
  int m_eventNr = 0;               
  int m_volumeID = 0;              
  int m_layerID = 0;               
  int m_surfaceID = 0;             
  float m_x = 0;                   
  float m_y = 0;                   
  float m_z = 0;                   
  float m_t = 0;                   
  float m_lx = 0;                  
  float m_ly = 0;                  
  float m_cov_lx = 0;              
  float m_cov_ly = 0;              
  std::vector<int> m_cell_IDx;     
  std::vector<int> m_cell_IDy;     
  std::vector<float> m_cell_lx;    
  std::vector<float> m_cell_ly;    
  std::vector<float> m_cell_data;  

  // (optional) the truth position
  std::vector<float> m_t_gx;          
  std::vector<float> m_t_gy;          
  std::vector<float> m_t_gz;          
  std::vector<float> m_t_gt;          
  std::vector<float> m_t_lx;          
  std::vector<float> m_t_ly;          
  std::vector<uint64_t> m_t_barcode;  

  ReadDataHandle<SimHitContainer> m_inputSimHits{this, "InputSimHits"};
};

}  // namespace ActsExamples