double EvtDecayBase::getProbMax( double prob ) {

  int i;

  //diagnostics
  sum_prob+=prob;
  if (prob>max_prob) max_prob=prob;


  if ( defaultprobmax && ntimes_prob<=500 ) { 
    //We are building up probmax with this iteration
     ntimes_prob += 1;
     if ( prob > probmax ) { probmax = prob;}
     if (ntimes_prob==500) { 
       probmax*=1.2;
     }
     return 1000000.0*prob;
  }

  if ( prob> probmax*1.0001) {

    report(INFO,"EvtGen") << "prob > probmax:("<<prob<<">"<<probmax<<")";
    report(INFO,"") << "("<<_modelname.c_str()<<") ";
    report(INFO,"") << EvtPDL::name(_parent).c_str()<<" -> ";
    for(i=0;i<_ndaug;i++){
       report(INFO,"") << EvtPDL::name(_daug[i]).c_str() << " ";
    }
    report(INFO,"") << endl;

    if (defaultprobmax) probmax = prob;

  }

  ntimes_prob += 1;


  return probmax;

} //getProbMax