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