SUSY::CrossSectionDB Class Reference

#include <SUSYCrossSection.h>

Collaboration diagram for SUSY::CrossSectionDB:

Classes
class	Process
class	Key

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.

{ return m_xsectDB.begin(); }

efficiency() [1/2]

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

inline

Definition at line 104 of file SUSYCrossSection.h.

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

efficiency() [2/2]

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

inline

Definition at line 103 of file SUSYCrossSection.h.

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

end()

iterator SUSY::CrossSectionDB::end ( ) const

inline

Definition at line 121 of file SUSYCrossSection.h.

{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.

{ 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.

{ 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.

{ 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.

{ 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.

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

rawxsect() [2/2]

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

inline

Definition at line 99 of file SUSYCrossSection.h.

{ 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.

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

rel_uncertainty() [2/2]

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

inline

Definition at line 105 of file SUSYCrossSection.h.

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

setExtended()

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

inline

Definition at line 85 of file SUSYCrossSection.h.

{ m_extended = isExtended; };

setUsePMGTool()

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

inline

Definition at line 86 of file SUSYCrossSection.h.

{ 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.

{ 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.

{ 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.

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

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:

• SUSYCrossSection.h
• SUSYCrossSection.cxx