PythiaChiGen.h

Go to the documentation of this file.

/*
 * @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 */