LikelihoodPdfProduct Class Reference

#include <LikelihoodPdfProduct.h>

Public Member Functions
void	addPdf (const char *specname, const char *name, bool splitPdf=false)
	add a pdf for a species, splitted or not More...
void	addSpecies (const char *name, float priorWeight=1.)
	add a species (hypothesis) to the likelihood, with a priori probability More...
float	getRatio (const char *specName, const std::vector< float > &measurements, const std::string &) const
	get the likelihood ratio p(a priori) * L(specName) / L_tot More...
void	initFromDB (const ElectronLikelihoodCalibration *calibration)
	initialize the PDFs from CondDB More...
	LikelihoodPdfProduct (const char *name, int ecalsubdet, int ptbin)
void	setSplitFrac (const char *specname, const char *catName, float frac=1.0)
	set the fraction of one category for a given species More...
	~LikelihoodPdfProduct ()

Private Member Functions
float	getSpeciesProb (const char *specName, const std::vector< float > &measurements, const std::string &gsfClass) const

Private Attributes
const ElectronLikelihoodCalibration *	_calibration
int	_ecalsubdet
std::string	_name
std::vector< float >	_priorList
int	_ptbin
std::vector< LikelihoodSpecies * >	_specList

Detailed Description

Definition at line 12 of file LikelihoodPdfProduct.h.

Constructor & Destructor Documentation

LikelihoodPdfProduct::LikelihoodPdfProduct	(	const char *	name,
		int	ecalsubdet,
		int	ptbin
	)

Definition at line 7 of file LikelihoodPdfProduct.cc.

References _ecalsubdet, _name, _ptbin, and AlCaHLTBitMon_QueryRunRegistry::string.

{
  _name=std::string(name);
  _ecalsubdet = ecalsubdet;
  _ptbin = ptbin;
}

LikelihoodPdfProduct::~LikelihoodPdfProduct

(

)

Definition at line 17 of file LikelihoodPdfProduct.cc.

References _specList.

{
  // do the delete's 
  std::vector<LikelihoodSpecies*>::iterator specItr;
  for(specItr=_specList.begin(); specItr!=_specList.end(); specItr++) {
    delete *specItr;
  }
}

Member Function Documentation

void LikelihoodPdfProduct::addPdf	(	const char *	specname,
		const char *	name,
		bool	splitPdf = false
	)

add a pdf for a species, splitted or not

Definition at line 46 of file LikelihoodPdfProduct.cc.

References _calibration, _ecalsubdet, _ptbin, _specList, LikelihoodSpecies::addPdf(), LikelihoodSpecies::getName(), LikelihoodSpecies::getSplitFractions(), LikelihoodPdf::initFromDB(), and LikelihoodPdf::split().

Referenced by ElectronLikelihood::Setup().

{
  std::vector<LikelihoodSpecies*>::const_iterator specItr;
  for(specItr=_specList.begin();specItr!=_specList.end();specItr++) {
    LikelihoodSpecies* species = *specItr;
    if(strcmp(species->getName(),specname)==0) {
      LikelihoodPdf *pdf = new LikelihoodPdf(name,species->getName(),_ecalsubdet,_ptbin);
      pdf->split(species->getSplitFractions(),splitPdf);
      pdf->initFromDB(_calibration);
      species->addPdf(pdf);
    }
  }
}

void LikelihoodPdfProduct::addSpecies	(	const char *	name,
		float	priorWeight = 1.
	)

add a species (hypothesis) to the likelihood, with a priori probability

Definition at line 37 of file LikelihoodPdfProduct.cc.

References _specList.

Referenced by ElectronLikelihood::Setup().

{
  LikelihoodSpecies* species = new LikelihoodSpecies(name,priorWeight);
  _specList.push_back(species);
}

float LikelihoodPdfProduct::getRatio	(	const char *	specName,
		const std::vector< float > &	measurements,
		const std::string &	gsfClass
	)		const

get the likelihood ratio p(a priori) * L(specName) / L_tot

if the pdf is not splitted, assign the split fraction = 1

Definition at line 80 of file LikelihoodPdfProduct.cc.

References _specList, LikelihoodSpecies::getName(), getSpeciesProb(), and LikelihoodSpecies::getSplitFractions().

Referenced by ElectronLikelihood::result(), and ElectronLikelihood::resultLog().

{
  float sigProb=0, bkgProb=0;
  std::vector<LikelihoodSpecies*>::const_iterator specItr;
  for(specItr=_specList.begin();specItr!=_specList.end();specItr++) {
    const LikelihoodSpecies* species = *specItr;
    std::map<std::string,float> const& splitFractions = species->getSplitFractions();
    std::map<std::string,float>::const_iterator iter = splitFractions.find(gsfClass);
    float splitFr= (splitFractions.size()==0) ? 1. : iter->second;
    if(strcmp(species->getName(),specname)==0) {
      sigProb=splitFr*getSpeciesProb(specname,measurements,gsfClass);
    }
    else {
      bkgProb+=splitFr*getSpeciesProb(species->getName(),measurements,gsfClass);
    } 
  }
  if(sigProb+bkgProb>0)
    return sigProb/(sigProb+bkgProb);
  else
    return -1.;
}

float LikelihoodPdfProduct::getSpeciesProb ( const char *  specName, const std::vector< float > &  measurements, const std::string &  gsfClass  ) const

private

Definition at line 108 of file LikelihoodPdfProduct.cc.

References _specList, LikelihoodSpecies::getListOfPdfs(), LikelihoodSpecies::getName(), and LikelihoodSpecies::getPrior().

Referenced by getRatio().

{
  float bareProb=1.;
  float priorWeight=1.;
  std::vector<LikelihoodSpecies*>::const_iterator specItr;
  for(specItr=_specList.begin();specItr!=_specList.end();specItr++) {
    const LikelihoodSpecies* species = *specItr;
    if(strcmp(species->getName(),specName)==0) {
      for(unsigned int ipdf=0; ipdf< species->getListOfPdfs().size(); ipdf++) {
        bareProb*=species->getListOfPdfs()[ipdf]->getVal(measurements.at(ipdf),gsfClass);
      }
      priorWeight=species->getPrior();
      break;
    }
  }
  edm::LogInfo("LikelihoodPdfProduct") << "Species: " << specName
                       << " bare probability = " << bareProb
                       << " with a priori probability = " << priorWeight;
  return priorWeight*bareProb;
}

void LikelihoodPdfProduct::initFromDB

(

const ElectronLikelihoodCalibration *

calibration

)

initialize the PDFs from CondDB

Definition at line 29 of file LikelihoodPdfProduct.cc.

References _calibration.

Referenced by ElectronLikelihood::Setup().

{
  _calibration=calibration;
}

void LikelihoodPdfProduct::setSplitFrac	(	const char *	specname,
		const char *	catName,
		float	frac = 1.0
	)

set the fraction of one category for a given species

Definition at line 65 of file LikelihoodPdfProduct.cc.

References _specList, LikelihoodSpecies::getName(), LikelihoodSpecies::setSplitFraction(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by ElectronLikelihood::Setup().

                                                    {
  std::vector<LikelihoodSpecies*>::const_iterator specItr;
  for(specItr=_specList.begin();specItr!=_specList.end();specItr++) {
    LikelihoodSpecies* species = *specItr;
    if(strcmp(species->getName(),specname)==0) {
      species->setSplitFraction( make_pair(std::string(catName), frac) );
      break;
    }
  }
}

Member Data Documentation

const ElectronLikelihoodCalibration* LikelihoodPdfProduct::_calibration

private

Definition at line 36 of file LikelihoodPdfProduct.h.

Referenced by addPdf(), and initFromDB().

int LikelihoodPdfProduct::_ecalsubdet

private

Definition at line 39 of file LikelihoodPdfProduct.h.

Referenced by addPdf(), and LikelihoodPdfProduct().

std::string LikelihoodPdfProduct::_name

private

Definition at line 35 of file LikelihoodPdfProduct.h.

Referenced by plotting.Subtract::__str__(), plotting.FakeDuplicate::__str__(), plotting.AggregateBins::__str__(), plotting.AggregateHistos::__str__(), plotting.ROC::__str__(), plotting.Plot::_createOne(), plotting.PlotGroup::_drawSeparate(), plotting.PlotGroup::_save(), plotting.Subtract::create(), plotting.FakeDuplicate::create(), plotting.AggregateBins::create(), plotting.AggregateHistos::create(), plotting.PlotGroup::draw(), plotting.Plot::getName(), plotting.PlotterFolder::getName(), plotting.PlotterFolder::getSelectionNameIterator(), LikelihoodPdfProduct(), validation.SimpleSample::name(), and plotting.PlotterItem::readDirs().

std::vector<float> LikelihoodPdfProduct::_priorList

private

Definition at line 38 of file LikelihoodPdfProduct.h.

int LikelihoodPdfProduct::_ptbin

private

Definition at line 40 of file LikelihoodPdfProduct.h.

Referenced by addPdf(), and LikelihoodPdfProduct().

std::vector<LikelihoodSpecies*> LikelihoodPdfProduct::_specList

private

Definition at line 37 of file LikelihoodPdfProduct.h.

Referenced by addPdf(), addSpecies(), getRatio(), getSpeciesProb(), setSplitFrac(), and ~LikelihoodPdfProduct().