iceCube.uhe.analysis
public class ModelParaPvalueCalculatorByPoissonBinnedLH extends java.lang.Object
PoissonBinnedLikelihoodCalculator calBG (for the atmospheric background) PoissonBinnedLikelihoodCalculator calSignal (for the neutrino model such as GZK with the parameters to be tested) PoissonBinnedLikelihoodCalculator calNuisanceSingal (for the nuisance signal - like E^-2 against the GZK test)are generated by the correponding binned data from the files and then executes the methods provied by ModelTestByPoissonBinnedLikelihoodFactory. The signal hypothesis with all the other values of parameters are stored in the list container calSignalSetsList. The example for using this class - The constraints on the source evolution parameters (m,Zmax). The likeliihood ratio
alternative (calBG + model stored in calSignalSetsList (+ nuisance ) with parameters to maximize Likelihood)/ null hypothesis(calBG + calSignal with a given neutrino parameter (+ nuisance to maximize likelihood))is used for the test-statstic
Modifier and Type | Field and Description |
---|---|
PoissonBinnedLikelihoodCalculator |
calBG
the poisson bin data container for BG
|
PoissonBinnedLikelihoodCalculator |
calNuisanceSignal
container for Nuisance Signal for the profole likelihood
|
PoissonBinnedLikelihoodCalculator |
calSignal
the container for Signal with the set of parameter values to be tested
|
protected java.util.List |
calSignalSetsList
List to store the PoissonBinnedLikelihoodCalculator calSignal with various model parameters
|
protected PoissonBinnedLikelihoodCalculator |
calSignalToMaxLH
container of the signal with the model parameters that maximized likelihood
|
protected boolean |
includeNuisance |
static int |
runTimesDefault
default number of runs in the replica experiment
to calculate the pdf of the likelihood ratio
|
protected ModelTestByPoissonBinnedLikelihoodFactory |
testFactory
The test factory
|
Constructor and Description |
---|
ModelParaPvalueCalculatorByPoissonBinnedLH(PoissonBinnedLikelihoodCalculator calBG,
PoissonBinnedLikelihoodCalculator calSignal,
java.util.List calSignalSetsList)
constructor for no nuislance case:
|
ModelParaPvalueCalculatorByPoissonBinnedLH(PoissonBinnedLikelihoodCalculator calBG,
PoissonBinnedLikelihoodCalculator calSignal,
java.util.Map calSignalMap)
constructor for no nuislance case:
|
ModelParaPvalueCalculatorByPoissonBinnedLH(PoissonBinnedLikelihoodCalculator calBG,
PoissonBinnedLikelihoodCalculator calSignal,
PoissonBinnedLikelihoodCalculator calNuisanceSignal,
java.util.List calSignalSetsList)
constructor:
|
ModelParaPvalueCalculatorByPoissonBinnedLH(PoissonBinnedLikelihoodCalculator calBG,
PoissonBinnedLikelihoodCalculator calSignal,
PoissonBinnedLikelihoodCalculator calNuisanceSignal,
java.util.Map calSignalMap)
constructor:
|
Modifier and Type | Method and Description |
---|---|
void |
buildCalSignalSetsListFromMap(java.util.Map calSignalMap)
build the list to store PoissonBinnedLikelihoodCalculator calSignal with various model parameters
from Map(Key=PoissonBinnedLikelihoodCalculator).
|
double |
getLLHRatio(boolean runReplicaExperiment)
Scan all the PoissonBinnedLikelihoodCalculators stored in List calSignalSetsList
to calculate the LLH.
|
java.util.ListIterator |
getllhRatioIterator()
return the list iterator to access the list container of the llhRatio data
obtained by running replica experiments with the method makeCollectionOfLogLikelihoodRatio()
|
PoissonBinnedLikelihoodCalculator |
getModelToMaximizeLH()
Return the PoissonBinnedLikelihoodCalculator to contain the model that maximized the poisson binned likelihood.
|
double |
getPvalue(double llhRatioObserved)
Calculate the p-value of the model parameters of calSignal.
|
double |
getPvalue(double llhRatioObserved,
int trialTimes)
Calculate the p-value of the model parameters of calSignal.
|
java.util.ListIterator |
getSignalSetsListIterator()
return the list iterator to access the signal poisson binned data
with various model parameters.
|
static void |
main(java.lang.String[] args)
A simple main method for debugging
|
void |
makeCollectionOfLogLikelihoodRatio(int runTimes)
Calculate collections of the log likelihhood ratio
|
void |
useTheFullSearchToMaximizeLLH()
Use the full search for the parameters to maximize LLH in the method getLLHRatio() or getPvalue().
|
void |
useTheQuickSearchToMaximizeLLH()
Use the quick search for the parameters to maximize LLH in the method getLLHRatio() or getPvalue().
|
public PoissonBinnedLikelihoodCalculator calBG
public PoissonBinnedLikelihoodCalculator calSignal
public PoissonBinnedLikelihoodCalculator calNuisanceSignal
protected PoissonBinnedLikelihoodCalculator calSignalToMaxLH
protected ModelTestByPoissonBinnedLikelihoodFactory testFactory
public static final int runTimesDefault
protected boolean includeNuisance
protected java.util.List calSignalSetsList
public ModelParaPvalueCalculatorByPoissonBinnedLH(PoissonBinnedLikelihoodCalculator calBG, PoissonBinnedLikelihoodCalculator calSignal, PoissonBinnedLikelihoodCalculator calNuisanceSignal, java.util.List calSignalSetsList)
PoissonBinnedLikelihoodCalculator calBG :for the atmospheric background PoissonBinnedLikelihoodCalculator calSignal :for the neutrino model with the parmeters to be tested PoissonBinnedLikelihoodCalculator calNuisanceSingal:for the nuisance signal - like E^-2 against the GZK test List calSignalSetsList : A container to store PoissonBinnedLikelihoodCalculator objects of the signal hypothesis but with all the parameter values allowd in their range. It will build the alternative hypothesis.
public ModelParaPvalueCalculatorByPoissonBinnedLH(PoissonBinnedLikelihoodCalculator calBG, PoissonBinnedLikelihoodCalculator calSignal, java.util.List calSignalSetsList)
PoissonBinnedLikelihoodCalculator calBG :for the atmospheric background PoissonBinnedLikelihoodCalculator calSignal :for the neutrino model with the parmeters to be tested List calSignalSetsList : A container to store PoissonBinnedLikelihoodCalculator objects of the signal hypothesis but with all the parameter values allowd in their range. It will build the alternative hypothesis.
public ModelParaPvalueCalculatorByPoissonBinnedLH(PoissonBinnedLikelihoodCalculator calBG, PoissonBinnedLikelihoodCalculator calSignal, PoissonBinnedLikelihoodCalculator calNuisanceSignal, java.util.Map calSignalMap)
PoissonBinnedLikelihoodCalculator calBG :for the atmospheric background PoissonBinnedLikelihoodCalculator calSignal :for the neutrino model with the parmeters to be tested PoissonBinnedLikelihoodCalculator calNuisanceSingal:for the nuisance signal - like E^-2 against the GZK test Map calSignalMap : A container to store PoissonBinnedLikelihoodCalculator objects of the signal hypothesis but with all the parameter values allowd in their range. It will build the alternative hypothesis. The PoissonBinnedLikelihoodCalculator objects should be stored as the key (not value) member of the map.
public ModelParaPvalueCalculatorByPoissonBinnedLH(PoissonBinnedLikelihoodCalculator calBG, PoissonBinnedLikelihoodCalculator calSignal, java.util.Map calSignalMap)
PoissonBinnedLikelihoodCalculator calBG :for the atmospheric background PoissonBinnedLikelihoodCalculator calSignal :for the neutrino model with the parmeters to be tested Map calSignalMap : A container to store PoissonBinnedLikelihoodCalculator objects of the signal hypothesis but with all the parameter values allowd in their range. It will build the alternative hypothesis. The PoissonBinnedLikelihoodCalculator objects should be stored as the key (not value) member of the map.
public void useTheQuickSearchToMaximizeLLH()
public void useTheFullSearchToMaximizeLLH()
public double getLLHRatio(boolean runReplicaExperiment)
alternative (calBG + model stored in calSignalSetsList (+ nuisance ) with parameters to maximize Likelihood)/ null hypothesis(calBG + calSignal with a given neutrino parameter (+ nuisance to maximize likelihood))is returned.
public PoissonBinnedLikelihoodCalculator getModelToMaximizeLH()
public void makeCollectionOfLogLikelihoodRatio(int runTimes)
alternative (calBG + model stored in calSignalSetsList (+ nuisance ) with parameters to maximize Likelihood)/ null hypothesis(calBG + calSignal with a given neutrino parameter (+ nuisance to maximize likelihood))by running replica experiments.
public double getPvalue(double llhRatioObserved, int trialTimes)
alternative (calBG + model stored in calSignalSetsList (+ nuisance ) with parameters to maximize Likelihood)/ null hypothesis(calBG + calSignal with a given neutrino parameter (+ nuisance to maximize likelihood))The p.d.f. of the test statstic is simulated by calling the method makeCollectionOfLogLikelihoodRatio(runTimesDefault) internally. int trialTimes is the trail number to run the pseudo-experiments in the p-value calculation. The acuracy of p-value will be O(1/trialTimes).
public double getPvalue(double llhRatioObserved)
alternative (calBG + model stored in calSignalSetsList (+ nuisance ) with parameters to maximize Likelihood)/ null hypothesis(calBG + calSignal with a given neutrino parameter (+ nuisance to maximize likelihood))The p.d.f. of the test statstic is simulated by calling the method makeCollectionOfLogLikelihoodRatio(runTimesDefault) internally. int runTimes is the trail number to run the pseudo-experiments in the p-value calculation. The default runTime number is used in the p-value calculation.
public java.util.ListIterator getSignalSetsListIterator()
public java.util.ListIterator getllhRatioIterator()
public void buildCalSignalSetsListFromMap(java.util.Map calSignalMap)
public static void main(java.lang.String[] args) throws java.io.IOException
java.io.IOException