//  Exemplo de "goodness of fit", Cowan, pg. 57
//  Joga-se uma moeda N vezes. Dados nh (caras) e N-nh (coroas), elas sao
//  consistentes com a hipotese de moeda "sem armacao" ?
//  Suponha N = 20 e nh = 17 
//  Pval = Soma f(n,0.5) , n=0,1,2,3,17,18,19,20
//
//  aqui se supoe sempre o mesmo numero N de moedas jogadas em um exp.
//                                 jtmn, Set. 2003

 Double_t binom(Int_t nh, Int_t N, Double_t p)
   {
     Double_t fval = TMath::Factorial(N)/TMath::Factorial(nh);
     fval = fval/TMath::Factorial(N-nh);
     fval = fval*TMath::Power(p,nh)*TMath::Power(p,N-nh);
     return fval;
   }

void moeda1(){

 gROOT->Reset();
 
 Int_t N = 20;
 Double_t Pval=0.0;
 Double_t p = 0.5;
 for (Int_t i=0; i<=3; i++){
   // para p=1/2 os dois binom abaixo sao iguais, poderia mult. por dois
   // ao inves de somar
   Pval = Pval + binom(i,N,p) + binom(N-i,N,p);
 }
 cout << "N = " << N << ", p = " << p << " Pval = " << Pval << endl;
}