LRAVertexPositioner.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration.
*/
 
#ifndef BEAMEFFECTS_LRAVERTEXPOSITIONER_H
#define BEAMEFFECTS_LRAVERTEXPOSITIONER_H 1
 
#include <tuple>
#include <vector>
 
#include <TFile.h>
#include <TH3.h>
 
#include "AthenaBaseComps/AthAlgTool.h"
#include "HepMC_Interfaces/ILorentzVectorGenerator.h"
 
#include "AthenaKernel/IAthRNGSvc.h"
 
namespace Simulation
{
  class LRAVertexPositioner : public extends<AthAlgTool, ILorentzVectorGenerator>
  {
  public:
    LRAVertexPositioner(const std::string &t, const std::string &n, const IInterface *p);

    virtual ~LRAVertexPositioner() = default;
 
    //---

    virtual StatusCode initialize() override final;

    virtual StatusCode finalize() override final;

    virtual CLHEP::HepLorentzVector *generate(const EventContext &ctx) const override final;
 
 
  private:
    Gaudi::Property<std::string> m_FileName{this, "FileName", "<<Unset>>", "LRA input file name."};      
    Gaudi::Property<std::string> m_HistName{this, "HistName", "<<Unset>>", "LRA input histogram name."}; 
 
    std::unique_ptr<TFile> m_LRAFile; 
    const TH3F *m_LRAHist = nullptr;  
 
    //---
 
    ServiceHandle<IAthRNGSvc> m_RNGService{this, "RNGService", "AthRNGSvc"};                                                
    Gaudi::Property<std::string> m_RNGStream{this, "RNGStream", "LRAVertexPositioner", "Stream name for the RNG service."}; 
    ATHRNG::RNGWrapper *m_RNGEngine ATLAS_THREAD_SAFE{};                                                                    
 
    //---
 
    using IntegralTuple = std::tuple<Int_t, Int_t, Int_t, Double_t>; 
    std::vector<IntegralTuple> m_Integral;                           
 
    const TAxis *m_xAxis = nullptr; 
    const TAxis *m_yAxis = nullptr; 
    const TAxis *m_zAxis = nullptr; 
  };
};
 
#endif //> !BEAMEFFECTS_LRAVERTEXPOSITIONER_H