numRecipes

Class PoissonBinnedLikelihoodRatioTest

java.lang.Object

numRecipes.PoissonBinnedLikelihoodRatioTest

public class PoissonBinnedLikelihoodRatioTest
extends java.lang.Object

Perform test based upon the ratio of the maximum likelihood ratio built from the binned Poisson. The null hypothesis (=background) and the alternative (=signal+background) hypothesis are describied, respectively, by the class PoissonBinnedLikelihoodCalculator. The class PoissonBinnedLikelihoodRatioFactory takes care of the relevant calculations

    The null hypothesis        :  bg only
    The alternative hypothesis :  bg + lambda * signal
                                (lambda>=0)

    

Written originally by S. Yoshida for the IceCube EHE analysis. 2013/6/24

Field Summary

Fields

Modifier and Type	Field and Description
protected PoissonBinnedLikelihoodRatioFactory	llhFactory PoissonBinnedLikelihoodRatioFactory object

Constructor Summary

Constructors

Constructor and Description
PoissonBinnedLikelihoodRatioTest(PoissonBinnedLikelihoodCalculator calNull, PoissonBinnedLikelihoodCalculator calAlter) The default constructor.

Method Summary

Methods

Modifier and Type	Method and Description
java.util.ListIterator	getListIterator()
double	getMaxLogLikelihoodRatio(double scaleFactor, boolean runReplicaBGexperiment) return the maximized log likelihood ratio log(L/L0).
double	getPvalue(double scaleFactor, double precision, double observed_llhratio) calculate the p-value for a given log-likelihood ratio (llh) value.
static void	main(java.lang.String[] args) A simple main method
void	makeCollectionOfLogLikelihoodRatio(double scaleFactor, int runTimes) make a list collection of the log-likelihood ratio obtained by running the replica experiment by runTimes under the null (background-only) hypothesis.
void	makeCollectionOfLogLikelihoodRatioAlternative(double scaleFactor, int runTimes) make a list collection of the log-likelihood ratio obtained by running the replica experiment by runTimes under the algernative (signal) hypothesis.

Methods inherited from class java.lang.Object

clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Field Detail

llhFactory

protected PoissonBinnedLikelihoodRatioFactory llhFactory

PoissonBinnedLikelihoodRatioFactory object

Constructor Detail

PoissonBinnedLikelihoodRatioTest

public PoissonBinnedLikelihoodRatioTest(PoissonBinnedLikelihoodCalculator calNull,
                                PoissonBinnedLikelihoodCalculator calAlter)

The default constructor. Sets the likelihhood calculator of the null and alternative hypotheses.

Method Detail

getMaxLogLikelihoodRatio

public double getMaxLogLikelihoodRatio(double scaleFactor,
                              boolean runReplicaBGexperiment)

return the maximized log likelihood ratio log(L/L0). if runReplicaBGexperiment sets true, run the replica MC experiment with the H0 (=background only) hypothesis, and calculate the ratio aganst the results of the replica experiment. This method may be used for evaluating p-value for a given likelihood ratio (which can be provided by this method also setting runReplicaBGexperiment=false).

 
        double scaleFactor:  scaleFactor which is multiplied to mu_bg_i; if set to 1.0, no fudge factor is introduced.
        boolean runReplicaBGexperiment:    true - run Replica experiment with null hypotheis
        

This method calls the same function of the PoissonBinnedLikelihoodRatioFactory class.

makeCollectionOfLogLikelihoodRatio

public void makeCollectionOfLogLikelihoodRatio(double scaleFactor,
                                      int runTimes)

make a list collection of the log-likelihood ratio obtained by running the replica experiment by runTimes under the null (background-only) hypothesis.

makeCollectionOfLogLikelihoodRatioAlternative

public void makeCollectionOfLogLikelihoodRatioAlternative(double scaleFactor,
                                                 int runTimes)

make a list collection of the log-likelihood ratio obtained by running the replica experiment by runTimes under the algernative (signal) hypothesis.

getListIterator

public java.util.ListIterator getListIterator()

getPvalue

public double getPvalue(double scaleFactor,
               double precision,
               double observed_llhratio)

calculate the p-value for a given log-likelihood ratio (llh) value. Estimated by the collection of the llh obtained by calling makeCollectionOfLogLikelihoodRatio(int runTimes), i.e., running many replica experiment under the null hypothesis.

main

public static void main(java.lang.String[] args)
                 throws java.io.IOException

A simple main method

Throws:

java.io.IOException