Rapgap.hh

Go to the documentation of this file.

#ifndef AGILE_RAPGAP_HH 
#define AGILE_RAPGAP_HH 

#include "AGILe/Generator.hh"
#include "HepMC/IO_HEPEVT.h"

namespace AGILe {


  class Rapgap : public Generator {
  public:
    Rapgap();
    
    virtual ~Rapgap() { }
    
    void setGenSpecificInitialState(PdgCode p1 = PROTON, double e1 = 6000.0, 
                                    PdgCode p2 = PROTON, double e2 = 7000.0);

    using Generator::setParam;

    virtual bool setParam(const string& name, const string& value);
       
    virtual const string getName() const{
      return "Rapgap";
    }

    void makeEvent(HepMC::GenEvent& evt);


    //virtual void setSeed(const int value) {
      //setParam("ISEED", value);
    //}


    //string getPDFSet(PdgCode pid);
    //int getPDFMember(PdgCode pid);
    //string getPDFScheme(PdgCode pid)const;
    
    virtual const double getCrossSection();
    
    //virtual const string getVersion();
    
    void finalize();

  protected:
    void initialize();

    void fillEvent(HepMC::GenEvent& evt);

    string _particleName1, _particleName2;
    double _e1, _e2;

    int _nevt;

    HepMC::IO_HEPEVT _hepevt;
  };
  

}

#endif