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);
    };
}

#endif  /* PYTHIAHADRONIZATIONMODEL_H */