// File:  getSignificance.cc
// Glen Cowan
// RHUL Physics

// Returns significance Z with which hypothesized value of mu is rejected; 
// Z = Phi^-1(1 - p), p = p-value of mu, Phi^-1 = quantile of standard Gaussian. 
// 
// Inputs:     
// mu    signal strength parameter (0 = background only, 1 = nominal).
// n     number of events seen in data, e.g. generate from Poisson(mu*s+b)
// s     expected number of signal events
// m     vector of numbers of events seen in subsidiary measurements.
// tau   defined by m_i ~ Poisson(tau_i*b_i), e.g., ratio of MC lumonisity
//        to that of data sample.

#include <iostream>
#include <fstream>
#include <sstream>
#include <vector>
#include <cstdlib>
#include <cmath>

#include "SigCalc.h"
#include "getSignificance.h"

double getSignificance(double mu, double n, double s,
                       vector<double> m, vector<double> tau){

  SigCalc* sc = new SigCalc (n, s, m, tau, 0);
  double qmu = sc->qmu(mu);
  double Z = sqrt(qmu);

  delete sc;
  return Z;

}