AlpgenPar Class Reference

#include <GeneratorInput.h>

Public Member Functions
	AlpgenPar (Pythia8::Info *infoPtrIn=NULL)
void	extractRunParam (string line)
double	getParam (const string &paramIn)
int	getParamAsInt (const string &paramIn)
bool	haveParam (const string &paramIn)
bool	parse (const string paramStr)
void	printParams ()

Private Member Functions
void	warnParamOverwrite (const string &paramIn, double val)

Static Private Member Functions
static string	trim (string s)

Private Attributes
Pythia8::Info *	infoPtr
map< string, double >	params

Static Private Attributes
static const double	ZEROTHRESHOLD = 1e-10

Detailed Description

Definition at line 33 of file GeneratorInput.h.

Constructor & Destructor Documentation

AlpgenPar::AlpgenPar ( Pythia8::Info *  infoPtrIn = NULL )

inline

Definition at line 38 of file GeneratorInput.h.

: infoPtr(infoPtrIn) {}

Member Function Documentation

void AlpgenPar::extractRunParam

(

string

line

)

Definition at line 62 of file GeneratorInput.cc.

References customizeTrackingMonitorSeedNumber::idx, and nevent.

                                           {

  // Extract information to the right of the final '!' character
  size_t idx = line.rfind("!");
  if (idx == string::npos) return;
  string paramName = trim(line.substr(idx + 1));
  string paramVal  = trim(line.substr(0, idx));
  istringstream iss(paramVal);

  // Special case: 'hard process code' - single integer input
  double val;
  if (paramName == "hard process code") {
    iss >> val;
    warnParamOverwrite("hpc", val);
    params["hpc"] = val;

  // Special case: 'Crosssection +- error (pb)' - two double values
  } else if (paramName.find("Crosssection") == 0) {
    double xerrup;
    iss >> val >> xerrup;
    warnParamOverwrite("xsecup", val);
    warnParamOverwrite("xerrup", val);
    params["xsecup"] = val;
    params["xerrup"] = xerrup;

  // Special case: 'unwtd events, lum (pb-1)' - integer and double values
  } else if (paramName.find("unwtd events") == 0) {
    int nevent;
    iss >> nevent >> val;
    warnParamOverwrite("nevent", val);
    warnParamOverwrite("lum", val);
    params["nevent"] = nevent;
    params["lum"]    = val;

  // Special case: 'mc,mb,...' - split on ',' for name and ' ' for values
  } else if (paramName.find(",") != string::npos) {
    
    // Simple tokeniser
    string        paramNameNow;
    istringstream issName(paramName);
    while (getline(issName, paramNameNow, ',')) {
      iss >> val;
      warnParamOverwrite(paramNameNow, val);
      params[paramNameNow] = val;
    }

  // Default case: assume integer and double on the left
  } else {
    int paramIdx;
    iss >> paramIdx >> val;
    warnParamOverwrite(paramName, val);
    params[paramName] = val;
  }
}

double AlpgenPar::getParam ( const string &  paramIn )

inline

Definition at line 52 of file GeneratorInput.h.

References haveParam(), and params.

Referenced by AlpgenHooks::initAfterBeams(), and LHAupAlpgen::setInit().

                                         {
    return (haveParam(paramIn)) ? params[paramIn] : 0.; }

int AlpgenPar::getParamAsInt ( const string &  paramIn )

inline

Definition at line 54 of file GeneratorInput.h.

References haveParam(), and params.

Referenced by AlpgenHooks::initAfterBeams(), and LHAupAlpgen::setInit().

                                              {
    return (haveParam(paramIn)) ? int(params[paramIn]) : 0.; }

bool AlpgenPar::haveParam ( const string &  paramIn )

inline

Definition at line 47 of file GeneratorInput.h.

References params.

Referenced by getParam(), getParamAsInt(), AlpgenHooks::initAfterBeams(), and LHAupAlpgen::setInit().

                                        {
    return (params.find(paramIn) == params.end()) ? false : true; }

bool AlpgenPar::parse

(

const string

paramStr

)

Definition at line 24 of file GeneratorInput.cc.

References createPayload::block, and geometryCSVtoXML::line.

Referenced by AlpgenHooks::initAfterBeams(), and LHAupAlpgen::LHAupAlpgen().

                                           {

  // Read par file in blocks:
  //   0 - process information
  //   1 - run parameters
  //   2 - cross sections
  int block = 0;

  // Loop over incoming lines
  stringstream paramStream(paramStr);
  string line;
  while (getline(paramStream, line)) {

    // Change to 'run parameters' block
    if        (line.find("run parameters") != string::npos) {
      block = 1;

    // End of 'run parameters' block
    } else if (line.find("end parameters") != string::npos) {
      block = 2;

    // Do not extract anything from block 0 so far
    } else if (block == 0) {

    // Block 1 or 2: extract parameters
    } else {
      extractRunParam(line);

    }
  } // while (getline(paramStream, line))

  return true;
}

void AlpgenPar::printParams

(

)

Definition at line 121 of file GeneratorInput.cc.

References gather_cfg::cout.

Referenced by AlpgenHooks::initAfterBeams().

                            {

  // Loop over all stored parameters and print
  cout << fixed << setprecision(3) << endl
       << " *-------  Alpgen parameters  -------*" << endl;
  for (map < string, double >::iterator it = params.begin();
       it != params.end(); ++it)
    cout << " |  " << left << setw(13) << it->first
         << "  |  " << right << setw(13) << it->second
         << "  |" << endl;
  cout << " *-----------------------------------*" << endl;
}

string AlpgenPar::trim ( string  s )

staticprivate

Definition at line 152 of file GeneratorInput.cc.

References i, and alignCSCRings::s.

                               {

  // Remove whitespace in incoming string
  size_t i;
  if ((i = s.find_last_not_of(" \t\r\n")) != string::npos)
    s = s.substr(0, i + 1);
  if ((i = s.find_first_not_of(" \t\r\n")) != string::npos)
    s = s.substr(i);
  return s;
}

void AlpgenPar::warnParamOverwrite ( const string &  paramIn, double  val  )

private

Definition at line 138 of file GeneratorInput.cc.

References funct::abs().

                                                                    {

  // Check if present and if new value is different
  if (haveParam(paramIn) &&
      abs(getParam(paramIn) - val) > ZEROTHRESHOLD) {
    if (infoPtr) infoPtr->errorMsg("Warning in LHAupAlpgen::"
        "warnParamOverwrite: overwriting existing parameter", paramIn);
  }
}

Member Data Documentation

Pythia8::Info* AlpgenPar::infoPtr

private

Definition at line 72 of file GeneratorInput.h.

map<string,double> AlpgenPar::params

private

Definition at line 69 of file GeneratorInput.h.

Referenced by getParam(), getParamAsInt(), and haveParam().

const double AlpgenPar::ZEROTHRESHOLD = 1e-10

staticprivate

Definition at line 75 of file GeneratorInput.h.