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 35 of file ObjectResolutionCalc.h.

Constructor & Destructor Documentation

pat::ObjectResolutionCalc::ObjectResolutionCalc ( )

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

Definition at line 9 of file ObjectResolutionCalc.cc.

References b, etaBin(), etaBinVals_, fResVsEt_, crabWrapper::key, relativeConstraints::keys, 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_)
     edm::LogError("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 ( )

Definition at line 59 of file ObjectResolutionCalc.cc.

References resoFile_.

59 { delete resoFile_; }

pat::ObjectResolutionCalc::resoFile_

TFile * resoFile_

Definition: ObjectResolutionCalc.h:57

Member Function Documentation

int ObjectResolutionCalc::etaBin ( float eta )

Definition at line 69 of file ObjectResolutionCalc.cc.

References newFWLiteAna::bin, reco::LeafCandidate::et(), reco::LeafCandidate::eta(), etaBinVals_, mps_fire::i, network_, getGTfromDQMFile::obj, obsRes(), operator()(), useNN_, and findQualityFiles::v.

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;
   }
   return bin;
 }

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

Definition at line 61 of file ObjectResolutionCalc.cc.

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

Referenced by etaBin().

                                                              {
   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;
 }