/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_7/src/GeneratorInterface/ThePEGInterface/interface/HepMCConverter.h

Go to the documentation of this file.

// -*- C++ -*-
//
// HepMCConverter.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2007 Leif Lonnblad
//
// ThePEG is licenced under version 2 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 <ThePEG/Vectors/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::EventT GenEvent;
  typedef typename Traits::VertexT GenVertex;
  typedef typename Traits::PdfInfoT PdfInfo;
  typedef map<tcPPtr,GenParticle*> ParticleMap;
  typedef map<tcColinePtr,long> FlowMap;
  typedef map<tcPPtr,Vertex*> VertexMap;
  typedef map<const Vertex *, GenVertex*> 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 &);

private:

  GenParticle * createParticle(tcPPtr p) const;

  void join(tcPPtr parent, tcPPtr child);

  GenVertex * createVertex(Vertex * v); 

  void setPdfInfo(const Event & e);

  void sortTopologically(tcPVector & allv);
private:

  GenEvent * geneve;

  ParticleMap pmap;

  FlowMap flowmap;

  vector<Vertex> vertices;

  VertexMap prov;

  VertexMap decv;

  GenVertexMap vmap;

  Energy energyUnit;

  Length lengthUnit;

  tcEHPtr eh;
};

        template<> struct HepMCTraits<HepMC::GenEvent> :
                public HepMCTraitsBase<
                        HepMC::GenEvent, HepMC::GenParticle,
                        HepMC::GenVertex, HepMC::Polarization,
                        HepMC::PdfInfo> {};
}


#endif /* ThePEG_HepMCConverter_H */