SampleXsection Class Reference

#include <SampleXsection.h>

Collaboration diagram for SampleXsection:

Public Types
enum	showering {   pythia =0, herwig =1, sherpa =2, pythia8 =4,   herwigpp =5, sherpa21 =6, amcatnlopythia8 =7, herwigpp713 =8,   sherpa228 =9, sherpa2210 =10, herwigpp721 =11, sherpa2212 =12,   amcatnloherwig7 =13, unknown =14 }

Public Member Functions
	SampleXsection ()
	SampleXsection (const SampleXsection &xs)
	~SampleXsection ()
SampleXsection &	operator= (const SampleXsection &xs)
bool	readFromFile (const char *fName)
void	setSample (const int dsid, const double xSect, const double kFact, const showering shower, const double xSectDw, const double xSectUp)
double	getRawXsection (const int dsid) const
double	getKfactor (const int dsid) const
double	getXsection (const int dsid) const
double	getXsectionUp (const int dsid) const
double	getXsectionDown (const int dsid) const
std::pair< double, double >	getXsectionDownUp (const int dsid) const
showering	getShowering (const int dsid) const
std::string	getShoweringString (const int dsid) const
int	getShoweringIndex (const int dsid) const
void	setTranslator (const std::unordered_map< std::string, std::string > &map)

Private Member Functions
showering	applyTranslation (const showering shower) const
std::string	showerToString (const showering shower) const
showering	stringToShower (const std::string &name) const

Private Attributes
std::map< int, std::pair< double, double > >	m_Xsects
std::map< int, std::pair< double, double > >	m_Uncert
std::map< int, showering >	m_Showering
std::unordered_map< std::string, std::string >	m_translator

Detailed Description

Definition at line 19 of file SampleXsection.h.

Member Enumeration Documentation

showering

enum SampleXsection::showering

Enumerator
pythia
herwig
sherpa
pythia8
herwigpp
sherpa21
amcatnlopythia8
herwigpp713
sherpa228
sherpa2210
herwigpp721
sherpa2212
amcatnloherwig7
unknown

Definition at line 21 of file SampleXsection.h.

                 {
    pythia=0,
    herwig=1,
    sherpa=2,
    pythia8=4,
    herwigpp=5,
    sherpa21=6,
    amcatnlopythia8=7,
    herwigpp713=8,
    sherpa228=9,
    sherpa2210=10,
    herwigpp721=11,
    sherpa2212=12,
    amcatnloherwig7=13,
    unknown=14
  };

Constructor & Destructor Documentation

SampleXsection() [1/2]

SampleXsection::SampleXsection ( )

inline

Definition at line 38 of file SampleXsection.h.

: m_Xsects() {}

SampleXsection() [2/2]

SampleXsection::SampleXsection ( const SampleXsection &  xs )

inline

Definition at line 39 of file SampleXsection.h.

: m_Xsects(xs.m_Xsects) {}

~SampleXsection()

SampleXsection::~SampleXsection ( )

inline

Definition at line 40 of file SampleXsection.h.

{}

Member Function Documentation

applyTranslation()

SampleXsection::showering SampleXsection::applyTranslation ( const showering  shower ) const

private

Definition at line 183 of file SampleXsection.cxx.

                                                                                                   {
  const std::string& stringName = showerToString(shower);
 
  // check the map
  auto itr = m_translator.find(stringName);
  if (itr == m_translator.end()) {
    return shower;
  } else {
    // need to apply the translation
    return stringToShower(itr->second);
  }
}

getKfactor()

double SampleXsection::getKfactor

(

const int

dsid

)

const

Definition at line 95 of file SampleXsection.cxx.

                                                      {
  map<int, pair<double, double> >::const_iterator it = m_Xsects.find(dsid);
  if (it != m_Xsects.end()) return it->second.second;
 
  return -1;
}

getRawXsection()

double SampleXsection::getRawXsection

(

const int

dsid

)

const

Definition at line 88 of file SampleXsection.cxx.

                                                          {
  map<int, pair<double, double> >::const_iterator it = m_Xsects.find(dsid);
  if (it != m_Xsects.end()) return it->second.first;
 
  return -1;
}

getShowering()

SampleXsection::showering SampleXsection::getShowering

(

const int

dsid

)

const

Definition at line 131 of file SampleXsection.cxx.

                                                                         {
  map<int, showering >::const_iterator it = m_Showering.find(dsid);
  if (it != m_Showering.end()) {
 
    showering shower = it->second;
    // apply the translation of set
    shower = applyTranslation(shower); 
 
    return shower;
  }
 
  return unknown;
}

getShoweringIndex()

int SampleXsection::getShoweringIndex

(

const int

dsid

)

const

Definition at line 150 of file SampleXsection.cxx.

                                                          {
  // S.Korn: labelling taken from here:
  // https://twiki.cern.ch/twiki/bin/view/AtlasProtected/BTagRecommendationsRelease22
 
  showering shower = getShowering(dsid);
 
  if (shower == showering::pythia) return 0;
  else if (shower == showering::herwigpp || shower == showering::herwig) // use the herwig++ maps for the few fortran
                                                                         // herwig samples we may have
    return 1;
  else if (shower == showering::pythia8) return 2;
  else if (shower == showering::sherpa) return 3;
  else if (shower == showering::sherpa21) return 4;
  else if (shower == showering::amcatnlopythia8) return 5;
  else if (shower == showering::herwigpp713) return 6;
  else if (shower == showering::sherpa228) return 7;
  else if (shower == showering::sherpa2210) return 8;
  else if (shower == showering::herwigpp721) return 9;
  else if (shower == showering::sherpa2212) return 10;
  else if (shower == showering::amcatnloherwig7) return 11;
  else {
    std::cout <<
    "==========================================================================================================================================="
              << std::endl;
    std::cout << "ERROR:TopDataPreparation: The DSID " << dsid <<
    " does not have a showering algorithm defined!!! This is needed for MC/MC SF!!! ---> EXIT." << std::endl;
    std::cout << "\t" << std::endl;
    exit(1);
  }
 
  return 0;
}

getShoweringString()

std::string SampleXsection::getShoweringString

(

const int

dsid

)

const

Definition at line 145 of file SampleXsection.cxx.

{
  return showerToString(getShowering(dsid));
}

getXsection()

double SampleXsection::getXsection

(

const int

dsid

)

const

Definition at line 102 of file SampleXsection.cxx.

                                                       {
  map<int, pair<double, double> >::const_iterator it = m_Xsects.find(dsid);
  if (it != m_Xsects.end()) return (it->second.first) * (it->second.second);
 
  return -1;
}

getXsectionDown()

double SampleXsection::getXsectionDown

(

const int

dsid

)

const

Definition at line 109 of file SampleXsection.cxx.

                                                           {
  map<int, pair<double, double> >::const_iterator it = m_Uncert.find(dsid);
  if (it != m_Uncert.end()) return it->second.first;
 
  return -1;
}

getXsectionDownUp()

pair< double, double > SampleXsection::getXsectionDownUp

(

const int

dsid

)

const

Definition at line 123 of file SampleXsection.cxx.

                                                                           {
  map<int, pair<double, double> >::const_iterator it = m_Uncert.find(dsid);
  if (it != m_Uncert.end()) return it->second;
 
  pair<double, double> temp = pair<double, double>(-1, -1);
  return temp;
}

getXsectionUp()

double SampleXsection::getXsectionUp

(

const int

dsid

)

const

Definition at line 116 of file SampleXsection.cxx.

                                                         {
  map<int, pair<double, double> >::const_iterator it = m_Uncert.find(dsid);
  if (it != m_Uncert.end()) return it->second.second;
 
  return -1;
}

operator=()

SampleXsection & SampleXsection::operator=

(

const SampleXsection &

xs

)

Definition at line 15 of file SampleXsection.cxx.

                                                                    {
  if (this != &xs) {
    m_Xsects = xs.m_Xsects;
    m_Uncert = xs.m_Uncert;
  }
  return *this;
}

readFromFile()

bool SampleXsection::readFromFile

(

const char *

fName

)

Definition at line 23 of file SampleXsection.cxx.

                                                   {
  ifstream in(fName);
 
  if (!in) return false;
 
  for (; !in.eof(); ) {
    string line;
    if (!getline(in, line)) break;
    if (!line.empty() && line[0] != '#') {
      istringstream istr(line);
      int dsid = -1;
      double xSect, kFact, xSectDw, xSectUp;
      string s_shower, s_xSectDw, s_xSectUp;
      showering shower;
      istr >> dsid >> xSect >> kFact >> s_shower >> s_xSectDw >> s_xSectUp;
      if (!s_shower.empty() && s_shower[0] != '#') {
        shower = stringToShower(s_shower);
        if (shower == showering::unknown) {
          cerr <<
          "ERROR!! TopDataPreparation::SampleXsection::readFromFile: unknown showering (which is needed for btagging SF!!!) : "
               << s_shower << " for DSID= " << dsid << endl;
        }
 
        if (!s_xSectDw.empty() && s_xSectDw[0] != '#') {
          xSectDw = atof(s_xSectDw.c_str());
          if (!s_xSectUp.empty()) {
            xSectUp = atof(s_xSectUp.c_str());
          } else {
            cerr << "ERROR!! SampleXsection::readFromFile: xSectDw of " << dsid << " is defined ( " << s_xSectDw <<
            " ), but not xSectUp" << endl;
            xSectUp = -1.;
          }
        } else {
          xSectDw = -1.;
          xSectUp = -1.;
        }
      }// !s_showering.empty()
      else {
        xSectDw = -1.;
        xSectUp = -1.;
        shower = unknown;
      }
      setSample(dsid, xSect, kFact, shower, xSectDw, xSectUp);
    }
  }
 
  in.close();
  return true;
}

setSample()

void SampleXsection::setSample	(	const int	dsid,
		const double	xSect,
		const double	kFact,
		const showering	shower,
		const double	xSectDw,
		const double	xSectUp
	)

Definition at line 73 of file SampleXsection.cxx.

                                                                           {
  if (dsid < 0) return;
 
  map<int, pair<double, double> >::const_iterator it = m_Xsects.find(dsid);
  if (it != m_Xsects.end()) {
    cerr << "ERROR!! SampleXsection::setSample: duplicate entry ! DSID= " << dsid << endl;
    //cerr << "ERROR!! SampleXsection::setSample:                   xSect=" << xSect << "\told xSect=" <<
    // it->second.first <<endl;
  }
  m_Xsects[dsid] = pair<double, double>(xSect, kFact);
  m_Uncert[dsid] = pair<double, double>(xSectDw, xSectUp);
  m_Showering[dsid] = shower;
}

setTranslator()

void SampleXsection::setTranslator ( const std::unordered_map< std::string, std::string > &  map )

inline

Definition at line 57 of file SampleXsection.h.

{m_translator = map;}

showerToString()

std::string SampleXsection::showerToString ( const showering  shower ) const

private

Definition at line 196 of file SampleXsection.cxx.

                                                                                   {
  if (shower == showering::pythia) return "pythia";
  else if (shower == showering::herwig) return "herwig";
  else if (shower == showering::herwigpp) return "herwig";
  else if (shower == showering::sherpa) return "sherpa";
  else if (shower == showering::sherpa21) return "sherpa21";
  else if (shower == showering::pythia8) return "pythia8";
  else if (shower == showering::amcatnlopythia8) return "amcatnlopythia8";
  else if (shower == showering::herwigpp713) return "herwigpp713";
  else if (shower == showering::sherpa228) return "sherpa228";
  else if (shower == showering::sherpa2210) return "sherpa2210";
  else if (shower == showering::herwigpp721) return "herwigpp721";
  else if (shower == showering::sherpa2212) return "sherpa2212";
  else if (shower == showering::amcatnloherwig7) return "amcatnloherwig7";
 
  return "";
}

stringToShower()

SampleXsection::showering SampleXsection::stringToShower ( const std::string &  name ) const

private

Definition at line 214 of file SampleXsection.cxx.

                                                                                  {
  if (name == "pythia") return showering::pythia;
  else if (name == "herwig") return showering::herwig;
  else if (name == "sherpa") return showering::sherpa;
  else if (name == "sherpa21") return showering::sherpa21;
  else if (name == "pythia8") return showering::pythia8;
  else if (name == "herwigpp") return showering::herwigpp;
  else if (name == "amcatnlopythia8") return showering::amcatnlopythia8;
  else if (name == "herwigpp713") return showering::herwigpp713;
  else if (name == "sherpa228") return showering::sherpa228;
  else if (name == "sherpa2210") return showering::sherpa2210;
  else if (name == "herwigpp721") return showering::herwigpp721;
  else if (name == "sherpa2212") return showering::sherpa2212;
  else if (name == "amcatnloherwig7") return showering::amcatnloherwig7;
 
  return showering::unknown;
}

Member Data Documentation

m_Showering

std::map<int, showering > SampleXsection::m_Showering

private

Definition at line 62 of file SampleXsection.h.

m_translator

std::unordered_map<std::string, std::string> SampleXsection::m_translator

private

Definition at line 63 of file SampleXsection.h.

m_Uncert

std::map<int, std::pair<double, double> > SampleXsection::m_Uncert

private

Definition at line 61 of file SampleXsection.h.

m_Xsects

std::map<int, std::pair<double, double> > SampleXsection::m_Xsects

private

Definition at line 60 of file SampleXsection.h.

---

The documentation for this class was generated from the following files:

• SampleXsection.h
• SampleXsection.cxx