TauolaPP.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
 
#ifndef TAUOLAPP_H
#define TAUOLAPP_H
 
#include "AthenaBaseComps/AthAlgorithm.h"
#include "AthenaKernel/IAthRNGSvc.h"
#include "CxxUtils/checker_macros.h"
#include "GaudiKernel/ServiceHandle.h"
 
#include "AtlasHepMC/GenEvent.h"
#include "AtlasHepMC/GenVertex.h"
 
#include "CLHEP/Random/RandomEngine.h"
 

class ATLAS_NOT_THREAD_SAFE TauolaPP : public AthAlgorithm {
public:

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

  virtual StatusCode initialize() override;

  virtual StatusCode execute(const EventContext& ctx) override;
 
  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;

  StringProperty m_key{this, "McEventKey", "GEN_EVENT"};
 
  // 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?



  IntegerProperty m_decay_particle{this, "decay_particle", 15};

  IntegerProperty m_decay_mode_same{this, "decay_mode_same", 1};

  IntegerProperty m_decay_mode_opp{this, "decay_mode_opposite", 2};

  DoubleProperty m_tau_mass{this, "tau_mass", 1.77684};

  BooleanProperty m_spin_correlation{this, "spin_correlation", true};

  BooleanProperty m_setRadiation{this, "setRadiation", true};

  DoubleProperty m_setRadiationCutOff{this, "setRadiationCutOff", 0.01};
 
};
 
#endif