23 int main(
int argc,
char **argv)
28 JFitK40Parameters fitk40;
43 JParser<> zap(
"Auxiliary program to draw K40 fit function.");
51 catch(
const exception &error) {
52 FATAL(error.what() << endl);
60 for (JModule::const_iterator i = module.begin(); i != module.end(); ++i) {
61 for (JModule::const_iterator
j = module.begin();
j != i; ++
j) {
63 if (
getDot(i->getDirection(),
j->getDirection()) > ct) {
64 ct =
getDot(i->getDirection(),
j->getDirection());
69 cout <<
"Cosine minimal angle between PMT axes " <<
FIXED(5,3) << ct << endl;
70 cout <<
"Default rate [Hz] " <<
FIXED(6,3) << JFitK40Parameters().getValue(ct) <<
'/' <<
FIXED(6,3) << JFitK40Parameters().getValue(+1.0) << endl;
71 cout <<
"Modified rate [Hz] " <<
FIXED(6,3) << fitk40 .getValue(ct) <<
'/' <<
FIXED(6,3) << fitk40 .getValue(+1.0) << endl;
75 TH1D h0(
"h0", NULL, 5000, -1.0, +1.0);
77 for (Int_t i = 1; i <= h0.GetXaxis()->GetNbins(); ++i) {
79 const Double_t x = h0.GetBinCenter(i);
81 h0.SetBinContent(i, fitk40.getValue(x));