22 int main(
int argc,
char **argv)
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.
Optical properties of Antares deep-sea site.
double getAbsorptionLength(const double lambda)
Absorption length.
double getScatteringProbability(const double x)
Function to describe light scattering in water.
Optical properties of KM3NeT deep-sea site.
#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)
General purpose messaging.
double getAbsorptionLength(const double lambda)
Absoption length.
Utility class to parse command line options.
double getScatteringLength(const double lambda)
Scattering length.
double f4(const double x)
Model specific function to describe light scattering in water (f4)
int main(int argc, char *argv[])