Analysis::CalibrationDataInterfaceROOT Class Reference

This tool provides an interface to flavour tagging performance estimates. More...

#include <CalibrationDataInterfaceROOT.h>

Inheritance diagram for Analysis::CalibrationDataInterfaceROOT:

Classes
class	HadronisationReferenceHelper

Public Types
enum	variableType { kEta , kAbsEta , kPt }
	known variable types that can be used as function arguments More...

Public Member Functions
	CalibrationDataInterfaceROOT (const std::string &taggerName, const std::string &configname="BTagCalibration.env", const std::string &pathname="")
	main constructor for "stand-alone" use (with information fed in from a .env configuration file read by TEnv)
	CalibrationDataInterfaceROOT (const std::string &taggerName, const char *fileSF, const char *fileEff, const std::vector< std::string > &jetAliases, const std::map< std::string, std::string > &SFNames, const std::map< std::string, std::vector< std::string > > &EffNames, const std::map< std::string, std::vector< std::string > > &excludeFromEV, const std::map< std::string, Analysis::EVReductionStrategy > &EVReductions, bool useEV=true, Uncertainty strat=SFEigen, bool useMCMCSF=true, bool useTopologyRescaling=false, bool useRecommendedEVExclusions=false, bool verbose=true, std::vector< std::string > flavours={"B", "C", "Light", "T"})
	alternative constructor passing configuration information explicitly (so that no .env file is needed)
	CalibrationDataInterfaceROOT ()
	default constructor for PROOF object retrieval
	CalibrationDataInterfaceROOT (const CalibrationDataInterfaceROOT &other)
	copy constructor
virtual	~CalibrationDataInterfaceROOT ()
	default destructor
CalibResult	getScaleFactor (const CalibrationDataVariables &variables, const std::string &label, const std::string &OP, Uncertainty unc, unsigned int numVariation=0, unsigned int mapIndex=0)
	efficiency scale factor retrieval by name.
CalibResult	getEfficiency (const CalibrationDataVariables &variables, const std::string &label, const std::string &OP, Uncertainty unc, const std::string &flavour, unsigned int numVariation=0, unsigned int mapIndex=0)
	efficiency retrieval by name
CalibResult	getInefficiencyScaleFactor (const CalibrationDataVariables &variables, const std::string &label, const std::string &OP, Uncertainty unc, unsigned int numVariation=0, unsigned int mapIndex=0)
	"MC" inefficiency scale factor retrieval by name
CalibResult	getInefficiency (const CalibrationDataVariables &variables, const std::string &label, const std::string &OP, Uncertainty unc, unsigned int numVariation=0, unsigned int mapIndex=0)
	inefficiency retrieval by name
CalibResult	getMCEfficiency (const CalibrationDataVariables &variables, const std::string &label, const std::string &OP, Uncertainty unc=None, unsigned int mapIndex=0)
	"MC" efficiency retrieval by name
CalibResult	getMCInefficiency (const CalibrationDataVariables &variables, const std::string &label, const std::string &OP, Uncertainty unc=None, unsigned int mapIndex=0)
	"MC" inefficiency retrieval by name
std::vector< std::string >	listScaleFactorUncertainties (const std::string &author, const std::string &label, const std::string &OP, bool named=false)
	retrieve the list of "uncertainties" relevant to the calibration object.
unsigned int	getNumVariations (const std::string &author, const std::string &label, const std::string &OP, Uncertainty unc)
	retrieve the number of variations relevant to the calibration object.
bool	retrieveCalibrationIndex (const std::string &label, const std::string &OP, const std::string &author, bool isSF, unsigned int &index, unsigned int mapIndex=0)
	Retrieve the index of the calibration object (container) starting from the label and operating point.
std::string	nameFromIndex (unsigned int index) const
	Retrieve the name of the calibration object (container) given its index.
CalibResult	getScaleFactor (const CalibrationDataVariables &variables, unsigned int indexSF, unsigned int indexEff, Uncertainty unc, const std::string &flavour, unsigned int numVariation=0)
	efficiency scale factor retrieval by index #2
CalibResult	getEfficiency (const CalibrationDataVariables &variables, unsigned int indexSF, unsigned int indexEff, Uncertainty unc, const std::string &flavour, unsigned int numVariation=0)
	efficiency retrieval by index
CalibResult	getInefficiencyScaleFactor (const CalibrationDataVariables &variables, unsigned int indexSF, unsigned int indexEff, Uncertainty unc, const std::string &flavour, unsigned int numVariation=0)
	"MC" inefficiency scale factor retrieval by index
CalibResult	getInefficiency (const CalibrationDataVariables &variables, unsigned int indexSF, unsigned int indexEff, Uncertainty unc, const std::string &flavour, unsigned int numVariation=0)
	inefficiency retrieval by index
CalibResult	getMCEfficiency (const CalibrationDataVariables &variables, unsigned int index, Uncertainty unc=None)
	"MC" efficiency retrieval by index
CalibResult	getMCInefficiency (const CalibrationDataVariables &variables, unsigned int index, Uncertainty unc=None)
	"MC" inefficiency retrieval by index
double	getMCMCScaleFactor (const CalibrationDataVariables &variables, unsigned indexSF, unsigned int indexEff) const
	MC/MC scale factor retrieval.
std::vector< std::string >	listScaleFactorUncertainties (unsigned int index, const std::string &flavour, bool named=false)
	retrieve the list of "uncertainties" relevant to the calibration object.
unsigned int	getNumVariations (unsigned int index, Uncertainty unc, const std::string &flavour)
	retrieve the number of variations relevant to the calibration object.
std::string	fullName (const std::string &author, const std::string &OP, const std::string &label, bool isSF, unsigned mapIndex=0) const
	@ brief construct the full object pathname from its individual components
CalibrationStatus	getScaleFactor (const CalibrationDataVariables &variables, unsigned int indexSF, unsigned int indexEff, Uncertainty unc, unsigned int numVariation, CalibResult &result, const std::string &flavour)
	efficiency scale factor retrieval by index #3
CalibrationStatus	getEfficiency (const CalibrationDataVariables &variables, unsigned int indexSF, unsigned int indexEff, Uncertainty unc, unsigned int numVariation, CalibResult &result, const std::string &flavour)
	efficiency retrieval by index
CalibrationStatus	getInefficiencyScaleFactor (const CalibrationDataVariables &variables, unsigned int indexSF, unsigned int indexEff, Uncertainty unc, unsigned int numVariation, CalibResult &result, const std::string &flavour)
	"MC" inefficiency scale factor retrieval by index
CalibrationStatus	getInefficiency (const CalibrationDataVariables &variables, unsigned int indexSF, unsigned int indexEff, Uncertainty unc, unsigned int numVariation, CalibResult &result, const std::string &flavour)
	inefficiency retrieval by index
CalibrationStatus	getMCEfficiency (const CalibrationDataVariables &variables, unsigned int index, Uncertainty unc, CalibResult &result)
	"MC" efficiency retrieval by index
CalibResult	getWeightScaleFactor (const CalibrationDataVariables &variables, const std::string &label, Uncertainty unc, unsigned int numVariation=0, unsigned int mapIndex=0)
	efficiency scale factor retrieval by name
CalibResult	getWeightScaleFactor (const CalibrationDataVariables &variables, unsigned int indexSF, unsigned int indexEff, Uncertainty unc, unsigned int numVariation=0)
	efficiency scale factor retrieval by index
CalibrationStatus	getWeightScaleFactor (const CalibrationDataVariables &variables, unsigned int indexSF, unsigned int indexEff, Uncertainty unc, unsigned int numVariation, CalibResult &result)
	efficiency scale factor retrieval by index, with different signature
CalibrationStatus	runEigenVectorRecomposition (const std::string &author, const std::string &label, const std::string &OP, unsigned int mapindex=0)
	run EigenVector Recomposition method
CalibrationStatus	runEigenVectorRecomposition (const std::string &label, unsigned int mapindex=0)
std::map< std::string, std::map< std::string, float > >	getEigenVectorRecompositionCoefficientMap ()
	Get Eigenvector recomposition map after running runEigenVectorRecomposition()
const TObject *	getMCEfficiencyObject (const std::string &author, const std::string &label, const std::string &OP, unsigned int mapIndex=0)
	retrieve the MC efficiency (central values) object for the given flavour label and operating point.
const TH1 *	getBinnedScaleFactors (const std::string &author, const std::string &label, const std::string &OP)
	retrieve the binned calibration object for the given flavour label and operating point.
const TH1 *	getShiftedScaleFactors (const std::string &author, const std::string &label, const std::string &OP, const std::string &unc, double sigmas)
	retrieve the binned calibration object for the given flavour label and operating point, with the result shifted by the given number of standard deviations for the given systematic uncertainty.
TMatrixDSym	getScaleFactorCovarianceMatrix (const std::string &author, const std::string &label, const std::string &OP, const std::string &unc="all")
	retrieve the named covariance matrix element corresponding to the binned calibration object.
void	initialize (const std::string &jetauthor, const std::string &OP, Uncertainty unc)
	initialization for PROOF usage
CalibrationDataContainer *	retrieveContainer (const std::string &label, const std::string &OP, const std::string &author, const std::string &cntname, bool isSF, bool doPrint=true)
	utility function taking care of object retrieval
const std::string &	EffCalibrationName (const std::string &flavour, unsigned int mapIndex=0) const
	Main interface methods accessing the flavour tagging performance information.
void	setEffCalibrationNames (const std::map< std::string, std::vector< std::string > > &names)
const std::string &	SFCalibrationName (const std::string &flavour) const
void	setSFCalibrationNames (const std::map< std::string, std::string > &names)

Protected Member Functions
std::string	getContainername (const std::string &flavour, bool SF, unsigned int mapIndex=0) const
	auxiliary function for string concatenation
std::string	getBasename (const std::string &name) const
	auxiliary function for retrieval of name within the directory
double	combinedUncertainty (double stat, const std::pair< double, double > &syst) const
	utility function for combination of statistical and (a priori asymmetric) systematic uncertainty.

Protected Attributes
std::string	m_taggerName
	tagging algorithm name

Private Member Functions
std::string	getAlias (const std::string &author) const
	associated alias retrieval method
bool	checkAbsEta (const CalibrationDataVariables &variables, unsigned int index)
void	increaseCounter (unsigned int index, OutOfBoundsType oob=Main)
void	checkWeightScaleFactors (unsigned int indexSF, unsigned int indexEff)

Private Attributes
TFile *	m_fileEff {}
	pointer to the TFile object providing access to the calibrations
TFile *	m_fileSF {}
	Do not attempt to persistify (PROOF)
std::map< std::string, std::string >	m_aliases
	Do not attempt to persistify (PROOF)
std::vector< CalibrationDataContainer * >	m_objects
	cache the objects themselves (so that the user will not have to delete them after each call etc.).
std::map< std::string, unsigned int >	m_objectIndices
std::string	m_filenameSF
	in addition, store also the filenames themselves (needed for the copy constructor)
std::string	m_filenameEff
std::vector< std::string >	m_flavours
std::map< const CalibrationDataContainer *, std::shared_ptr< CalibrationDataEigenVariations > >	m_eigenVariationsMap
	store the eigenvector class and associate to its CalibrationDataContainer
bool	m_runEigenVectorMethod {}
	decide whether to run the eigenvector method or not
Uncertainty	m_EVStrategy {}
std::map< std::string, Analysis::EVReductionStrategy >	m_EVReductions
	Eigenvector reduction strategy (per flavour)
std::map< std::string, std::vector< std::string > >	m_excludeFromCovMatrix
	store the uncertainties which should be excluded from building the full covariance matrix
bool	m_useRecommendedEVExclusions {}
	if true, exclude pre-recommended lists of uncertainties from the covariance matrix building, in addition to the above user specified lists
bool	m_verbose {}
	if true, allow also for some informational (and not only error/warning) messages
bool	m_useMCMCSF {}
	specify whether or not to use MC/MC (hadronisation) scale factors (the fact that this is steerable is intended to be temporary only)
bool	m_useTopologyRescaling {}
	specify whether or not to use MC/MC (topology) scale factors (also this steering option may be removed)
std::map< std::string, HadronisationReferenceHelper * >	m_refMap
	the following maps (one for each directory) specify the name of the container serving as the 'hadronisation' reference for each object
std::vector< int >	m_hadronisationReference
	store the 'hadronisation' reference for each object (-1 means no reference found)
std::map< std::string, std::map< std::string, float > >	m_coefficientMap
double	m_maxAbsEta {}
	|eta| bounds and strategy for dealing with out-of-bounds conditions
OutOfBoundsStrategy	m_absEtaStrategy {}
OutOfBoundsStrategy	m_otherStrategy {}
std::vector< unsigned int >	m_etaCounters
	counters for flagging out-of-bound cases
std::vector< unsigned int >	m_mainCounters
std::vector< unsigned int >	m_extrapolatedCounters
std::vector< std::pair< unsigned int, unsigned int > >	m_checkedWeightScaleFactors
double	m_maxTagWeight {}
std::map< std::string, std::vector< std::string > >	m_calibrationEffNames
	this simply collects the per-flavour properties.
std::map< std::string, std::string >	m_calibrationSFNames

Detailed Description

This tool provides an interface to flavour tagging performance estimates.

A separate instance should be used for each different tagging algorithm. For each instance, all appropriate jet collections and tagger operating points need to be specified.

The model:

• b-jets: data-MC scale factor (factorised 2D function of eta, pt) MC reference the product is the data efficiency; alternatively, the scale factor may be used
• c-jets: as for b-jets, but with a different MC reference
• light-flavour jets: data-MC scale factor (factorised 2D function of eta, pt) MC reference Besides the results, it is also possible to retrieve associated uncertainties. This need not be configured, and a choice as to the uncertainty component can be made on a case-by-case basis.

The idea is to use the same physical ROOT file that is also accessed through COOL, but to do so in a stand-alone fashion, so that there is no COOL or Athena dependence. Apart from this, the same infrastructure and limitations as with COOL access (

Definition at line 87 of file CalibrationDataInterfaceROOT.h.

Member Enumeration Documentation

variableType

enum Analysis::CalibrationDataInterfaceBase::variableType

inherited

known variable types that can be used as function arguments

Enumerator
kEta
kAbsEta
kPt

Definition at line 67 of file CalibrationDataInterfaceBase.h.

{ kEta, kAbsEta, kPt };

Constructor & Destructor Documentation

CalibrationDataInterfaceROOT() [1/4]

Analysis::CalibrationDataInterfaceROOT::CalibrationDataInterfaceROOT	(	const std::string &	taggerName,
		const std::string &	configname = "BTagCalibration.env",
		const std::string &	pathname = "" )

main constructor for "stand-alone" use (with information fed in from a .env configuration file read by TEnv)

CalibrationDataInterfaceROOT() [2/4]

Analysis::CalibrationDataInterfaceROOT::CalibrationDataInterfaceROOT	(	const std::string &	taggerName,
		const char *	fileSF,
		const char *	fileEff,
		const std::vector< std::string > &	jetAliases,
		const std::map< std::string, std::string > &	SFNames,
		const std::map< std::string, std::vector< std::string > > &	EffNames,
		const std::map< std::string, std::vector< std::string > > &	excludeFromEV,
		const std::map< std::string, Analysis::EVReductionStrategy > &	EVReductions,
		bool	useEV = true,
		Uncertainty	strat = SFEigen,
		bool	useMCMCSF = true,
		bool	useTopologyRescaling = false,
		bool	useRecommendedEVExclusions = false,
		bool	verbose = true,
		std::vector< std::string >	flavours = {"B", "C", "Light", "T"} )

alternative constructor passing configuration information explicitly (so that no .env file is needed)

Definition at line 481 of file CalibrationDataInterfaceROOT.cxx.

                                                    :
  m_filenameSF(fileSF), m_filenameEff(""), m_flavours(std::move(flavours)),
  m_runEigenVectorMethod(useEV), m_EVStrategy(strat), m_EVReductions(EVReductions),
  m_useRecommendedEVExclusions(useRecommendedEEVExclusions), m_verbose(verbose),
  m_useMCMCSF(useMCMCSF), m_useTopologyRescaling(useTopologyRescaling),
  m_maxAbsEta(2.5), m_absEtaStrategy(GiveUp),
  m_otherStrategy(Flag), m_maxTagWeight(10.0)
{
  // Normal constructor avoiding the need for a .env file.
  //
  //     taggerName:     this should correspond to the tagger name as used in the calibration ROOT file
  //     fileSF:         full path of the calibration ROOT file containing the calibration scale factors
  //     fileEff:        optional full path name of a ROOT file containing additional MC efficiency maps
  //                     (use a null pointer to disable the use of such additional file)
  //     flavours;       This should correspond to the list of flavour labels that's used by a tagger, and
  //                     which corresponds to the labels used in internal maps
  //     jetAliases:     this can be used to convert jet collection names to the corresponding names in the
  //                     calibration ROOT file (this may be useful as e.g. the collection names in the
  //                     calibration ROOT file have the JVF criterion attached as a suffix).
  //                     Each alias is specified as
  //                                        nameOrig->nameTarget,
  //                     where nameOrig and nameTarget are the names of the input jet collection and the
  //                     jet collection name as used in the calibration ROOT file, respectively.
  //     SFNames:        map specifying for each of the calibration flavours ("B", "C", "T", "Light") the
  //                     name of the scale factor calibration object
  //     EffNames:       map specifying for each of the calibration flavours ("B", "C", "T", "Light") the
  //                     names of the possibly relevant efficiency calibration objects
  //     excludeFromEV:  map specifying for each of the calibration flavours ("B", "C", "T", "Light") the
  //                     systematic uncertainties to be excluded from the Eigenvector variation treatment
  //                     (this is used only if Eigenvector variations are used to begin with)
  //     EVReductions:   Eigenvector variation reduction strategies for "B", "C", "Light" jets (again,
  //                     this is only relevant if Eigenvector variations are used to begin with)
  //     useEV:          switch specifying if Eigenvector variations will be used or not
  //     useMCMCSF:      switch specifying if generator-dependent scale factors are to be applied or not
 
  // Note: at present, the means to change the strategies and maximum values initialized above do not exist
  // when using this constructor
 
  if (m_verbose) {
    cout << "=== CalibrationDataInterfaceROOT::CalibrationDataInterfaceROOT ===" << endl;
    cout << " taggerName           : " << taggerName.c_str() << endl;
    cout << " Systematic strategy : ";
    if (m_EVStrategy == Analysis::Uncertainty::SFEigen){
      cout << "SFEigen" << endl;
    } else if (m_EVStrategy == Analysis::Uncertainty::SFGlobalEigen){
      cout << "SFGlobalEigen" << endl;
    } else {
      cout << " Other" << endl;
    }
    if (fileEff) cout << " Efficiency file name : " << fileEff << endl;
    cout << " SF         file name : " << fileSF << endl
     << endl;
  }
 
  m_taggerName = taggerName;
 
  m_fileSF = TFile::Open(fileSF, "READ");
  if (fileEff && strcmp(fileSF, fileEff) != 0) {
    m_filenameEff = string(fileEff);
    m_fileEff = TFile::Open(fileEff, "READ");
  } else
    m_fileEff = m_fileSF;
 
  if (m_verbose) {
    TObjString* s;
    m_fileSF->GetObject("VersionInfo/BuildNumber", s);
    if (s) cout << " CDI file build number: " << s->GetName() << endl;
    cout << endl;
  }
  
  for (unsigned int i = 0; i < jetAliases.size(); ++i) {
    // Each alias specification uses an arrow ("->"). Forget about entries
    // not properly following this specification.
    string::size_type arrow = jetAliases[i].find("->");
    if (arrow == string::npos) continue;
    m_aliases[jetAliases[i].substr(0,arrow)] = jetAliases[i].substr(arrow+2);
  }
 
  setEffCalibrationNames(EffNames);
  setSFCalibrationNames(SFNames);
 
  if (m_runEigenVectorMethod) { 
    // if we want to run EV, then decide which one
    // The following should hold for both eigenvector decomposition methods (SFEigen and SFGlobalEigen)
    // The global one simply adapts itself to using the m_excludeFromCovMatrix to perform the same task
    
    m_excludeFromCovMatrix = excludeFromEV;
    unsigned int n_excluded = 0;
    for (auto const& flavour : m_flavours) {
      n_excluded += m_excludeFromCovMatrix[flavour].size();
    }
    if (m_verbose) {
      cout << " List of uncertainties to exclude:";
      if (n_excluded == 0) cout << " none";
      for (auto const& flavour : m_flavours) {
        if (m_excludeFromCovMatrix[flavour].size() > 0) {
          cout << "\n\t" << flavour << ":\t";
          for (unsigned int i = 0; i < m_excludeFromCovMatrix[flavour].size(); ++i) {
            cout << m_excludeFromCovMatrix[flavour].at(i);
            if (i+1 == m_excludeFromCovMatrix[flavour].size()) cout << ";  ";
          }
          cout << endl;
        }
            }
      cout << endl;
    }
 
  }
  
  if (m_verbose) cout << "======= end of CalibrationDataInterfaceROOT instantiation ========" << endl;
}

CalibrationDataInterfaceROOT() [3/4]

Analysis::CalibrationDataInterfaceROOT::CalibrationDataInterfaceROOT

(

)

default constructor for PROOF object retrieval

Definition at line 603 of file CalibrationDataInterfaceROOT.cxx.

{ 
  // Default constructor for PROOF purposes
 
  m_fileEff=0; 
  m_fileSF=0; 
} 

CalibrationDataInterfaceROOT() [4/4]

Analysis::CalibrationDataInterfaceROOT::CalibrationDataInterfaceROOT

(

const CalibrationDataInterfaceROOT &

other

)

copy constructor

Definition at line 612 of file CalibrationDataInterfaceROOT.cxx.

                                                                                                                    :
  Analysis::CalibrationDataInterfaceBase(other), m_aliases(other.m_aliases), m_objects(), m_objectIndices(),
  m_filenameSF(other.m_filenameSF), m_filenameEff(other.m_filenameEff),
  m_eigenVariationsMap(), m_runEigenVectorMethod(other.m_runEigenVectorMethod), m_EVStrategy(other.m_EVStrategy),
  m_excludeFromCovMatrix(other.m_excludeFromCovMatrix), m_useMCMCSF(other.m_useMCMCSF), m_useTopologyRescaling(other.m_useTopologyRescaling),
  m_refMap(), m_hadronisationReference(),
  m_maxAbsEta(other.m_maxAbsEta), m_absEtaStrategy(other.m_absEtaStrategy), m_otherStrategy(other.m_otherStrategy),
  m_etaCounters(other.m_etaCounters), m_mainCounters(other.m_mainCounters), m_extrapolatedCounters(other.m_extrapolatedCounters),
  m_checkedWeightScaleFactors(other.m_checkedWeightScaleFactors), m_maxTagWeight(other.m_maxTagWeight)
{
  // Copy constructor. Note that the "cacheable" items aren't copied (they will be re-created if needed)
 
  // The TFile objects cannot be copied. Therefore, create duplicate objects starting from the filenames
  m_fileSF = TFile::Open(m_filenameSF.c_str(), "READ");
  if (m_filenameEff == m_filenameSF)
    m_fileSF = m_fileEff;
  else
    m_fileEff = TFile::Open(m_filenameEff.c_str(), "READ");
}

~CalibrationDataInterfaceROOT()

Analysis::CalibrationDataInterfaceROOT::~CalibrationDataInterfaceROOT ( )

virtual

default destructor

Definition at line 633 of file CalibrationDataInterfaceROOT.cxx.

{
  // Destructor
  if ((m_fileEff!=0) && (m_fileSF!=0)) { 
    if (m_fileEff == m_fileSF) { 
      m_fileEff->Close(); 
      delete m_fileEff; m_fileEff = 0; 
    } else { 
      m_fileEff->Close(); 
      m_fileSF->Close(); 
      delete m_fileEff; m_fileEff = 0; 
      delete m_fileSF;  m_fileSF = 0; 
    } 
  }
  // delete also the stored objects (these are owned by us)
  for (std::vector<CalibrationDataContainer*>::iterator it = m_objects.begin(); it != m_objects.end(); ++it) {
    if (*it) {
      delete *it; *it = 0;
    }
  }
 
  for (std::map<std::string, HadronisationReferenceHelper*>::iterator it = m_refMap.begin();
       it != m_refMap.end(); ++it) {
    if(it->second)
      { delete it->second; it->second=nullptr; }
  }
 
  // Print summary output on out-of-bounds issues
  if (m_absEtaStrategy == Flag && m_verbose) {
    bool found = false;
    cout << "\t\tCalibrationDataInterfaceROOT |eta| out-of-bounds summary:" << endl;
    for (unsigned int index = 0; index < m_mainCounters.size(); ++index)
      if (m_etaCounters[index] > 0) {
        found = true;
        cout << "\t\t\t" << nameFromIndex(index) << ": " << m_etaCounters[index] << endl;
      }
    if (!found) cout << "\t\t\tNo issues found" << endl;
  }
  if (m_otherStrategy == Flag && m_verbose) {
    bool found = false;
    cout << "\t\tCalibrationDataInterfaceROOT object out-of-bounds summary:" << endl;
    for (unsigned int index = 0; index < m_mainCounters.size(); ++index)
      if (m_mainCounters[index] + m_extrapolatedCounters[index] > 0) {
        found = true;
        cout << "\t\t\t" << nameFromIndex(index)
            << " general: "       << m_mainCounters[index]
            << ", extrapolated: " << m_extrapolatedCounters[index]
              << endl;
      }
    if (!found) cout << "\t\t\tNo issues found" << endl;
  }
}

Member Function Documentation

checkAbsEta()

bool Analysis::CalibrationDataInterfaceROOT::checkAbsEta ( const CalibrationDataVariables & variables, unsigned int index )

nodiscardprivate

Definition at line 1916 of file CalibrationDataInterfaceROOT.cxx.

{
  // Check whether the jet eta value is outside the range of validity, subject to the strategy
  // specified in the configuration file.
  bool pass = true;
  if (m_absEtaStrategy == Ignore) return pass; 
 
  switch (m_absEtaStrategy) {
  case GiveUp:
    if (std::fabs(variables.jetEta) > m_maxAbsEta) {
      pass = false;  
    }
    break;
  case Flag:
  default:
    if (std::fabs(variables.jetEta) > m_maxAbsEta) {
      increaseCounter(index, Eta);
    }
 
  }
  return pass;
}

checkWeightScaleFactors()

void Analysis::CalibrationDataInterfaceROOT::checkWeightScaleFactors ( unsigned int indexSF, unsigned int indexEff )

private

Definition at line 1763 of file CalibrationDataInterfaceROOT.cxx.

{
  // Check the tag weight scale factors that would result from the combination of
  // the provided scale factor and MC tag weight objects.
  // The way this is done is by determining the binning that would apply to the
  // combination of the two individual inputs, and then by explicitly computing
  // the scale factors in each of these resulting bins.
 
  std::vector<std::pair<unsigned int, unsigned int> >::const_iterator it = std::find(m_checkedWeightScaleFactors.begin(), m_checkedWeightScaleFactors.end(), std::make_pair(indexSF, indexEff));
  if (it != m_checkedWeightScaleFactors.end()) return;
 
 
  // Assume that only histogram containers are involved here (this should be the case
  // as at least a strict tag weight binning should be applied).
  CalibrationDataHistogramContainer* container = dynamic_cast<CalibrationDataHistogramContainer*>(m_objects[indexSF]);
  if (! container) {
    cerr << "CalibrationDataInterfaceROOT::checkWeightScaleFactors: error: container for object " << nameFromIndex(indexSF) << " not found!" << endl;
    return;
  } else if (! container->GetValue("MCreference")) {
    cerr << "CalibrationDataInterfaceROOT::checkWeightScaleFactors: error: no MCreference histogram for object " << nameFromIndex(indexSF) << "!" << endl;
    return;
  }
  CalibrationDataHistogramContainer* effContainer = dynamic_cast<CalibrationDataHistogramContainer*>(m_objects[indexEff]);
  if (! effContainer) {
    cerr << "CalibrationDataInterfaceROOT::checkWeightScaleFactors: error: container for object " << nameFromIndex(indexEff) << " not found!" << endl;
    return;
  }
 
  // Retrieve the variable types and corresponding bin boundaries
  std::vector<unsigned int> vars    = container->getVariableTypes();
  std::vector<unsigned int> effVars = effContainer->getVariableTypes();
  // Retrieve the corresponding bin boundaries
  std::map<unsigned int, std::vector<double> > boundaries, effBoundaries, mergedBoundaries;
  for (unsigned int t = 0; t < vars.size(); ++t)
    boundaries[vars[t]] = container->getBinBoundaries(vars[t]);
  for (unsigned int t = 0; t < effVars.size(); ++t)
    effBoundaries[effVars[t]] = effContainer->getBinBoundaries(effVars[t]);
 
  // Special case: handle |eta| versus eta differences, by transforming to the latter
  if (boundaries.find(CalibrationDataContainer::kEta) == boundaries.end() && boundaries.find(CalibrationDataContainer::kAbsEta) != boundaries.end()) {
    boundaries[CalibrationDataContainer::kEta] = boundaries[CalibrationDataContainer::kAbsEta];
    boundaries.erase(CalibrationDataContainer::kAbsEta);
  }
  if (effBoundaries.find(CalibrationDataContainer::kEta) == effBoundaries.end() && effBoundaries.find(CalibrationDataContainer::kAbsEta) != effBoundaries.end()) {
    effBoundaries[CalibrationDataContainer::kEta] = effBoundaries[CalibrationDataContainer::kAbsEta];
    effBoundaries.erase(CalibrationDataContainer::kAbsEta);
  }
  if (boundaries.find(CalibrationDataContainer::kEta) != boundaries.end() && effBoundaries.find(CalibrationDataContainer::kEta) != effBoundaries.end()) {
    std::vector<double>& v = boundaries[CalibrationDataContainer::kEta];
    std::vector<double>& vEff = effBoundaries[CalibrationDataContainer::kEta];
    if (v[0] < 0 && vEff[0] >= 0) {
      // in this case, supplement the positive entries in vEff with their negative analogues
      std::vector<double> vtmp(vEff);
      for (std::vector<double>::iterator it = vtmp.begin(); it != vtmp.end(); ++it)
        if (*it > 0) vEff.insert(vEff.begin(), -(*it));
    } else if (v[0] >= 0 && vEff[0] < 0) {
      // in this case, supplement the positive entries in v with their negative analogues
      std::vector<double> vtmp(v);
      for (std::vector<double>::iterator it = vtmp.begin(); it != vtmp.end(); ++it)
        if (*it > 0) v.insert(v.begin(), -(*it));
    }
  }
 
  // Now that the individual sets of boundaries have been determined, merge these
  for (unsigned int t = 0; t < vars.size(); ++t) {
    if (effBoundaries.find(vars[t]) == effBoundaries.end())
      // Variables not present in the efficiency object can go in unmodified
      mergedBoundaries[vars[t]] = boundaries[vars[t]];
    else {
      // Merge the boundaries for variables existing in both objects.
      // Take the MC array as a starting point, as it's likely to be the longest.
      mergedBoundaries[vars[t]] = effBoundaries[vars[t]];
      
      for (std::vector<double>::iterator it = boundaries[vars[t]].begin(); it != boundaries[vars[t]].end(); ++it) {
        std::vector<double>::iterator itcmp = mergedBoundaries[vars[t]].begin();
        // Iterate until we've found a value in the target array equal to
        // or larger than the given element
        while (itcmp != mergedBoundaries[vars[t]].end() &&
               (! CalibrationDataContainer::isNearlyEqual(*itcmp, *it)) &&
               *itcmp < *it) ++itcmp;
        // Nothing needs to be done if the values are "nearly identical"
        // (or if we don't find such an element).
        if (itcmp == mergedBoundaries[vars[t]].end() || CalibrationDataContainer::isNearlyEqual(*itcmp, *it)) continue;
        // Otherwise insert the given element (this can mean adding to the end)
        mergedBoundaries[vars[t]].insert(itcmp, *it);
      }
    }
  }
  // Variables not present in the scale factor object still need to go in
  for (unsigned int t = 0; t < effVars.size(); ++t)
    if (boundaries.find(effVars[t]) == boundaries.end()) 
      mergedBoundaries[effVars[t]] = effBoundaries[effVars[t]];
  
  // Carry out a rudimentary cross-check of the tag weight bin
  // (the binning used for the scale factor and MC objects should be identical).
  if (boundaries.find(CalibrationDataContainer::kTagWeight) == boundaries.end()) {
    cerr << "CalibrationDataInterfaceROOT::checkWeightScaleFactors: " << "no tag weight axis found for object " << nameFromIndex(indexSF) << endl;
  } else if (effBoundaries.find(CalibrationDataContainer::kTagWeight) == effBoundaries.end()) {
    cerr << "CalibrationDataInterfaceROOT::checkWeightScaleFactors: " << "no tag weight axis found for object " << nameFromIndex(indexEff) << endl;
  } else if (boundaries[CalibrationDataContainer::kTagWeight].size() != effBoundaries[CalibrationDataContainer::kTagWeight].size()) {
    cerr << "CalibrationDataInterfaceROOT::checkWeightScaleFactors: " << "different tag weight binning for objects " << nameFromIndex(indexSF) << " (";
    std::vector<double>& v = boundaries[CalibrationDataContainer::kTagWeight];
    for (unsigned int ib = 0; ib < v.size()-1; ++ib) cerr << v[ib] << ",";
    cerr << v[v.size()-1] << ") and " << nameFromIndex(indexEff) << " (";
    v = effBoundaries[CalibrationDataContainer::kTagWeight];
    for (unsigned int ib = 0; ib < v.size()-1; ++ib) cerr << v[ib] << ",";
    cerr << v[v.size()-1] << ") do not match!" << endl;
  } else {
    // Make sure that (possibly) dummy vectors exist for _all_ known variables
    // (this is a mere technicality allowing to loop over all variables explicitly).
    mergedBoundaries.try_emplace(CalibrationDataContainer::kPt, std::vector<double>{20.,300.});
    mergedBoundaries.try_emplace(CalibrationDataContainer::kEta, std::vector<double>{-2.5, 2.5});
   
    // Finally, carry out the cross-check that all this is about: recompute the scale factor
    // in each pseudo-bin
    if (m_verbose){
      cout << "CalibrationDataInterfaceROOT::checkWeightScaleFactors: cross-checking scale factors for objects " << nameFromIndex(indexSF) << " and " << nameFromIndex(indexEff) << "\n" << std::setfill('-') << std::setw(100) << "-" << endl;
      cout << std::setfill(' ');
    }
    CalibrationDataVariables x;
    std::vector<double>& vPt = mergedBoundaries[CalibrationDataContainer::kPt], vEta = mergedBoundaries[CalibrationDataContainer::kEta], vTagWeight = mergedBoundaries[CalibrationDataContainer::kTagWeight];
    for (unsigned int ipt = 0; ipt < vPt.size()-1; ++ipt) {
      x.jetPt = (vPt[ipt] + vPt[ipt+1]) * 500.; // account for MeV -> GeV conversion
      for (unsigned int ieta = 0; ieta < vEta.size()-1; ++ieta) {
          x.jetEta = (vEta[ieta] + vEta[ieta+1]) / 2.;
          for (unsigned int iwt = 0; iwt < vTagWeight.size()-1; ++iwt) {
            x.jetTagWeight = (vTagWeight[iwt] + vTagWeight[iwt+1]) / 2.;
          // Retrieve the central scale factor value and the old and new MC tag weight fractions
          double value;
          container->getResult(x, value);
          Analysis::UncertaintyResult uncertaintyResult(0,0);
          container->getUncertainty("MCreference", x, uncertaintyResult);
          double fracMCref = uncertaintyResult.first;
          double fracMCnew;
          effContainer->getResult(x, fracMCnew);
          // Compute the new scale factor value
          if (!(fracMCnew > 0.)) {
              cout << "\tfor (pt=" << x.jetPt << ",eta=" << x.jetEta << ",tagweight=" << x.jetTagWeight << "): invalid new MC fraction: " << fracMCnew << endl;
            } else {
              double newvalue = 1.0 + (value - 1.0) * fracMCref/fracMCnew;
              if (newvalue <= 0 || newvalue > m_maxTagWeight) cout << "\tfor (pt=" << x.jetPt << ",eta=" << x.jetEta << ",tagweight=" << x.jetTagWeight << "): old (value=" << value << ",MC=" << fracMCref << "), new (value=" << newvalue << ",MC=" << fracMCnew << ")" << endl;
            }
          }
      }
    }
  }
 
  m_checkedWeightScaleFactors.push_back(std::make_pair(indexSF, indexEff));
}

combinedUncertainty()

double Analysis::CalibrationDataInterfaceBase::combinedUncertainty ( double stat, const std::pair< double, double > & syst ) const

protectedinherited

utility function for combination of statistical and (a priori asymmetric) systematic uncertainty.

NB perhaps this should be in its own

Definition at line 147 of file CalibrationDataInterfaceBase.cxx.

{
  // Return the total (combined statistical and systematic) uncertainty started from
  // its individual components. The result is symmetrised by using only the larger of
  // the (a priori asymmetric) up- and downward systematic uncertainties.
 
  // The systematic uncertainty is (a priori) asymmetric, but this interface doesn't
  // at present allow for asymmetric uncertainties.
  // Address this by taking the larger (absolute) value of the two.
  double largest = syst.first;
  if (TMath::Abs(syst.second) > TMath::Abs(largest)) largest = syst.second;
 
  return TMath::Sqrt(stat*stat + largest*largest);
}

EffCalibrationName()

const std::string & Analysis::CalibrationDataInterfaceBase::EffCalibrationName ( const std::string & flavour, unsigned int mapIndex = 0 ) const

inherited

Main interface methods accessing the flavour tagging performance information.

Note that for both of the following, the label is assumed to adhere to the TruthInfo conventions (see package PhysicsAnalysis/JetTagging/JetTagInfo).

Definition at line 47 of file CalibrationDataInterfaceBase.cxx.

{
  // Return the MC efficiency name for the given flavour.
  // Note that no check is performed on the validity of the flavour.
  
  try {
    return m_calibrationEffNames.at(flavour)[mapIndex];
  }
  catch (const std::out_of_range& e) {
    std::cerr << "EffCalibrationName: flavour '" << flavour << "' is not known." << std::endl;
    throw e;
  }
}

fullName()

string Analysis::CalibrationDataInterfaceROOT::fullName	(	const std::string &	author,
		const std::string &	OP,
		const std::string &	label,
		bool	isSF,
		unsigned	mapIndex = 0 ) const

@ brief construct the full object pathname from its individual components

Definition at line 2718 of file CalibrationDataInterfaceROOT.cxx.

{
  // Construct the full calibration object's pathname within the calibration ROOT file.
  //
  //    author:   jet collection name
  //    OP:       tagger working point
  //    label:    jet flavour label
  //    isSF:     set to true (false) for scale factors (MC efficiencies)
  //    mapIndex: index in the list of MC efficiency calibration objects
 
  string flavour = (label == "N/A") ? "Light" : label;
  string full(m_taggerName + "/" + getAlias(author) + "/" + OP + "/" + flavour + "/");
  full += getContainername(flavour, isSF, mapIndex);
  // full += getAlias(author); full += "/";
  // string name = (isSF) ?
  //   getBasename(OP, label, "_SF", true) :
  //   getBasename(OP, label, "_Eff", false, mapIndex);
  // full += name;
  return full;
}

getAlias()

string Analysis::CalibrationDataInterfaceROOT::getAlias ( const std::string & author ) const

private

associated alias retrieval method

Definition at line 2707 of file CalibrationDataInterfaceROOT.cxx.

{
  // Return the alias for the given jet collection name, if an alias exists.
  // If this is not the case, the return value will simply equal the input jet collection name.
 
  std::map<string,string>::const_iterator it = m_aliases.find(author);
  return (it == m_aliases.end()) ? author : it->second;
}

getBasename()

std::string Analysis::CalibrationDataInterfaceBase::getBasename ( const std::string & name ) const

protectedinherited

auxiliary function for retrieval of name within the directory

Definition at line 138 of file CalibrationDataInterfaceBase.cxx.

{
  // Retrieve the name within the directory starting from the full name
 
  return name.substr(name.find_last_of('/')+1, std::string::npos);
}

getBinnedScaleFactors()

const TH1 * Analysis::CalibrationDataInterfaceROOT::getBinnedScaleFactors	(	const std::string &	author,
		const std::string &	label,
		const std::string &	OP )

retrieve the binned calibration object for the given flavour label and operating point.

A null result will be returned in case of error (e.g. if the calibration object isn't binned to begin with).

Definition at line 2105 of file CalibrationDataInterfaceROOT.cxx.

{
  // Retrieve the actual histogrammed calibration scale factors, identifying the object by name.
  //
  //    author:  jet collection name
  //    label:   jet flavour label
  //    OP:      tagger working point
 
  unsigned int index;
  if (! retrieveCalibrationIndex (label, OP, author, true, index)) {
    // Return a dummy result if the object is not found
    cerr << "getBinnedScaleFactors: unable to find SF calibration for object " << fullName(author, OP, label, true) << endl;
    return 0;
  }
  CalibrationDataHistogramContainer* container = dynamic_cast<CalibrationDataHistogramContainer*>(m_objects[index]);
  return (container) ? dynamic_cast<TH1*>(container->GetValue("result")) : 0;
}

getContainername()

std::string Analysis::CalibrationDataInterfaceBase::getContainername ( const std::string & flavour, bool SF, unsigned int mapIndex = 0 ) const

protectedinherited

auxiliary function for string concatenation

auxiliary function for retrieval of container name

Definition at line 118 of file CalibrationDataInterfaceBase.cxx.

{
  // Construct the full pathname corresponding to the container indicated by the combination
  // of tagging operating point, jet flavour, and a possible extra extension. The calibration
  // container name (stored internally) is also attached.
 
  const std::vector<std::string>& effNames = m_calibrationEffNames.at(flavour);
  if (!SF && mapIndex >= effNames.size()) {
    std::cerr << "getContainername: given mapIndex=" << mapIndex << " incompatible with array size "
          << effNames.size() << "; resetting to 0" << std::endl;
    mapIndex = 0;
  }
  std::string name = SF ? m_calibrationSFNames.at(flavour) : effNames[mapIndex];
  name += SF ? "_SF" : "_Eff";
 
  return name;
}

getEfficiency() [1/3]

Analysis::CalibResult Analysis::CalibrationDataInterfaceROOT::getEfficiency	(	const CalibrationDataVariables &	variables,
		const std::string &	label,
		const std::string &	OP,
		Uncertainty	unc,
		const std::string &	flavour,
		unsigned int	numVariation = 0,
		unsigned int	mapIndex = 0 )

efficiency retrieval by name

Definition at line 1071 of file CalibrationDataInterfaceROOT.cxx.

{
  // Data efficiency retrieval identifying the requested calibration objects by name.
  // The data efficiency is computed as the product of MC efficiency and data/MC efficiency scale factor.
  // The return value is either a (value, uncertainty) or an (up, down) variation pair, as documented
  // above, and will be a dummy value in case an error occurs.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    label:         jet flavour label
  //    OP:            tagger operating point
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
  //    mapIndex:      index to the efficiency map to be used
 
  unsigned int indexSF, indexEff;
  if (! (retrieveCalibrationIndex (label, OP, variables.jetAuthor, false, indexEff, mapIndex) &&
     retrieveCalibrationIndex (label, OP, variables.jetAuthor, true, indexSF))) {
    cerr << "getEfficiency: unable to find Eff calibration for object " << fullName(variables.jetAuthor, OP, label, false, mapIndex) << " or SF calibration for object " << fullName(variables.jetAuthor, OP, label, true) << endl;
    // Return a dummy result if the object is not found
    return Analysis::dummyResult;
  }
  
  Analysis::CalibResult result;
  return (getEfficiency(variables, indexSF, indexEff, unc, numVariation, result, flavour) == Analysis::kError) ? Analysis::dummyResult : result;
}

getEfficiency() [2/3]

Analysis::CalibResult Analysis::CalibrationDataInterfaceROOT::getEfficiency	(	const CalibrationDataVariables &	variables,
		unsigned int	indexSF,
		unsigned int	indexEff,
		Uncertainty	unc,
		const std::string &	flavour,
		unsigned int	numVariation = 0 )

efficiency retrieval by index

Definition at line 1103 of file CalibrationDataInterfaceROOT.cxx.

{
  // Data efficiency retrieval identifying the requested calibration objects by index.
  // The data efficiency is computed as the product of MC efficiency and data/MC efficiency scale factor.
  // The return value is either a (value, uncertainty) or an (up, down) variation pair, as documented
  // above, and will be a dummy value in case an error occurs.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    indexSF:       index to scale factor calibration object
  //    indexEff:      index to MC efficiency object
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
 
  Analysis::CalibResult result;
  return (getEfficiency(variables, indexSF, indexEff, unc, numVariation, result, flavour) == Analysis::kError) ?
    Analysis::dummyResult : result;
}

getEfficiency() [3/3]

Analysis::CalibrationStatus Analysis::CalibrationDataInterfaceROOT::getEfficiency	(	const CalibrationDataVariables &	variables,
		unsigned int	indexSF,
		unsigned int	indexEff,
		Uncertainty	unc,
		unsigned int	numVariation,
		Analysis::CalibResult &	result,
		const std::string &	flavour )

efficiency retrieval by index

Definition at line 1126 of file CalibrationDataInterfaceROOT.cxx.

{
  // Data efficiency retrieval identifying the requested calibration objects by index.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    indexSF:       index to scale factor calibration object
  //    indexEff:      index to MC efficiency object
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
  //    result:        (value, uncertainty) or (up, down) variation pair, depending on the unc value.
  //                   A dummy value will be returned in case of an error.
 
  Analysis::CalibResult sfResult;
  Analysis::CalibrationStatus sfStatus = getScaleFactor(variables, indexSF, indexEff, unc, numVariation, sfResult, flavour);
  if (sfStatus == Analysis::kError) return sfStatus;
  Analysis::CalibResult effResult;
  Analysis::CalibrationStatus effStatus= getMCEfficiency(variables, indexEff, unc, effResult);
  if (effStatus == Analysis::kError) return effStatus;
 
  double relative = 0;
  double value = effResult.first;
  if (TMath::Abs(sfResult.first) > Analysis::CalibZERO) {
    value = std::min(effResult.first*sfResult.first, 1.);
 
    // Treat the scale factor variation cases separately since the contents of the CalibResult are different
    // (e.g. 'value' above contains the upward variation)
    if (unc == SFEigen || unc == SFNamed) {
      double valueDown = effResult.first*sfResult.second;
      result.first  = value; // up/down variataions of data-efficiency
      result.second = valueDown;
      return sfStatus;
    }
    if (value > 0.) {
      relative = effResult.second/effResult.first;
      double sfRelative = sfResult.second/sfResult.first;
      /*
    cout << "sferr=" << sfResult.second
    << "btag Calib relative=" << relative << " sfRelative=" << sfRelative << endl;
      */
      relative = TMath::Sqrt(sfRelative*sfRelative + relative*relative);
    }
  } else {
    // now never happens due to protection of SF return value:
    cerr << "ERROR: CalibrationDataInterfaceROOT::getEfficiency: SF null result, SF=" << sfResult.first << " MC eff=" << effResult.first  << "; setting SF=1." << endl;
    relative = Analysis::dummyValue;
  }
 
  result.first  = value;
  result.second = value*relative;
  // "Select" the status code for the actual calibration (it is subject to more constraints)
  return sfStatus;
}

getEigenVectorRecompositionCoefficientMap()

std::map< std::string, std::map< std::string, float > > Analysis::CalibrationDataInterfaceROOT::getEigenVectorRecompositionCoefficientMap

(

)

Get Eigenvector recomposition map after running runEigenVectorRecomposition()

Definition at line 2264 of file CalibrationDataInterfaceROOT.cxx.

                                                                               {
  if(m_coefficientMap.empty())
    cerr << "getCoefficientMap: Call runEigenVectorRecomposition() before retrieving coefficient map! " <<endl;
  return m_coefficientMap;
}

getInefficiency() [1/3]

Analysis::CalibResult Analysis::CalibrationDataInterfaceROOT::getInefficiency	(	const CalibrationDataVariables &	variables,
		const std::string &	label,
		const std::string &	OP,
		Uncertainty	unc,
		unsigned int	numVariation = 0,
		unsigned int	mapIndex = 0 )

inefficiency retrieval by name

Definition at line 1307 of file CalibrationDataInterfaceROOT.cxx.

{
  // Data inefficiency retrieval identifying the requested calibration objects by name.
  // The data efficiency is computed as the product of MC efficiency and data/MC efficiency scale factor;
  // the inefficiency is then computed as the 1 minus the efficiency.
  // The return value is either a (value, uncertainty) or an (up, down) variation pair, as documented
  // above, and will be a dummy value in case an error occurs.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    label:         jet flavour label
  //    OP:            tagger operating point
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
  //    mapIndex:      index to the efficiency map to be used
 
  unsigned int indexSF, indexEff;
  if (! (retrieveCalibrationIndex (label, OP, variables.jetAuthor, false, indexEff, mapIndex) &&
     retrieveCalibrationIndex (label, OP, variables.jetAuthor, true, indexSF))) {
    cerr << "getInefficiency: unable to find Eff calibration for object "
     << fullName(variables.jetAuthor, OP, label, false, mapIndex)
     << " or SF calibration for object "
     << fullName(variables.jetAuthor, OP, label, true) << endl;
    // Return a dummy result if the object is not found
    return Analysis::dummyResult;
  }
  
  Analysis::CalibResult result;
  return (getInefficiency(variables, indexSF, indexEff, unc, numVariation, result, label) == Analysis::kError) ?
    Analysis::dummyResult : result;
}

getInefficiency() [2/3]

Analysis::CalibResult Analysis::CalibrationDataInterfaceROOT::getInefficiency	(	const CalibrationDataVariables &	variables,
		unsigned int	indexSF,
		unsigned int	indexEff,
		Uncertainty	unc,
		const std::string &	flavour,
		unsigned int	numVariation = 0 )

inefficiency retrieval by index

Definition at line 1344 of file CalibrationDataInterfaceROOT.cxx.

{
  // Data inefficiency retrieval identifying the requested calibration objects by index.
  // The data efficiency is computed as the product of MC efficiency and data/MC efficiency scale factor;
  // the inefficiency is then computed as the 1 minus the efficiency.
  // The return value is either a (value, uncertainty) or an (up, down) variation pair, as documented
  // above, and will be a dummy value in case an error occurs.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    indexSF:       index to scale factor calibration object
  //    indexEff:      index to MC efficiency object
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
 
  Analysis::CalibResult result;
  return (getInefficiency(variables, indexSF, indexEff, unc, numVariation, result, flavour) == Analysis::kError) ?
    Analysis::dummyResult : result;
}

getInefficiency() [3/3]

Analysis::CalibrationStatus Analysis::CalibrationDataInterfaceROOT::getInefficiency	(	const CalibrationDataVariables &	variables,
		unsigned int	indexSF,
		unsigned int	indexEff,
		Uncertainty	unc,
		unsigned int	numVariation,
		Analysis::CalibResult &	result,
		const std::string &	flavour )

inefficiency retrieval by index

Definition at line 1368 of file CalibrationDataInterfaceROOT.cxx.

{
  // Data inefficiency retrieval identifying the requested calibration objects by index.
  // The data efficiency is computed as the product of MC efficiency and data/MC efficiency scale factor;
  // the inefficiency is then computed as the 1 minus the efficiency.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    indexSF:       index to scale factor calibration object
  //    indexEff:      index to MC efficiency object
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
  //    result:        (value, uncertainty) or (up, down) variation pair, depending on the unc value.
  //                   A dummy value will be returned in case of an error.
 
  Analysis::CalibResult sfResult;
  Analysis::CalibrationStatus sfStatus = getScaleFactor(variables, indexSF, indexEff, unc, numVariation, sfResult, flavour);
  if (sfStatus == Analysis::kError) return sfStatus;
  Analysis::CalibResult effResult;
  Analysis::CalibrationStatus effStatus= getMCEfficiency(variables, indexEff, unc, effResult);
  if (effStatus == Analysis::kError) return effStatus;
  
  double val = std::max(0., 1. - effResult.first * sfResult.first);
  double err = 0.; // Analysis::dummyValue;
 
  // Bail out here if not both results are strictly positive
  if (effResult.first <= 0. || sfResult.first <= 0.) return Analysis::kError;
 
  // Treat the scale factor variation cases separately since the contents of the CalibResult are different
  // (e.g. 'val' above contains the upward variation)
  if (unc == SFEigen || unc == SFNamed) {
    double valDown = std::max(0., 1. - effResult.first*sfResult.second);
 
    result.first = val;
    result.second = valDown;
  } else {
    // safer than pow(x, 2):
    err = effResult.second/effResult.first*effResult.second/effResult.first
      + sfResult.second/sfResult.first*sfResult.second/sfResult.first;
    err = val*TMath::Sqrt(err);
 
    result.first  = std::max(0., std::min(1., val));
    result.second = err;
  }
 
  // "Select" the status code for the actual calibration (it is subject to more constraints)
  return sfStatus;
}

getInefficiencyScaleFactor() [1/3]

Analysis::CalibResult Analysis::CalibrationDataInterfaceROOT::getInefficiencyScaleFactor	(	const CalibrationDataVariables &	variables,
		const std::string &	label,
		const std::string &	OP,
		Uncertainty	unc,
		unsigned int	numVariation = 0,
		unsigned int	mapIndex = 0 )

"MC" inefficiency scale factor retrieval by name

Definition at line 1186 of file CalibrationDataInterfaceROOT.cxx.

{
  // Inefficiency scale factor retrieval identifying the requested calibration objects by name.
  // The data efficiency is computed as the product of MC efficiency and data/MC efficiency scale factor;
  // the inefficiency scale factor is then computed as the ratio of data to MC inefficiencies.
  // The return value is either a (value, uncertainty) or an (up, down) variation pair, as documented
  // above, and will be a dummy value in case an error occurs.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    label:         jet flavour label
  //    OP:            tagger operating point
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
  //    mapIndex:      index to the efficiency map to be used
 
  unsigned int indexSF, indexEff;
  if (! (retrieveCalibrationIndex (label, OP, variables.jetAuthor, false, indexEff, mapIndex) &&
     retrieveCalibrationIndex (label, OP, variables.jetAuthor, true, indexSF))) {
    cerr << "getInefficiencyScaleFactor: unable to find Eff calibration for object "
     << fullName(variables.jetAuthor, OP, label, false, mapIndex)
     << " or SF calibration for object "
     << fullName(variables.jetAuthor, OP, label, true) << endl;
    // Return a dummy result if the object is not found
    return Analysis::dummyResult;
  }
  
  Analysis::CalibResult result;
  return (getInefficiencyScaleFactor(variables, indexSF, indexEff, unc, numVariation, result, label) == Analysis::kError) ?
    Analysis::dummyResult : result;
}

getInefficiencyScaleFactor() [2/3]

Analysis::CalibResult Analysis::CalibrationDataInterfaceROOT::getInefficiencyScaleFactor	(	const CalibrationDataVariables &	variables,
		unsigned int	indexSF,
		unsigned int	indexEff,
		Uncertainty	unc,
		const std::string &	flavour,
		unsigned int	numVariation = 0 )

"MC" inefficiency scale factor retrieval by index

Definition at line 1223 of file CalibrationDataInterfaceROOT.cxx.

{
  // Inefficiency scale factor retrieval identifying the requested calibration objects by index.
  // The data efficiency is computed as the product of MC efficiency and data/MC efficiency scale factor;
  // the inefficiency scale factor is then computed as the ratio of data to MC inefficiencies.
  // The return value is either a (value, uncertainty) or an (up, down) variation pair, as documented
  // above, and will be a dummy value in case an error occurs.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    indexSF:       index to scale factor calibration object
  //    indexEff:      index to MC efficiency object
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
 
  Analysis::CalibResult result;
  return (getInefficiencyScaleFactor(variables, indexSF, indexEff, unc, numVariation, result, flavour) == Analysis::kError) ?
    Analysis::dummyResult : result;
}

getInefficiencyScaleFactor() [3/3]

Analysis::CalibrationStatus Analysis::CalibrationDataInterfaceROOT::getInefficiencyScaleFactor	(	const CalibrationDataVariables &	variables,
		unsigned int	indexSF,
		unsigned int	indexEff,
		Uncertainty	unc,
		unsigned int	numVariation,
		Analysis::CalibResult &	result,
		const std::string &	flavour )

"MC" inefficiency scale factor retrieval by index

Definition at line 1247 of file CalibrationDataInterfaceROOT.cxx.

{
  // Inefficiency scale factor retrieval identifying the requested calibration objects by index.
  // The data efficiency is computed as the product of MC efficiency and data/MC efficiency scale factor;
  // the inefficiency scale factor is then computed as the ratio of data to MC inefficiencies.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    indexSF:       index to scale factor calibration object
  //    indexEff:      index to MC efficiency object
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
  //    result:        (value, uncertainty) or (up, down) variation pair, depending on the unc value.
  //                   A dummy value will be returned in case of an error.
 
  Analysis::CalibResult sfResult;
  Analysis::CalibrationStatus sfStatus = getScaleFactor(variables, indexSF, indexEff, unc, numVariation, sfResult, flavour);
  if (sfStatus == Analysis::kError) return sfStatus;
  Analysis::CalibResult effResult;
  Analysis::CalibrationStatus effStatus= getMCEfficiency(variables, indexEff, unc, effResult);
  if (effStatus == Analysis::kError) return effStatus;
 
  double eff = std::min(effResult.first, 1.);
  //  double efferr = effResult.second; // not needed as (per the code change indicated below) we are not doing anything with MC statistical uncertainties
  double sf = sfResult.first; 
  double sferr = sfResult.second; 
 
  double val = 0.; // Analysis::dummyValue;
  double err = 0.; // Analysis::dummyValue;
  if (1. - eff > CalibZERO) {
    // Protect against negative scale factors
    val = std::max((1. - eff*sf), CalibZERO) / (1. - eff);
    // Treat the scale factor variation cases separately since the contents of the CalibResult are different
    // ('sf' and 'sferr' above contain the upward and downward variations, respectively).
    if (unc == SFEigen || unc == SFNamed) {
      double valDown = std::max((1. - eff*sferr), CalibZERO) / (1. - eff);
      result.first  = val;
      result.second = valDown;
      return sfStatus;
    }
    // When using eigenvector (or named) variations (as above), only scale factor variations are considered.
    // For the sake of consistency, it has been decided (see https://its.cern.ch/jira/browse/AFT-350) to remove them also when EV variations aren't used
    //err = pow((1. - sf) / (1. - eff) * efferr, 2) + pow(eff*sferr, 2);
    err = pow(eff*sferr, 2);
    if (err > 0.)
      err = 1./(1. - eff) * TMath::Sqrt(err);
    // cout << "btag Calib Ineff err=" << err << endl;
  }
 
  result.first  = std::max(CalibZERO, val);
  result.second = err;
  // "Select" the status code for the actual calibration (it is subject to more constraints)
  return sfStatus;
}

getMCEfficiency() [1/3]

Analysis::CalibResult Analysis::CalibrationDataInterfaceROOT::getMCEfficiency	(	const CalibrationDataVariables &	variables,
		const std::string &	label,
		const std::string &	OP,
		Uncertainty	unc = None,
		unsigned int	mapIndex = 0 )

"MC" efficiency retrieval by name

Definition at line 963 of file CalibrationDataInterfaceROOT.cxx.

{
  // MC efficiency retrieval identifying the requested calibration object by name.
  // The return value is a (value, uncertainty) pair, as documented above, and will
  // be a dummy value in case an error occurs.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    label:         jet flavour label
  //    OP:            tagger operating point
  //    unc:           keyword indicating what uncertainties to evaluate
  //    mapIndex:      index to the efficiency map to be used
 
  unsigned int index;
  if (! retrieveCalibrationIndex (label, OP, variables.jetAuthor, false, index, mapIndex)) {
    cerr << "getMCEfficiency: unable to find Eff calibration for object " << fullName(variables.jetAuthor, OP, label, false, mapIndex) << endl;
    // Return a dummy result if the object is not found
    return Analysis::dummyResult;
  }
  
  Analysis::CalibResult result;
  return (getMCEfficiency(variables, index, unc, result) == Analysis::kError) ?
    Analysis::dummyResult : result;
}

getMCEfficiency() [2/3]

Analysis::CalibrationStatus Analysis::CalibrationDataInterfaceROOT::getMCEfficiency	(	const CalibrationDataVariables &	variables,
		unsigned int	index,
		Uncertainty	unc,
		Analysis::CalibResult &	result )

"MC" efficiency retrieval by index

Definition at line 1009 of file CalibrationDataInterfaceROOT.cxx.

{
  // MC efficiency retrieval identifying the requested calibration object by index.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    index:         index to calibration object
  //    unc:           keyword indicating what uncertainties to evaluate
  //    result:        (value, uncertainty) variation pair.
  //                   A dummy value will be returned in case of an error.
 
  CalibrationDataContainer* container = m_objects[index];
  if (! container) return Analysis::kError;
  
  // perform out-of-bound check of jet eta
  if (!checkAbsEta(variables, index)) {
    if (m_verbose)
      cerr << "Jet |eta| is outside of the boundary!" << endl;
    return Analysis::kRange;
  }
 
 
  // always retrieve the result itself
  double value;
  Analysis::CalibrationStatus status = container->getResult(variables, value);
  if (status == Analysis::kError) return status;
  if (m_otherStrategy == GiveUp)
    assert (status != Analysis::kRange); // no need to test also statDown
  else if (m_otherStrategy == Flag)
    if (status == Analysis::kRange)
      this->increaseCounter(index);
 
  // retrieve the statistical uncertainty if desired
  double stat(0);
  if (unc == Total || unc == Statistical) {
    if (container->getStatUncertainty(variables, stat) == Analysis::kError) {
      cerr << "getMCEfficiency: error retrieving MC efficiency parameter covariance matrix!" << endl;
    return Analysis::kError;
    }
  }
 
  // Temporary(?) hack: comment this out since the present MC results don't have "systematics" contributions
  // Analysis::UncertaintyResult resSyst(0,0);
  // if (unc == Total || unc == Systematic) {
  //   if (container->getSystUncertainty(variables, resSyst) == Analysis::kError)
  //     cerr << "getScaleFactor: error retrieving Scale factor parameter covariance matrix!"
  //       << endl;
  // }
 
  // since there is no combination of stat/syst uncertainties to be made, comment this out too
  double uncertainty = stat; // combinedUncertainty(stat, resSyst);
  result.first  = std::max(0., std::min(1., value));
  result.second = uncertainty;
 
  return status;
}

getMCEfficiency() [3/3]

Analysis::CalibResult Analysis::CalibrationDataInterfaceROOT::getMCEfficiency	(	const CalibrationDataVariables &	variables,
		unsigned int	index,
		Uncertainty	unc = None )

"MC" efficiency retrieval by index

Definition at line 991 of file CalibrationDataInterfaceROOT.cxx.

{
  // MC efficiency retrieval identifying the requested calibration object by index.
  // The return value is a (value, uncertainty) pair, as documented above, and will
  // be a dummy value in case an error occurs.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    index:         index to calibration object
  //    unc:           keyword indicating what uncertainties to evaluate
 
  Analysis::CalibResult result;
  return (getMCEfficiency(variables, index, unc, result) == Analysis::kError) ?
    Analysis::dummyResult : result;
}

getMCEfficiencyObject()

const TObject * Analysis::CalibrationDataInterfaceROOT::getMCEfficiencyObject	(	const std::string &	author,
		const std::string &	label,
		const std::string &	OP,
		unsigned int	mapIndex = 0 )

retrieve the MC efficiency (central values) object for the given flavour label and operating point.

A null result will be returned in case of error (e.g. if the calibration object isn't binned to begin with). It is the user's responsibility to verify whether the object derives from a TH1 or a TF1.

Definition at line 2127 of file CalibrationDataInterfaceROOT.cxx.

{
  // Retrieve the actual central values object for the MC efficiences, identifying the object by name.
  // The object returned can be either a TH1 or a TF1; it is up to the user to determine which.
  //
  //    author:   jet collection name
  //    label:    jet flavour label
  //    OP:       tagger working point
  //    mapIndex: index to the efficiency map to be used
 
  unsigned int index;
  if (! retrieveCalibrationIndex (label, OP, author, false, index, mapIndex)) {
    // Return a dummy result if the object is not found
    cerr << "getMCEfficiencyObject: unable to find efficiency calibration for object "
     << fullName(author, OP, label, false, mapIndex) << endl;
    return 0;
  }
  CalibrationDataContainer* container = m_objects[index];
  return (container) ? container->GetValue("result") : 0;
}

getMCInefficiency() [1/2]

Analysis::CalibResult Analysis::CalibrationDataInterfaceROOT::getMCInefficiency	(	const CalibrationDataVariables &	variables,
		const std::string &	label,
		const std::string &	OP,
		Uncertainty	unc = None,
		unsigned int	mapIndex = 0 )

"MC" inefficiency retrieval by name

Definition at line 1422 of file CalibrationDataInterfaceROOT.cxx.

{
  // Data inefficiency retrieval identifying the requested calibration objects by name.
  // The inefficiency is computed as the 1 minus the efficiency.
  // The return value is a (value, uncertainty), as documented above, and will be a dummy value
  // in case an error occurs.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    label:         jet flavour label
  //    OP:            tagger operating point
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
  //    mapIndex:      index to the efficiency map to be used
 
  Analysis::CalibResult effResult = getMCEfficiency(variables, label, OP, unc, mapIndex);
  return std::make_pair(std::max(0., 1. - effResult.first), effResult.second);
}

getMCInefficiency() [2/2]

Analysis::CalibResult Analysis::CalibrationDataInterfaceROOT::getMCInefficiency	(	const CalibrationDataVariables &	variables,
		unsigned int	index,
		Uncertainty	unc = None )

"MC" inefficiency retrieval by index

Definition at line 1445 of file CalibrationDataInterfaceROOT.cxx.

{
  // MC inefficiency retrieval identifying the requested calibration object by index.
  // The inefficiency is computed as the 1 minus the efficiency.
  // The return value is a (value, uncertainty), as documented above, and will be a dummy value
  // in case an error occurs.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    index:         index to MC efficiency object
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
 
  Analysis::CalibResult effResult = getMCEfficiency(variables, index, unc);
  return std::make_pair(std::max(0., 1. - effResult.first), effResult.second);
}

getMCMCScaleFactor()

double Analysis::CalibrationDataInterfaceROOT::getMCMCScaleFactor	(	const CalibrationDataVariables &	variables,
		unsigned	indexSF,
		unsigned int	indexEff ) const

MC/MC scale factor retrieval.

Normally this is to be used only internally; however, since this information may be of interest it is made public anyway.

Definition at line 1465 of file CalibrationDataInterfaceROOT.cxx.

{
  // Retrieve the MC/MC scale factor given the set of scale factor and efficiency indices.
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    indexSF:       index to scale factor calibration object
  //    indexEff:      index to MC efficiency object
  
  // If either reference doesn't exist, or if they are the same, nothing can / needs to be done.
  int indexSFRef = m_hadronisationReference[indexSF], indexEffRef = m_hadronisationReference[indexEff];
  if (indexSFRef < 0 || indexEffRef < 0 || indexSFRef == indexEffRef) return 1;
 
  // Verify also that the individual efficiencies are physically meaningful.
  double effSFRef;  m_objects[indexSFRef]->getResult(variables, effSFRef);
  double effEffRef; m_objects[indexEffRef]->getResult(variables, effEffRef);
  return (effSFRef > 0 && effEffRef > 0) ? effSFRef/effEffRef : 1;
}

getNumVariations() [1/2]

unsigned int Analysis::CalibrationDataInterfaceROOT::getNumVariations	(	const std::string &	author,
		const std::string &	label,
		const std::string &	OP,
		Uncertainty	unc )

retrieve the number of variations relevant to the calibration object.

The Uncertainty enum is used to specify the category.

Definition at line 2060 of file CalibrationDataInterfaceROOT.cxx.

{
  // Retrieve the number of eigenvector variations or named variations relevant for
  // the given scale factor calibration object, identifying the object by name.
  //
  //    author:  jet collection name
  //    label:   jet flavour label
  //    OP:      tagger working point
  //    unc:     should be set to SFEigen or SFNamed for the cases of
  //             eigenvector variations or named variations, respectively
 
  unsigned int index;
 
  if (! retrieveCalibrationIndex (label, OP, author, true, index)) return 0;
  return getNumVariations(index, unc, label);
}

getNumVariations() [2/2]

unsigned int Analysis::CalibrationDataInterfaceROOT::getNumVariations	(	unsigned int	index,
		Uncertainty	unc,
		const std::string &	flavour )

retrieve the number of variations relevant to the calibration object.

The Uncertainty enum is used to specify the category.

Definition at line 2082 of file CalibrationDataInterfaceROOT.cxx.

{
  // Retrieve the number of eigenvector variations or named variations relevant for
  // the given scale factor calibration object, identifying the object by index.
  //
  //    index:   index to calibration scale factor object
  //    unc:     should be set to SFEigen or SFNamed for the cases of
  //             eigenvector variations or named variations, respectively
 
  if (! (unc == SFEigen || unc == SFNamed || unc == SFGlobalEigen)) return 0;
  CalibrationDataContainer* container = m_objects[index];
  if (! container) return 0;
  std::shared_ptr<CalibrationDataEigenVariations> eigenVariation=m_eigenVariationsMap.at(container);
  if (unc == SFGlobalEigen){
    std::shared_ptr<CalibrationDataGlobalEigenVariations> GEV = std::dynamic_pointer_cast<CalibrationDataGlobalEigenVariations>(eigenVariation);
    return GEV->getNumberOfEigenVariations(flavour);
  }
  return (unc == SFEigen) ? eigenVariation->getNumberOfEigenVariations() : eigenVariation->getNumberOfNamedVariations();
}

getScaleFactor() [1/3]

Analysis::CalibResult Analysis::CalibrationDataInterfaceROOT::getScaleFactor	(	const CalibrationDataVariables &	variables,
		const std::string &	label,
		const std::string &	OP,
		Uncertainty	unc,
		unsigned int	numVariation = 0,
		unsigned int	mapIndex = 0 )

efficiency scale factor retrieval by name.

#1

Definition at line 736 of file CalibrationDataInterfaceROOT.cxx.

{
  // Scale factor retrieval identifying the requested calibration object by name.
  // The return value is either a (value, uncertainty) or an (up, down) variation pair, as documented
  // above, and will be a dummy value in case an error occurs.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    label:         jet flavour label
  //    OP:            tagger operating point
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
  //    mapIndex:      index to the efficiency map to be used (this is needed for MC/MC scale factor
  //                   application)
  unsigned int indexEff, indexSF;
  if (! (retrieveCalibrationIndex (label, OP, variables.jetAuthor, false, indexEff, mapIndex) && retrieveCalibrationIndex (label, OP, variables.jetAuthor, true, indexSF))) {
    cerr << "getScaleFactor: unable to find SF calibration for object " << fullName(variables.jetAuthor, OP, label, false, mapIndex) << " or SF calibration for object " << fullName(variables.jetAuthor, OP, label, true) << endl;
    // Return a dummy result if the object is not found
    return Analysis::dummyResult;
  }
 
  Analysis::CalibResult result; // the following is SF #3
  return (getScaleFactor(variables, indexSF, indexEff, unc, numVariation, result, label) == Analysis::kError) ?
    Analysis::dummyResult : result;
}

getScaleFactor() [2/3]

Analysis::CalibResult Analysis::CalibrationDataInterfaceROOT::getScaleFactor	(	const CalibrationDataVariables &	variables,
		unsigned int	indexSF,
		unsigned int	indexEff,
		Uncertainty	unc,
		const std::string &	flavour,
		unsigned int	numVariation = 0 )

efficiency scale factor retrieval by index #2

Definition at line 767 of file CalibrationDataInterfaceROOT.cxx.

{
  // Scale factor retrieval identifying the requested calibration object by index.
  // The return value is either a (value, uncertainty) or an (up, down) variation pair, as documented
  // above, and will be a dummy value in case an error occurs.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    indexSF:       index to scale factor calibration object
  //    indexEff:      index to MC efficiency object
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
  Analysis::CalibResult result; // the following is SF #3
  return (getScaleFactor(variables, indexSF, indexEff, unc, numVariation, result, flavour) == Analysis::kError) ?
    Analysis::dummyResult : result;
}

getScaleFactor() [3/3]

Analysis::CalibrationStatus Analysis::CalibrationDataInterfaceROOT::getScaleFactor	(	const CalibrationDataVariables &	variables,
		unsigned int	indexSF,
		unsigned int	indexEff,
		Uncertainty	unc,
		unsigned int	numVariation,
		Analysis::CalibResult &	result,
		const std::string &	flavour )

efficiency scale factor retrieval by index #3

Definition at line 788 of file CalibrationDataInterfaceROOT.cxx.

{
  // Scale factor retrieval identifying the requested calibration object by index.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    indexSF:       index to scale factor calibration object
  //    indexEff:      index to MC efficiency object
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
  //    result:        (value, uncertainty) or (up, down) variation pair, depending on the unc value.
  //                   A dummy value will be returned in case of an error.
 
  CalibrationDataContainer* container = m_objects[indexSF];
  if (! container) {
    cerr << "getScaleFactor: error retrieving container!" << endl;
    return Analysis::kError;
  }
 
  // perform out-of-bound check of jet eta
  if (!checkAbsEta(variables, indexSF)) {
    if (m_verbose)
      cerr << "Jet |eta| is outside of the boundary!" << endl;
    return Analysis::kRange;
  }
 
  // retrieve the MC/MC scale factor
  double MCMCSF = m_useMCMCSF ? getMCMCScaleFactor(variables, indexSF, indexEff) : 1; // if we don't want to switch generator, MCMCSF = 1, as it should  be
 
  if (!m_runEigenVectorMethod && (unc == SFEigen || unc == SFNamed || unc == SFGlobalEigen))
  {
     cerr << " ERROR. Trying to call eigenvector method but initialization not switched on in b-tagging configuration." << endl;
     cerr << " Please correct your configuration first. Nominal uncertainties used. " << endl;
  }
 
  // Procede with eigenvariations methods i.e. return the SF variations
  if (unc == SFEigen || unc == SFNamed || unc==SFGlobalEigen) {
    std::shared_ptr<CalibrationDataEigenVariations> eigenVariation;
    try {
      eigenVariation=m_eigenVariationsMap.at(container);
    } catch (const std::out_of_range&) {
      cerr << " Could not retrieve eigenvector variation, while it should have been there." << endl;
      return Analysis::kError;
    }
    TH1* up=0;
    TH1* down=0;
    bool extrapolate = false; // store if the numVariation is the extrapolation named uncertainty index
    if (unc == SFEigen || unc==SFNamed){
      unsigned int maxVariations = (unc == SFEigen) ? eigenVariation->getNumberOfEigenVariations() : eigenVariation->getNumberOfNamedVariations();
      if (numVariation > maxVariations-1) {
        cerr << "Asked for " << ((unc == SFEigen) ? "eigenvariation" : "named variation") << " number: " << numVariation << " but overall number of available variations is: " << maxVariations << endl;
        return Analysis::kError;
      }
      bool isOK = eigenVariation->getEigenvectorVariation(numVariation,up,down);
      if (!isOK) {
        cerr << "Eigenvector object is there but cannot retrieve up and down uncertainty histograms." << endl;
        return Analysis::kError;
      }
      // the 'extrapolation' uncertainty (always a named one) needs a somewhat special treatment
      extrapolate = (unc == SFNamed) ? eigenVariation->isExtrapolationVariation(numVariation) : false;
 
    } else if (unc == SFGlobalEigen) {
      std::shared_ptr<CalibrationDataGlobalEigenVariations> GEV = std::dynamic_pointer_cast<CalibrationDataGlobalEigenVariations>(eigenVariation); //dynamic_cast<std::shared_ptr<CalibrationDataGlobalEigenVariations> >(eigenVariation);
      if (not GEV){
        cerr << "Analysis::CalibrationDataInterfaceROOT::getScaleFactor: dynamic cast failed\n";
        return Analysis::kError;
      }
      unsigned int maxVariations = GEV->getNumberOfEigenVariations(flavour); // <----- This gets the number of variations of the flavour
      if (numVariation > maxVariations-1) {
        cerr << "Asked for global eigenvariation number: " << numVariation << " but overall number of available variations is: " << maxVariations << endl;
        return Analysis::kError;
      }
      bool isOK = GEV->getEigenvectorVariation(flavour, numVariation,up,down);
      if (!isOK) {
        cerr << "Eigenvector object is there but cannot retrieve up and down uncertainty histograms." << endl;
        return Analysis::kError;
      }
      // the 'extrapolation' uncertainty (always a named one) needs a somewhat special treatment
      extrapolate = GEV->isExtrapolationVariation(numVariation, flavour);
    } else {
      std::cerr << "ERROR: you requested " << unc << " but that isn't in the set of (SFEigen, SFGlobalEigen, SFNamed)  for eigenvariations. " << std::endl;
      return Analysis::kError;
    }
    
    double valueUp;
    double valueDown;
    Analysis::CalibrationStatus statUp   = container->getResult(variables, valueUp,  up,   extrapolate); // This is what actually retrieves results from the container
    Analysis::CalibrationStatus statDown = container->getResult(variables, valueDown,down, extrapolate);
 
    if (statUp == Analysis::kError || statDown == Analysis::kError)
      return Analysis::kError;
    if (m_otherStrategy == GiveUp)
      assert (statUp != Analysis::kRange); // no need to test also statDown
    else if (m_otherStrategy == GiveUpExtrapolated)
      assert (statUp != Analysis::kExtrapolatedRange); // no need to test also statDown
    else if (m_otherStrategy == Flag) {
      if (statUp == Analysis::kRange)
        increaseCounter(indexSF);
      else if (statUp == Analysis::kExtrapolatedRange)
        increaseCounter(indexSF, Extrapolated);
    }
 
    result.first  = MCMCSF*valueUp;
    result.second = MCMCSF*valueDown;
 
    // Prevent negative return values. Should the comparison be against a strict 0?
    result.first  = std::max(Analysis::CalibZERO, result.first);
    result.second = std::max(Analysis::CalibZERO, result.second);
 
    return statUp; // end of getScaleFactor if SFEigen, SFGlobalEigen, or SFNamed is set
 
 
  } // The above returns the up/down varied scale factor
  //Proceed with no-eigenvector result
 
  // always retrieve the result itself
  double value;
  Analysis::CalibrationStatus status = container->getResult(variables, value);
  if (status == Analysis::kError) {
    cerr << "getScaleFactor: error retrieving result in non-EV context!" << endl;
    return status;
  }
  if (m_otherStrategy == GiveUp){
    assert (status != Analysis::kRange);
  } else if (m_otherStrategy == GiveUpExtrapolated) {
    assert (status != Analysis::kExtrapolatedRange);
  } else if (m_otherStrategy == Flag) {
    if (status == Analysis::kRange){
      increaseCounter(indexSF);
    } else if (status == Analysis::kExtrapolatedRange) {
      increaseCounter(indexSF, Extrapolated);
    }
  }
 
  // retrieve the statistical uncertainty if desired
  double stat(0);
  if (unc == Total || unc == Statistical) {
    if (container->getStatUncertainty(variables, stat) == Analysis::kError) {
      cerr << "getScaleFactor: error retrieving Scale factor parameter covariance matrix!" << endl;
      return Analysis::kError;
    }
  }
 
  Analysis::UncertaintyResult resSyst(0,0);
  if (unc == Total || unc == Systematic) {
    if (container->getSystUncertainty(variables, resSyst) == Analysis::kError) {
      cerr << "getScaleFactor: error retrieving Scale factor parameter systematic uncertainty!" << endl;
      return Analysis::kError;
    }
  } else if (unc == Extrapolation) {
    // this uncertainty is special, since it is not normally to be combined into the overall systematic uncertainty
    if (container->getUncertainty("extrapolation", variables, resSyst) == Analysis::kError)
      cerr << "getScaleFactor: error retrieving Scale factor parameter extrapolation uncertainty!" << endl;
  } else if (unc == TauExtrapolation) {
    // also this uncertainty is special, since it it singles out an uncertainty relevant only for tau "jets",
    // and some care has to be taken not to duplicate or omit uncertainties
    if (container->getUncertainty("extrapolation from charm", variables, resSyst) == Analysis::kError)
      cerr << "getScaleFactor: error retrieving Scale factor parameter extrapolation uncertainty!" << endl;
  }
 
  double uncertainty = combinedUncertainty(stat, resSyst);
  result.first  = MCMCSF*value;
  result.second = MCMCSF*uncertainty;
 
  // Prevent negative return values. Should the comparison be against a strict 0?
  result.first = std::max(Analysis::CalibZERO, result.first);
  return status;
 
}

getScaleFactorCovarianceMatrix()

TMatrixDSym Analysis::CalibrationDataInterfaceROOT::getScaleFactorCovarianceMatrix	(	const std::string &	author,
		const std::string &	label,
		const std::string &	OP,
		const std::string &	unc = "all" )

retrieve the named covariance matrix element corresponding to the binned calibration object.

The unc argument should correspond to a given source of statistical or systematic uncertainty, or "all" (in case the full covariance matrix is required) For 2D and 3D histograms, the bin numbering follows the "global" bin number as defined by class TH1.

Definition at line 2412 of file CalibrationDataInterfaceROOT.cxx.

{
  // Return the scale factor covariance matrix for the given calibration object.
  // This function is deprecated since its functionality is duplicated in the
  // CalibrationDataEigenVariations class.
  //
  //    author:  jet collection name
  //    label:   jet flavour label
  //    OP:      tagger working point
  //    unc:     source of uncertainty to consider
  // Catch issues with the specified input as early as possible
  TMatrixDSym dummy;
  if (unc == "comment" || unc == "result" || unc == "combined") return dummy;
 
  unsigned int index;
  if (! retrieveCalibrationIndex (label, OP, author, true, index)) {
    // Return a dummy result if the object is not found
    cerr << "getScaleFactorCovarianceMatrix: unable to find SF calibration for object " << fullName(author, OP, label, true) << endl;
    return dummy;
  }
  CalibrationDataHistogramContainer* container = dynamic_cast<CalibrationDataHistogramContainer*>(m_objects[index]);
  if (!container) return dummy;
 
  // retrieve the central calibration and its axes
  TH1* result = dynamic_cast<TH1*>(container->GetValue("result"));
  if (! result) return dummy;
  // "normal" case: single source of uncertainty
  if (unc != "all") {
    if (unc == "statistics") {
      return getStatCovarianceMatrix(result);
    } else {
      TH1* hunc = dynamic_cast<TH1*>(container->GetValue(unc.c_str()));
      if (! hunc) {
    cout << "getScaleFactorCovarianceMatrix: no uncertainty object found "
         << "corresponding to name " << unc << endl;
    return dummy;
      }
      return getSystCovarianceMatrix(result, hunc, container->isBinCorrelated(unc), unc, container->getTagWeightAxis());
    }
  }
 
  // special case: complete covariance matrix. This is to be constructed
  // as the sum over all individual contributions.
  // First, treat the statistics separately (as above)
  TMatrixDSym cov = getStatCovarianceMatrix(result);
 
  // Then loop through the list of (other) uncertainties
  std::vector<string> uncs = container->listUncertainties();
  for (unsigned int t = 0; t < uncs.size(); ++t) {
    if (uncs[t] == "comment" || uncs[t] == "result" || uncs[t] == "combined" || 
        uncs[t] == "statistics" || uncs[t]=="extrapolation" || uncs[t]=="MChadronisation" || 
        uncs[t]=="ReducedSets" || uncs[t]=="systematics") continue;
    TH1* hunc = dynamic_cast<TH1*>(container->GetValue(uncs[t].c_str()));
    if (not hunc) {
      std::cerr<<"Analysis::CalibrationDataInterfaceROOT::getScaleFactorCovarianceMatrix : dynamic cast failed\n";
      continue;
    }
    TMatrixDSym syst_cov = getSystCovarianceMatrix(result, hunc, container->isBinCorrelated(uncs[t]), uncs[t], container->getTagWeightAxis());
    cov += syst_cov;
  }
 
  return cov;
}

getShiftedScaleFactors()

const TH1 * Analysis::CalibrationDataInterfaceROOT::getShiftedScaleFactors	(	const std::string &	author,
		const std::string &	label,
		const std::string &	OP,
		const std::string &	unc,
		double	sigmas )

retrieve the binned calibration object for the given flavour label and operating point, with the result shifted by the given number of standard deviations for the given systematic uncertainty.

A null result will be returned in case of error (e.g. if the calibration object isn't binned to begin with, or if the uncertainty asked for isn't fully correlated from bin to bin).

Definition at line 2155 of file CalibrationDataInterfaceROOT.cxx.

{
  // Retrieve the actual histogrammed calibration scale factors, identifying the object by name
  // and with the scale factors shifted by the uncertainties due to the given source of uncertainty
  // (where bin-to-bin correlations are accounted for, i.e., shifts may be either positive or negative).
  //
  //    author:  jet collection name
  //    label:   jet flavour label
  //    OP:      tagger working point
  //    unc:     source of uncertainty to consider
  //    sigmas:  number of standard deviations by which to shift the scale factor central values
 
  // quick sanity check
  if (unc == "comment" || unc == "result" || unc == "combined" || unc == "statistics") return 0;
 
  unsigned int index;
  if (! retrieveCalibrationIndex (label, OP, author, true, index)) {
    // Return a null result if the object is not found
    cerr << "getShiftedScaleFactors: unable to find SF calibration for object " << fullName(author, OP, label, true) << endl;
    return nullptr;
  }
  CalibrationDataHistogramContainer* container = dynamic_cast<CalibrationDataHistogramContainer*>(m_objects[index]);
  if (! container) return nullptr;
 
  TH1* result = dynamic_cast<TH1*>(container->GetValue("result"));
  TH1* hunc = dynamic_cast<TH1*>(container->GetValue(unc.c_str()));
  // another sanity check...
  if ((! hunc) || (! result)) return nullptr;
  if (hunc->GetDimension() != result->GetDimension() || hunc->GetNbinsX() != result->GetNbinsX() ||
      hunc->GetNbinsX() != result->GetNbinsX() || hunc->GetNbinsX() != result->GetNbinsX())
    return nullptr;
  // also check that the uncertainty is to be treated as correlated from bin to bin
  // (for the variation is applied coherently, which isn't appropriate for uncertainties
  // that aren't correlated from bin to bin)
  if (! container->isBinCorrelated(unc)) return 0;
 
  // if everything is consistent, the actual operation simply consists of adding histograms...
  std::string name(container->GetName()); name += "_"; name += unc; name += "_";
  TH1* shifted = dynamic_cast<TH1*>(result->Clone(name.c_str()));
  if (not shifted) return nullptr;
  shifted->Add(hunc, sigmas);
  return shifted;
}

getWeightScaleFactor() [1/3]

Analysis::CalibResult Analysis::CalibrationDataInterfaceROOT::getWeightScaleFactor	(	const CalibrationDataVariables &	variables,
		const std::string &	label,
		Uncertainty	unc,
		unsigned int	numVariation = 0,
		unsigned int	mapIndex = 0 )

efficiency scale factor retrieval by name

Definition at line 1485 of file CalibrationDataInterfaceROOT.cxx.

{
    // #1
  // Tag weight fraction scale factor retrieval identifying the requested calibration object by name.
  // The return value is either a (value, uncertainty) or (if eigenvector or named variations are specified)
  // an (up, down) variation pair, and will be a dummy value in case an error occurs.
  // Note that in contrast to the "regular" (non-continuous) case, the computation of the scale factor in
  // general needs the (selection- or even process-specific) MC tag weight fractions, in order to rescale
  // scale factors. This is used to ensure that the tag weight fractions (both in data and in MC) sum up to
  // unity for each given kinematic bin.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    label:         jet flavour label
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
  //    mapIndex:      index to the MC efficiency map to be used for scale factor rescaling
 
  static const string cont("Continuous");
 
  unsigned int indexSF, indexEff;
  if (! (retrieveCalibrationIndex (label, cont, variables.jetAuthor, false, indexEff, mapIndex) &&
     retrieveCalibrationIndex (label, cont, variables.jetAuthor, true, indexSF))) {
    cerr << "getWeightScaleFactor: unable to find Eff calibration for object "
     << fullName(variables.jetAuthor, cont, label, false, mapIndex)
     << " or SF calibration for object "
     << fullName(variables.jetAuthor, cont, label, true) << endl;
    return Analysis::dummyResult;
  }
 
  Analysis::CalibResult result;
  return (getWeightScaleFactor(variables, indexSF, indexEff, unc, numVariation, result) == Analysis::kError) ? Analysis::dummyResult : result;
}

getWeightScaleFactor() [2/3]

Analysis::CalibrationStatus Analysis::CalibrationDataInterfaceROOT::getWeightScaleFactor	(	const CalibrationDataVariables &	variables,
		unsigned int	indexSF,
		unsigned int	indexEff,
		Uncertainty	unc,
		unsigned int	numVariation,
		Analysis::CalibResult &	result )

efficiency scale factor retrieval by index, with different signature

Definition at line 1550 of file CalibrationDataInterfaceROOT.cxx.

{
    // #3
  // Tag weight fraction scale factor retrieval identifying the requested calibration object by index.
  // Note that in contrast to the "regular" (non-continuous) case, the computation of the scale factor in
  // general needs the (selection- or even process-specific) MC tag weight fractions, in order to rescale
  // scale factors. This is used to ensure that the tag weight fractions (both in data and in MC) sum up to
  // unity for each given kinematic bin.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    indexSF:       index to calibration object
  //    indexEff:      index to MC tag weight
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
  //    result:        (value, uncertainty) or (up, down) variation pair, depending on the unc value.
  //                   A dummy value will be returned in case of an error.
  CalibrationDataContainer* container = m_objects[indexSF];
  if (! container) return Analysis::kError;
  CalibrationDataContainer* effContainer = m_objects[indexEff];
  if (! effContainer) return Analysis::kError;
 
  // the first time this combination of scale factor and "efficiency" objects is given, check on the
  // scale factors that will result from their combination (where the computations reproduce those
  // shown below)
  checkWeightScaleFactors(indexSF, indexEff);
 
  // perform out-of-bound check of jet eta
  if (!checkAbsEta(variables, indexSF)) {
    if (m_verbose)
      cerr << "Jet |eta| is outside of the boundary!" << endl;
    return Analysis::kRange;
  }
 
  // Always retrieve the result itself 
  double value;
  Analysis::CalibrationStatus status = container->getResult(variables, value);
  if (status == Analysis::kError) return status;
  if (m_otherStrategy == GiveUp) assert (status != Analysis::kRange);
  else if (m_otherStrategy == GiveUpExtrapolated) assert (status != Analysis::kExtrapolatedRange);
  else if (m_otherStrategy == Flag) {
    if (status == Analysis::kRange)
      increaseCounter(indexSF);
    else if (status == Analysis::kExtrapolatedRange)
      increaseCounter(indexSF, Extrapolated);
  }
 
  // Retrieve the reference MC tag weight fraction (corresponding to the calibration scale factors)
  Analysis::UncertaintyResult refMCResult(0,0);
  if (container->getUncertainty("MCreference", variables, refMCResult) == Analysis::kError)
    return Analysis::kError;
  double fracMCref = refMCResult.first;
  // Retrieve the MC reference information, if requested (the initialisation below is to make sure
  // that no exceptions in the code will be needed)
  double fracSFref = fracMCref, fracEffref = fracMCref;
  if (m_useMCMCSF) {
    int indexSFref = m_hadronisationReference[indexSF], indexEffref = m_hadronisationReference[indexEff];
    if (indexSFref < 0 || indexEffref < 0) {
      cerr << "getWeightScaleFactor: error: generator-specific corrections requested but necessary reference containers lacking " << endl;
      return Analysis::kError;
    } else {
      m_objects[indexSFref]->getResult(variables, fracSFref);
      m_objects[indexEffref]->getResult(variables, fracEffref);
      if (! (fracSFref > 0. && fracEffref > 0.)) {
        cerr << "getWeightScaleFactor: error: invalid reference tag weight fraction " <<fracSFref <<" " <<fracEffref << std::endl;
        return Analysis::kError;
      }
    }
  }
 
  // Retrieve the MC tag weight fraction for the sample we need to reweight to
  double fracMCnew;
  Analysis::CalibrationStatus effStatus = effContainer->getResult(variables, fracMCnew);
  if (effStatus == Analysis::kError) return effStatus;
  if (m_otherStrategy == GiveUp) assert (effStatus != Analysis::kRange);
  else if (m_otherStrategy == Flag)
    if (effStatus == Analysis::kRange) increaseCounter(indexEff);
  // since we need to divide by this quantity, check that it is well-defined
  if (!(fracMCnew > 0.) and m_useTopologyRescaling) {// but we only care if using topology rescaling
    cerr << "getWeightScaleFactor: error: null fracMCnew would lead to invalid operation" << endl;
    return Analysis::kError;
  }
 
  if (!m_runEigenVectorMethod && (unc == SFEigen || unc == SFNamed)) {
     cerr << "getWeightScaleFactor: ERROR. Trying to call eigenvector method but initialization not switched on in b-tagging .env config file." << endl;
     cerr << "                             Please correct your .env config file first. Nominal uncertainties used. " << endl;
  }
 
  if (unc == SFEigen || unc == SFNamed) {
    std::shared_ptr<CalibrationDataEigenVariations> eigenVariation;
    try {
      eigenVariation = m_eigenVariationsMap.at(container);
    } catch (const std::out_of_range&) {
      cerr << "getWeightScaleFactor: could not retrieve eigenvector variation, while it should have been there." << endl;
      return Analysis::kError;
    }
    unsigned int maxVariations = (unc == SFEigen) ? eigenVariation->getNumberOfEigenVariations() : eigenVariation->getNumberOfNamedVariations();
    if (numVariation > maxVariations-1) {
      cerr << "getWeightScaleFactor: asked for " << ((unc == SFEigen) ? "eigenvariation" : "named variation") << " number: " << numVariation << " but overall number of available variations is: " << maxVariations << endl;
      return Analysis::kError;
    }
    TH1* up=0;
    TH1* down=0;
    bool isOK = (unc == SFEigen) ? eigenVariation->getEigenvectorVariation(numVariation,up,down) : eigenVariation->getNamedVariation(numVariation,up,down);
    if (!isOK) {
      cerr << "getWeightScaleFactor: Eigenvector object is there but cannot retrieve up and down uncertainty histograms." << endl;
      return Analysis::kError;
    }
    // the 'extrapolation' uncertainty (always a named one) needs a somewhat special treatment
    bool extrapolate = ( unc == SFNamed ) ? eigenVariation->isExtrapolationVariation(numVariation) : false;
    
    double valueUp;
    double valueDown;
    Analysis::CalibrationStatus statusUp   = container->getResult(variables, valueUp,  up,   extrapolate);
    Analysis::CalibrationStatus statusDown = container->getResult(variables, valueDown,down, extrapolate);
    if (statusUp == Analysis::kError || statusDown == Analysis::kError) return Analysis::kError;
 
    // now carry out the rescaling. Protect against unphysical or suspiciously large scale factors
    double variationUp = valueUp - value;
    double variationDown = valueDown - value;
    // First step: from the calibration sample to its reference sample
    if (m_useTopologyRescaling) value = 1.0 + (value - 1.0) * (fracMCref / fracSFref);
    // Second step: from the calibration reference sample to the MC object's reference sample
    if (m_useMCMCSF) value *= (fracSFref / fracEffref);
    // Third step: from the MC object's reference sample to the MC sample itself
    if (m_useTopologyRescaling) value = 1.0 + (value - 1.0) * (fracEffref / fracMCnew);
    // Since all transformations of the scale factor itself are linear, the transformation of the variations is simpler.
    if (m_useTopologyRescaling) {
      double f = (fracMCref / fracMCnew);
      variationUp   *= f;
      variationDown *= f;
    } else if (m_useMCMCSF) {
      double f = (fracSFref / fracEffref);
      variationUp   *= f;
      variationDown *= f;
    }
    valueUp   = value + variationUp;
    valueDown = value + variationDown;
    if (valueUp < 0) {
      valueUp = 0; increaseCounter(indexSF, TagWeight);
    } else if (valueUp > m_maxTagWeight) {
      valueUp = m_maxTagWeight; increaseCounter(indexSF, TagWeight);
    }
    if (valueDown < 0) {
      valueDown = 0; increaseCounter(indexSF, TagWeight);
    } else if (valueDown > m_maxTagWeight) {
      valueDown = m_maxTagWeight; increaseCounter(indexSF, TagWeight);
    }
 
    result.first  = valueUp;
    result.second = valueDown;
    return statusUp;
  } //end eigenvector method
 
  //Proceed with no-eigenvector result
    
  // retrieve the statistical uncertainty if desired
  double stat(0);
  if (unc == Total || unc == Statistical) {
    if (container->getStatUncertainty(variables, stat) == Analysis::kError) {
      cerr << "getWeightScaleFactor: error retrieving Scale factor parameter covariance matrix!" << endl;
      return Analysis::kError;
    }
  }
  Analysis::UncertaintyResult uncertaintyResult(0,0);
  if (unc == Total || unc == Systematic) {
    if (container->getSystUncertainty(variables, uncertaintyResult) == Analysis::kError) {
      cerr << "getWeightScaleFactor: error retrieving Scale factor parameter systematic uncertainty!" << endl;
      return Analysis::kError;
    }
  } else if (unc == Extrapolation) {
    // this uncertainty is special, since it is not normally to be combined into the overall systematic uncertainty
    if (container->getUncertainty("extrapolation", variables, uncertaintyResult) == Analysis::kError)
      cerr << "getWeightScaleFactor: error retrieving Scale factor parameter extrapolation uncertainty!" << endl;
  } else if (unc == TauExtrapolation) {
    // also this uncertainty is special, since it it singles out an uncertainty relevant only for tau "jets",
    // and some care has to be taken not to duplicate or omit uncertainties
    if (container->getUncertainty("extrapolation from charm", variables, uncertaintyResult) == Analysis::kError)
      cerr << "getWeightScaleFactor: error retrieving Scale factor parameter extrapolation uncertainty!" << endl;
  }  
 
  double uncertainty = combinedUncertainty(stat, uncertaintyResult);
 
  // Now carry out the rescaling. Again protect against unphysical or suspiciously large scale factors
  // First step: from the calibration sample to its reference sample
  if (m_useTopologyRescaling) value = 1.0 + (value - 1.0) * (fracMCref / fracSFref);
  // Second step: from the calibration reference sample to the MC object's reference sample
  if (m_useMCMCSF) value *= (fracSFref / fracEffref);
  // Third step: from the MC object's reference sample to the MC sample itself
  if (m_useTopologyRescaling) value = 1.0 + (value - 1.0) * (fracEffref / fracMCnew);
  if (value < 0) {
    value = 0; increaseCounter(indexSF, TagWeight);
  } else if (value > m_maxTagWeight) {
    value = m_maxTagWeight; increaseCounter(indexSF, TagWeight);
  }
  // Since all transformations of the scale factor itself are linear, the transformation of the uncertainty is simpler.
  if (m_useTopologyRescaling) {
    uncertainty *= (fracMCref / fracMCnew);
  } else if (m_useMCMCSF) {
    uncertainty *= (fracSFref / fracEffref);
  }
 
  result.first  = std::max(0., value);
  result.second = uncertainty;
  // "Select" the status code for the actual calibration object (it is subject to more constraints)
  return status;
}

getWeightScaleFactor() [3/3]

Analysis::CalibResult Analysis::CalibrationDataInterfaceROOT::getWeightScaleFactor	(	const CalibrationDataVariables &	variables,
		unsigned int	indexSF,
		unsigned int	indexEff,
		Uncertainty	unc,
		unsigned int	numVariation = 0 )

efficiency scale factor retrieval by index

Definition at line 1523 of file CalibrationDataInterfaceROOT.cxx.

{
    // #2
  // Tag weight fraction scale factor retrieval identifying the requested calibration object by index.
  // The return value is either a (value, uncertainty) or (if eigenvector or named variations are specified)
  // an (up, down) variation pair, and will be a dummy value in case an error occurs.
  // Note that in contrast to the "regular" (non-continuous) case, the computation of the scale factor in
  // general needs the (selection- or even process-specific) MC tag weight fractions, in order to rescale
  // scale factors. This is used to ensure that the tag weight fractions (both in data and in MC) sum up to
  // unity for each given kinematic bin.
  //
  //    variables:     object holding kinematic (and other) information needed to compute the result
  //    indexSF:       index to calibration object
  //    indexEff:      index to MC tag weight
  //    unc:           keyword indicating what uncertainties to evaluate (or whether eigenvector or
  //                   named variations are to be computed)
  //    numVariation:  variation index (in case of eigenvector or named variations)
 
  Analysis::CalibResult result;
  return (getWeightScaleFactor(variables, indexSF, indexEff, unc, numVariation, result) == Analysis::kError) ?
    Analysis::dummyResult : result;
}

increaseCounter()

void Analysis::CalibrationDataInterfaceROOT::increaseCounter ( unsigned int index, OutOfBoundsType oob = Main )

private

Definition at line 1956 of file CalibrationDataInterfaceROOT.cxx.

{
  // Internal method bumping the relevant counter out-of-bounds counter for the specified object.
  //
  //   oob:    further classification of out-of-bounds case
  //   index:  object index
 
  // make sure the vectors are appropriately dimensioned
  if (index >= m_mainCounters.size()) {
    unsigned int minsize = (index == 0) ? 2 : 2*index;
    m_mainCounters.resize(minsize, 0);
    m_etaCounters.resize(minsize, 0);
    m_extrapolatedCounters.resize(minsize, 0);
  }
  switch (oob) {
  case Main:
    m_mainCounters[index]++; break;
  case Eta:
    m_etaCounters[index]++; break;
  case Extrapolated:
  default:
    m_extrapolatedCounters[index]++;
  }
}

initialize()

void Analysis::CalibrationDataInterfaceROOT::initialize	(	const std::string &	jetauthor,
		const std::string &	OP,
		Uncertainty	unc )

initialization for PROOF usage

Definition at line 2481 of file CalibrationDataInterfaceROOT.cxx.

{ 
  // Preload objects necessary so that the input calibration file can be closed. 
  // This functionality is only needed when using PROOF. 
 
  if((!m_fileEff)||(!m_fileSF)) { 
    cerr << "initialize can only be called once per CalibrationDataInterfaceROOT object" << endl; 
    return; 
  } else { 
    cout << "initializing BTagCalibrationDataInterfaceROOT for PROOF with jetAuthor = " << jetauthor  << ", tagger = " << m_taggerName << ", operating point = " << OP << ", uncertainty = " << unc << endl; 
  } 
         
  CalibrationDataVariables BTagVars; 
  BTagVars.jetAuthor = jetauthor; 
  BTagVars.jetPt  = 100000.; //Irrelevant, just has to be valid to retrieve objects 
  BTagVars.jetEta = 1.5; //Irrelevant, just has to be valid to retrieve objects 
         
  for(const auto& flavour : m_flavours){
    std::pair<double, double> BTagCalibResult; 
    BTagCalibResult = getScaleFactor(BTagVars, flavour, OP, unc);
    std::cout << "CalibrationDataInterfaceROOT->initialize : BTagCalibResult " << std::endl;
         
    std::pair<double, double> BTagCalibMCEff; 
    BTagCalibMCEff = getMCEfficiency(BTagVars, flavour, OP, unc); 
    std::cout << "CalibrationDataInterfaceROOT->initialize : BTagCalibMCEff " << std::endl;
  }      
         
  if (m_fileEff != m_fileSF) { 
    m_fileEff->Close(); 
    delete m_fileEff; 
  } 
  m_fileSF->Close(); 
  delete m_fileSF; 
  m_fileEff = 0; //prevents repeat deletion in destructor 
  m_fileSF = 0; //prevents repeat deletion in destructor 
} 

listScaleFactorUncertainties() [1/2]

std::vector< string > Analysis::CalibrationDataInterfaceROOT::listScaleFactorUncertainties	(	const std::string &	author,
		const std::string &	label,
		const std::string &	OP,
		bool	named = false )

retrieve the list of "uncertainties" relevant to the calibration object.

A few uncertainty names are predetermined: "result", "comment", "statistics", "systematics". Other sources of systematic uncertainty may be added. Note that the "systematics" source does not give access to correlations between bins. If the 'named' argument is true, the list does not include all uncertainties but only those excluded from the eigenvector construction (this option is only relevant if eigenvector use has been switched on to begin with). In this case the order of the uncertainties listed is important since it indicates the index by which the given named uncertainty is identified.

Definition at line 1984 of file CalibrationDataInterfaceROOT.cxx.

{
  // Retrieve the sources of uncertainty relevant for the given scale factor calibration object,
  // identifying the object by name.
  //
  //    author:  jet collection name
  //    label:   jet flavour label
  //    OP:      tagger working point
  //    named:   if false, an unsorted list of sources of uncertainties will be returned.
  //             if true, only 'named' uncertainties will be returned, and the position in
  //             the vector that is the return value determines the 'numVariation' index
  //             that is to be used if named variations are to be retrieved.
 
  unsigned int index;
  if (! retrieveCalibrationIndex (label, OP, author, true, index)) {
    // Return a dummy result if the object is not found
    cerr << "listScaleFactorUncertainties: unable to find SF calibration for object " << fullName(author, OP, label, true) << endl;
    std::vector<string> dummy;
    return dummy;
  }
  return listScaleFactorUncertainties(index, label, named);
}

listScaleFactorUncertainties() [2/2]

std::vector< string > Analysis::CalibrationDataInterfaceROOT::listScaleFactorUncertainties	(	unsigned int	index,
		const std::string &	flavour,
		bool	named = false )

retrieve the list of "uncertainties" relevant to the calibration object.

A few uncertainty names are predetermined: "result", "comment", "statistics", "systematics". Other sources of systematic uncertainty may be added. Note that the "systematics" source does not give access to correlations between bins. If the 'named' argument is true, the list does not include all uncertainties but only those excluded from the eigenvector construction (this option is only relevant if eigenvector use has been switched on to begin with). In this case the order of the uncertainties listed is important since it indicates the index by which the given named uncertainty is identified.

Definition at line 2012 of file CalibrationDataInterfaceROOT.cxx.

{
  // Note: this method already works on a per-flavour basis, so passing flavour in is a simple addition
  // Method is called primarily from within the BTaggingEfficiencyTool - W.L.
  // Retrieve the sources of uncertainty relevant for the given scale factor calibration object,
  // identifying the object by index.
  //
  //    index:   index to scale factor calibration object
  //    named:   if false, an unsorted list of sources of uncertainties will be returned.
  //             if true, only 'named' uncertainties will be returned, and the position in
  //             the vector that is the return value determines the 'numVariation' index
  //             that is to be used if named variations are to be retrieved.
 
  std::vector<string> dummy;
  CalibrationDataContainer* container = m_objects[index];
 
  if (container) {
    if (named) {
      // Find out which uncertainties are excluded from eigenvector construction
      if (! m_runEigenVectorMethod) return dummy;
      std::shared_ptr<CalibrationDataEigenVariations> eigenVariation=m_eigenVariationsMap.at(container);
      if (m_EVStrategy == Analysis::Uncertainty::SFEigen){
        std::vector<string> unordered = eigenVariation->listNamedVariations(); // this is for the regular EV
        std::vector<string> ordered(unordered.size());
        for (unsigned int i = 0; i < unordered.size(); ++i) {
          ordered[eigenVariation->getNamedVariationIndex(unordered[i])] = unordered[i];
        }
        return ordered;
      } else if (m_EVStrategy == Analysis::Uncertainty::SFGlobalEigen){
        // here we want to get the named uncertainties from the global eigenvariations flavour container specifically...
        std::shared_ptr<CalibrationDataGlobalEigenVariations> GEV = std::dynamic_pointer_cast<CalibrationDataGlobalEigenVariations>(eigenVariation);
        std::vector<std::string> unordered = GEV->listNamedVariations(flavour);
        std::vector<std::string> ordered(unordered.size()); // ordered by the NAMED VARIATION (internal) ORDERING
        for (unsigned int i = 0; i < unordered.size(); ++i) {
          ordered[GEV->getNamedVariationIndex(unordered[i], flavour)] = unordered[i]; 
        }
        return ordered;
      }
    }
    return container->listUncertainties(); // return this if not named 
  }
 
  return dummy;
}

nameFromIndex()

std::string Analysis::CalibrationDataInterfaceROOT::nameFromIndex

(

unsigned int

index

)

const

Retrieve the name of the calibration object (container) given its index.

Definition at line 1942 of file CalibrationDataInterfaceROOT.cxx.

{
  // Return the object name corresponding to the given index.
 
  for (std::map<std::string, unsigned int>::const_iterator it = m_objectIndices.begin();
       it != m_objectIndices.end(); ++it)
    if (it->second == index) return it->first;
 
  // This should never happen..
  return string("");
}

retrieveCalibrationIndex()

bool Analysis::CalibrationDataInterfaceROOT::retrieveCalibrationIndex	(	const std::string &	label,
		const std::string &	OP,
		const std::string &	author,
		bool	isSF,
		unsigned int &	index,
		unsigned int	mapIndex = 0 )

Retrieve the index of the calibration object (container) starting from the label and operating point.

The return value will be false if the requested object cannot be found.

Definition at line 688 of file CalibrationDataInterfaceROOT.cxx.

{
  // Retrieve the integer index corresponding to a given combination of 
  // flavour label / tagger / working point / jet collection name, and separately
  // for calibration scale factors and MC efficiencies (all these ingredients are needed
  // to specify fully the calibration object).
  // In fact this method will also trigger the retrieval of the object itself, if not already
  // done, and will cache it internally. The absence of the requested calibration object will
  // be flagged by a false return value.
  // This method is used internally but should also be called by users in order to exploit the
  // "code speed-up" features documented above.
  //
  //    label:     jet flavour label
  //    OP:        tagger working point
  //    author:    jet collection name
  //    isSF:      set to true (false) for scale factors (MC efficiencies)
  //    index:     resulting index (meaningful only for a 'true' function return value)
  //    mapIndex:  index to the MC efficiency map to be used
 
  index = 0;
 
  // construct the full name from the label, operating point, SF/Eff choice;
  // then look up this full name
  string name = fullName(author, OP, label, isSF, mapIndex);
  std::map<string, unsigned int>::const_iterator it = m_objectIndices.find(name);
  if (it == m_objectIndices.end()) {
    // If no container is found, attempt to retrieve it here (this is so that users won't
    // have to call the named scale factor etc. methods once just to retrieve the container).
    string flavour = (label == "N/A") ? "Light" : label;
    string cntname = getContainername(flavour, isSF, mapIndex);
    if (m_verbose) std::cout << "CalibrationDataInterfaceROOT->retrieveCalibrationIndex :  container name is " << cntname  << std::endl;
    retrieveContainer(flavour, OP, author, cntname, isSF, m_verbose); // Only call this if you want to retrieve a currently not available container
    it = m_objectIndices.find(name);
    if (it == m_objectIndices.end()) return false;
  } else {
    if (m_verbose) std::cout << "CalibrationDataInterfaceROOT->retrieveCalibrationIndex : container " << name << " already cached! " << std::endl;
  }
 
  index = it->second;
  return true;
}

retrieveContainer()

CalibrationDataContainer * Analysis::CalibrationDataInterfaceROOT::retrieveContainer	(	const std::string &	label,
		const std::string &	OP,
		const std::string &	author,
		const std::string &	cntname,
		bool	isSF,
		bool	doPrint = true )

utility function taking care of object retrieval

Definition at line 2520 of file CalibrationDataInterfaceROOT.cxx.

{
  // Attempt to retrieve the given container from file. Note that also the corresponding
  // "hadronisation" reference is retrieved (if possible and not yet done).
  //
  //     dir:     name of the directory containing the requested container
  //     cntname: name of the requested container itself (not including the full path)
  //     isSF:    set to false (true) if the object is to be retrieved from the MC efficiencies
  //              file (the calibration scale factor file). Note that it is assumed that scale
  //              factor objects will always be retrieved from the calibration scale factor file.
  //     doPrint: if true, print out some basic information about the successfully retrieved container
  //              (note that this is typically steered by the m_verbose setting;
  //               only for the retrieval of the maps used for MC/MC SF calculations, this printout is always switched off)
 
  string dir = m_taggerName + "/" + getAlias(author) + "/" + OP + "/" + label;
  // construct the full object name
  string name = dir + "/" + cntname;
 
  // If the object cannot be found, then each call will result in a new attempt to
  // retrieve the object from the ROOT file. Hopefully this will not happen too often...
  unsigned int idx = m_objectIndices[name] = m_objects.size();
  // CalibrationDataContainer* cnt =
  //   dynamic_cast<CalibrationDataContainer*>((isSF ? m_fileSF : m_fileEff) ->Get(name.c_str()));
  CalibrationDataContainer* cnt;
  (isSF ? m_fileSF : m_fileEff)->GetObject(name.c_str(), cnt);
  // If the requested object is a MC efficiency container and is not found, make a second attempt
  // to retrieve it from the calibration scale factor file. This will avoid the need to duplicate
  // efficiency containers so that the MC efficiency file needs to store only those containers
  // not already present in the calibration scale factor file. Of course this is meaningful only
  // if separate files are used to begin with.
  if (!isSF && !cnt && m_fileSF != m_fileEff) m_fileSF->GetObject(name.c_str(), cnt);
  m_objects.push_back(cnt);
  if (!cnt) {
    cerr << "btag Calib: retrieveContainer: failed to retrieve container named " << name.c_str() << " from file" << endl;
    return 0;
  }
 
  // For successfully retrieved containers, also print some more information (implemented on user request)
  if (doPrint) {
    cout << "CalibrationDataInterface: retrieved container " << name << " (with comment: '" << cnt->getComment() << "' and hadronisation setting '" << cnt->getHadronisation() << "')" << endl;
  }
 
 
  // If the requested object is a MC efficiency container, make sure to retrieve the corresponding
  // calibration scale factor container first (a feature first thought to be necessary, erroneously,
  // but left in since this ordering should not hurt in any case).
  if (m_refMap.find(dir) == m_refMap.end()) {
    if (isSF) {
      // Retrieve the mapping objects from both files and merge their information using the 'helper' class.
      // The map resulting from this is used to retrieve the information required to compute MC/MC scale factors.
      string hadronisationRefs(dir + "/MChadronisation_ref");
      TMap* mapSF = 0; m_fileSF->GetObject(hadronisationRefs.c_str(), mapSF);
      TMap* mapEff = 0; if (m_fileEff != m_fileSF) m_fileEff->GetObject(hadronisationRefs.c_str(), mapEff);
      m_refMap[dir] = new HadronisationReferenceHelper(mapSF, mapEff);
      delete mapSF;
      delete mapEff;
    } else {
      string SFCalibName = getContainername(getBasename(dir), true);
      if (m_objectIndices.find(SFCalibName) == m_objectIndices.end()) retrieveContainer(label, OP, author, SFCalibName, true, doPrint);
    }
  }
 
  // Attempt to find the corresponding hadronisation reference container needed for the application of
  // MC/MC scale factors.
  if (idx+1 > m_hadronisationReference.size()) m_hadronisationReference.resize(idx+1, -1);
  m_hadronisationReference[idx] = -1;
  string spec = cnt->getHadronisation();
  if (spec != "") {
    std::map<string, HadronisationReferenceHelper*>::const_iterator mapit = m_refMap.find(dir);
    if (mapit != m_refMap.end()) {
      string ref;
      if (mapit->second->getReference(spec, ref)) {
        // Retrieve the hadronisation reference if not already done. Note that the "isSF" is left unchanged:
        // this allows to retrieve the reference from the same file as the scale factor object. An exception
        // is the reference for the calibration scale factor object, which should always be obtained from
        // the scale factor file. 
        // An efficiency container can be its own hadronisation reference (this is not "protected" against).
        string refname(dir + "/" + ref);
        std::map<string, unsigned int>::const_iterator it = m_objectIndices.find(refname);
        // If the reference cannot be found, assume that it hasn't yet been retrieved so attempt it now.
        if (it == m_objectIndices.end()) {
          // Omit the printout of container information here (the idea being that showing MC/MC SF information would confuse rather than help)
          retrieveContainer(label, OP, author, ref, isSF, false); it = m_objectIndices.find(refname);
        }
        m_hadronisationReference[idx] = it->second;
      }
    } else if (m_useMCMCSF) {
      cerr << "btag Calib: retrieveContainer: MC hadronisation reference map not found -- this should not happen!" << endl;
    }
  }
  if (m_hadronisationReference[idx] == -1 || ! m_objects[m_hadronisationReference[idx]]){
    // Not being able to construct the MC/MC scale factors will lead to a potential bias.
    // However, this is not considered sufficiently severe that we will flag it as an error.
    if (m_useMCMCSF){
      cerr << "btag Calib: retrieveContainer: warning: unable to apply MC/MC scale factors for container " << name << " with hadronisation reference = '" << spec << "'" << endl;
    }
  }
 
  // Initialize the Eigenvector variation object corresponding to this object, if applicable. Notes:
  // - the dual use of "isSF" (both referring to the file and to the object, see above) requires another protection here
  // - the constructor's second argument is used to determine whether to exclude a pre-determined set of uncertainties from the EV decomposition
  //
  // We also want to separate behavior between SFEigen and SFGlobalEigen systematic strategies
  // The former requires a CalibrationDataEigenVariations object to be made per flavour.
  // The latter combines all corresponding flavours, so once it's been made for a single flavour, it's cached under all the corresponding "flavour containers"
  // simulataneously in m_eigenVariationsMap, and is checked for on each subsequent call to this method.
  if (m_runEigenVectorMethod && isSF && name.find("_SF") != string::npos) {
    CalibrationDataHistogramContainer* histoContainer=dynamic_cast<CalibrationDataHistogramContainer*>(cnt);
    if (histoContainer==0) {
      cerr << "Could not cast Container to a HistogramContainer. " << endl;
      return 0;
    }
    if (m_EVStrategy == Analysis::Uncertainty::SFEigen){
      std::shared_ptr<CalibrationDataEigenVariations> newEigenVariation(new CalibrationDataEigenVariations(m_filenameSF, m_taggerName, OP, author, histoContainer, m_useRecommendedEVExclusions));
      newEigenVariation->setVerbose(m_verbose);
 
      // At this point we may also want to reduce the number of eigenvector variations.
      // The choices are stored with the container object; but first we need to know what flavour we are dealing with.
      string flavour = dir.substr(dir.find_last_of("/")+1);
 
      for (const auto & entry : m_excludeFromCovMatrix[flavour]) {
        newEigenVariation->excludeNamedUncertainty(entry, cnt);
      }
      newEigenVariation->initialize();
      int to_retain = histoContainer->getEigenvectorReduction(m_EVReductions[flavour]); // returns the number of eigenvariations to retain as per the EV reduction strategy
      if (to_retain > -1) {
        if (m_verbose) cout << "btag Calib: reducing number of eigenvector variations for flavour " << flavour << " to " << to_retain << endl;
        // The merged variations will end up as the first entry in the specified list, i.e., as the last of the variations to be "retained"
        newEigenVariation->mergeVariationsFrom(size_t(to_retain-1)); // All variations stored with indices larger than this are merged
      } else if (m_EVReductions[flavour] != Loose) {
        cerr << "btag Calib: unable to retrieve eigenvector reduction information for flavour " << flavour << " and scheme " << m_EVReductions[flavour] << "; not applying any reduction" << endl;
      }
      m_eigenVariationsMap[cnt]=std::move(newEigenVariation);

    } else if (m_EVStrategy == Analysis::Uncertainty::SFGlobalEigen) {
      std::map<const CalibrationDataContainer*, std::shared_ptr<CalibrationDataEigenVariations> >::iterator evit = m_eigenVariationsMap.find(cnt);
      // The global implementation internally combines all the "flavour containers" (containers that correspond to each other, only with different flavours)
      // But the CalibrationDataInterfaceROOT object doesn't need to know that, so we want to get all the flavour containers in one go here
      // and map them (with m_eigenVariationsMap) to the same CalibrationDataGlobalEigenVariations pointer.
      // Then, in methods like getScaleFactor, we call the virtual methods which will give the proper result e.g. if you want the SF for a b-jet, it'll call the 
 
      if (evit == m_eigenVariationsMap.end()){
        // now to see if it's completely empty or not
        if (m_eigenVariationsMap.empty()){
          std::shared_ptr<CalibrationDataGlobalEigenVariations> newEigenVariation(new CalibrationDataGlobalEigenVariations(m_filenameSF, m_taggerName, OP, author, m_flavours, histoContainer, m_useRecommendedEVExclusions));
          for (const auto & entry : m_excludeFromCovMatrix[label]) {
            newEigenVariation->excludeNamedUncertainty(entry, label); // <---- custom exclude named uncertainties method for global variations
          }
 
          newEigenVariation->initialize();
          
          // flavour loop to get the flavour reduction schemes and apply them
          for (std::string& flavour : m_flavours){
            int to_retain = histoContainer->getEigenvectorReduction(m_EVReductions[flavour]); // returns the number of eigenvariations to retain as per the EV reduction strategy
            if (to_retain > -1) {
              if (m_verbose) cout << "btag Calib: reducing number of eigenvector variations for flavour " << flavour << " to " << to_retain << endl;
              // The merged variations will end up as the first entry in the specified list, i.e., as the last of the variations to be "retained"
              newEigenVariation->mergeVariationsFrom(size_t(to_retain-1), flavour); // All variations stored with indices larger than this are merged
            } else if (m_EVReductions[flavour] != Loose) {
              cerr << "btag Calib: unable to retrieve eigenvector reduction information for flavour " << flavour << " and scheme " << m_EVReductions[flavour] << "; not applying any reduction" << endl;
            }
          }
 
          m_eigenVariationsMap.insert({cnt, newEigenVariation});
        } else {
          // Need to point to the CDGEV object four times in the m_eigenVariationsMap to appease the CDIROOT backend design...
          // Ok, turns out I can't retrieve the containers from CDGEV and insert them directly, because I'd have to use the containers directly instead..
          // So the strategy is to just take the CGEV objects that are already in the map, and mpa the present container to it
          std::shared_ptr<CalibrationDataEigenVariations> previous_eigenvariation = m_eigenVariationsMap.begin()->second;
          m_eigenVariationsMap.insert({cnt, previous_eigenvariation});
        }
 
      } else {
        std::cout << "CalibrationDataInterfaceROOT->retrieveContainer : the CDGEV object for " << name << " already exists! " << std::endl;
      }
    }
  }
 
  return cnt;
}

runEigenVectorRecomposition() [1/2]

Analysis::CalibrationStatus Analysis::CalibrationDataInterfaceROOT::runEigenVectorRecomposition	(	const std::string &	author,
		const std::string &	label,
		const std::string &	OP,
		unsigned int	mapindex = 0 )

run EigenVector Recomposition method

Definition at line 2204 of file CalibrationDataInterfaceROOT.cxx.

                                                           {
  // run eigen vector recomposition method. If success, stored the retrieved coefficient map
  // in m_coefficientMap and return success. Otherwise return error and keep m_coefficientMap
  // untouched.
  //    author:  jet collection name
  //    label:   jet flavour label
  //    OP:      tagger working point
  //    mapIndex:     index to the MC efficiency map to be used. Should be 0?
  // Todo: What is mapindex?
  // Todo: Check the way xAODBTaggingTool initialize CDI. Check if that is the as how we are initialize CDI. 
  if(!m_runEigenVectorMethod) {
    cerr << "runEigenVectorRecomposition: Recomposition need to be ran with CalibrationDataInterfaceRoot initialized in eigenvector mode" << endl;
    return Analysis::kError;
  }
  
  unsigned int indexSF;
  if (! retrieveCalibrationIndex (label, OP, author, true, indexSF, mapIndex)) {
    cerr << "runEigenVectorRecomposition: unable to find SF calibration for object "
     << fullName(author, OP, label, true) << endl;
    return Analysis::kError;
  }
 
  return runEigenVectorRecomposition (label, indexSF);
}

runEigenVectorRecomposition() [2/2]

Analysis::CalibrationStatus Analysis::CalibrationDataInterfaceROOT::runEigenVectorRecomposition	(	const std::string &	label,
		unsigned int	mapindex = 0 )

Definition at line 2233 of file CalibrationDataInterfaceROOT.cxx.

                                                          {
  // run eigen vector recomposition method. If success, stored the retrieved coefficient map
  // in m_coefficientMap and return success. Otherwise return error and keep m_coefficientMap
  // untouched.
  //    label:   jet flavour label  
  //    indexSF:       index to scale factor calibration object
  CalibrationDataContainer* container = m_objects[indexSF];
  if (! container) {
    cerr << "runEigenVectorRecomposition: error retrieving container!" << endl;
    return Analysis::kError;
  }
 
  // Retrieve eigenvariation
  std::shared_ptr<CalibrationDataEigenVariations> eigenVariation;
  try {
    eigenVariation = m_eigenVariationsMap.at(container);
  } catch (const std::out_of_range&) {
    cerr << "runEigenVectorRecomposition: Could not retrieve eigenvector variation, while it should have been there." << endl;
    return Analysis::kError;
  }
  // Doing eigenvector recomposition
  std::map<std::string, std::map<std::string, float>> coefficientMap;
  if(!eigenVariation->EigenVectorRecomposition(label, coefficientMap))
    return Analysis::kError;
 
  m_coefficientMap = std::move(coefficientMap);
  return Analysis::kSuccess;
}

setEffCalibrationNames()

void Analysis::CalibrationDataInterfaceBase::setEffCalibrationNames ( const std::map< std::string, std::vector< std::string > > & names )

inherited

Definition at line 63 of file CalibrationDataInterfaceBase.cxx.

{
  // Set the MC efficiency names.
 
  m_calibrationEffNames = names;
}

setSFCalibrationNames()

void Analysis::CalibrationDataInterfaceBase::setSFCalibrationNames ( const std::map< std::string, std::string > & names )

inherited

Definition at line 87 of file CalibrationDataInterfaceBase.cxx.

{
  // Set the efficiency scale factor calibration names.
 
  m_calibrationSFNames = names;
}

SFCalibrationName()

const std::string & Analysis::CalibrationDataInterfaceBase::SFCalibrationName ( const std::string & flavour ) const

inherited

Definition at line 72 of file CalibrationDataInterfaceBase.cxx.

{
  // Return the efficiency scale factor calibration name for the given flavour.
  // Note that no check is performed on the validity of the flavour.
 
  try {
    return m_calibrationSFNames.at(flavour);
  }
  catch (const std::out_of_range& e) {
    std::cerr << "SFCalibrationName: flavour '" << flavour << "' is not known." << std::endl;
    throw e;
  }
}

Member Data Documentation

m_absEtaStrategy

OutOfBoundsStrategy Analysis::CalibrationDataInterfaceROOT::m_absEtaStrategy {}

private

Definition at line 474 of file CalibrationDataInterfaceROOT.h.

{};

m_aliases

std::map<std::string, std::string> Analysis::CalibrationDataInterfaceROOT::m_aliases

private

Do not attempt to persistify (PROOF)

jet author aliases (there is no single CalibrationBroker object here to take care of this, so we do it in this class)

Definition at line 405 of file CalibrationDataInterfaceROOT.h.

m_calibrationEffNames

std::map<std::string, std::vector<std::string> > Analysis::CalibrationDataInterfaceBase::m_calibrationEffNames

privateinherited

this simply collects the per-flavour properties.

Definition at line 72 of file CalibrationDataInterfaceBase.h.

m_calibrationSFNames

std::map<std::string, std::string> Analysis::CalibrationDataInterfaceBase::m_calibrationSFNames

privateinherited

Definition at line 73 of file CalibrationDataInterfaceBase.h.

m_checkedWeightScaleFactors

std::vector<std::pair<unsigned int, unsigned int> > Analysis::CalibrationDataInterfaceROOT::m_checkedWeightScaleFactors

private

Definition at line 488 of file CalibrationDataInterfaceROOT.h.

m_coefficientMap

std::map<std::string, std::map<std::string, float> > Analysis::CalibrationDataInterfaceROOT::m_coefficientMap

private

Definition at line 466 of file CalibrationDataInterfaceROOT.h.

m_eigenVariationsMap

std::map<const CalibrationDataContainer*, std::shared_ptr<CalibrationDataEigenVariations> > Analysis::CalibrationDataInterfaceROOT::m_eigenVariationsMap

private

store the eigenvector class and associate to its CalibrationDataContainer

Definition at line 426 of file CalibrationDataInterfaceROOT.h.

m_etaCounters

std::vector<unsigned int> Analysis::CalibrationDataInterfaceROOT::m_etaCounters

private

counters for flagging out-of-bound cases

Definition at line 480 of file CalibrationDataInterfaceROOT.h.

m_EVReductions

std::map<std::string, Analysis::EVReductionStrategy> Analysis::CalibrationDataInterfaceROOT::m_EVReductions

private

Eigenvector reduction strategy (per flavour)

Definition at line 433 of file CalibrationDataInterfaceROOT.h.

m_EVStrategy

Uncertainty Analysis::CalibrationDataInterfaceROOT::m_EVStrategy {}

private

Definition at line 430 of file CalibrationDataInterfaceROOT.h.

{}; // <--- In addition, specify whether to use the global eigenvariations method

m_excludeFromCovMatrix

std::map<std::string, std::vector<std::string> > Analysis::CalibrationDataInterfaceROOT::m_excludeFromCovMatrix

private

store the uncertainties which should be excluded from building the full covariance matrix

Definition at line 436 of file CalibrationDataInterfaceROOT.h.

m_extrapolatedCounters

std::vector<unsigned int> Analysis::CalibrationDataInterfaceROOT::m_extrapolatedCounters

private

Definition at line 482 of file CalibrationDataInterfaceROOT.h.

m_fileEff

TFile* Analysis::CalibrationDataInterfaceROOT::m_fileEff {}

private

pointer to the TFile object providing access to the calibrations

Definition at line 400 of file CalibrationDataInterfaceROOT.h.

{}; 

m_filenameEff

std::string Analysis::CalibrationDataInterfaceROOT::m_filenameEff

private

Definition at line 418 of file CalibrationDataInterfaceROOT.h.

m_filenameSF

std::string Analysis::CalibrationDataInterfaceROOT::m_filenameSF

private

in addition, store also the filenames themselves (needed for the copy constructor)

Definition at line 417 of file CalibrationDataInterfaceROOT.h.

m_fileSF

TFile* Analysis::CalibrationDataInterfaceROOT::m_fileSF {}

private

Do not attempt to persistify (PROOF)

Definition at line 401 of file CalibrationDataInterfaceROOT.h.

{}; 

m_flavours

std::vector<std::string> Analysis::CalibrationDataInterfaceROOT::m_flavours

private

Definition at line 419 of file CalibrationDataInterfaceROOT.h.

m_hadronisationReference

std::vector<int> Analysis::CalibrationDataInterfaceROOT::m_hadronisationReference

private

store the 'hadronisation' reference for each object (-1 means no reference found)

Definition at line 461 of file CalibrationDataInterfaceROOT.h.

m_mainCounters

std::vector<unsigned int> Analysis::CalibrationDataInterfaceROOT::m_mainCounters

private

Definition at line 481 of file CalibrationDataInterfaceROOT.h.

m_maxAbsEta

double Analysis::CalibrationDataInterfaceROOT::m_maxAbsEta {}

private

|eta| bounds and strategy for dealing with out-of-bounds conditions

Definition at line 473 of file CalibrationDataInterfaceROOT.h.

{};

m_maxTagWeight

double Analysis::CalibrationDataInterfaceROOT::m_maxTagWeight {}

private

Definition at line 489 of file CalibrationDataInterfaceROOT.h.

{};

m_objectIndices

std::map<std::string, unsigned int> Analysis::CalibrationDataInterfaceROOT::m_objectIndices

private

Definition at line 411 of file CalibrationDataInterfaceROOT.h.

m_objects

std::vector<CalibrationDataContainer*> Analysis::CalibrationDataInterfaceROOT::m_objects

private

cache the objects themselves (so that the user will not have to delete them after each call etc.).

The caching is done so that objects can be retrieved by number as well as by (OP, flavour, calibration name) combination.

Definition at line 410 of file CalibrationDataInterfaceROOT.h.

m_otherStrategy

OutOfBoundsStrategy Analysis::CalibrationDataInterfaceROOT::m_otherStrategy {}

private

Definition at line 475 of file CalibrationDataInterfaceROOT.h.

{};

m_refMap

std::map<std::string, HadronisationReferenceHelper*> Analysis::CalibrationDataInterfaceROOT::m_refMap

private

the following maps (one for each directory) specify the name of the container serving as the 'hadronisation' reference for each object

Definition at line 459 of file CalibrationDataInterfaceROOT.h.

m_runEigenVectorMethod

bool Analysis::CalibrationDataInterfaceROOT::m_runEigenVectorMethod {}

private

decide whether to run the eigenvector method or not

Definition at line 429 of file CalibrationDataInterfaceROOT.h.

{};

m_taggerName

std::string Analysis::CalibrationDataInterfaceBase::m_taggerName

protectedinherited

tagging algorithm name

Definition at line 94 of file CalibrationDataInterfaceBase.h.

m_useMCMCSF

bool Analysis::CalibrationDataInterfaceROOT::m_useMCMCSF {}

private

specify whether or not to use MC/MC (hadronisation) scale factors (the fact that this is steerable is intended to be temporary only)

Definition at line 452 of file CalibrationDataInterfaceROOT.h.

{};

m_useRecommendedEVExclusions

bool Analysis::CalibrationDataInterfaceROOT::m_useRecommendedEVExclusions {}

private

if true, exclude pre-recommended lists of uncertainties from the covariance matrix building, in addition to the above user specified lists

Definition at line 441 of file CalibrationDataInterfaceROOT.h.

{};

m_useTopologyRescaling

bool Analysis::CalibrationDataInterfaceROOT::m_useTopologyRescaling {}

private

specify whether or not to use MC/MC (topology) scale factors (also this steering option may be removed)

Definition at line 455 of file CalibrationDataInterfaceROOT.h.

{};

m_verbose

bool Analysis::CalibrationDataInterfaceROOT::m_verbose {}

private

if true, allow also for some informational (and not only error/warning) messages

Definition at line 444 of file CalibrationDataInterfaceROOT.h.

{};

---

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

• CalibrationDataInterfaceROOT.h
• CalibrationDataInterfaceROOT.cxx