iceCube.uhe.analysis
public class ModelTestByPoissonBinnedLikelihoodFactory extends java.lang.Object
PoissonBinnedLikelihoodCalculator calBG (for the atmospheric background) PoissonBinnedLikelihoodCalculator calSinal (for the neutrino model such as GZK to be tested) PoissonBinnedLikelihoodCalculator calNuisanceSingal (for the nuisance signal - like E^-2 against the GZK test) PossionBinnedLikelihoodCalculator calNull Null Hypothesis PossionBinnedLikelihoodCalculator calAlter Alternative Hypothesis boolean includeNuisance true for including the nuisance signal contribution
Modifier and Type | Field and 説明 |
---|---|
PoissonBinnedLikelihoodCalculator |
calAlter |
PoissonBinnedLikelihoodCalculator |
calBG |
PoissonBinnedLikelihoodCalculator |
calHybrid |
PoissonBinnedLikelihoodCalculator |
calNuisanceSignal |
PoissonBinnedLikelihoodCalculator |
calNull |
PoissonBinnedLikelihoodCalculator |
calSignal |
protected boolean |
includeNuisance |
protected double |
llhAlter |
protected double |
llhNull |
protected double |
maximizedNuisanceFactor |
protected double |
maximizedSignalFactor |
Constructor and 説明 |
---|
ModelTestByPoissonBinnedLikelihoodFactory(PoissonBinnedLikelihoodCalculator calBG,
PoissonBinnedLikelihoodCalculator calSignal)
The default constructor.
|
ModelTestByPoissonBinnedLikelihoodFactory(PoissonBinnedLikelihoodCalculator calBG,
PoissonBinnedLikelihoodCalculator calSignal,
PoissonBinnedLikelihoodCalculator calNuisanceSignal)
The default constructor.
|
Modifier and Type | Method and 説明 |
---|---|
double |
buildLikelihoodForAlternativeHypothesis(boolean saturatedLikelihood,
boolean useTheResultsOfReplicaExperiment)
Build the alternaive hypothesis :
saturatedLikelihood == true: Build the saturated likelihood from the null hypothesis calNull.
|
double |
buildLikelihoodForAlternativeHypothesis(boolean saturatedLikelihood,
boolean nuisanceAlternative,
boolean useTheResultsOfReplicaExperiment)
Build the alternaive hypothesis :
saturatedLikelihood == true: Build the saturated likelihood from the null hypothesis calNull.
|
double |
buildLikelihoodForHybridHypothesis(boolean useTheResultsOfReplicaExperiment)
Build the "hybrid" hypothesis - signal with floating normalization + nuisance with floating nomalization
to maximize the likelihood
|
double |
buildLikelihoodForNullHypothesis(boolean runReplicaExperiment)
Build the null hypothesis : the (GZK or astro) signal + background.
|
double |
buildLikelihoodForNullHypothesis(boolean backgroundONLY,
boolean runReplicaExperiment)
Build the null hypothesis : the (GZK or astro) signal + background.
|
void |
doNotIncludeNuisanceSignal()
Do not let the likelihood calculation include the nuisance signal
|
java.util.ListIterator |
getllhRatioIterator() |
double |
getModelNormalizationToMaximizeLikelihood()
Return the normalization of the model to maximize the binned poisson pdf.
|
double |
getNuisanceNormalizationToMaximizeLikelihood()
Return the nuisance normalization of the model to maximize the binned poisson pdf.
|
static java.util.LinkedList |
inputList(java.io.InputStream in) |
boolean |
isTheNuisanceSignalIncluded() |
void |
makeCollectionOfLogLikelihoodRatio(int type,
int runTimes)
make a list collection of the log-likelihood -log(binned Poisson pdf) obtained by running the replica experiment
by runTimes under the null hypothesis of bg + neutrino model (+ nuisance signal)
The likeliihood ratio type is following:
1 saturated likelihood/null hypothesis(bg + model signal (+ nuisance))
2 alternative hypothesis(bg + model signal with floated normalization)/null hypothesis(bg + model signal (+ nuisance))
3 saturated likelihood/alternative hypothesis (bg + model signal with floated normalization)
4 (bg + nuisance model signal with floated normalization)/(bg + model signal with floated normalization)
5 alternative hypothesis(bg + either signal or nuisance model with floated normalization,
select the one giving larger likelihood)/null hypothesis(bg ONLY)
6 (bg + nuisance model signal with floated normalization + signal with floated norm)/
(bg + model signal with floated normalization)
|
static void |
outputLikelihhoRatioList(ModelTestByPoissonBinnedLikelihoodFactory testFac,
java.io.OutputStream out) |
static void |
outputMaximizedSignalFactorList(ModelTestByPoissonBinnedLikelihoodFactory testFac,
java.io.OutputStream out) |
void |
setAnotherNeutrinoSignalModel(PoissonBinnedLikelihoodCalculator calAnotherSignal)
Set the PoissonBinnedLikelihoodCalculator that stores the binned poisson data
for the signal model to be tested.
|
void |
setNuisanceSignal(PoissonBinnedLikelihoodCalculator calNuisanceSignal)
Set the PoissonBinnedLikelihoodCalculator that stores the binned poisson data
for the nuisance signal
|
public PoissonBinnedLikelihoodCalculator calBG
public PoissonBinnedLikelihoodCalculator calSignal
public PoissonBinnedLikelihoodCalculator calNuisanceSignal
public PoissonBinnedLikelihoodCalculator calNull
public PoissonBinnedLikelihoodCalculator calAlter
public PoissonBinnedLikelihoodCalculator calHybrid
protected boolean includeNuisance
protected double llhNull
protected double llhAlter
protected double maximizedSignalFactor
protected double maximizedNuisanceFactor
public ModelTestByPoissonBinnedLikelihoodFactory(PoissonBinnedLikelihoodCalculator calBG, PoissonBinnedLikelihoodCalculator calSignal, PoissonBinnedLikelihoodCalculator calNuisanceSignal)
public ModelTestByPoissonBinnedLikelihoodFactory(PoissonBinnedLikelihoodCalculator calBG, PoissonBinnedLikelihoodCalculator calSignal)
public void setNuisanceSignal(PoissonBinnedLikelihoodCalculator calNuisanceSignal)
public void setAnotherNeutrinoSignalModel(PoissonBinnedLikelihoodCalculator calAnotherSignal)
public void doNotIncludeNuisanceSignal()
public boolean isTheNuisanceSignalIncluded()
public double buildLikelihoodForNullHypothesis(boolean backgroundONLY, boolean runReplicaExperiment)
public double buildLikelihoodForNullHypothesis(boolean runReplicaExperiment)
public double buildLikelihoodForAlternativeHypothesis(boolean saturatedLikelihood, boolean nuisanceAlternative, boolean useTheResultsOfReplicaExperiment)
public double buildLikelihoodForAlternativeHypothesis(boolean saturatedLikelihood, boolean useTheResultsOfReplicaExperiment)
public double buildLikelihoodForHybridHypothesis(boolean useTheResultsOfReplicaExperiment)
public double getModelNormalizationToMaximizeLikelihood()
public double getNuisanceNormalizationToMaximizeLikelihood()
public void makeCollectionOfLogLikelihoodRatio(int type, int runTimes)
1 saturated likelihood/null hypothesis(bg + model signal (+ nuisance)) 2 alternative hypothesis(bg + model signal with floated normalization)/null hypothesis(bg + model signal (+ nuisance)) 3 saturated likelihood/alternative hypothesis (bg + model signal with floated normalization) 4 (bg + nuisance model signal with floated normalization)/(bg + model signal with floated normalization) 5 alternative hypothesis(bg + either signal or nuisance model with floated normalization, select the one giving larger likelihood)/null hypothesis(bg ONLY) 6 (bg + nuisance model signal with floated normalization + signal with floated norm)/ (bg + model signal with floated normalization)
public java.util.ListIterator getllhRatioIterator()
public static void outputLikelihhoRatioList(ModelTestByPoissonBinnedLikelihoodFactory testFac, java.io.OutputStream out) throws java.io.IOException
java.io.IOException
public static void outputMaximizedSignalFactorList(ModelTestByPoissonBinnedLikelihoodFactory testFac, java.io.OutputStream out) throws java.io.IOException
java.io.IOException
public static java.util.LinkedList inputList(java.io.InputStream in) throws java.io.IOException
java.io.IOException