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

#ifndef _JERPROVIDER_
#define _JERPROVIDER_

#include "TNamed.h"
#include "TFile.h"
#include "TF1.h"
#include <TRandom.h>
#include <iostream>
//#include <stream.h>
#include <cmath>
#include <map>

using std::cout;
using std::endl;
using std::map;

class JERProvider : public TNamed
{

 public: 

  // Constructor, destructor
  JERProvider() { }
  JERProvider(std::string CollectionName, std::string MethodName = "Truth", std::string FileName="JERProviderPlots.root");
  virtual ~JERProvider();
  
  //initialization
  void init();

  // Read the parametrization
  TF1* getParam(double y);

  // Read sigma
  float getSigma(double pt, double y);

  // Read the uncertainty
  float getUncert(double pt, double y);

  // inputCollection 
  void setInputCollection(std::string CollectionName, std::string MethodName);
  
#ifdef JER_STANDALONE
        ClassDef(JERProvider,1);
#endif

 private:
  
  static const int m_nParam = 6;
  //static const double m_GeV = 1000.;
  static const int m_nY = 6;

  // strings
  std::string m_collectionName;
  std::string m_methodName;
  std::string m_fileName;
  
  //initialization flag
  bool m_isInit;

  // Input file
  TFile* m_inputFile; 

  // Pointer to the correct graph
  map<int,TF1*> m_jerFunc;

  // Parameters and their errors
  map<int,float> m_param;

  // Uncertainty from Data/MC comparison
  float m_uncert;

  // Offset from Data/MC comparison
  float m_offset;

};

#endif

---