MultiplicityInfo.h

// -*- C++ -*-
//
// MultiplicityInfo.h is a part of Herwig - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2019 The Herwig Collaboration
//
// Herwig is licenced under version 3 of the GPL, see COPYING for details.
// Please respe
#include "Herwig/Utilities/Statistic.h"
 
namespace Herwig {
using namespace ThePEG;
 
enum ParticleSpecies {
  lightMeson=0,strangeMeson,lightBaryon,other
};
 
struct MultiplicityInfo {
 
  MultiplicityInfo(double mult=0.,double error=0.,
                   ParticleSpecies type=other)
    : obsMultiplicity(mult), obsError(error), type(type) {};
 
  double obsMultiplicity;
 
  double obsError;
 
  ParticleSpecies type;
 
  Statistic count;
 
  double simMultiplicity() { return count.mean();}
 
  double simError() { return count.mean_stdDev();}
 
  double nSigma() {
    return obsMultiplicity == 0.0 ? 0.0 :
      (simMultiplicity() - obsMultiplicity) 
      / sqrt(sqr(simError()) + sqr(obsError));
  }
 
  string bargraph() {
    if (obsMultiplicity == 0.0) return "     ?     ";
    else if (nSigma() >= 6.0)   return "-----|---->";
    else if (nSigma() >= 5.0)   return "-----|----*";
    else if (nSigma() >= 4.0)   return "-----|---*-";
    else if (nSigma() >= 3.0)   return "-----|--*--";
    else if (nSigma() >= 2.0)   return "-----|-*---";
    else if (nSigma() >= 1.0)   return "-----|*----";
    else if (nSigma() > -1.0)   return "-----*-----";
    else if (nSigma() > -2.0)   return "----*|-----";
    else if (nSigma() > -3.0)   return "---*-|-----";
    else if (nSigma() > -4.0)   return "--*--|-----";
    else if (nSigma() > -5.0)   return "-*---|-----";
    else if (nSigma() > -6.0)   return "*----|-----";
    else                        return "<----|-----";
  }
};
}