Log of /trunk/hist/hist/src/TLimit.cxx
Parent Directory
Revision
48992 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Thu Mar 28 15:26:26 2013 UTC (21 months, 3 weeks ago) by
rdm
File length: 17997 byte(s)
Diff to
previous 44507
From Lifeng Sun:
The attached patchset fixes a bunch of typo in the source:
0001-succes-success.patch
0002-preceed-preced.patch
0003-informations-information.patch
0004-childs-children.patch
0005-avaliable-available.patch
0006-writeable-writable.patch
0007-comand-command.patch
0008-unkown-unknown.patch
0009-wierd-weird.patch
0010-wheter-whether.patch
0011-unecessary-unnecessary.patch
0012-splitted-split.patch
0013-registerd-registered.patch
0014-recieve-receive.patch
0015-processsing-processing.patch
0016-ouput-output.patch
0017-mutiple-multiple.patch
0018-lenght-length.patch
0019-interupted-interrupted.patch
0020-independant-independent.patch
0021-inconsistant-inconsistent.patch
0022-expresion-expression.patch
0023-explicitely-explicitly.patch
0024-enviroment-environment.patch
0025-deafult-default.patch
0026-continous-continuous.patch
0027-completly-completely.patch
0028-commited-committed.patch
0029-choosen-chosen.patch
0030-backgroud-background.patch
0031-auxilliary-auxiliary.patch
0032-authentification-authentication.patch
0033-appropiate-appropriate.patch
0034-an-other-another.patch
0035-environement-environment.patch
0036-targetting-targeting.patch
0037-suppported-supported.patch
0038-paramater-parameter.patch
Revision
44507 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Mon Jun 4 12:30:41 2012 UTC (2 years, 7 months ago) by
axel
File length: 17998 byte(s)
Diff to
previous 26221
Remove
using namespace std;
from Riostream.h, which has huge consequences for all of ROOT.
Riostream.h is now a simple wrapper for fstream, iostream, iomanip for backward compatibility; Riosfwd.h simply wraps iosfwd.
Because of templates and their inline functions, Riostream.h needed to be included in headers, too (e.g. TParameter.h), which violated the assumption that Riostream.h is not exposing its using namespace std to headers.
ROOT now requires R__ANSISTREAM, R__SSTREAM, which does not change the set of supported compilers.
Without "using namespace std", several identifiers are now prefixed by std::; e.g. roofit/* source files now have a using namespace std to keep their coding style.
TFile::MakeProject() now generates "using namespace std" to convert the CINT-style class names into C++ ones.
Revision
26221 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Mon Nov 17 08:10:26 2008 UTC (6 years, 2 months ago) by
brun
File length: 17918 byte(s)
Diff to
previous 22694
From Antje Huettmann huettmaa@mail.desy.de
Fix a memory leak when calling ComputeLimit many times.The initialization of serrf and berrf is much more critical in cases where you have to retoss many times (i.e. when you have large errors).
Revision
16530 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Sun Oct 15 06:23:48 2006 UTC (8 years, 3 months ago) by
brun
Original Path:
trunk/hist/src/TLimit.cxx
File length: 17825 byte(s)
Diff to
previous 15181
From Christophe Delaere:
new version of TLimit. It implements Tom Junk's solution for a
bias appearing sometimes when introducing systematics.
Bellow is the detailed information.
Following an exchange of mails between Laurent Duflot and Tom Junk,
Tom wrote the following:
"Laurent Duflot has given me an example which showed me I was wrong with
how I reweight s+b pseudoexperiments to get background-only
pseudoexperiments and vice versa. I am sending this along to Alex Read
and Bill Murray in case they might say "Yes, you should have known that
all along, Tom!". It certainly was a surprise to me.
We do "importance sampling" of s+b and b pseudoexperiments and reweight
them to look like the other kind, and the weight factor is the likelihood
ratio, which coincidentally happens to look like our test statistic.
But the test statistic uses central values of s and b so as not to split
outcomes, and the weights for reweighting pseudoexperiments I was using
were the likelihood ratios using systematically varied s and b on each
pseudoexperiment, where I pick smeared values on each pseudoexperiment.
The problem is that this is biased. Pseudoexperiments in the
background-hypothesis distribution that are more signal-like
preferentially have their background expectations systematically varied
upwards. Pseudoexperiments in the s+b distribution that are more
background-like preferentially have their background expectations
systematically varied downwards. These two kinds of outcomes may
populate the same bin of -2lnQ. Reweighting the bg pseudoexperiments
to look like corresponding s+b pseudoexperiments should take into
account that the corresponding backgrounds are different.
I've approximated this by re-fluctuating the signal and the background
on each pseudoexperiment to get an unbiased target distribution. So the
weight ratio is
Ratio = P(pd|s1+b1)/P(pd|b2)
where pd=pseudodata, s1, b1 are a randomly re-fluctuated signal and
background, and b2 is the fluctuated background which gave me that
pseudodata, and this weight should turn a biased background
pseudoexperiment into a representative s+b pseudoexperiment.
The problem I see now is that the s+b pseudoexperiment should also be
biased, but I do not know how to compute this without running a lot of
pseudoexperiments. It also gets tough as the number of variable
systematic errors gets large.
Nonetheless, the approximation above helps make the reweighted quantities
match more their unweighted counterparts in the large-pseudoexperiment
limit.
Revision
10219 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Tue Oct 5 17:41:08 2004 UTC (10 years, 3 months ago) by
brun
Original Path:
trunk/hist/src/TLimit.cxx
File length: 14726 byte(s)
Diff to
previous 8396
From Christophe Delaere
I changed the indexing so that (NbinsX+2) bins are used for each histogram.
When I loop on histograms, I then go from 0 to and including (NbinsX+1).
I also removed the unused nbins variable.
The old version was only partially wrong: the overflow bin was not used,
but the memory corruption problem arrised only in some rare cases were
maxbins was wrongly computed because of a missing "+1" at line 127.
Revision
8049 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Tue Jan 27 13:41:53 2004 UTC (10 years, 11 months ago) by
brun
Original Path:
trunk/hist/src/TLimit.cxx
File length: 14477 byte(s)
Diff to
previous 6337
From Christophe Delaere:
TLimit has been modified to take into account statistical errors from the MC.
Tutorial limit.C modified accordingly.
In TMultiLayerPerceptron, add the possibilite to export the NN
in Python format.
Corrected a problem with TEventlist when 2 NN are in the same directory.
Update the example mlpHiggs.C
The web page describing this work has also been updated. see:
http://www.fynu.ucl.ac.be/users/c.delaere/MLP/
Revision
5225 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Fri Sep 6 20:24:18 2002 UTC (12 years, 4 months ago) by
brun
Original Path:
trunk/hist/src/TLimit.cxx
File length: 12937 byte(s)
Diff to
previous 5222
Add the following comments in the class documentation:
// class TLimit
// ------------
//
// Algorithm to compute 95% C.L. limits using the Likelihood ratio
// semi-bayesian method.
// It takes signal, background and data histograms wrapped in a
// TLimitDataSource as input and runs a set of Monte Carlo experiments in
// order to compute the limits. If needed, inputs are fluctuated according
// to systematics. The output is a TConfidenceLevel.
//
// class TLimitDataSource
// ----------------------
//
// Takes the signal, background and data histograms as well as different
// systematics sources to form the TLimit input.
//
// class TConfidenceLevel
// ----------------------
//
// Final result of the TLimit algorithm. It is created just after the
// time-consuming part and can be stored in a TFile for further processing.
// It contains light methods to return CLs, CLb and other interesting
// quantities.
This form allows you to request diffs between any two revisions of this file.
For each of the two "sides" of the diff,
enter a numeric revision.
