// File:  getExclusionSignificance.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.
// S1    sum of weights, w
// S2    sum of squares of w
// tau   defined by m_i ~ Poisson(tau_i*b/omega), where omega = E[w]

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

#include "SigCalc.h"

double getExclusionSignificance(double mu, double n, double s,
       vector<double> m, vector<double> tau,
       vector<double> S1, vector<double> S2){

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

  delete sc;
  return Z;

}