26 int main(
int argc,
char **argv)
46 JParser<> zap(
"Auxiliary program to draw K40 fit function.");
54 catch(
const exception &error) {
55 FATAL(error.what() << endl);
59 const JModule module = getModule<JKM3NeT_t>(1001);
63 for (JModule::const_iterator
i = module.begin();
i != module.end(); ++
i) {
64 for (JModule::const_iterator
j = module.begin();
j !=
i; ++
j) {
66 if (
getDot(
i->getDirection(),
j->getDirection()) > ct) {
67 ct =
getDot(
i->getDirection(),
j->getDirection());
72 const JCombinatorics_t combinatorics(module);
74 cout <<
"Cosine minimal angle between PMT axes " <<
FIXED(5,3) << ct << endl;
76 cout <<
"Modified rate [Hz] " <<
FIXED(6,3) << K40 .getValue(ct) <<
'/' <<
FIXED(6,3) << K40 .getValue(+1.0) << endl;
80 const Int_t nx = combinatorics.getNumberOfPairs();
81 const Double_t
xmin = -0.5;
82 const Double_t
xmax = nx - 0.5;
84 TH1D h0(
"h0", NULL, 5000, -1.0, +1.0);
87 for (Int_t
i = 1;
i <= h0.GetXaxis()->GetNbins(); ++
i) {
89 const Double_t
x = h0.GetBinCenter(
i);
91 h0.SetBinContent(
i, K40.getValue(x));
94 for (Int_t
i = 1;
i <= h1.GetXaxis()->GetNbins(); ++
i) {
96 const int index = h1.GetBinCenter(
i);
100 const double x =
getDot(module[pair.first].getDirection(), module[pair.second].getDirection());
102 h1.SetBinContent(
i, K40.getValue(x));
Utility class to parse command line options.
JCombinatorics::pair_type pair_type
int main(int argc, char *argv[])
#define gmake_property(A)
macros to convert (template) parameter to JPropertiesElement object
double getDot(const JNeutrinoDirection &first, const JNeutrinoDirection &second)
Dot product.
Utility class to parse parameter values.
Empty structure for specification of parser element that is initialised (i.e. does not require input)...
Auxiliary data structure for floating point format specification.
Utility class to parse parameter values.
I/O formatting auxiliaries.
T & getInstance(const T &object)
Get static instance from temporary object.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
General purpose messaging.
Utility class to parse command line options.
Data structure for optical module.