Photospp_i.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef GENERATOR_PHOTOSPP_H
#define GENERATOR_PHOTOSPP_H
 
#include "AthenaBaseComps/AthAlgorithm.h"
#include "GaudiKernel/ServiceHandle.h"
 
#include "AthenaKernel/IAthRNGSvc.h"
#include "CLHEP/Random/RandomEngine.h"
 
extern "C" double phoranc_(int *idum);
 
class Photospp_i : public AthAlgorithm{
 
public:

  Photospp_i(const std::string &name, ISvcLocator *pSvcLocator);

  StatusCode initialize();

  StatusCode execute();

  void setupPhotos();

  friend double ::phoranc_(int *idum);
 
  static CLHEP::HepRandomEngine* p_rndmEngine;
 
private:


  void reseedRandomEngine(const std::string& streamName, const EventContext& ctx);
  CLHEP::HepRandomEngine* getRandomEngine(const std::string& streamName, unsigned long int randomSeedOffset, const EventContext& ctx) const;
  CLHEP::HepRandomEngine* getRandomEngineDuringInitialize(const std::string& streamName, unsigned long int randomSeedOffset, unsigned int conditionsRun=1, unsigned int lbn=1) const;
 
  // Random number service
  ServiceHandle<IAthRNGSvc> m_rndmSvc{this, "RndmSvc", "AthRNGSvc"};
 
  //Gen_tf run args.
  IntegerProperty m_dsid{this, "Dsid", 999999};

  IntegerProperty m_randomSeed{this, "RandomSeed", 1234567, "Random seed for the built-in random engine"}; // FIXME make this into an unsigned long int?

  StringProperty m_genEventKey{this, "MCEventKey", "GEN_EVENT"};

  BooleanProperty m_exponentiation{this, "ExponentiationMode", true};

  BooleanProperty m_createHistory{this, "CreateHistory", false};

  BooleanProperty m_stopCritical{this, "StopCriticalErrors", false};

  BooleanProperty m_delayInitialisation{this, "DelayInitialisation", false};

  BooleanProperty m_ZMECorrection{this, "ZMECorrection", false};

  BooleanProperty m_WMECorrection{this, "WMECorrection", false};

  BooleanProperty m_photonSplitting{this, "PhotonSplitting", false};

  DoubleProperty m_infraRedCutOff{this, "InfraRedCutOff", -1./*, 1.e-07, 0.01/91.187*/};
  DoubleProperty m_maxWtInterference{this, "WtInterference", 3.};

  DoubleProperty m_alphaQED{this, "AlphaQED", 0.00729735039};
 
};
#endif