#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().

 {
   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.

 {
   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().

 {
   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.

 {
   _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.

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