EmissionVetoHook1.h

Go to the documentation of this file.

#include "GeneratorInterface/Pythia8Interface/plugins/PowhegHooksBB4L.h"

class EmissionVetoHook1 : public Pythia8::UserHooks {
public:
  // Constructor and destructor.
  EmissionVetoHook1(int nFinalIn,
                    bool vetoOnIn,
                    int vetoCountIn,
                    int pThardModeIn,
                    int pTemtModeIn,
                    int emittedModeIn,
                    int pTdefModeIn,
                    bool MPIvetoOnIn,
                    int QEDvetoModeIn,
                    int nFinalModeIn,
                    int VerbosityIn)
      : nFinalExt(nFinalIn),
        vetoOn(vetoOnIn),
        vetoCount(vetoCountIn),
        pThardMode(pThardModeIn),
        pTemtMode(pTemtModeIn),
        emittedMode(emittedModeIn),
        pTdefMode(pTdefModeIn),
        MPIvetoOn(MPIvetoOnIn),
        QEDvetoMode(QEDvetoModeIn),
        nFinalMode(nFinalModeIn),
        nISRveto(0),
        nFSRveto(0),
        Verbosity(VerbosityIn) {}
  ~EmissionVetoHook1() override {
    std::cout << "Number of ISR vetoed = " << nISRveto << std::endl;
    std::cout << "Number of FSR vetoed = " << nFSRveto << std::endl;
  }

  //--------------------------------------------------------------------------

  bool canVetoMPIStep() override { return true; }
  int numberVetoMPIStep() override { return 1; }
  bool doVetoMPIStep(int nMPI, const Pythia8::Event &e) override;

  bool canVetoISREmission() override { return vetoOn; }
  bool doVetoISREmission(int, const Pythia8::Event &e, int iSys) override;

  bool canVetoFSREmission() override { return vetoOn; }
  bool doVetoFSREmission(int, const Pythia8::Event &e, int iSys, bool) override;

  bool canVetoMPIEmission() override { return MPIvetoOn; }
  bool doVetoMPIEmission(int, const Pythia8::Event &e) override;

  void fatalEmissionVeto(std::string message);

  double pTpythia(const Pythia8::Event &e, int RadAfterBranch, int EmtAfterBranch, int RecAfterBranch, bool FSR);

  double pTpowheg(const Pythia8::Event &e, int i, int j, bool FSR);

  double pTcalc(const Pythia8::Event &e, int i, int j, int k, int r, int xSRin);

  //--------------------------------------------------------------------------

private:
  int nFinalExt, vetoOn, vetoCount, pThardMode, pTemtMode, emittedMode, pTdefMode, MPIvetoOn, QEDvetoMode, nFinalMode;
  int nFinal;
  double pThard, pTMPI;
  bool accepted, isEmt;
  // The number of accepted emissions (in a row)
  int nAcceptSeq;
  // Statistics on vetos
  unsigned long int nISRveto, nFSRveto;
  int Verbosity;
};