Pythia6Interface.h

Go to the documentation of this file.

/*
 *  CepGen: a central exclusive processes event generator
 *  Copyright (C) 2013-2024  Laurent Forthomme
 *
 *  This program is free software: you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation, either version 3 of the License, or
 *  any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with this program.  If not, see <http://www.gnu.org/licenses/>.
 */
 
#ifndef CepGenAddOns_Pythia6Wrapper_Pythia6Interface_h
#define CepGenAddOns_Pythia6Wrapper_Pythia6Interface_h
 
#include <array>
#include <string>
#include <vector>
 
extern "C" {
extern struct {
  int n;  
  int npad;
  int k[5][4000];     
  double p[5][4000];  
  double v[5][4000];  
} pyjets_;
extern struct {
  std::array<int, 200> mstu;
  std::array<double, 200> paru;
  std::array<int, 200> mstj;
  std::array<double, 200> parj;
} pydat1_;
}
 
namespace cepgen {
  namespace pythia6 {
    int pythia6Status(int);
    int cepgenStatus(int);
    void checkPDGid(int);
 
    double pyalem(double);
    double pyalps(double);
    void pyckbd();
    void pyexec();
    void pygive(const std::string&);
    void pyjoin(std::vector<int> join);
    int pyk(int id, int qty);
    void pylist(int mlist);
    double pymass(int pdgid_);
    std::string pyname(int pdgid);
    double pyp(int id, int qty);
    inline static int& mstu(size_t i) { return pydat1_.mstu.at(i - 1); }
    inline static double& paru(size_t i) { return pydat1_.paru.at(i - 1); }
  }  // namespace pythia6
}  // namespace cepgen
 
#endif