/afs/hep.man.ac.uk/u/markowen/ATLAS/SFrameProof/2011EPS/topUtils/topUtils/CalibrationDataInterfaceBase.h

// CalibrationDataInterfaceTool.h, (c) ATLAS Detector software

// 
// for some details on weighting event using the SF, see
// https://twiki.cern.ch/twiki/bin/view/AtlasProtected/TopCommonObjects#B_tagging
// http://indico.cern.ch/getFile.py/access?contribId=4&resId=0&materialId=slides&confId=106205
// http://cdsweb.cern.ch/record/1269912
// 

#ifndef ANALYSISCALIBRATIONDATAINTERFACEBASE_H
#define ANALYSISCALIBRATIONDATAINTERFACEBASE_H

#include "Rtypes.h"
#include "TMatrixDSym.h"
#include <vector>
#include <string>
#include <map>
#include <utility>
#include <algorithm>

class TF1;

namespace Analysis 
{

  const double CalibZERO = 1.e-6;
  const double dummyValue = -1.;
  typedef std::pair<double,double> CalibResult;
  // The following result is returned upon error
  const CalibResult dummyResult(dummyValue,dummyValue);
  //  const CalibResult dummyResult(std::make_pair(dummyValue,dummyValue));

  class CalibrationDataInterfaceBase {
    public:
      CalibrationDataInterfaceBase();

      virtual ~CalibrationDataInterfaceBase() = 0;
      
      const std::string& EffCalibrationName(const std::string& flavour) const;
      void setEffCalibrationNames(const std::map<std::string, std::string>& names);
      
      const std::string& SFCalibrationName(const std::string& flavour) const;
      void setSFCalibrationNames(const std::map<std::string, std::string>& names);
      
      enum variableType { kEta, kAbsEta, kPt };

      const std::vector<int>& functionArguments(const std::string& object) const;

      bool existFunctionArguments(const std::string& object) const;

    private:

      mutable std::map<std::string, std::vector<int> > m_functionArguments;

      mutable std::map<std::string, std::string> m_calibrationEffNames;
      mutable std::map<std::string, std::string> m_calibrationSFNames;
      
      double parametricVariance(TF1* fct, const TMatrixDSym* cov, const Double_t* x) const;

    protected:

      std::string getBasename (const std::string& OP, const std::string& flavour,
                               const std::string& extra = "", bool SF = true) const;

      CalibResult getResults (const Double_t* x, const Double_t* xsyst,
                              TF1* fct, TF1* fctsys,
                              const TMatrixDSym* cov) const;

      bool computeFunctionArguments(const std::string& object, const char* ctitle) const;

      std::string m_taggerName;

      std::vector<std::string> m_operatingPoints;

    }; 
} // end of namespace

#endif 

---