pat::ObjectResolutionCalc Class Reference

Class to calculate MC resolutions for pat objects. More...

#include "PhysicsTools/PatUtils/interface/ObjectResolutionCalc.h"

Public Member Functions
int	etaBin (float eta)
	ObjectResolutionCalc ()
	ObjectResolutionCalc (const TString &resopath, bool useNN)
float	obsRes (int obs, int eta, float eT)
template<typename T >
void	operator() (T &obj)
	~ObjectResolutionCalc ()

Private Attributes
std::vector< float >	etaBinVals_
TF1	fResVsEt_ [10][10]
TMultiLayerPerceptron *	network_ [10]
TFile *	resoFile_
bool	useNN_

Detailed Description

Class to calculate MC resolutions for pat objects.

Author

Jan Heyninck, Petra Van Mulders, Christophe Delaere

Version

Id

ObjectResolutionCalc.h,v 1.5 2008/10/08 19:19:25 gpetrucc Exp

Definition at line 28 of file ObjectResolutionCalc.h.

Constructor & Destructor Documentation

ObjectResolutionCalc() [1/2]

pat::ObjectResolutionCalc::ObjectResolutionCalc

(

)

ObjectResolutionCalc() [2/2]

ObjectResolutionCalc::ObjectResolutionCalc	(	const TString &	resopath,
		bool	useNN = false
	)

Definition at line 15 of file ObjectResolutionCalc.cc.

References b, etaBin(), etaBinVals_, Exception, fResVsEt_, submitPVResolutionJobs::key, relativeConstraints::keys, edm::errors::LogicError, Skims_PA_cff::name, network_, AlCaHLTBitMon_ParallelJobs::p, resoFile_, createJobs::tmp, and useNN_.

                                                                                      : useNN_(useNN) {
  edm::LogVerbatim("ObjectResolutionCalc")
      << ("ObjectResolutionCalc") << "=== Constructing a TopObjectResolutionCalc...";
  resoFile_ = new TFile(resopath);
  if (!resoFile_)
    throw edm::Exception(edm::errors::LogicError)
        << "ObjectResolutionCalc: no resolutions fits for this file available: " << resopath << "...";

  TString resObsName[8] = {"_ares", "_bres", "_cres", "_dres", "_thres", "_phres", "_etres", "_etares"};
  TList* keys = resoFile_->GetListOfKeys();
  TIter nextitem(keys);
  TKey* key = nullptr;
  while ((key = (TKey*)nextitem())) {
    TString name = key->GetName();
    if (useNN_) {
      for (Int_t ro = 0; ro < 8; ro++) {
        TString obsName = resObsName[ro];
        obsName += "_NN";
        if (name.Contains(obsName)) {
          network_[ro] = (TMultiLayerPerceptron*)resoFile_->GetKey(name)->ReadObj();
        }
      }
    } else {
      if (name.Contains("etabin") && (!name.Contains("etbin"))) {
        for (int p = 0; p < 8; p++) {
          if (name.Contains(resObsName[p])) {
            TString etabin = name;
            etabin.Remove(0, etabin.Index("_") + 1);
            etabin.Remove(0, etabin.Index("_") + 7);
            int etaBin = etabin.Atoi();
            TH1F* tmp = (TH1F*)(resoFile_->GetKey(name)->ReadObj());
            fResVsEt_[p][etaBin] = (TF1)(*(tmp->GetFunction("F_" + name)));
          }
        }
      }
    }
  }
  // find etabin values
  TH1F* tmpEta = (TH1F*)(resoFile_->GetKey("hEtaBins")->ReadObj());
  for (int b = 1; b <= tmpEta->GetNbinsX(); b++)
    etaBinVals_.push_back(tmpEta->GetXaxis()->GetBinLowEdge(b));
  etaBinVals_.push_back(tmpEta->GetXaxis()->GetBinUpEdge(tmpEta->GetNbinsX()));
  edm::LogVerbatim("ObjectResolutionCalc") << "Found " << etaBinVals_.size() - 1 << " eta-bins with edges: ( ";
  for (size_t u = 0; u < etaBinVals_.size(); u++)
    edm::LogVerbatim("ObjectResolutionCalc") << etaBinVals_[u] << ", ";
  edm::LogVerbatim("ObjectResolutionCalc") << "\b\b )" << std::endl;

  edm::LogVerbatim("ObjectResolutionCalc") << "=== done." << std::endl;
}

~ObjectResolutionCalc()

ObjectResolutionCalc::~ObjectResolutionCalc

(

)

Definition at line 66 of file ObjectResolutionCalc.cc.

References resoFile_.

{ delete resoFile_; }

Member Function Documentation

etaBin()

int ObjectResolutionCalc::etaBin

(

float

eta

)

Definition at line 76 of file ObjectResolutionCalc.cc.

References newFWLiteAna::bin, PVValHelper::eta, etaBinVals_, and mps_fire::i.

Referenced by ObjectResolutionCalc().

                                          {
  int nrEtaBins = etaBinVals_.size() - 1;
  int bin = nrEtaBins - 1;
  for (int i = 0; i < nrEtaBins; i++) {
    if (fabs(eta) > etaBinVals_[i] && fabs(eta) < etaBinVals_[i + 1]) {
      bin = i;
      break;
    }
  }
  return bin;
}

obsRes()

float ObjectResolutionCalc::obsRes	(	int	obs,
		int	eta,
		float	eT
	)

Definition at line 68 of file ObjectResolutionCalc.cc.

References PVValHelper::eta, Exception, fResVsEt_, edm::errors::LogicError, and useNN_.

                                                             {
  if (useNN_)
    throw edm::Exception(edm::errors::LogicError,
                         "TopObjectResolutionCalc::obsRes should never be called when using a NN for resolutions.");
  float res = fResVsEt_[obs][eta].Eval(eT);
  return res;
}

operator()()

template<typename T >

void pat::ObjectResolutionCalc::operator() ( T &  obj )

inline

Definition at line 46 of file ObjectResolutionCalc.h.

{}

Member Data Documentation

etaBinVals_

std::vector<float> pat::ObjectResolutionCalc::etaBinVals_

private

Definition at line 51 of file ObjectResolutionCalc.h.

Referenced by etaBin(), and ObjectResolutionCalc().

fResVsEt_

TF1 pat::ObjectResolutionCalc::fResVsEt_[10][10]

private

Definition at line 52 of file ObjectResolutionCalc.h.

Referenced by ObjectResolutionCalc(), and obsRes().

network_

TMultiLayerPerceptron* pat::ObjectResolutionCalc::network_[10]

private

Definition at line 53 of file ObjectResolutionCalc.h.

Referenced by ObjectResolutionCalc().

resoFile_

TFile* pat::ObjectResolutionCalc::resoFile_

private

Definition at line 50 of file ObjectResolutionCalc.h.

Referenced by ObjectResolutionCalc(), and ~ObjectResolutionCalc().

useNN_

bool pat::ObjectResolutionCalc::useNN_

private

Definition at line 54 of file ObjectResolutionCalc.h.

Referenced by ObjectResolutionCalc(), and obsRes().