PythiaChiGen.h

Go to the documentation of this file.

//****************************************************************************
//*                   This file is part of PandaRoot.                        *
//*                                                                          *
//*            PandaRoot is distributed under the terms of the               *
//*              GNU General Public License (GPL) version 3,                 *
//*                 copied verbatim in the file "LICENSE".                   *
//*                                                                          *
//*  Copyright (C) 2006 - 2024 FAIR GmbH and copyright holders of PandaRoot  *
//*     The copyright holders are listed in the file "COPYRIGHTHOLDERS".     *
//*               The authors are listed in the file "AUTHORS".              *
//****************************************************************************

/*
 * @author Alexey Luchinsky
 * @author Stanislav Poslavsky (stvlpos (at) mail.ru)
 */

#ifndef PYTHIAHADRONIZATIONMODEL_H
#define PYTHIAHADRONIZATIONMODEL_H

#include "ChiGenContext.h"
#include "PartonicModel.h"
#include "ColorSingletModels.h"
#include "EvtGenBase/EvtStdHep.hh"

using namespace Pythia8;

namespace chigen {

class PythiaChiGen {
 public:
  PythiaChiGen(PartonicModel &partonicModel, int maxTries = 30, bool doHadronization = true);

  PartonicModel &getPartonicModel() const { return partonicModel; }

  virtual ~PythiaChiGen();
  Event *next();

  ClassDef(PythiaChiGen, 1);

 private:
  double eCM;
  PartonicModel &partonicModel;
  int maxTries;
  double protonPz;
  bool doHadronization;
  Vec4 uMomentum, dMomentum;
  Vec4 ubarMomentum, dbarMomentum;
  Vec4 mesonMomentum, gluonMomentum;
  Vec4 boost;
  EvtStdHep evtstdhep;
  static int COLORS[];
  void calculateKinematics();
  void calculateHadronRemnants(double x, double hadronPz, Vec4 &u, Vec4 &d);

  void decay(Event &event);
  void fillColorEvent(Event &event);
};
} // namespace chigen

#endif /* PYTHIAHADRONIZATIONMODEL_H */