AlpgenHooks Class Reference

#include <GeneratorInput.h>

Inheritance diagram for AlpgenHooks:

Public Member Functions
	AlpgenHooks (Pythia8::Pythia &pythia)
bool	initAfterBeams ()
	~AlpgenHooks ()

Private Attributes
LHAupAlpgen *	LHAagPtr

Detailed Description

Definition at line 134 of file GeneratorInput.h.

Constructor & Destructor Documentation

AlpgenHooks::AlpgenHooks

(

Pythia8::Pythia &

pythia

)

Definition at line 853 of file GeneratorInput.cc.

References LHAagPtr.

                                       : LHAagPtr(NULL) {

  // If LHAupAlpgen needed, construct and pass to Pythia
  string agFile = pythia.settings.word("Alpgen:file");
  if (agFile != "void") {
    LHAagPtr = new LHAupAlpgen(agFile.c_str(), &pythia.info);
    pythia.settings.mode("Beams:frameType", 5);
    pythia.setLHAupPtr(LHAagPtr);
  }
}

AlpgenHooks::~AlpgenHooks ( )

inline

Definition at line 140 of file GeneratorInput.h.

References LHAagPtr.

{ if (LHAagPtr) delete LHAagPtr; }

Member Function Documentation

bool AlpgenHooks::initAfterBeams

(

)

Definition at line 872 of file GeneratorInput.cc.

References AlpgenPar::getParam(), AlpgenPar::getParamAsInt(), AlpgenPar::haveParam(), max(), AlpgenPar::parse(), and AlpgenPar::printParams().

Referenced by JetMatchingAlpgenInputAlpgen::initAfterBeams(), and JetMatchingMadgraphInputAlpgen::initAfterBeams().

                                 {

  // Read in ALPGEN specific configuration variables
  bool setMasses = settingsPtr->flag("Alpgen:setMasses");
  bool setNjet   = settingsPtr->flag("Alpgen:setNjet");
  bool setMLM    = settingsPtr->flag("Alpgen:setMLM");

  // If ALPGEN parameters are present, then parse in AlpgenPar object
  AlpgenPar par(infoPtr);
  string parStr = infoPtr->header("AlpgenPar");
  if (!parStr.empty()) {
    par.parse(parStr);
    par.printParams();
  }

  // Set masses if requested
  if (setMasses) {
    if (par.haveParam("mc")) particleDataPtr->m0(4,  par.getParam("mc"));
    if (par.haveParam("mb")) particleDataPtr->m0(5,  par.getParam("mb"));
    if (par.haveParam("mt")) particleDataPtr->m0(6,  par.getParam("mt"));
    if (par.haveParam("mz")) particleDataPtr->m0(23, par.getParam("mz"));
    if (par.haveParam("mw")) particleDataPtr->m0(24, par.getParam("mw"));
    if (par.haveParam("mh")) particleDataPtr->m0(25, par.getParam("mh"));
  }

  // Set MLM:nJets if requested
  if (setNjet) {
    if (par.haveParam("njets"))
      settingsPtr->mode("JetMatching:nJet", par.getParamAsInt("njets"));
    else
      infoPtr->errorMsg("Warning in AlpgenHooks:init: "
          "no ALPGEN nJet parameter found");
  }

  // Set MLM merging parameters if requested
  if (setMLM) {
    if (par.haveParam("ptjmin") && par.haveParam("drjmin") &&
        par.haveParam("etajmax")) {
      double ptjmin = par.getParam("ptjmin");
      ptjmin = max(ptjmin + 5., 1.2 * ptjmin);
      settingsPtr->parm("JetMatching:eTjetMin",   ptjmin);
      settingsPtr->parm("JetMatching:coneRadius", par.getParam("drjmin"));
      settingsPtr->parm("JetMatching:etaJetMax",  par.getParam("etajmax"));

    // Warn if setMLM requested, but parameters not present
    } else {
      infoPtr->errorMsg("Warning in AlpgenHooks:init: "
          "no ALPGEN merging parameters found");
    }
  }

  // Initialisation complete.
  return true;
}

Member Data Documentation

LHAupAlpgen* AlpgenHooks::LHAagPtr

private

Definition at line 148 of file GeneratorInput.h.

Referenced by AlpgenHooks(), and ~AlpgenHooks().