31 JParser<> zap(
"Auxiliary program histogram plot absorption and scattering lengths as well as dispersion.");
38 catch(
const exception &error) {
39 FATAL(error.what() << endl);
47 TH1D h0(
"abs[Antares]", NULL, 400, 280.0, 680.0);
48 TH1D h1(
"abs[KM3NeT]", NULL, 400, 280.0, 680.0);
49 TH1D h2(
"scat[Antares]", NULL, 400, 280.0, 680.0);
50 TH1D h3(
"scat[KM3NeT]", NULL, 400, 280.0, 680.0);
52 TH1D h4(
"vp[Antares]", NULL, 400, 280.0, 680.0);
53 TH1D h5(
"vp[KM3NeT]", NULL, 400, 280.0, 680.0);
54 TH1D h6(
"vg[Antares]", NULL, 400, 280.0, 680.0);
55 TH1D h7(
"vg[KM3NeT]", NULL, 400, 280.0, 680.0);
57 TH1D h8(
"Ps[Antares]", NULL, 800, -1.0, +1.0);
58 TH1D h9(
"Ps[KM3NeT]", NULL, 800, -1.0, +1.0);
59 TH1D ha(
"f4", NULL, 800, -1.0, +1.0);
60 TH1D hb(
"p00075", NULL, 800, -1.0, +1.0);
62 JDispersion antares(250.0);
63 JDispersion km3net (350.0);
65 for (
int i = 1; i <= h0.GetNbinsX(); ++i) {
67 const double x = h0.GetBinCenter (i);
74 h4.SetBinContent(i, antares.getIndexOfRefractionPhase(x));
75 h5.SetBinContent(i, km3net .getIndexOfRefractionPhase(x));
76 h6.SetBinContent(i, antares.getIndexOfRefractionGroup(x));
77 h7.SetBinContent(i, km3net .getIndexOfRefractionGroup(x));
80 for (
int i = 1; i <= h8.GetNbinsX(); ++i) {
82 const double x = h8.GetBinCenter (i);
Utility class to parse command line options.
double getScatteringLength(const double lambda)
Scattering length.
double getScatteringProbability(const double x)
Function to describe light scattering in water.
double getAbsorptionLength(const double lambda)
Absorption length.
double getScatteringProbability(const double x)
Function to describe light scattering in water.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
double p00075(const double x)
Model specific function to describe light scattering in water (p00075)
double getAbsorptionLength(const double lambda)
Absoption length.
double getScatteringLength(const double lambda)
Scattering length.
double f4(const double x)
Model specific function to describe light scattering in water (f4)