---

m_bar.C

Go to the documentation of this file.

#include "stdinc.h"

#ifdef TOOLBOX

main(int argc, char *argv[])
{
    real ml = -1.0, mu = -1.0, x = -2.35, factor = 1;
    bool verbose = 0;

    check_help();

    extern char *poptarg;
    int c;
    char* param_string = "e:f:l:L:u:U:vx:";

    while ((c = pgetopt(argc, argv, param_string)) != -1)
        switch(c) {
              case 'e':
              case 'x': x = atof(poptarg);
                        break;
              case 'f': factor = atof(poptarg);
                        break;
              case 'l':
              case 'L': ml = atof(poptarg);
                        break;
              case 'u':
              case 'U': mu = atof(poptarg);
                        break;
              case 'v': verbose = true;
                        break;
              case '?': params_to_usage(cerr, argv[0], param_string);
                        get_help();
                        exit(1);
          }

    if (ml <= 0 || mu <= 0) exit(2);

    if (mu < ml) {
        real tmp = ml;
        ml = mu;
        mu = tmp;
    }

    real mbar = 1;

    if (ml < mu) {
        if (x != -1) {
            mbar = ml * (1+x) * (pow(mu/ml, 2+x)-1)
                                  / (pow(mu/ml, 1+x)-1) / (2+x);
        } else {
            mbar = (mu - ml) / log(mu/ml);
        }
    }

    if (verbose) cout << "ml = " << ml << ", mu = " << mu << ", <m> = ";

    cout << mbar*factor << endl;
}

#endif

---