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

// copied from Reconstruction/Jet/JetUncertainties/tags/JetUncertainties-00-03-04
#ifndef _JESUNCERTAINTYPROVIDER_
#define _JESUNCERTAINTYPROVIDER_

#include "TNamed.h"
#include "TFile.h"
#include "TH2D.h"
#include <iostream>
#include <string>
#include <sstream>
#include <cmath>
#include <map>
using std::cout;
using std::map;

class JESUncertaintyProvider : public TNamed
{

public:

  enum Components { 
    // Basic components
    CALORIMETER = 1,                  // Calorimeter uncertainty from single particle propagation
    NOISETHRESHOLDS = 2,              // Use topocluster noise thresholds from data
    PERUGIATUNE = 4,                  // Perugia 2010 Pythia tune 
    ALPGENHERWIGJIMMY = 8,            // Alpgen+Herwig+Jimmy
    ETAINTERCALIBRATION = 16,         // Uncertainty due to intercalibration (endcap wrt central region)
    CLOSURE = 32,                     // Non-closure of numerical inversion constants
    PILEUP = 64,                      // Uncertainty due to in-time pile-up
    
    // Suggested combinations 
    NOCLOSURE = 95,                   // Allow users to include their own non-closure
    NOPILEUP = 63,                    // For data with no pile-up
    NOPILEUPNOCLOSURE = 31,           // Data with no pile-up, no non-closure
    ALL = 127                         // Everything
  }; 
  

  // Constructors, destructor
  JESUncertaintyProvider(std::string CollectionName="AntiKt6TopoJetsEM", std::string FileName="JESUncertainty.root");
   ~JESUncertaintyProvider();
  
  // Read the uncertainties off of the graph - By default, no Pile-Up term!
  double getRelUncert(double pT, double Eta, Components UncertComps = JESUncertaintyProvider::NOPILEUP, unsigned int nVtx=1);
  double getAbsUncert(double pT, double Eta, Components UncertComps = JESUncertaintyProvider::NOPILEUP, unsigned int nVtx=1);
  
  // Pass a pointer to the correct graph
  //TH2D* getUncGraphCopy(Components UncertComps = JESUncertaintyProvider::NOPILEUP, unsigned int nVtx=1);

  // Initialize the provider
  virtual void init();

  // Open input file
  virtual TFile* openInputFile(std::string FileName = "JESUncertainty.root");

        #ifdef JES_STANDALONE
        ClassDef(JESUncertaintyProvider,1);
        #endif //JES_STANDALONE

 protected:

  // Change the inputCollection
  bool setInputCollection(std::string CollectionName);
    
  // Maximum number of components
  static const unsigned int m_nUncertainties_EMJES = 7;
  static const unsigned int m_nUncertainties_inSitu = 2;
  
  // Maximum number of vertices available
  static const unsigned int m_nVertices = 7;
  
  // For standalone compatibility
  const double m_GeV;
  
  // Pointers to the uncertainty histograms 
  map<int,TH2D*> m_uncGraph;

  // Pointers to the pileup histograms
  map<int,TH2D*> m_pileupUncGraph;

  // Input File
  TFile* m_inputFile; 

  // Helper function to check requested components
  double getComponents(int currentBin, Components UncertComps, unsigned int nVtx = 1);

  // Name of jet collection
  std::string m_collectionName;

  // Name of file
  std::string m_fileName; 
  
  // Initialization flag
  bool m_isInit;
  
  // Deal with in-situ uncertainties (temporary)
  bool m_doInSitu;

  
};


#endif

---