1 #ifndef __JCOMPAREHISTOGRAMS__JTESTCHI2_T__
2 #define __JCOMPAREHISTOGRAMS__JTESTCHI2_T__
14 namespace JCOMPAREHISTOGRAMS {
50 if(h1 -> GetNbinsX() != h2 -> GetNbinsX())
51 ERROR(
"Histograms with different bining. The objects: " << h1 -> GetName() <<
" can not be compared." << endl);
53 double chi2 = h1 -> Chi2Test (h2 , options.c_str());
55 double M = h1->Integral();
56 double N = h2->Integral();
58 TH2D* h3 = (TH2D*)h1->Clone(h1->GetName()==h2->GetName() ?
64 for (
int i=1 ; i < h1->GetNbinsX() ; ++i){
66 double m = h1->GetBinContent(i);
67 double n = h2->GetBinContent(i);
70 double c = (
M*n -
N*m)/sqrt((n+m)*(
N*
M));
71 h3->SetBinContent(i,c);
77 if (options.find(
"CHI2") != std::string::npos) {
78 (chi2 > threshold ? passed =
false : passed =
true);
80 (chi2 < threshold ? passed =
false : passed =
true);
83 JResultTitle title(testName, parameterName, passed , chi2);
85 h3->SetTitle(title.
getTitle().c_str());
90 parameterName, chi2, threshold, h3, passed);
124 if(slice ==
'x' || slice ==
'X'){
126 int nSlices1 = h1->GetNbinsX();
127 int nSlices2 = h2->GetNbinsX();
129 TH1* h3 = h1->ProjectionX(h1->GetName()==h2->GetName() ?
133 if(nSlices1 != nSlices2)
134 ERROR(
"Histograms with different binning. The objects: " << h1->GetName() <<
" and " << h2->GetName() <<
" can not be compared." << endl);
136 for (
int i=1 ; i<=nSlices1 ; ++i){
140 TH1D* s1 = h1->ProjectionY (sliceName.c_str(),i,i);
141 TH1D* s2 = h2->ProjectionY (sliceName.c_str(),i,i);
143 double chi2 = s1 -> Chi2Test (s2 , options.c_str());
147 if (options.find(
"CHI2") != std::string::npos) {
148 (chi2 > threshold ? passed =
false : passed =
true);
150 (chi2 < threshold ? passed =
false : passed =
true);
153 if (!passed) nFailures++;
155 h3->SetBinContent(i,chi2);
161 (nFailures/nSlices1 > failuresThreshold ? passed =
false : passed =
true);
163 JResultTitle title(testName, parameterName, passed , nFailures);
165 h3->SetTitle(title.
getTitle().c_str());
170 parameterName, nFailures, failuresThreshold, h3, passed);
172 }
else if (slice ==
'y' || slice ==
'Y'){
174 int nSlices1 = h1->GetNbinsY();
175 int nSlices2 = h2->GetNbinsY();
177 TH1* h3 = h1->ProjectionY(h1->GetName()==h2->GetName() ?
181 if(nSlices1 != nSlices2)
182 ERROR(
"Histograms with different binning. The objects: " << h1->GetName() <<
" can not be compared." << endl);
184 for (
int i=1 ; i<=nSlices1 ; ++i){
188 TH1D* s1 = h1->ProjectionX (sliceName.c_str(),i,i);
189 TH1D* s2 = h2->ProjectionX (sliceName.c_str(),i,i);
191 double chi2 = s1 -> Chi2Test (s2 , options.c_str());
195 if (options.find(
"CHI2") != std::string::npos) {
196 (chi2 > threshold ? passed =
false : passed =
true);
198 (chi2 < threshold ? passed =
false : passed =
true);
201 if (!passed) nFailures++;
203 h3->SetBinContent(i,chi2);
208 (nFailures/nSlices1 > failuresThreshold ? passed =
false : passed =
true);
210 JResultTitle title(testName, parameterName, passed , nFailures);
212 h3->SetTitle(title.
getTitle().c_str());
217 parameterName, nFailures, failuresThreshold, h3, passed);
245 int nx1 = h1->GetNbinsX();
246 int nx2 = h2->GetNbinsX();
247 int ny1 = h1->GetNbinsY();
248 int ny2 = h2->GetNbinsY();
250 double M = h1->Integral();
251 double N = h2->Integral();
253 if(nx1 != nx2 || ny1 != ny2 ||
M == 0 ||
N == 0)
254 ERROR(
"Histograms with different binning. The objects: " << h1->GetName() <<
" can not be compared." << endl);
256 TH2D* h3 = (TH2D*)h1->Clone(h1->GetName()==h2->GetName() ?
264 for (
int i=1 ; i<nx1 ; ++i){
265 for (
int j=1 ;
j<ny1 ; ++
j){
267 double m = h1 -> GetBinContent(i,
j);
268 double n = h2 -> GetBinContent(i,
j);
269 double chi2 = (n-m*
N/
M)/sqrt(m*
N/
M);
270 (fabs(chi2) > chi2Threshold ? outliers+=1./(nx1*ny1) : outliers+=0 );
271 h3->SetBinContent(i,
j,chi2);
277 (outliers > outliersThreshold ? passed =
false : passed =
true);
279 JResultTitle title(testName, parameterName , passed , 100*outliers);
281 h3->SetTitle(title.
getTitle().c_str());
286 parameterName, 100*outliers, 100*outliersThreshold, h3, passed);
JTestResult JChi2TestSlice(TH2 *h1, TH2 *h2, double threshold, double failuresThreshold, std::string testName, std::string parameterName, std::string options, char slice)
Chi2 test for sliced 2D histograms.
JTestResult JChi2Test(TH1 *h1, TH1 *h2, double threshold, std::string testName, std::string parameterName, std::string options)
Chi2 test for 1D histograms.
Class dedicated to standardize the title of the graphical objects produced by the JTest_t() derived c...
JTestResult JChi2TestBin_2D(TH2 *h1, TH2 *h2, double outliersThreshold, double chi2Threshold, std::string testName, std::string parameterName)
Bin-by-Bin Chi2 comparison of 2D histograms.
then JShowerPostfit f $INPUT_FILE o $OUTPUT_FILE N
#define MAKE_CSTRING(A)
Make C-string.
then for HISTOGRAM in h0 h1
Auxiliary class to handle file name, ROOT directory and object name.
#define MAKE_STRING(A)
Make string.
JTestChi2_t()
Default constructor.
then usage $script[port]< option > nPossible options
std::string to_string(const T &value)
Convert value to string.
Structure containing the result of the test.
$WORKDIR ev_configure_domsimulator txt echo process $DOM_SIMULATOR $i $SOURCE_HOST[$index] csh c(setenv ROOTSYS $ROOTSYS &&source $JPP_DIR/setenv.csh $JPP_DIR &&($DOM_SIMULATOR\-u\$NAME\$\-H\$SERVER\$\-M\$LOGGER\$\-d $DEBUG</dev/null > &/dev/null &))'
std::string getTitle()
Returns a standard string to be used as title of a graphical root object.
Implementation of the different Chi2-related tests.