void printSignificance(TH1F sig, TH1F bkg) { TFile *f = TFile::Open("significance1.root"); TH1F *S = (TH1F *)f->Get("sig"); TH1F *B = (TH1F *)f->Get("bkg"); for (int i=0; i<=S->GetNbinsX(); i++){ for (int j=0;i<=B->GetNbinsX(); j++){ Double_t S = histosignal->Integral(i,j); Double_t B = histobackground->Integral(i,j); //Double_t Sig1 = S/sqrt(B); Double_t Sig2 = S/sqrt(S+B); effsig=sig->Integral(i,j)/sig->integral(); effbkg=bkg->Integral(i,j)/bkg->integral(); if( B != 0. && (S+B) != 0. ) { // Sig = S/sqrt(B); Sig = S/sqrt(S+B); } std::cout << " effsig " << effsig() << " effbkg " << effbkg << " Sig" << Sig << std::endl; } }