ArrayBM.h

Go to the documentation of this file.

/*  -*- C++ -*- */
 
/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef PILEUPCOMPS_ARRAYBM
#define PILEUPCOMPS_ARRAYBM 1
#include "GaudiKernel/ContextSpecificPtr.h"
#include "GaudiKernel/ServiceHandle.h"
#include "Gaudi/Property.h"
 
#include "PileUpTools/IBeamIntensity.h"
#include "CxxUtils/checker_macros.h"
#include "AthenaBaseComps/AthService.h"
#include "AthenaKernel/IAthRNGSvc.h"
 
#include <random>
 
namespace CLHEP
{
  class RandGeneral;
}
 
class ArrayBM : public extends<AthService, IBeamIntensity>
{
public:

  ArrayBM(const std::string& name,ISvcLocator* svc);
  virtual ~ArrayBM();

  virtual StatusCode initialize() override final;

  virtual float normFactor(int iXing) const override final;
  virtual float largestElementInPattern() const override final { return m_largestElementInPattern; }
  virtual void selectT0(unsigned int run, unsigned long long event) override final;
  virtual unsigned int getCurrentT0BunchCrossing() const override final { return m_t0Offset; }
  virtual unsigned int getBeamPatternLength() const override final { return m_ipLength; }
private:
  unsigned int m_maxBunchCrossingPerOrbit;
  Gaudi::Hive::ContextSpecificData<unsigned int> m_t0Offset;
  Gaudi::Property<std::uint64_t> m_seed{this, "Seed", 0, "Seed for FastReseededPRNG. Zero seed switches to AthRNGSvc mode."};
  Gaudi::Property<std::vector<float>> m_intensityPatternProp;
  unsigned int m_ipLength;
  double* m_intensityPattern;
  CLHEP::RandGeneral* m_biRandom;
  std::unique_ptr<std::discrete_distribution<unsigned int>> m_t0Dist{nullptr};
  ServiceHandle<IAthRNGSvc> m_randomSvc{this, "RandomSvc", "AthRNGSvc","The random number service that will be used."};
  ATHRNG::RNGWrapper*           m_rngWrapper ATLAS_THREAD_SAFE{};
  float m_largestElementInPattern;
  int m_emptyBunches;
  double* m_signalPattern;
};
#endif