#include <SUSYCrossSection.h>

Collaboration diagram for SUSY::CrossSectionDB:

Classes
class	Key
class	Process

Public Types
typedef xsDB_t::const_iterator	iterator

Public Member Functions
	CrossSectionDB (const std::string &txtfilename="dev/PMGTools/PMGxsecDB_mc16.txt", bool usePathResolver=true, bool isExtended=false, bool usePMGTool=true)
void	loadFile (const std::string &)
void	extend (const std::string &)
void	setExtended (bool isExtended=true)
void	setUsePMGTool (bool usePMGTool=true)
Process	process (int id, int proc=0) const
Process	process (int id, int pdgId1, int pdgId2) const
float	xsectTimesEff (int id, int proc=0) const
float	xsectTimesEff (int id, int pdgId1, int pdgId2) const
float	rawxsect (int id, int proc=0) const
float	rawxsect (int id, int pdgId1, int pdgId2) const
float	kfactor (int id, int proc=0) const
float	kfactor (int id, int pdgId1, int pdgId2) const
float	efficiency (int id, int proc=0) const
float	efficiency (int id, int pdgId1, int pdgId2) const
float	rel_uncertainty (int id, int proc=0) const
float	rel_uncertainty (int id, int pdgId1, int pdgId2) const
float	sumweight (int id, int proc=0) const
float	sumweight (int id, int pdgId1, int pdgId2) const
std::string	name (int id) const
iterator	begin () const
iterator	end () const

Private Types
typedef std::map< Key, Process >	xsDB_t

Private Member Functions
xsDB_t::iterator	my_find (const int proc)

Private Attributes
bool	m_extended
bool	m_usePMGTool
asg::AnaToolHandle< PMGTools::IPMGCrossSectionTool >	m_pmgxs
xsDB_t	m_xsectDB
xsDB_t	m_cache

Detailed Description

Definition at line 28 of file SUSYCrossSection.h.

Member Typedef Documentation

◆ iterator

typedef xsDB_t::const_iterator SUSY::CrossSectionDB::iterator

Definition at line 119 of file SUSYCrossSection.h.

◆ xsDB_t

typedef std::map<Key, Process> SUSY::CrossSectionDB::xsDB_t

private

Definition at line 115 of file SUSYCrossSection.h.

Constructor & Destructor Documentation

◆ CrossSectionDB()

SUSY::CrossSectionDB::CrossSectionDB	(	const std::string &	txtfilename = "dev/PMGTools/PMGxsecDB_mc16.txt",
		bool	usePathResolver = true,
		bool	isExtended = false,
		bool	usePMGTool = true )

Definition at line 113 of file SUSYCrossSection.cxx.

  : m_pmgxs("")
{
  setExtended(isExtended);
  setUsePMGTool(usePMGTool);
 
  // configuring the PMG tool... 
  m_pmgxs.setTypeAndName("PMGTools::PMGCrossSectionTool/PMGCrossSectionTool");
  m_pmgxs.retrieve().ignore(); // Ignore the status code
  std::vector<std::string> inFiles={};
  std::string fullPath = usePathResolver ? PathResolverFindCalibFile(txtfilename) : txtfilename;
  std::ifstream in(fullPath.c_str());
  if (!in) return;
  inFiles.push_back(fullPath);
  m_pmgxs->readInfosFromFiles( inFiles );
 
}

Member Function Documentation

◆ begin()

iterator SUSY::CrossSectionDB::begin ( ) const

inline

Definition at line 120 of file SUSYCrossSection.h.

120{ return m_xsectDB.begin(); }

SUSY::CrossSectionDB::m_xsectDB

xsDB_t m_xsectDB

Definition SUSYCrossSection.h:128

◆ efficiency() [1/2]

float SUSY::CrossSectionDB::efficiency	(	int	id,
		int	pdgId1,
		int	pdgId2 ) const

inline

Definition at line 104 of file SUSYCrossSection.h.

104{ return efficiency(id, finalState(pdgId1, pdgId2)); }

SUSY::CrossSectionDB::efficiency

float efficiency(int id, int proc=0) const

Definition SUSYCrossSection.h:103

SUSY::finalState

unsigned int finalState(const int SUSY_Spart1_pdgId, const int SUSY_Spart2_pdgId)

Definition SUSYCrossSection.cxx:227

◆ efficiency() [2/2]

float SUSY::CrossSectionDB::efficiency	(	int	id,
		int	proc = 0 ) const

inline

Definition at line 103 of file SUSYCrossSection.h.

103{ return process(id, proc).efficiency(); }

process

const std::string process

Definition fbtTestBasics.cxx:76

◆ end()

iterator SUSY::CrossSectionDB::end ( ) const

inline

Definition at line 121 of file SUSYCrossSection.h.

121{return m_xsectDB.end(); }

◆ extend()

void SUSY::CrossSectionDB::extend ( const std::string & txtfilename )

Definition at line 167 of file SUSYCrossSection.cxx.

                                                           {
  // Just like the above function, but with more functionality
  std::string line;
 
  std::ifstream in(txtfilename.c_str());
  if (!in) return;
  while ( getline(in, line) )
    {
      // skip leading blanks (in case there are some in front of a comment)
      if ( !line.empty() ){
          while ( line[0] == ' ' ) line.erase(0, 1);
      }
      // skip lines that do not start with a number, they are comments
      if ( !line.empty() && isdigit(line[0]) ){
          std::stringstream is(line);
          int id;
          std::string name;
          float xsect, kfactor, efficiency, relunc;
          float sumweight = -1, stat = -1;
          is >> id;
          auto my_it = my_find( id );
          if (my_it==m_xsectDB.end()){
            is >> name >> xsect >> kfactor >> efficiency >> relunc;
            if (m_extended == true){
              is >> sumweight >> stat;
            }
            m_xsectDB[Key(id, name)] = Process(id, name, xsect, kfactor, efficiency, relunc, sumweight, stat);
          } else {
            // Now we have extended records
            if (!m_extended) m_extended=true;
            is >> sumweight >> stat;
            my_it->second.sumweight(sumweight);
            my_it->second.stat(stat);
          }
        }
    }
 
}

◆ kfactor() [1/2]

float SUSY::CrossSectionDB::kfactor	(	int	id,
		int	pdgId1,
		int	pdgId2 ) const

inline

Definition at line 102 of file SUSYCrossSection.h.

102{ return kfactor(id, finalState(pdgId1, pdgId2)); }

◆ kfactor() [2/2]

float SUSY::CrossSectionDB::kfactor	(	int	id,
		int	proc = 0 ) const

inline

Definition at line 101 of file SUSYCrossSection.h.

101{ return process(id, proc).kfactor(); }

◆ loadFile()

void SUSY::CrossSectionDB::loadFile ( const std::string & txtfilename )

Definition at line 131 of file SUSYCrossSection.cxx.

                                                             {
 
  std::string line;
  
  std::ifstream in(txtfilename.c_str());
  if (!in) return;
  while ( getline(in, line) ){
    // skip leading blanks (in case there are some in front of a comment)
    if ( !line.empty() ){
      while ( line[0] == ' ' ) line.erase(0, 1);
    }
    // skip lines that do not start with a number, they are comments
    if ( !line.empty() && isdigit(line[0]) ){
      std::stringstream is(line);
      int id;
      std::string name;
      float xsect, kfactor, efficiency, relunc;
      float sumweight = -1, stat = -1;
      is >> id >> name >> xsect >> kfactor >> efficiency >> relunc;
      if (m_extended == true){
          is >> sumweight >> stat;
      }
      m_xsectDB[Key(id, name)] = Process(id, name, xsect, kfactor, efficiency, relunc, sumweight, stat);
    }
  }
}

◆ my_find()

SUSY::CrossSectionDB::xsDB_t::iterator SUSY::CrossSectionDB::my_find ( const int proc )

private

Definition at line 159 of file SUSYCrossSection.cxx.

                                                                             {
  for (SUSY::CrossSectionDB::xsDB_t::iterator it = m_xsectDB.begin(); it!=m_xsectDB.end();++it){
    if (it->second.ID()==proc) return it;
  }
  return m_xsectDB.end();
}

◆ name()

std::string SUSY::CrossSectionDB::name ( int id ) const

inline

Definition at line 111 of file SUSYCrossSection.h.

111{ return process(id, 0).name(); }

◆ process() [1/2]

Process SUSY::CrossSectionDB::process	(	int	id,
		int	pdgId1,
		int	pdgId2 ) const

inline

Definition at line 91 of file SUSYCrossSection.h.

91{ return process(id, finalState(pdgId1, pdgId2)); }

◆ process() [2/2]

SUSY::CrossSectionDB::Process SUSY::CrossSectionDB::process	(	int	id,
		int	proc = 0 ) const

Definition at line 206 of file SUSYCrossSection.cxx.

{
  // for background x-sections, use the PMG tool
  if(proc==0 && m_usePMGTool) {
    return Process( id, m_pmgxs->getSampleName(id), m_pmgxs->getAMIXsection(id), m_pmgxs->getKfactor(id), m_pmgxs->getFilterEff(id), m_pmgxs->getXsectionUncertainty(id), -1, -1 );  
  } else {
    const Key k(id, proc);
    xsDB_t::const_iterator pos = m_cache.find(k);
    if (pos != m_cache.end()) {
      return pos->second;
    } else {
      pos = m_xsectDB.find(k);
      if (pos != m_xsectDB.end()) {
    return pos->second;
      }
    }
  }
  return Process();
}

◆ rawxsect() [1/2]

float SUSY::CrossSectionDB::rawxsect	(	int	id,
		int	pdgId1,
		int	pdgId2 ) const

inline

Definition at line 100 of file SUSYCrossSection.h.

100{ return rawxsect(id, finalState(pdgId1, pdgId2)); }

SUSY::CrossSectionDB::rawxsect

float rawxsect(int id, int proc=0) const

Definition SUSYCrossSection.h:99

◆ rawxsect() [2/2]

float SUSY::CrossSectionDB::rawxsect	(	int	id,
		int	proc = 0 ) const

inline

Definition at line 99 of file SUSYCrossSection.h.

99{ return process(id, proc).xsect(); }

◆ rel_uncertainty() [1/2]

float SUSY::CrossSectionDB::rel_uncertainty	(	int	id,
		int	pdgId1,
		int	pdgId2 ) const

inline

Definition at line 106 of file SUSYCrossSection.h.

106{ return rel_uncertainty(id, finalState(pdgId1, pdgId2)); }

SUSY::CrossSectionDB::rel_uncertainty

float rel_uncertainty(int id, int proc=0) const

Definition SUSYCrossSection.h:105

◆ rel_uncertainty() [2/2]

float SUSY::CrossSectionDB::rel_uncertainty	(	int	id,
		int	proc = 0 ) const

inline

Definition at line 105 of file SUSYCrossSection.h.

105{ return process(id, proc).relunc(); }

◆ setExtended()

void SUSY::CrossSectionDB::setExtended ( bool isExtended = true )

inline

Definition at line 85 of file SUSYCrossSection.h.

85{ m_extended = isExtended; };

◆ setUsePMGTool()

void SUSY::CrossSectionDB::setUsePMGTool ( bool usePMGTool = true )

inline

Definition at line 86 of file SUSYCrossSection.h.

86{ m_usePMGTool = usePMGTool; };

◆ sumweight() [1/2]

float SUSY::CrossSectionDB::sumweight	(	int	id,
		int	pdgId1,
		int	pdgId2 ) const

inline

Definition at line 108 of file SUSYCrossSection.h.

108{ return sumweight(id, finalState(pdgId1, pdgId2)); }

◆ sumweight() [2/2]

float SUSY::CrossSectionDB::sumweight	(	int	id,
		int	proc = 0 ) const

inline

Definition at line 107 of file SUSYCrossSection.h.

107{ return process(id, proc).sumweight(); }

◆ xsectTimesEff() [1/2]

float SUSY::CrossSectionDB::xsectTimesEff	(	int	id,
		int	pdgId1,
		int	pdgId2 ) const

inline

Definition at line 98 of file SUSYCrossSection.h.

98{ return xsectTimesEff(id, finalState(pdgId1, pdgId2)); }

SUSY::CrossSectionDB::xsectTimesEff

float xsectTimesEff(int id, int proc=0) const

Definition SUSYCrossSection.h:94

◆ xsectTimesEff() [2/2]

float SUSY::CrossSectionDB::xsectTimesEff	(	int	id,
		int	proc = 0 ) const

inline

Definition at line 94 of file SUSYCrossSection.h.

                                                  {
    Process p = process(id, proc);
    return p.xsect() * p.kfactor() * p.efficiency();
  }

Member Data Documentation

◆ m_cache

xsDB_t SUSY::CrossSectionDB::m_cache

private

Definition at line 129 of file SUSYCrossSection.h.

◆ m_extended

bool SUSY::CrossSectionDB::m_extended

private

Definition at line 116 of file SUSYCrossSection.h.

◆ m_pmgxs

asg::AnaToolHandle<PMGTools::IPMGCrossSectionTool> SUSY::CrossSectionDB::m_pmgxs

private

Definition at line 125 of file SUSYCrossSection.h.

◆ m_usePMGTool

bool SUSY::CrossSectionDB::m_usePMGTool

private

Definition at line 117 of file SUSYCrossSection.h.

◆ m_xsectDB

xsDB_t SUSY::CrossSectionDB::m_xsectDB

private

Definition at line 128 of file SUSYCrossSection.h.

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

Classes

Public Types

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ iterator

◆ xsDB_t

Constructor & Destructor Documentation

◆ CrossSectionDB()

Member Function Documentation

◆ begin()

◆ efficiency() [1/2]

◆ efficiency() [2/2]

◆ end()

◆ extend()

◆ kfactor() [1/2]

◆ kfactor() [2/2]

◆ loadFile()

◆ my_find()

◆ name()

◆ process() [1/2]

◆ process() [2/2]

◆ rawxsect() [1/2]

◆ rawxsect() [2/2]

◆ rel_uncertainty() [1/2]

◆ rel_uncertainty() [2/2]

◆ setExtended()

◆ setUsePMGTool()

◆ sumweight() [1/2]

◆ sumweight() [2/2]

◆ xsectTimesEff() [1/2]

◆ xsectTimesEff() [2/2]

Member Data Documentation

◆ m_cache

◆ m_extended

◆ m_pmgxs

◆ m_usePMGTool

◆ m_xsectDB