new way of calculation Nsigma correction

[u/mrichter/AliRoot.git] / PWG0 / CorrectionMatrix2D.cxx

diff --git a/PWG0/CorrectionMatrix2D.cxx b/PWG0/CorrectionMatrix2D.cxx
index fdfdd79f140a46af7c99bd3e98dae7fbd98a41ff..59f1dde84e63380b39a76aa489643496a0484de1 100644 (file)
--- a/PWG0/CorrectionMatrix2D.cxx
+++ b/PWG0/CorrectionMatrix2D.cxx
@@ -39,6 +39,10 @@ CorrectionMatrix2D::CorrectionMatrix2D(Char_t* name, Char_t* title,
   fhMeas  = new TH2F(Form("meas_%s",name), Form("meas_%s",title),  nBinX, Xmin, Xmax, nBinY, Ymin, Ymax);
   fhGene  = new TH2F(Form("gene_%s",name), Form("gene_%s",title),  nBinX, Xmin, Xmax, nBinY, Ymin, Ymax);
   fhCorr  = new TH2F(Form("corr_%s",name), Form("corr_%s",title),  nBinX, Xmin, Xmax, nBinY, Ymin, Ymax);
+
+  fhMeas->Sumw2();
+  fhGene->Sumw2();
+  fhCorr->Sumw2();
 }
 
 //____________________________________________________________________
@@ -53,6 +57,10 @@ CorrectionMatrix2D::CorrectionMatrix2D(Char_t* name,Char_t* title,
   fhMeas  = new TH2F(Form("meas_%s",name), Form("meas_%s",title),  nBinX, X, nBinY, Y);
   fhGene  = new TH2F(Form("gene_%s",name), Form("gene_%s",title),  nBinX, X, nBinY, Y);
   fhCorr  = new TH2F(Form("corr_%s",name), Form("corr_%s",title),  nBinX, X, nBinY, Y);
+
+  fhMeas->Sumw2();
+  fhGene->Sumw2();
+  fhCorr->Sumw2();
 }
 
 
@@ -83,8 +91,8 @@ TH1F* CorrectionMatrix2D::Get1DCorrection(Char_t* opt) {
   // integrate the correction over one variable 
   // 
 
-  TH1D* meas1D; 
-  TH1D* gene1D; 
+  TH1D* meas1D = 0;
+  TH1D* gene1D = 0; 
 
   if (strcmp(opt,"x")==0) {
     meas1D = fhMeas->ProjectionX();