/Users/jmonk/Physics/ForIA/ForIA/AnalysisTools/JERProvider.hh

Go to the documentation of this file.

#ifndef _JERPROVIDER_
#define _JERPROVIDER_

#include "TNamed.h"
#include "TFile.h"
#include "TF1.h"
#include "TGraphErrors.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();
  
  // filename -- only allowed before initialisation
  void setFileName(std::string fileName);

  //initialization
  void init();

  // Read the parametrization
  TF1* getParam(double y);
 
  // Read the offset from data/MC
  TGraphErrors* getOffset(double y);
  
  // Read the error from data/MC
  TGraphErrors* getSyst(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;

  // Pointer to the correct graph
  map<int,TGraphErrors*> m_jerOffset;
  
  // Pointer to the correct graph
  map<int,TGraphErrors*> m_jerSyst;
  
  // Uncertainty from Data/MC comparison
  float m_uncert;

  // Offset from Data/MC comparison
  float m_offset;

  // Systematics
  float m_syst;


};

#endif