FixedArrayBM.h

Go to the documentation of this file.

/*  -*- C++ -*- */
 
/*
  Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef PILEUPCOMPS_FIXEDARRAYBM
#define PILEUPCOMPS_FIXEDARRAYBM 1
#include "Gaudi/Property.h"
 
#include "PileUpTools/IBeamIntensity.h"
#include "AthenaBaseComps/AthService.h"
 
class FixedArrayBM : public extends<AthService, IBeamIntensity>
{
public:

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

  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(const EventContext& /*ctx*/) override final { }
  virtual unsigned int getCurrentT0BunchCrossing() const override final { return m_t0Offset; }
  virtual unsigned int getBeamPatternLength() const override final { return m_ipLength; }
private:
  virtual unsigned int getBCID(int iXing) const
  {
    //FIXME to be completely safe this should should probably depend on the bunch spacing too. Perhaps that concept should be deprecated though?
    return static_cast<unsigned int>((((iXing + static_cast<int>(m_t0Offset)) % static_cast<int>(m_maxBunchCrossingPerOrbit)) + static_cast<int>(m_maxBunchCrossingPerOrbit) )  % static_cast<int>(m_maxBunchCrossingPerOrbit));
  }
 
  Gaudi::Property<unsigned int> m_maxBunchCrossingPerOrbit{this, "MaxBunchCrossingPerOrbit", 3564,
    "The number of slots in each LHC beam. Default: 3564."};
 
  Gaudi::Property<unsigned int> m_t0Offset{this, "T0Offset", 0,
    "Offset of the T0 w.r.t. our intensity pattern"};
 
  Gaudi::Property<bool> m_allowEmptyT0BunchCrossing{this, "AllowEmptyT0BunchCrossing", false,
    "Allow the offset of the T0 to sit in an empty bunch crossing."};
 
  Gaudi::Property<std::vector<float>> m_intensityPatternProp{this, "IntensityPattern", {},
    "An array of floats containing the beam intensity distribution as a function of "
    "time in bins of 25ns. FixedArrayBM normalizes the distribution and uses it as a "
    "stencil to determine the relative intensity at each beam xing in the simulated range"};

  unsigned int m_ipLength{1};
  double* m_intensityPattern;
  float m_largestElementInPattern{1.0};
};
#endif