StepArrayBM.h

Go to the documentation of this file.

/*  -*- C++ -*- */
 
/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef PILEUPCOMPS_STEPARRAYBM
#define PILEUPCOMPS_STEPARRAYBM 1
 
#include "GaudiKernel/ContextSpecificPtr.h"
#include "Gaudi/Property.h"
 
#include "PileUpTools/IBeamIntensity.h"
#include "AthenaBaseComps/AthService.h"
 
 
class StepArrayBM : public extends<AthService, IBeamIntensity>
{
public:
 
  StepArrayBM(const std::string& name,ISvcLocator* svc);
  virtual ~StepArrayBM();
 
  virtual StatusCode initialize() override final;
 
  virtual float normFactor(int iXing) const override final
  {
    unsigned int index = static_cast<unsigned int>((((iXing + static_cast<int>(m_t0Offset)) % static_cast<int>(m_ipLength)) + static_cast<int>(m_ipLength) )  % static_cast<int>(m_ipLength));
    //The array itself has max-value 1.0, but we want it to have mean value 1.0 (for filled bunch crossings), so we multiple by the m_largestElementInPattern.
    ATH_MSG_VERBOSE("normFactor for BCID " << iXing
                    << " (offset " << m_t0Offset
                    << " index " << index
                    << ") is = " <<  m_largestElementInPattern*m_intensityPattern[ index ]);
    return m_largestElementInPattern*m_intensityPattern[ index ];
  }
  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::Hive::ContextSpecificData<unsigned int> m_signalOffset;
  Gaudi::Property<std::vector<float>> m_intensityPatternProp;
  Gaudi::Property<std::vector<float>> m_signalPatternProp;
  unsigned int m_ipLength;
  unsigned int m_spLength;
  double* m_intensityPattern;
  bool* m_signalPattern;
  float m_largestElementInPattern;
};
#endif