PythiaGun.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file PythiaGun.h

/*******************************************************************************
 * Copyright (c) The JETSCAPE Collaboration, 2018
 *
 * Modular, task-based framework for simulating all aspects of heavy-ion collisions
 *
 * For the list of contributors see AUTHORS.
 *
 * Report issues at https://github.com/JETSCAPE/JETSCAPE/issues
 *
 * or via email to bugs.jetscape@gmail.com
 *
 * Distributed under the GNU General Public License 3.0 (GPLv3 or later).
 * See COPYING for details.
 ******************************************************************************/

// Create a pythia collision at a specified point and return the two inital hard partons

#ifndef PYTHIAGUN_H
#define PYTHIAGUN_H

#include "HardProcess.h"
#include "JetScapeLogger.h"
#include "Pythia8/Pythia.h"

using namespace Jetscape;

class PythiaGun : public HardProcess, public Pythia8::Pythia {

private:
  double pTHatMin;
  double pTHatMax;
  double eCM;
  double vir_factor;
  double softMomentumCutoff;
  bool FSR_on;
  bool softQCD;

  // Allows the registration of the module so that it is available to be used by the Jetscape framework.
  static RegisterJetScapeModule<PythiaGun> reg;

public:
  PythiaGun(string xmlDir = "DONTUSETHIS", bool printBanner = false)
      : Pythia8::Pythia(xmlDir, printBanner), HardProcess() {
    SetId("UninitializedPythiaGun");
  }

  ~PythiaGun();

  void InitTask();
  void Exec();

  // Getters
  double GetpTHatMin() const { return pTHatMin; }
  double GetpTHatMax() const { return pTHatMax; }

  // Cross-section information in mb and event weight.
  double GetSigmaGen() { return info.sigmaGen(); };
  double GetSigmaErr() { return info.sigmaErr(); };
  double GetPtHat() { return info.pTHat(); };
  double GetEventWeight() { return info.weight(); };
};

#endif // PYTHIAGUN_H