public class ConfidenceLevel extends Object
For discoveries, 1-CLb indicates the probability that the background fluctuates to produce a distribution of candidates at least as signal-like as those observed in the data. For discovery, 1-CLb is required to be no more than 2.87x10-7, or twice that, depending on how one interprets what is meant by “five sigma,” including just one side of a Gaussian tail or both. A “three sigma” excess is defined to be 1-CLb = 1.3x10-3 or twice that. But forming discovery p-values, we must compute 1-CLb values of the order of 10-7. This computation involves generating of the order of 10^8 pseudoexperiments, just to be on the safe side.
Read Reference: HEP-EX/9902006. see: Tom Junk,NIM A434, p. 435-443,
Constructor and Description |
---|
ConfidenceLevel()
Default constructor.
|
ConfidenceLevel(int mc)
Construct ConfLevel
|
ConfidenceLevel(int mc,
boolean onesided)
Build confidence level.
|
Modifier and Type | Method and Description |
---|---|
void |
doc()
Show online documentation.
|
double |
get3sProbability()
Get 3s probability.
|
double |
get5sProbability()
Get 5s probability.
|
double |
getAverageCLs()
Get average CLs.
|
double |
getAverageCLsb()
Get average CLsb.
|
double |
getCLb()
Get the Confidence Level for the background only.
|
double |
getCLb(boolean use_sMC)
Get the Confidence Level for the background only.
|
double |
getCLs()
Get the Confidence Level defined by CLs = CLsb/CLb.
|
double |
getCLs(boolean use_sMC)
Get the Confidence Level defined by CLs = CLsb/CLb.
|
double |
getCLsb()
Get the Confidence Level for the signal plus background hypothesis
The confidence level for excluding the possibility of simultaneous presence
of new particle production and background (the s + b hypothesis)
|
double |
getCLsb(boolean use_sMC)
Get the Confidence Level for the signal plus background hypothesis
|
double |
getExpectedCLb_b()
Get the expected Confidence Level for the background only if there is
only background.These are indications of how well an experiment would
do on average in excluding a signal if the signal truly is not present, and are
the important figures of merit when optimizing an analysis for exclusion.
|
double |
getExpectedCLb_b(int sigma)
Get the expected Confidence Level for the background only if there is
only background.
|
double |
getExpectedCLb_sb()
Get the expected Confidence Level for the background only if there is
signal and background.
|
double |
getExpectedCLb_sb(int sigma)
Get the expected Confidence Level for the background only if there is
signal and background.
|
double |
getExpectedCLs_b()
Get getExpectedCLsb_b/getExpectedCLb_b.
|
double |
getExpectedCLs_b(int sigma)
Get getExpectedCLsb_b/getExpectedCLb_b
|
double |
getExpectedCLsb_b()
Get the expected Confidence Level for the signal plus background
hypothesis if there is only background.
|
double |
getExpectedCLsb_b(int sigma)
Get the expected Confidence Level for the signal plus background
hypothesis if there is only background.
|
double |
getExpectedStatistic_b() |
double |
getExpectedStatistic_b(int sigma)
Get the expected statistic value in the background only hypothesis
|
double |
getExpectedStatistic_sb(int sigma)
Get the expected statistic value in the signal plus background hypothesis
|
H1D |
getLNQb(int bins,
double min,
double max)
Get a histogram of a canonical -2lnQ plot for
background hypothesis (full)
|
H1D |
getLNQsb(int bins,
double min,
double max)
Get a histogram of a canonical -2lnQ plot for
for signal and background hypothesis
|
ArrayList<H1D> |
getResults(String Option)
Display sort of a "canonical" -2lnQ plot.
|
double |
getStatistic() |
void |
setBtot(double in) |
void |
setDtot(int in) |
void |
setLRB(double[] in) |
void |
setLRS(double[] in) |
void |
setStot(double in) |
void |
setTSB(double[] in) |
void |
setTSD(double in) |
void |
setTSS(double[] in) |
public ConfidenceLevel()
public ConfidenceLevel(int mc)
mc
- number of MonteCarlo experimentspublic ConfidenceLevel(int mc, boolean onesided)
mc
- is the number of Monte Carlo experimentsonesided
- specifies if the intervals are one-sided or not.public double getExpectedStatistic_b(int sigma)
sigma
- between -2 and 2public double getExpectedStatistic_sb(int sigma)
sigma
- public double getCLb()
public double getCLb(boolean use_sMC)
But forming discovery p-values, we must compute 1-CLb values of the order of 10-7. This computation involves generating of the order of 10^8 pseudoexperiments, just to be on the safe side.
use_sMC
- public double getCLsb()
public double getCLsb(boolean use_sMC)
use_sMC
- public double getCLs()
This hypothesis is excluded at the 95% CL if CLs = 0.05, and at more than the 95% CL if CLs < 0.05, assuming that signal is present.
public double getCLs(boolean use_sMC)
This hypothesis is excluded at the 95% CL if CLs = 0.05, and at more than the 95% CL if CLs < 0.05, assuming that signal is present.
use_sMC
- use or not MC.public double getExpectedCLsb_b(int sigma)
sigma
- public double getExpectedCLb_sb(int sigma)
sigma
- public double getExpectedCLb_b(int sigma)
sigma
- public double getAverageCLsb()
public double getAverageCLs()
public double get3sProbability()
public double get5sProbability()
public ArrayList<H1D> getResults(String Option)
Option
- public void setTSD(double in)
public void setLRS(double[] in)
public void setLRB(double[] in)
public void setBtot(double in)
public void setStot(double in)
public void setDtot(int in)
public double getStatistic()
public void setTSB(double[] in)
public void setTSS(double[] in)
public double getExpectedStatistic_b()
public double getExpectedCLb_sb()
public double getExpectedCLs_b(int sigma)
sigma
- public double getExpectedCLs_b()
sigma
- public double getExpectedCLb_b()
public double getExpectedCLsb_b()
public H1D getLNQb(int bins, double min, double max)
bins
- number of binsmin
- min valuemax
- max valuepublic H1D getLNQsb(int bins, double min, double max)
bins
- number of binsmin
- min valuemax
- max valuepublic void doc()
Jas4pp 1.5 © Java Analysis Studio for Particle Physics