d2/dd8/ObjectResolutionCalc_8h_source.html

//

//


#ifndef PhysicsTools_PatUtils_ObjectResolutionCalc_h

#define PhysicsTools_PatUtils_ObjectResolutionCalc_h


#include "FWCore/Framework/interface/EventSetup.h"

#include "FWCore/Framework/interface/Event.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "FWCore/Utilities/interface/Exception.h"


#include "DataFormats/PatCandidates/interface/Electron.h"

#include "DataFormats/PatCandidates/interface/Muon.h"

#include "DataFormats/PatCandidates/interface/Tau.h"

#include "DataFormats/PatCandidates/interface/MET.h"

#include "DataFormats/PatCandidates/interface/Jet.h"


#include "TF1.h"

#include "TH1.h"

#include "TFile.h"

#include "TKey.h"

#include "TString.h"

#include "TMultiLayerPerceptron.h"


namespace pat {


  class ObjectResolutionCalc {

  public:

    ObjectResolutionCalc();

    ObjectResolutionCalc(const TString& resopath, bool useNN);

    ~ObjectResolutionCalc();


    float obsRes(int obs, int eta, float eT);

    int etaBin(float eta);


#ifdef OBSOLETE

    void operator()(Electron& obj);

    void operator()(Muon& obj);

    void operator()(Tau& obj);

    void operator()(Jet& obj);

    void operator()(MET& obj);

#else

    // WORKAROUND

    template <typename T>

    void operator()(T& obj) {}

#endif


  private:

    TFile* resoFile_;

    std::vector<float> etaBinVals_;

    TF1 fResVsEt_[10][10];

    TMultiLayerPerceptron* network_[10];

    bool useNN_;

  };


}  // namespace pat


#endif