/data/refman/pasoursint/CMSSW_5_3_9_patch3/src/RecoEgamma/ElectronIdentification/interface/LikelihoodPdf.h

Go to the documentation of this file.

#ifndef LikelihoodPdf_H
#define LikelihoodPdf_H

#include "CondFormats/PhysicsToolsObjects/interface/Histogram.h"
#include "CondFormats/EgammaObjects/interface/ElectronLikelihoodCalibration.h"
#include <TH1F.h>
#include <TDirectory.h>
#include <string>
#include <vector>
#include <map>

class LikelihoodPdf {

 public:
  LikelihoodPdf() {};
  LikelihoodPdf(const LikelihoodPdf& pdf) {}; 
  LikelihoodPdf(const char* name, const char* species, int ecalsubdet, int ptbin);
  virtual ~LikelihoodPdf();
  
  void initFromDB(const ElectronLikelihoodCalibration *calibration);

  void split(std::map<std::string,float> splitFractions, bool splitPdf = false);

  float getVal(float x, std::string catName="NOSPLIT", bool normalized = true);

  std::string getName() { return _name; }

  std::string getSpecies() { return _species; }



 private:

  float normalization(const PhysicsTools::Calibration::HistogramF *thePdf);
  
  std::string _name;
  std::string _species;
  int _ecalsubdet;
  int _ptbin;

  std::map<std::string,const PhysicsTools::Calibration::HistogramF*> _splitPdf;
  std::map<std::string,std::string> _splitRule;

};

#endif