GoSamAmplitude.h

// -*- C++ -*-
//
// GoSamAmplitude.h is a part of Herwig - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2019 The Herwig Collaboration
//
// Herwig is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef Herwig_GoSamAmplitude_H
#define Herwig_GoSamAmplitude_H
//
// This is the declaration of the GoSamAmplitude class.
//
 
#include "Herwig/MatrixElement/Matchbox/Base/MatchboxOLPME.h"
#include "ThePEG/Utilities/DynamicLoader.h"
 
 
namespace Herwig {
 
using namespace ThePEG;
 
 
class gosamprocinfo{
 
        public:
                gosamprocinfo(){};
                gosamprocinfo(int HID,int GID, string procstr,string typestr):
                        theHOlpId(HID),theGOlpId(GID),theProcstr(procstr),theTypestr(typestr){
                }
                ~gosamprocinfo(){}
                int HID() const {return theHOlpId;}
                int GID() const {return theGOlpId;}
                string Pstr() const {return theProcstr;}
                string Tstr() const {return theTypestr;}
                void setGID(int g){theGOlpId=g;}
                void setOAs(int i){ orderAlphas=i;}
                int orderAs(){return orderAlphas;}
                void setOAew(int j){ orderAlphaew=j;}
                int orderAew(){return orderAlphaew;}
        private:
                int theHOlpId;
                int theGOlpId;
                string theProcstr;
                string theTypestr;
                int orderAlphas;
                int orderAlphaew;
        public:
                void persistentOutput(PersistentOStream & os) const{os<<theHOlpId<<theGOlpId<<theProcstr<<theTypestr<<orderAlphas<<orderAlphaew;}
                void persistentInput(PersistentIStream &is) {is>>theHOlpId>>theGOlpId>>theProcstr>>theTypestr>>orderAlphas>>orderAlphaew;}
};
 
 
class GoSamAmplitude: public MatchboxOLPME {
 
public:
 
  GoSamAmplitude();
 
public:
 
  virtual void fillOrderFile(const map<pair<Process,int>,int>& procs, string OrderFileName);
  
  virtual bool isCS() const { return false; }
  virtual bool isExpanded() const { return true; }
  virtual bool isBDK() const { return false; }
  virtual bool isDR() const { return isitDR; }
  virtual bool isDRbar() const {return false;}
 
  virtual void signOLP(const string&, const string&);
 
  virtual bool checkOLPContract(string contractFileName);
 
  // virtual bool hasFinalStateSymmetry() const { return true; }
  virtual bool hasFinalStateSymmetry() const { return false; }
  
  virtual bool buildGoSam();
 
  virtual void startOLP(const string&, int& status);
 
  // virtual Energy2 mu2() const { return lastSHat(); }
 
  virtual bool startOLP(const map<pair<Process,int>,int>& procs);
 
  void evalSubProcess() const;
 
  virtual void evalColourCorrelator(pair<int,int> ij) const;
 
  virtual LorentzVector<Complex> plusPolarization(const Lorentz5Momentum& p,
                                                  const Lorentz5Momentum& n, int inc) const;
 
  virtual void evalSpinColourCorrelator(pair<int,int> ij) const;
 
  void getids() const;
 
  int accuracyTargetNegExp() const { return theAccuracyTarget; }
 
public:
 
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
 
  static void Init();
 
 
protected:
 
 
  virtual void doinit();
 
  virtual void doinitrun();
 
 
protected:
 
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
 
 
// If needed, insert declarations of virtual function defined in the
// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).
 
 
private:
 
  GoSamAmplitude & operator=(const GoSamAmplitude &) = delete;
 
  mutable vector<double> colourCorrelatorResults;
 
  mutable vector<double> spinColourCorrelatorResults;
 
  mutable vector<int> idpair;
 
  mutable vector<string> idtypepair;
 
  map<int , gosamprocinfo > processmap;
 
  mutable string gosamPathInterface;
  mutable string gosamSetupInFileNameInterface;
  mutable string gosamBuildScript;
 
  mutable string gosamPath;
  mutable string gosamSourcePath;
  mutable string gosamInstallPath;
 
  mutable string gosamSetupInFileName;
 
  mutable string orderFileTitle;
  mutable string contractFileTitle;
  mutable string contractFileName;
  mutable string orderFileName;
 
  mutable string accuracyFileTitle;
  mutable string accuracyFile;
  int theAccuracyTarget;
 
  bool theCodeExists;
  bool theFormOpt;
  bool theNinja;
  bool theHiggsEff;
  bool theMassiveLeptons;
 
  int theLoopInducedOption;
 
  bool isitDR;
 
  mutable bool doneGoSamInit;
  mutable bool doneGoSamInitRun;
  
  vector<int> massiveParticles; //theMassiveParticles;
 
  void setupGoSamIn(string setupGoSamInFile);
 
protected:
 
  string bindir_;
 
  string pkgdatadir_;
 
  string GoSamPrefix_;
 
}; // end "class GoSamAmplitude: public MatchboxOLPME"
 
  inline PersistentOStream& operator<<(PersistentOStream& os, const gosamprocinfo& p) {
    p.persistentOutput(os); return os;
  }
 
  inline PersistentIStream& operator>>(PersistentIStream& is, gosamprocinfo& p) {
    p.persistentInput(is); return is;
  }
 
} // end "namespace Herwig"
 
#endif /* Herwig_GoSamAmplitude_H */