HepMCConverter.h

// -*- C++ -*-
//
// HepMCConverter.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2019 Leif Lonnblad
//
// ThePEG is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef ThePEG_HepMCConverter_H
#define ThePEG_HepMCConverter_H
// This is the declaration of the HepMCConverter class.

#include "ThePEG/Config/ThePEG.h"
#include "ThePEG/EventRecord/Event.h"
#include "HepMCTraits.h"

namespace ThePEG {


template <typename HepMCEventT, typename Traits = HepMCTraits<HepMCEventT> >
class HepMCConverter  {

public:

  struct HepMCConverterException: public Exception {};
  struct Vertex {
    tcParticleSet in;
    tcParticleSet out;
  };

  typedef typename Traits::ParticleT GenParticle;
  typedef typename Traits::ParticlePtrT GenParticlePtrT;
  typedef typename Traits::EventT GenEvent;
  typedef typename Traits::VertexT GenVertex;
  typedef typename Traits::VertexPtrT GenVertexPtrT;
  typedef typename Traits::PdfInfoT PdfInfo;
  typedef map<tcPPtr,GenParticlePtrT> ParticleMap;
  typedef map<tcColinePtr,long> FlowMap;
  typedef map<tcPPtr,Vertex*> VertexMap;
  typedef map<const Vertex *, GenVertexPtrT> GenVertexMap;

public:

  static GenEvent * convert(const Event & ev, bool nocopies = false,
                            Energy eunit = Traits::defaultEnergyUnit(),
                            Length lunit = Traits::defaultLengthUnit());

  static void
  convert(const Event & ev, GenEvent & gev, bool nocopies,
          Energy eunit, Length lunit);

  static void
  convert(const Event & ev, GenEvent & gev, bool nocopies = false);

private:

  HepMCConverter(const Event & ev, bool nocopies, Energy eunit, Length lunit);

  HepMCConverter(const Event & ev, GenEvent & gev, bool nocopies,
                 Energy eunit, Length lunit);

  void init(const Event & ev, bool nocopies);

  HepMCConverter();

  HepMCConverter(const HepMCConverter &);

  HepMCConverter & operator=(const HepMCConverter &) = delete;

private:

  GenParticlePtrT createParticle(tcPPtr p) const;

  void join(tcPPtr parent, tcPPtr child);

  GenVertexPtrT createVertex(Vertex * v); 

  void setPdfInfo(const Event & e);

private:

  GenEvent * geneve;

  ParticleMap pmap;

  FlowMap flowmap;

  vector<Vertex> vertices;

  VertexMap prov;

  VertexMap decv;

  GenVertexMap vmap;

  Energy energyUnit;

  Length lengthUnit;

};

}

#include "HepMCConverter.tcc"

#endif /* ThePEG_HepMCConverter_H */