Hijing.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef GENERATORMODULESHIJING_H
#define GENERATORMODULESHIJING_H
 
#include "GeneratorModules/GenModule.h"
#include "CLHEP/Vector/LorentzVector.h"
 
#include "Hijing_i/HiParnt.h"
#include "Hijing_i/RanSeed.h"
#include "Hijing_i/HiMain1.h"
#include "Hijing_i/HiMain2.h"
#include "Hijing_i/HijJet1.h"
#include "Hijing_i/HijJet2.h"
#include "Hijing_i/HijJet4.h"
#include "Hijing_i/HiStrng.h"
#include "Hijing_i/HijCrdn.h"
 
typedef std::vector<std::string> CommandVector;
 
// new to store hijing event parameters
#include "GeneratorObjects/HijingEventParams.h"
 
class Hijing:public GenModule {
public:
  Hijing(const std::string& name, ISvcLocator* pSvcLocator);
  virtual ~Hijing() = default;
 
  virtual StatusCode genInitialize() override;
  virtual StatusCode callGenerator() override;
  virtual StatusCode genFinalize() override;
  virtual StatusCode fillEvt(HepMC::GenEvent* evt) override;
 
protected:
  virtual CLHEP::HepLorentzVector randomizeVertex(CLHEP::HepRandomEngine* engine);
 
  //Gen_tf run args.
  IntegerProperty m_dsid{this, "Dsid", 999999};
 
  // inputs to HIJSET (Hijing initialization) ...
  float m_efrm{0.};
  std::string         m_frame;
  std::string         m_proj;
  std::string         m_targ;
  int         m_iap{0};
  int         m_iat{0};
  int         m_izp{0};
  int         m_izt{0};
  // ... and HIJING (event generation)
  float m_bmin{0.};
  float m_bmax{0.};
 
  // Vertex shifting
 
  float m_x{0.};
  float m_y{0.};
  float m_z{0.};
  bool m_sel{false}; // Vertex shifting on or off
  bool m_spec{true}; // TRUE  will keep specator information (default)
 
  //Set internal randomization of vertices for beam gas gen (on or off).
  bool m_rand{false};
  bool m_wide{false}; // True allows particles with x,y distributions
  bool m_prand{false}; //BPK mirror event randomly
  bool m_keepAllDecayVertices{true};
 
  // Limit settings
  double m_partonStoreMinPt{5.0};
  double m_vertexOffsetCut{1.0E-7};
 
  // Random numbers seed
  int  m_randomseed{0};
  std::vector<long int> m_seeds;
 
  // event counter
  int m_events{0};
 
  // Commands to setup hijing
  CommandVector m_InitializeVector;
 
  // Accessor to HIJING Options and parameters COMMON
  HiParnt m_hiparnt;
 
  // Accessor to HIJING Random number seed COMMON
  RanSeed m_ranseed;
 
  // Accessors to HIJING Event information COMMONS
  HiMain1 m_himain1;
  HiMain2 m_himain2;
  HijJet1 m_hijjet1;
  HijJet2 m_hijjet2;
  HijJet4 m_hijjet4;
  HiStrng m_histrng;
  HijCrdn m_hijcrdn;
 
  void        set_user_params(void);
 
  SG::WriteHandleKey<HijingEventParams> m_event_paramsKey{this, "HijingOutputKey", "Hijing_event_params"};
};
 
#endif